★ p値のみかた。 ★

 382 p値のみかた。  ひで  2001/10/10 (水) 13:53
  383 Re: p値のみかた。  青木繁伸  2001/10/10 (水) 16:03
   384 Re^2: p値のみかた。  青木繁伸  2001/10/10 (水) 16:21
    400 Re^3: p値のみかた。  ともも  2001/10/16 (火) 21:37
     402 Re^4: p値のみかた。  青木繁伸  2001/10/16 (火) 22:14
      410 Re^5: p値のみかた。  ともも  2001/10/17 (水) 01:38
    387 Re^3: p値のみかた。(単なる感想)  DISIR  2001/10/12 (金) 11:09
     388 Re^4: p値のみかた。(単なる感想)  青木繁伸  2001/10/12 (金) 15:45
      390 Re^5: p値のみかた。(単なる感想)教えて下さい。  sb812109  2001/10/13 (土) 07:48
       391 Re^6: p値のみかた。(単なる感想)教えて下さい。  青木繁伸  2001/10/13 (土) 20:23
        392 Re^7: p値のみかた。(単なる感想)教えて下さい。  sb812109  2001/10/14 (日) 10:19
         396 Re^8: p値のみかた。(単なる感想)教えて下さい。  sb812109  2001/10/15 (月) 22:52
          397 Re^9: p値のみかた。(単なる感想)教えて下さい。  青木繁伸  2001/10/16 (火) 11:02


382. p値のみかた。  ひで  2001/10/10 (水) 13:53
STATA6@を使用していますが,単ロジスティック回帰分析にて有意性を見る場合 Prob > chi2 , P>|z|のどちらを見ればよいもんなんでしょうか。
2*2表なんですが,χの2乗のP値とも微妙にちがいます。
あとzとはどのようなもんなんでしょうか。
logistic bitem sex
 Prob > chi2     =     0.0028
 Pseudo R2       =     0.0089
   bitem | Odds Ratio   Std. Err.       z     P>|z|       
     sex |    .64        .0957243    -2.984  0.003

     [このページのトップへ]


383. Re: p値のみかた。  青木繁伸  2001/10/10 (水) 16:03
> STATA6@を使用していますが

マニュアルないのですか?

http://www.ats.ucla.edu/stat/stata/output/lognoframe.htm
でも見るといいのではないでしょうか?
出力中のそれぞれの数値が何を意味しているか説明してくれていますよ。

つづく

     [このページのトップへ]


384. Re^2: p値のみかた。  青木繁伸  2001/10/10 (水) 16:21
> 単ロジスティック回帰分析にて有意性を見る場合 Prob > chi2 ,P>|z|のどちらを見ればよいもんなんでしょうか。

それぞれ意味が違います。

> Prob > chi2 = 0.0028
モデルのあてはまり具合を示す,ピアソンのカイ二乗統計量の上側確率ですね。
Chi2 の値は出力されていませんか?
chi2(数字) = 数字
みたいな書式で出ているのではないかな?かっこ内の数値は自由度。
> bitem | Odds Ratio   Std. Err.   z         P>|z|
> sex   |    .64       .0957243    -2.984    0.003

ここにあるzは独立変数のオッズ比(Odds Ratio)の有意性の検定ですね(注)
Excel で計算してみると, P>|z| = (1-normsdist(abs(-2.984)))*2=0.002845206です。

(注)logistic のときでも,logit のときの Z が出力されているんだって!

     [このページのトップへ]


400. Re^3: p値のみかた。  ともも  2001/10/16 (火) 21:37
stata使っています。

> > Prob > chi2 = 0.0028
> モデルのあてはまり具合を示す,ピアソンのカイ二乗統計量の上側確率ですね。
> Chi2 の値は出力されていませんか?
> chi2(数字) = 数字
> みたいな書式で出ているのではないかな?かっこ内の数値は自由度。

これはLoglikelihood Ratioのchi2と検定で,
ピアソンのカイ二乗統計量はlfitコマンド(goodness of fit)によって表示されるのだと理解していますがどうでしょうか??

     [このページのトップへ]


402. Re^4: p値のみかた。  青木繁伸  2001/10/16 (火) 22:14
> > > Prob > chi2 = 0.0028
> > モデルのあてはまり具合を示す,ピアソンのカイ二乗統計量の上側確率ですね。
> > Chi2 の値は出力されていませんか?
> > chi2(数字) = 数字
> > みたいな書式で出ているのではないかな?かっこ内の数値は自由度。
>
> これはLoglikelihood Ratioのchi2と検定で,
> ピアソンのカイ二乗統計量はlfitコマンド(goodness of fit)によって表示されるのだと理解していますがどうでしょうか??

おっしゃる通り。
383 で
http://www.ats.ucla.edu/stat/stata/output/lognoframe.htm
を見てくださいといっていながら,
384 で,そこを確認しなかったのは小生のミス。

しかし,
「chi2(4) = 13.93」みたいな行があると思うが,そういう出力はないかと逆質問したんですね。

「LR chi2」とあれば,それは log-likelihood ratio を表していると分かるじゃないですか。

# 質問するなら,適当に要約しないで,出力全部見せてほしい X-P

# 私は stata なんか使ったことないんだから,本来は,
# 使っている人が自己責任で解決すべき問題なんです X-P

     [このページのトップへ]


410. Re^5: p値のみかた。  ともも  2001/10/17 (水) 01:38
> # 私は stata なんか使ったことないんだから,本来は,
> # 使っている人が自己責任で解決すべき問題なんです X-P

いやいや,何とか手助けしようという青木さんの姿勢は感心します。

統計分野に限らずありがちですが,もし,

>chi2といわれても何を意味するのか全くわかりませんが。

>つまり何がききたいのでしょうか。

的な回答だったら,無意味に質問者を恐縮させて,自分で調べる気も失せて,デメリットしかないと思うんです。

     [このページのトップへ]


387. Re^3: p値のみかた。(単なる感想)  DISIR  2001/10/12 (金) 11:09
http://www.ats.ucla.edu/stat/stata/faq/oratio.htm
の例題だと
z = (coef.)/(Std. Err)=1.694596/0.9759001=1.736
ってことなんですよね。パラメータの効果の検定を正規近似を使って検定してるんですね。
それはいいとして,coef.とodds ratio のz同じものっていうのは,出力見たときとまどうし,このソフトいいとは思えないな〜

     [このページのトップへ]


388. Re^4: p値のみかた。(単なる感想)  青木繁伸  2001/10/12 (金) 15:45
> http://www.ats.ucla.edu/stat/stata/faq/oratio.htm
> の例題だと
> z = (coef.)/(Std. Err)=1.694596/0.9759001=1.736
> ってことなんですよね。パラメータの効果の検定を正規近似を使って検定してるんですね。
> それはいいとして,coef.とodds ratio のz同じものっていうのは,出力見たときとまどうし,このソフトいいとは思えないな〜

そうなんですよ。私もさんざん悩んで,上の階の先生がstata 持っていたのを思い出して,マニュアル見せてもらって,「あ〜あ」と思ったのです。
5.444444/5.313234=1.024694941ですからね,ずいぶん結果は違います。
相対危険度については,信頼区間で表現するのが通常なので上のようなZ統計量は計算されないのでいいのですが,
 z    P>|z|  [95% Conf. Interval]
1.736 0.082  .8040183    36.86729
のように係数の検定結果とオッズ比の信頼区間を並べられるととまどいますね。
それに,見た目結果が矛盾することもあるのではないでしょうか?

     [このページのトップへ]


390. Re^5: p値のみかた。(単なる感想)教えて下さい。  sb812109  2001/10/13 (土) 07:48
> (DISIRさん)それはいいとして,coef.とodds ratio のz同じものっていうのは,
> (青木先生)相対危険度については,(中略),,,見た目結果が矛盾すること
> もあるのではないでしょうか?

zは,Wald検定のzであり,ただ単にOR==exp(coef),5.44444==exp(1.694596)ではないのでしょうか。

相対危険度(オッズ比?)で表現すれば,p>0.05であれば95%信頼区間に1を含むであろうし,coefで表現すれば,95%信頼区間に0を含むだけであり,矛盾を来たす心配はないと思います。

ORとcoefは全く同じものを違う形で表現しただけですので,誤解を生じることはないと思うのですが,,,。

つまり,何が問題なのかよく理解できないのですが,,,。

     [このページのトップへ]


391. Re^6: p値のみかた。(単なる感想)教えて下さい。  青木繁伸  2001/10/13 (土) 20:23
> zは,Wald検定のzであり,
これは確認してみます

> ただ単にOR==exp(coef),5.44444==exp(1.694596)ではないのでしょうか。
そうですね。確かにソウです。
admit  Odds Ratio Std. Err.
gender 5.444444   5.313234

では std.err.は何のstd.err.なんでしょう?
おしえてください。

     [このページのトップへ]


392. Re^7: p値のみかた。(単なる感想)教えて下さい。  sb812109  2001/10/14 (日) 10:19
> Admit  Odds Ratio Std. Err.
> gender 5.444444   5.313234
> では std.err.は何のstd.err.なんでしょう?
> おしえてください。

何を問題としているのかやっと判りました。
STATAはよく使っているのですが,今まで何の疑問も感じていませんでした。

logistic bitem sex /* 結果はOR表示 */
logistic bitem sex, coef /* 結果はcoef表示 */

両者ともz値は全く同じですので,基本的に回帰係数の検定は,Wald検定を行い,そのまま,OR表示の場合もそれを流用しているのだと思います。ご指摘の通り,全く,不自然です。

     [このページのトップへ]


396. Re^8: p値のみかた。(単なる感想)教えて下さい。  sb812109  2001/10/15 (月) 22:52
http://www.stata.com/support/faqs/stat/2deltameth.html

だそうです。

     [このページのトップへ]


397. Re^9: p値のみかた。(単なる感想)教えて下さい。  青木繁伸  2001/10/16 (火) 11:02
> http://www.stata.com/support/faqs/stat/2deltameth.html
>
> だそうです。

http://aoki2.si.gunma-u.ac.jp/lecture/mb-arc/arc010/004.html#008
と関係することですね。
私も,このことを忘れていました。

     [このページのトップへ]


● 「統計学関連なんでもあり」の過去ログ--- 015 の目次へジャンプ
● 「統計学関連なんでもあり」の目次へジャンプ
● 直前のページへ戻る