## ==== パラメータ設定 ====
set.seed(123)
p <- 0.4 # 母比率
n <- 6 # 標本サイズ
B <- 10000 # シミュレーション回数
## ==== シミュレーション ====
x <- rbinom(B, size = n, prob = p)
## ==== 集計（有病者数・標本比率・度数）====
k_vals <- 0:n
freq <- tabulate(x + 1, nbins = n + 1)
res <- data.frame(有病者数 = k_vals, 標本比率 = round(k_vals / n, 2), 度数 = as.integer(freq),
row.names = NULL, check.names = FALSE
)
res 有病者数 標本比率 度数
1 0 0.00 487
2 1 0.17 1835
3 2 0.33 3169
4 3 0.50 2766
5 4 0.67 1369
6 5 0.83 336
7 6 1.00 38