超基本的な分析方法のメモ.
2群の平均値の差の検定を行う. 対応のないt検定は以下の通り.
次は対応のあるt検定.
2群の平均値の差の検定を行う. 対応のないt検定は以下の通り.
> #以下のAとBのデータに2群の検定を行う.
> dat
A B
1 23.5 24.0
2 23.0 23.0
3 25.0 25.5
4 25.5 27.0
5 22.5 24.5
6 24.5 26.0
7 26.0 27.5
8 24.0 26.0
9 25.0 26.5
10 24.0 25.0
>
>
> #以下でt.test()関数を使って,対応のないt検定を実施する.
> #dat$Aの意味は,datのデータの中のA列のデータということ.paried = で対応の有無を指定.今回は無しなのでFALSE
> t.test(dat$A,dat$B, paired = FALSE)
Welch Two Sample t-test
data: dat$A and dat$B
t = -2.1287, df = 17.142, p-value = 0.04807
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-2.38859256 -0.01140744
sample estimates:
mean of x mean of y
24.3 25.5
> #tはt値で,dfは自由度.p-valueがp値.p<0.05で有意差ありとすることが多い.
>
次は対応のあるt検定.
> #以下のAとBのデータに2群の検定を行う.
> dat
A B
1 23.5 24.0
2 23.0 23.0
3 25.0 25.5
4 25.5 27.0
5 22.5 24.5
6 24.5 26.0
7 26.0 27.5
8 24.0 26.0
9 25.0 26.5
10 24.0 25.0
>
> #以下でt.test()関数を使って,対応のあるt検定を実施する.
> #dat$Aの意味は,datのデータの中のA列のデータということ.paried = で対応の有無を指定.今回は有りなのでTRUE
> t.test(dat$A,dat$B, paired = TRUE)
Paired t-test
data: dat$A and dat$B
t = -5.6223, df = 9, p-value = 0.0003249
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-1.6828291 -0.7171709
sample estimates:
mean of the differences
-1.2
> #tはt値で,dfは自由度.p-valueがp値.p<0.05で有意差ありとすることが多い.
>
コメント