10人の授業前後のテスト成績の exam1.csv というのを作った。
no,pre,post
1,50,51
2,42,54
3,57,53
4,52,58
5,61,64
6,42,60
7,45,65
8,52,57
9,53,58
10,45,62
1,50,51
2,42,54
3,57,53
4,52,58
5,61,64
6,42,60
7,45,65
8,52,57
9,53,58
10,45,62
こんなの。
で、
(wd <- getwd())
library(psych)
exam <- read.csv(paste0(wd,"/exam1.csv"))
with(exam, t.test(pre, post, paired = T))
ってやると
Paired t-test
data: pre and post
t = -3.2692, df = 9, p-value = 0.009694
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-14.043187 -2.556813
sample estimates:
mean of the differences
-8.3
となる。
> describe(exam[c("pre","post")])
とやると
vars n mean sd median trimmed mad min max range skew kurtosis se
pre 1 10 49.9 6.37 51 49.50 8.90 42 61 19 0.21 -1.37 2.01
post 2 10 58.2 4.66 58 58.25 5.93 51 65 14 -0.02 -1.47 1.47
となる。
