# 在 R 中,可以使用 pnorm 函数得到正态分布的分位数
# 先定义参数
alpha <- 0.05 # α
beta <- 0.1 # β
delta <- 0.08 # 率差
sd <- 0.14 # 标准差
d <- delta / sd # 效应量
# 找到正态分布的分位数
Z_alpha <- qnorm(1 - alpha / 2) # α的分位数
Z_beta <- qnorm(1 - beta) # β的分位数
# 计算样本量
n <- (Z_alpha + Z_beta)^2 * 2 * sd^2 / delta^2
n <- round(n) # 对样本量进行四舍五入
print(n)
基于上述条件和计算,每组的样本量应为64。这意味着你需要总共128个病人(64个在实验组,64个在对照组)来达到预设的统计效力(即,使得第一类错误的概率小于5%,第二类错误的概率小于10%)。