library(readr)

# データの読み込み
df <- read_csv("yamagata.csv", locale = locale(encoding = "SHIFT-JIS"))

df$`/項目` <- NULL
df$調査年 <- as.numeric(sapply(df$調査年, gsub, pattern="年度", replacement = ""
))
df <- df[order(df$調査年),]

# 2015年度の人口増減率を標本とする
sdata <- df[df$調査年 == 2015,][["#A05101_人口増減率【％】"]]
n <- length(sdata)

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

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

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

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

# t検定（母分散が未知）
# 右片側検定
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)