ユーザ用ツール

サイト用ツール


rでt検定

文書の過去の版を表示しています。


Rでt検定

独立した2群間の連続変数の平均値を比較する方法

1.t検定

2.Welch検定

3.Mann-Whitney U検定

2群間に対応がある場合(インスリン導入前後のHbA1cの値など)は、「対応のあるt検定」へ。

1.どのような場合にt検定を用いることができるのか?

2群が独立しており、

(1)各群のデータが正規分布に従う ⇒ 両群の分散が等しければ、t検定。分散が等しくなければ、Welch検定を用いる。

(2)各群のデータが正規分布に従わない、あるいは順序変数 ⇒ Mann-Whitney U検定

2.各群が正規分布に従うかどうかの判定”Kolmogorov-Smirnov検定(ks.test()関数)”

3.正規分布に従う2群の分散が等しいかどうかの判定”F検定(var.test()関数)”

F検定の後にt検定することについては、議論があるらしいが、気にしないことにする。

標本の分散が帰無仮説に従うかどうかを検定

F 検定

2018.01.12

https://stats.biopapyrus.jp/stats/f-test.html

p>0.05であれば、2群は分散が等しいと言えるので、2群の平均値を比較するときはt検定を用いる。

p>0.05であれば、2群は分散が等しいとは言えないので、2群の平均値を比較するときはWelch検定を用いることになる。

4.t検定"t.test(group1, group2, paired=F)"

p<0.05であれば、2群の平均値に有意差あり。

p>0.05であれば、2群の平均値に有意差無し。

5.Welch検定"t.test(group1, group2, var.equal=F, paired=F)"

t.test()関数の中で、var.equal=F と入れればWelch検定ができる。

p<0.05であれば、2群の平均値に有意差あり。

p>0.05であれば、2群の平均値に有意差無し。

6.Mann-Whitney U検定"wilcox.test()関数"???

http://cse.naro.affrc.go.jp/takezawa/r-tips/r/65.html

上記サイトでは、ウィルコクソンの順位和検定と、マン・ホイットニーの U 検定が同じであると記載されているが、本当かどうかは未確認。

memo

ks.test()関数を用いるときは、下記コードでは、y1[“sBP_00pre”]ではエラーになる。正しくは、y1$sBP_00preのように記載する。データフレームとリストの違い?

y1 <- subset(x2, Hb_Improve=="TRUE")
y1 <- y1$sBP_00pre
ks.test(x=y1,y="pnorm",mean=mean(y1),sd=sd(y1))

参考リンク

参考文献

EZRで行う場合は、下の本のp92以降が非常に分かりやすいです。

rでt検定.1517665680.txt.gz · 最終更新: 2018/10/07 (外部編集)

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki