# 母平均の仮説検定（t検定）

# 母集団データ
pdata <- na.omit(airquality$Wind)

# 標本データ
n <- 10
set.seed(1)
sdata <- sample(pdata, n)

# 標本平均
xbar <- mean(sdata)

# 不偏標準偏差
s <- sd(sdata)

print(sprintf("標本 : n = %d, xbar = %g, s = %g",
              n, xbar, s), quote = FALSE)

# 帰無仮説
mu0 <- 12
print(sprintf("帰無仮説 : μ = %g", mu0), quote = FALSE)

# 右片側検定
t.right <- t.test(sdata, alternative = "greater", mu = mu0)
# 左片側検定
t.left <- t.test(sdata, alternative = "less", mu = mu0)
# 両側検定
t.twosided <- t.test(sdata, alternative = "two.sided", mu = mu0)

print(sprintf("t = %g, 右片側 p = %g, 左片側 p = %g, 両側 p = %g",
              t.right$statistic, t.right$p.value, t.left$p.value, t.twosided$p.value),
      quote = FALSE)