第3章 統計理論の基礎

練習問題 3-2

mydata <- read.csv("02_practice_01.csv", header = F)

mu <- 8 # 帰無仮説のもとでの母平均
xbar <- mean(mydata$V1) # 標本平均
sig2 <- var(mydata$V1) # 不偏分散
sig <- sqrt(sig2) # 標準偏差
n <- length(mydata$V1) # サンプルサイズ
se <- sig / sqrt(n) # 標準誤差
  
tscore <- (xbar - mu) / se # t統計量
dof <- n - 1 # 自由度

pval <- 2 * pt(abs(tscore), df = dof, lower.tail = F) # p値
pval
## [1] 0.04911477

補足

t.test()関数を用いれば簡単にt検定を実施できます．

t.test(mydata$V1, mu = 8)
## 
##  One Sample t-test
## 
## data:  mydata$V1
## t = 2.102, df = 19, p-value = 0.04911
## alternative hypothesis: true mean is not equal to 8
## 95 percent confidence interval:
##   8.008586 11.991414
## sample estimates:
## mean of x 
##        10