理学療法士がまとめたノート

統計学備忘録(R言語のメモ)

since2016 ときどきTEXのメモ

t分布からの信頼区間

投稿: 2017-06-12
更新: 2022-1-16

例)p=0.025、自由度5の場合

# 分位点関数 quantile
(q <- qt(0.025, 5, lower.tail = F))
# 確率関数 probability
(p <- pt(q, 5, lower.tail = F))
# 密度関数 density
(d <- dt(q, 5))

サンプルサイズ6の95%信頼区間を求める

#データサンプル
DBP<-c(26,19,28,24,29,32)
vDBP<-var(DBP)

統計量T値の求め方

(t <- mean(DBP) /sqrt(vDBP/6))

95%信頼区間の求め方

mean(DBP)-qt(0.025,5,lower.tail = F)*sqrt(vDBP/6)
mean(DBP)+qt(0.025,5,lower.tail = F)*sqrt(vDBP/6)

Rの関数で確認

t.test(DBP)

One Sample t-test
data:  DBP
t = 14.328, df = 5, p-value = 2.985e-05
alternative hypothesis: true mean is not equal to 0
95 percent confidence interval:
 21.60893 31.05774
sample estimates:
mean of x 
 26.33333 

t分布のグラフ

x <- seq(-4, 4, 0.1)
plot(x, dt(x, 20), type="l")

for(i in c(2,3,4)){
for(i2 in 2:4){
curve(dt(x, 6-i), col=i, type="l",add=T, )
}
}
labels<-c("自由度","黒:20","赤:4","緑:3","青:2")
legend("topleft", legend = labels)

f:id:yoshida931:20220116223222p:plain:w500