私の関わりのある法人
ksbutton.png omemebuttan.png sowerbuttan.png
※このブログに書いていることは、私の関わりある法人の意見ではなく、
 あくまでも、私個人の意見です。

2023年09月20日

R で対応あるt検定



 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

こんなの。

で、

(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

となる。
posted by kingstone at 03:33| Comment(0) | TrackBack(0) | 統計 | このブログの読者になる | 更新情報をチェックする
この記事へのコメント
コメントを書く
お名前:

メールアドレス:

ホームページアドレス:

コメント:

※ブログオーナーが承認したコメントのみ表示されます。

この記事へのトラックバック