今回は,R プログラミングではなく,統計学の方でダメ出し
例題をやってみる その2 で,
> ここでは、抗生物質を使うとプラセボより感染率は小さくなるはずだ、という期待が込められているので片側検定を行なっていると思われる。
のに
> (ここではカイ二乗検定を行なっている。)
とは,何を考えているのだろうか?そのようなときには,prop.test を使う。
そして,alternarive="less" を指定する。そうすると,
>ここでp=0.02と返ってきているが、論文では0.01と言っている。
などと,寝言を言わなくても,ちゃんと p= 0.01043 になる。
> prop.test(c(32, 53), c(488, 484), alternative="less")
2-sample test for equality of proportions with continuity
correction
data: c(32, 53) out of c(488, 484)
X-squared = 5.3389, df = 1, p-value = 0.01043
alternative hypothesis: less
95 percent confidence interval:
-1.00000000 -0.01212703
sample estimates:
prop 1 prop 2
0.06557377 0.10950413
※コメント投稿者のブログIDはブログ作成者のみに通知されます