No.15171 適合度検定  【ラヂオ】 2011/08/18(Thu) 22:15

初めて質問致します。
小さな医学統計で,35病変が,どの領域に生じるかを調べました。
領域Aに21例,領域Bに8例,Cに3例,Dに3例となり,この観測値の不均衡を適合度検定で導き出しました(SPSS)。ところが,これのみでは,分布が不均衡であることは証明できても,領域Aにこの病変が生じやすいとまでは結論付けることができません。
この数値をみれば当たり前?と思ったりもしますが,そこは論文の科学的根拠としては薄いということで困っています。
ここからさらに突っ込んで,領域Aにこの病変が生じやすいと統計学的に示すにはどうしたらいいでしょうか?
そもそも35病変で統計を取ることが,統計学的にはおかしいという突っ込みはできれば避けたいです。
さらにこの病変を別の観点から分類したもの(2)があるのですが,これに関しては項目が細分化されている上に観測値が小さいためにさらに難しくなっています。適合度検定はいくつかをカテゴリ化して施行したのですが,
何とかお知恵を貸していただけないでしょうか。

No.15172 Re: 適合度検定  【ラヂオ】 2011/08/18(Thu) 22:17

すみません。ファイルがうまく添付できないみたいなので。
取り敢えず上記内容のみで

No.15173 Re: 適合度検定  【青木繁伸】 2011/08/18(Thu) 22:56

サンプルサイズ35程度のファイルをどうやって添付しようとして失敗したのやら,疑問ですが。

21, 8, 3, 3 が一様かという検定は,
> chisq.test(c(21, 8, 3, 3))

Chi-squared test for given probabilities

data: c(21, 8, 3, 3)
X-squared = 24.7714, df = 3, p-value = 1.724e-05

ということになり,そうではないよという結果になる。

では,

> 領域Aにこの病変が生じやすいとまでは結論付けることができません。

ということなら,理論的に,領域Aにこの病変が生じやすいという理論的確率分布を定義して検定すればよいのでは。

当然,観察されたものが理論分布だと言い張れば,カイ二乗検定は0になり,帰無仮説は棄却される確率なんて0になるわけで,最初から勝負は決まっている。例えば他のカテゴリーの2倍ということで,観察数の割合が 2:1:1:1 ということを検定するなら

> chisq.test(c(21, 8, 3, 3), p=c(2, 1, 1, 1)/5)

Chi-squared test for given probabilities

data: c(21, 8, 3, 3)
X-squared = 8.2143, df = 3, p-value = 0.04178

ということで,ギリギリではあるけど帰無仮説は棄却され,つまり,観察比率は 2:1:1:1 ではないということになる。
このようなことをあなたはどのように評価するか?

No.15174 Re: 適合度検定  【ラヂオ 】 2011/08/19(Fri) 00:53

ありがとうございます.
理論的確率分布をどう定義するか,よく考えてみます.
ありがとうございました.

(ファイルはエラーがでるばかりで添付できなかったのです.原因はよくわかりません)

No.15176 Re: 適合度検定  【青木繁伸】 2011/08/19(Fri) 10:08

> (ファイルはエラーがでるばかりで添付できなかったのです.原因はよくわかりません)

エクセルのファイルを添付しようとしました?
それはできませんね。上の「利用上の注意」に書いてあります。

35例くらいのデータなら,テキストで,書き込めば良いだろうと言うことです。

No.15190 Re: 適合度検定  【ラヂオ】 2011/08/20(Sat) 22:51

「利用上の注意」を未確認でした.申し訳ありません.
素人的な質問で申し訳ないのですが,この領域ABCDの4つに均一に分布する事を仮定して,今回の観測値の偏りを何か視覚的にわかりやすく表示する方法はないでしょうか?
確率で領域A 60%,領域B 22.86%,領域cおよびD 8.57%と表示しても,あまりインパクトがなくて.....
当たり前ですが,オッズ比は算出出来ませんし...

見当はずれの質問だったら申し訳ありません

No.15191 Re: 適合度検定 追記  【ラヂオ】 2011/08/20(Sat) 22:52

いくら探しても,適合度検定で偏りが存在すると判明した後の事について触れているものがなくて....
なにかご教授いただけましたら,ありがたいです.

No.15192 Re: 適合度検定  【ラヂオ】 2011/08/20(Sat) 23:04

細かくわけて申し訳ありません.
母集団比率の区間推定として,95%信頼区間を出すというのはおかしいでしょうか

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