「統計学関連なんでもあり」の過去ログ--- 043

No.12032　Rのdrop1()について　　【yama】　2010/02/11(Thu) 22:49

よろしくお願いします。
ロジスティック回帰でdrop1(モデル)を指定したときに示される，χ2とPは何を表していますでしょうか。
また，使われているのは尤度比検定であっていますでしょうか。

No.12039　Re: Rのdrop1()について　　【青木繁伸】　2010/02/12(Fri) 09:29

drop1, add1 のオンラインヘルプの test 引数の説明は以下のようになっていますよ

test

should the results include a test statistic relative to the original model? The F test is only appropriate for lm and aov models or perhaps for glm fits with estimated dispersion. The Chisq test can be an exact test (lm models with known scale) or a likelihood-ratio test or a test of the reduction in scaled deviance depending on the method.

No.12047　Re: Rのdrop1()について　　【yama】　2010/02/13(Sat) 17:53

青木先生，ありがとうございます。
ちょっと分からないところがあるので質問してもよろしいでしょうか。
複数の説明変数を含んだ一般化線形モデルで，二項分布を指定してdrop1()を行いました。
χ2検定の結果の各説明変数が自由度が1になっていて，各変数に対してp値が出ています。
これは，null model に対して各変数ひとつずつのモデルとを尤度比検定したと解釈していいのでしょうか。

No.12048　Re: Rのdrop1()について　　【青木繁伸】　2010/02/13(Sat) 18:19

H0:回帰係数=0 の検定では？

該当部分をコピー・ペーストしてくれるとはっきりする

No.12051　Re: Rのdrop1()について　　【yama】　2010/02/13(Sat) 19:52

青木先生，ありがとうございます。
該当部分は以下の通りです。
よろしくお願いします。
single term deletions

Model:
result ~ A+B+C
       Df Deviance    AIC     LRT   Pr(Chi)    
      27.525 35.525                      
A       1   27.863 33.863  0.3388 0.5605402    
B       1   39.155 45.155 11.6303 0.0006489 ***
C       1   27.628 33.628  0.1032 0.7480354  

No.12054　Re: Rのdrop1()について　　【青木繁伸】　2010/02/13(Sat) 21:09

どれか一つの変数を除いたときのモデルの適合度ですかね。上の例で言えば，Bを除くと有意な変化がある，つまり，Bは重要な変数であると。。。（結果を示すときに，その結果を導いたプログラムもあれば。データまでとは言わないので）

No.12056　Re: Rのdrop1()について　　【yama】　2010/02/13(Sat) 21:32

青木先生，ありがとうございます。
プログラムは以下の通りです。
よろしくお願いいたします。

> model<-glm(formula=result~A+B+C,binomial)
> drop1(model,test="Chisq")

No.12061　Re: Rのdrop1()について　　【青木繁伸】　2010/02/13(Sat) 21:49

オンラインヘルプは読んでいますか？

add1 / drop1 のオンラインヘルプに曰く

Description

Compute all the single terms in the scope argument that can be added to or dropped from the model, fit those models and compute a table of the changes in fit.

No.12183　Re: Rのdrop1()について　　【yama】　2010/02/26(Fri) 16:12

返信が遅くなってしまい申し訳ございません。
ありがとうございます。
ちゃんと読んでみます。