Question

我很难确定如何对混合效果模型进行功效分析。我想要的是一系列模拟，这些模拟显示了在不同样本大小（例如n = 50、100、200、300、500等）下设定效果大小的功效，最好以图形或表格形式显示。

我一直在关注一些在线资源，以了解有关功率仿真的信息。例如，我一直在遵循指南here，here，here和here。有些示例有效，有些则无效。更重要的是，我无法通过自己的学习设计使这些模拟用于我的目的。

具体来说，我想测试双向交互的能力，但是在此我简化了分析，仅出于此问题的目的而仅考虑单个预测变量。

我已经尝试了以下操作，使用了我拥有的一些试验数据的效果大小。

首先，一些数据：

set.seed(1234)
smaller.df <- data.frame(
  subj <- rep(1:106, 2),
  x = sample(1:6, 212, replace = T),
  m = sample(1:3, 212, replace = T),
  y = rnorm(212)
)

接下来，我可以拟合模型并将我观察到的效果插入我的试验数据中：

model.of.interest2 <- lmer(y ~ x + (1|subj), data=smaller.df)
fixef(model.of.interest2)['x'] <- .45

首先，结果只是告诉我我的功率为0或100％...我不确定为什么。我在网上找到的所有资源都没有清楚地说明我在这里应该做什么...

powerSim(model.of.interest2, nsim=100, test = fcompare(y~x))

第二，我无法使用以下代码生成（我认为应该是）不同样本量的功效图：

pc <- powerCurve(model.of.interest2 , along = "subj", nsim =  100)
plot(pc)

它仍然只显示我具有0，然后是100％的功率。

我还尝试使用“扩展”命令，我认为这会创建不同的样本量。例如，

model500 = extend(model.of.interest2, along = "subj", n = 500)

但是它只会返回一个观测值为0的模型。

我真的不确定自己在做错什么，甚至应该做些什么。

不同样本量下混合效应模型的功效仿真

0 个答案: