﻿rm(list = ls())  # ワークスペースをクリア
library(stats)  # ライブラリの読み込み
# サンプルデータの設定
success <- 5383 ; trials <- 5383 + 5125 ; p0 <- 51/100
# 有意水準
alpha <- 0.05
# 二項検定を実行
result <- binom.test(success, trials, p0, alternative = "two.sided")
# 検定結果の表示
if (result$p.value < alpha) {
  cat("p値", result$p.value, "が有意水準", alpha, 
      "未満であるから帰無仮説を棄却\n")
} else {
  cat("p値", result$p.value, "が有意水準", alpha, 
      "以上であるから帰無仮説を受容\n")
}