[統計]統計学でいう、検定の繰り返しと多重比較についての質問です。
2014/11/120:47:17
[統計]統計学でいう、検定の繰り返しと多重比較についての質問です。
例えば学校のクラスがA〜Bの3つあって、国語の平均点を比較して「クラスによって国語の成績が異なるか」を調べるため、分散分析を行ったとします。
クラスの効果が有意であったとしても、A⇔B、B⇔C、C⇔Aのどの組み合せ間に有意差があるのかは分からないわけですが、ここで全ての組み合わせについてt検定を繰り返すと第一種の過誤を犯す危険が増えるので、Tukey法などの多重比較を行うようにと教科書には書かれています。
「A=B=C」という帰無仮説を棄却するということは、「A=B」「B=C」「C=A」という3つの帰無仮説のうち少なくとも1つを棄却するということに等しく、3つのうち少なくとも1つが棄却される確率は、「A=B」という1つの帰無仮説が棄却される確率よりは高いだろうから、有意水準の調整が必要だという話はわかります。
しかし疑問があって、それはあくまで「A=B=C」という全体の帰無仮説との整合性にこだわっているから生じる問題であって、「全部同じかどうか」に興味はなく、単にA⇔B間に差があるかどうかを知りたいのだという人がいた場合は、普通にA⇔Bのt検定を行えば良いような気がします。
つまり、
・全体仮説と部分仮説があり、
・部分仮説が「少なくとも1つ」棄却されると全体仮説も自動的に棄却されるという関係があり、
・関心があるのはあくまで全体仮説のほうである
という場合にのみ、部分仮説の検定の繰り返しが問題になるのであって、全体仮説にそもそも関心がなければ繰り返しの問題は気にしなくていい……と思って、そういう内容をブログに書きました。
http://statsbeginner.hatenablog.com/entry/2014/11/01/140721
しかし、よく考えるとおかしいということに気づきました。
いまA群とB群のデータだけがあって、t検定したらA⇔B間に有意差がありましたというときに、後からC群のデータが追加された場合、有無を言わさず「A=B=C」を棄却したいところですが、多重比較のロジックだと、さっき棄却したばかりの「A=B」が今度は棄却されない可能性も出てきますよね?
これはどういうことなんでしょうか。
新たにC群が見つかって情報が増え、全体の誤差が増えたので、「A=B」の検定についても見なおさなければならないと理解すればいいのでしょうか?
ベストアンサーに選ばれた回答
2014/11/122:44:05
なんだか,揚げ足取りみたいで,申し訳ありませんが,ブログで
>伝統的にまず分散分析→有意だったら多重比較
というのは,伝統ではなく,きちんと方法を理解してない人の説明です。
もしかして,ご存知かもしれませんが,滋賀医大・小山氏の解説に,参考になる点があります。
http://www.shiga-med.ac.jp/~koyama/stat/com-ph.html
特に注意が必要なのは,小山氏も書いていますが,ANOVAと多重検定は、そもそも検定の目的が異なり,前者で有意でなくても,後者では有意となる場合があることです。
この点は,統計ソフト Prism の解説も参照
http://graphpad.com/guides/prism/6/statistics/index.htm?stat_result...
さらに,ブログでの誤解
>3群以上になると分散分析
2群以上,です。ここが,しばしば誤解される点です。
もちろん,2群で分散分析するとt検定の結果と同じになります。Rなどでやってみると分かります。だから,研究者によっては,t検定は不要,と言う人もいます。
さらに,検定法として問題があるのは
>2つのクラスの担任の先生の教育能力は等しい
に対して
>有意水準の調整が必要だ
これは多重検定ではなく,2要因混合モデルの分散分析の対象です。
要因A:クラス2水準
要因B:各生徒に対して,教科4水準
>2つのクラスの担任の先生の教育能力は等しい
を棄却するのに,多重検定は不要です。
質問に関して
>後からC群のデータが追加された場合
>有無を言わさず「A=B=C」を棄却したい
なぜ?
最初から,AB間だけ興味があれば,そこだけ検定すれば良いことになります。
これに関連して,別の質問での私の回答も参照
http://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q1011053757...
例えば,最初から,AB間だけに関心があれば,Cが加わったとしても,ABだけ検定すれば良い。
これを多重検定の点から言えば,Bonferroni法で1ヶ所だけ比較するので
t検定のp値 / 1 = t検定のp値そのもの
と考えれば良いと分かります。
以下の論文も参照
小野秀樹(2003)
多重比較 : 薬理データの統計解析における Bonferroni 補正法のメリット・デメリット
ト
-
2014/11/123:08:07
それから,質問では,t検定と分散分析の関係にのみ注目していますが,例えば,A,B,C の3クラスで,全平均(学年平均)と比べて,有意に異なるかどうか,という検定もできます。カイ二乗検定で残差分析をすることに似た方法です。
これを,Analysis of Means (ANOM) と言います。例えば,次のサイト参照
http://www.wludyka.com/anom.html
返信を取り消しますが
よろしいですか?
- 取り消す
- キャンセル
あわせて知りたい
- エクセルで、年齢が平均になるようにグループ分けする方法はありますか? 例えば...
- t検定のボンフェローニの補正について教えてください。 あるグループ、A、B...
- 統計で,t分布とq分布の違いが判りません.分母の√内で誤差分散か不偏誤差分散...
- セキュリティソフトに検知されないウイルスに感染してしまいました。 パソコンの...
- 二要因分散分析 エクセルでの計算 二要因分散分析をエクセルで計算しようと思った...
- スチレンが光(紫外線)に弱い理由を教えて下さい!ラジカル安定性が関係あるとか聞...
- 佐久間象山の子孫はいるのですか?
- スマートフォンでスクリーンショットを撮った場合、その写真に位置情報は含まれて...
- ネガフィルムとポジフィルムの違いをおしえてください。
- 「御手元」と「お手許」の違い、用法を教えて下さい。 よろしくお願い致します。
- 若乃花の父親は初代若乃花なんですか?
- マグニチュード7は震度で言うと何ですか?
- エクセルでAなら1,Bなら2、Cなら3となる関数を教えてください。
- サイズW×D×Hとはそれぞれどこの部分を指すのですか? ちなみに額縁です。
このカテゴリの回答受付中の質問
- 簿記の講義で総平均法の問題が出たので、解答と求め方をやさしく教えてください。...
- 数学のルールブックみたいなもので有名なものを教えて欲しいです!
- クイズの答えを教えて下さい。 叔父に出されましたが答えを教えてくれないまま… ...
- 横国理工志望です 青チャートのレベル4の例題はセンター終了後からやっても大丈夫...
- メロペネム0.25gの250mg+生食10mlという指示の場合生食を何mlで溶解すればよいの...
- 1.冷却水設備について質問です。 図の様な、水冷式圧縮設備があります。(簡...
- 点滴の滴下数の計算と実際の調節方法について教えて下さい。(まだ学校では習って...
- 【工学系】どなたかこちらの問題解いてください!!!お願いします!
- 射影と鏡像について、教えてください。例などを使ってもらえると助かります。
- r(t)=(a cost/1+sin^2t, a costsint/1+sin^2t) これが正則曲線であることを証明...
このカテゴリの投票受付中の質問
- 一次関数の変化の割合を求める簡単な方法わかる方がいましたら、教えてください。
- なぜ2進法表示の1000-1が111になるんですか?
- x^2-2xy+3y^2-2x-8y+13を満たす整数x、yの組を求めてください。 チップ25です
- これはどう解くのでしょうか?確率です。
- x乗+y2乗=25 x2乗−4=y2乗−9 の連立方程式の解を教えて下さい!!
- 高一ながら良問作っちゃったんだよなぁ 暇人は解いてみて 解法待ってるゾ
- 数aの問題です(><)教えてください
- 大至急です。立志舎高校の土曜コース単位認定試験について分かる方回答お願い致し...
- 数aの問題です(><)教えてください
- この問題の解説をお願いします!
専門家が解決した質問
初めてのいろいろ、心配がいっぱい
回答受付中の質問
カテゴリQ&Aランキング
- 戻る
- 次へ
総合Q&Aランキング
あくまで「貸した」だけなんですけど
みんなのアンテナ
Yahoo!知恵袋カテゴリ
お客様自身の責任と判断で、ご利用ください。