My Notes

統計学とかR(R言語)とかPython3の覚え書きとか走り書きとか。 座右の銘にしたい: All work and no play makes Jack a dull boy.

一標本t検定と一標本t検定の信頼区間の計算をR(R言語)で

Rコード

# 一標本t検定の計算
# t = (xbar - μ0) / (s / sqrt(n))
# 具体例
# 標本平均10, 標本標準偏差10, 標本数15
t <- (90 - 100) / (10 / sqrt(15))
t
# [1] -3.872983
# 自由度15 - 1 = 14 α = 0.05 上方棄却値2.145
# |-3.87|>2.145 で帰無仮説は棄却


# 一標本t検定の信頼区間を求める式
信頼区間の上界 <- 90 + (2.145)*(10 / sqrt(15))
信頼区間の上界
# 95.53837


信頼区間の下界 <- 90 - (2.145)*(10 / sqrt(15))
信頼区間の下界
# 84.46163

R Console

> # 一標本t検定の計算
> # t = (xbar - μ0) / (s / sqrt(n))
> # 具体例
> # 標本平均10, 標本標準偏差10, 標本数15
> t <- (90 - 100) / (10 / sqrt(15))
> t
[1] -3.872983
> # [1] -3.872983
> # 自由度15 - 1 = 14 α = 0.05 上方棄却値2.145
> # |-3.87|>2.145 で帰無仮説は棄却
> 
> 
> # 一標本t検定の信頼区間を求める式
> 信頼区間の上界 <- 90 + (2.145)*(10 / sqrt(15))
> 信頼区間の上界
[1] 95.53837
> # 95.53837
> 
> 
> 信頼区間の下界 <- 90 - (2.145)*(10 / sqrt(15))
> 信頼区間の下界
[1] 84.46163
> # 84.46163

参考文献

統計クイックリファレンス 第2版

統計クイックリファレンス 第2版