我正在尝试R为调查实验做一些功率计算。
我知道如何使用pwr命令来做到这一点。例如，以下命令显示我检测效果大小为 0.2 的最佳样本量是每组有 394 名受访者。
p.out <- power.t.test(
# sample size (leave blank!)
n = ,
# minimum detectable effect size
delta = 0.2,
# alpha and power (1-kappa)
sig.level = 0.05, power = 0.8,
alternative = "two.sided"
)
plot(p.out)
我想说明我有一个两阶段整群抽样策略的事实。例如，选择一组 30 个教室，然后在每个教室中随机抽取 10 名学生进行访谈。如何在功效分析中考虑这种抽样差异？
请注意，这不是类似于集群随机试验 (CRT)的场景。鉴于调查实验的性质，治疗不是在小组层面分配，而是在单位层面分配。换句话说，每个班级都会随机分配经过治疗和控制的学生。
在这种情况下，是否有一种简单直观的方法来运行功率分析？也许指定某种形式的类内相关？
非常感谢你的帮助