如何在power.t.test中设置增量范围

# 导入所需的库和数据
library(pwr)

# 定义两组样本的均值和标准差
mean1 <- 10
mean2 <- 12
sd1 <- 2
sd2 <- 2

# 设置显著性水平和统计功效
alpha <- 0.05
power <- 0.8

# 设置增量范围
delta_range <- seq(0.5, 2, by = 0.5)

# 遍历增量范围，进行假设检验
for (delta in delta_range) {
  result <- power.t.test(n = NULL, delta = delta, sd = sqrt(sd1^2 + sd2^2), 
                         sig.level = alpha, power = power, type = "two.sample", 
                         alternative = "two.sided")
  
  # 输出结果和相关信息
  cat("增量范围:", delta, "\n")
  cat("所需的样本大小:", ceiling(result$n), "\n")
  cat("实际统计功效:", result$power, "\n")
  cat("\n")
}