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

★ 分割表で多重比較は可能ですか？ ★

3083.　分割表で多重比較は可能ですか？　bonbon　2004/05/07 (金) 10:35
├3086.　Re: 分割表で多重比較は可能ですか？　勝　2004/05/07 (金) 13:36
│└3087.　Re^2: 分割表で多重比較は可能ですか？　bonbon　2004/05/07 (金) 16:39
│　└3089.　Re^3: 分割表で多重比較は可能ですか？　青木繁伸　2004/05/07 (金) 18:29
│　　└3103.　Re^4: 分割表で多重比較は可能ですか？　ぼんぼn　2004/05/10 (月) 08:13
│　　　└3104.　Re^5: 分割表で多重比較は可能ですか？　青木繁伸　2004/05/10 (月) 09:22
└3085.　Re: 分割表で多重比較は可能ですか？　青木繁伸　2004/05/07 (金) 10:58

3083.　分割表で多重比較は可能ですか？　bonbon　　2004/05/07 (金) 10:35

初めまして。よろしくお願いいたします。

データは患者一人に対して反応がプラスかマイナスのどちらかで計数して集計しました。
下記に集計結果を分割表で示します。

2000年(A) 2001年(B) 2002年(C) 2003年(D) プラス 122 106 95 64 マイナス 332 226 152 116
比較したいのは経時的に，A vs B，B vs CおよびC vs Dでの比較を行ってプラスの反応に有意差があるかどうかです。
各群間に対応はありません。

分割表から検討するに当たり，どういった多重比較検定法が適当かヒントをいただけると助かります。
単にt検定を複数回行うと検出率が悪くなるので不適当と判断しました。
分散分析ですと，ScheffeのF testやBonferroniの多重比較検定のどのプログラムを使うのがベストなのか？
ノンパラメトリック検定ではKruskal-Wallis検定の適用を考えておりますが，データが分割表の場合はχ2乗しか思い浮かびません。

集計前のデータは例として・・・
患者 2000年(A) 00001 ＋ 00002 ー 00003 ＋ 00004 ＋・・
といった並び方となります。
お手数ですがよろしくお願いいたします。

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

3086.　Re: 分割表で多重比較は可能ですか？　勝　　2004/05/07 (金) 13:36

> 比較したいのは経時的に，A vs B，B vs CおよびC vs Dでの比較を行ってプラスの反応に有意差があるかどうかです。

　Bonferroni式などで有意水準を調整して，比較したいペアのカイ2乗検定を繰り返し行う方法もあります。
　他に方法がないときや比較したいペア数が少ないときにはこの方法が有効であると思います。

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

3087.　Re^2: 分割表で多重比較は可能ですか？　bonbon　　2004/05/07 (金) 16:39

＜＜青木先生へ＞＞
早速アドバイスありがとうございます。

＞この中に，同じ患者は重複して数えられていないのですね。
そうです。重複はしておりません。
ちなみに一部重複している場合はどのように考えたらよろしいでしょうか？
重複している場合でも【対応している】と考えてはおりませんがいかがでしょうか？
教えて頂いたアドレスより対比較の検討を行いました。
独立性の検定を chi-sq-test-s.xlsにて算出したところ，p値は0.0103でした。このあとどの式を適用すればよろしいでしょうか？

また，傾向検定でコクラン・アーミテージ検定では独立変数が定義できないので検定に不向きなデータではないでしょうか？

お手数をお掛けしますが，よろしくお願いいたします。

＜＜勝様＞＞
早々にアドバイスをありがとうございます。
> Bonferroni式などで有意水準を調整して，比較したいペアのカイ2乗検定を繰り返し行う方法もあります。
早速試させていただきました。
上記の3種の組み合わせで比較をしたいのでそれぞれのP値は0.05/3＝0.017と算出されました。
A vs B，B vs CおよびC vs Dのカイ2乗検定を行ったところ，それぞれ
A vs BのP値は0.1229，
B vs CのP値は0.1023
C vs DのP値は0.5396
と算出され，上記の組み合わせでは有意差は無いと判断されました。

以上，途中経過ですが取急ぎ報告まで。

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

3089.　Re^3: 分割表で多重比較は可能ですか？　青木繁伸　　2004/05/07 (金) 18:29

> ちなみに一部重複している場合はどのように考えたらよろしいでしょうか？

対象としている事象が何なのかにもよるでしょうが，例えば病院の外来患者である病気にかかっているかどうかという場合なら，一番最初の来院時のデータだけを使うと言うことでしょう。
重複が少ない場合には無視してもよいでしょう。

> 独立性の検定を chi-sq-test-s.xlsにて算出したところ，
> p値は0.0103でした。このあとどの式を適用すればよろしいでしょうか？

ライアンの方法なら検定手順の3番以降，チューキーの方法なら検定手順の4番以降でしょうが，手計算をやるのは面倒でしょうから，各ページの最後に R によるプログラムへのリンクを用意しております。

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

3103.　Re^4: 分割表で多重比較は可能ですか？　ぼんぼn　　2004/05/10 (月) 08:13

青木先生へ
お忙しいところご返事ありがとうございます。

>各ページの最後に R によるプログラムへのリンクを用意しております。
上記のページに書かれているソースをどのように使用させていただければよろしいでしょうか？
自分でソースを元にプログラムを組み上げる必要があるのでしょうか？
excel組み上げる事が可能であれば，ポイントをご教授いただけると助かりますがいかがでしょうか？

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

3104.　Re^5: 分割表で多重比較は可能ですか？　青木繁伸　　2004/05/10 (月) 09:22

R を使うにはどうしたらいいかは，下記ページの一番上のリンクをクリック。
http://aoki2.si.gunma-u.ac.jp/R/index.html

Excel でもなんでも，プログラム言語でプログラムを組むのは可能です。
アルゴリズムは，先に示したページに書いてあるとおりです。

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

3085.　Re: 分割表で多重比較は可能ですか？　青木繁伸　　2004/05/07 (金) 10:58

この中に，同じ患者は重複して数えられていないのですね。

> 比較したいのは経時的に，A vs B，B vs CおよびC vs Dでの比較を行ってプラスの反応に有意差があるかどうかです。

以下を参照してみてください。ただし，対比較ですから全部の組み合わせを比較します。
http://aoki2.si.gunma-u.ac.jp/lecture/Hiritu/pmul.html

傾向検定なら，
http://aoki2.si.gunma-u.ac.jp/lecture/Hiritu/Armitage.html

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