1

我有以下模型：

ModelPower <- lmer(DV ~ GroupAbstract * Condition_Cat_Abs + (1|Participant) + (1 + GroupAbstract|Stimulus), data = Dataset)

该模型给出以下输出：

Random effects:
 Groups      Name                  Variance Std.Dev. Corr 
 Participant (Intercept)           377.401  19.427        
 Stimulus    (Intercept)            91.902   9.587        
             GroupAbstractOutgroup   2.003   1.415   -0.40
 Residual                          338.927  18.410        
Number of obs: 16512, groups:  Participant, 344; Stimulus, 32

Fixed effects:
                                                        Estimate Std. Error         df t value             Pr(>|t|)    
(Intercept)                                              65.8962     2.0239    59.6906  32.559 < 0.0000000000000002 ***
GroupAbstractOutgroup                                    -0.9287     0.5561   129.9242  -1.670               0.0973 .  
Condition_Cat_AbsSecondOrderIn                           -2.2584     0.4963 16103.9277  -4.550           0.00000539 ***
Condition_Cat_AbsSecondOrderOut                          -7.0821     0.4963 16103.9277 -14.270 < 0.0000000000000002 ***
GroupAbstractOutgroup:Condition_Cat_AbsSecondOrderIn     -3.0229     0.7019 16103.9277  -4.307           0.00001665 ***
GroupAbstractOutgroup:Condition_Cat_AbsSecondOrderOut     7.8765     0.7019 16103.9277  11.222 < 0.0000000000000002 ***

我对交互“GroupAbstractOutgroup:Condition_Cat_AbsSecondOrderIn”很感兴趣，我正在尝试使用 R 包 simr 估计样本量以检测至少 -2 的效应量。原来的斜率是-3.02，所以我指定了新的：

ModelPower@beta[names(fixef(ModelPower)) %in% "GroupAbstractOutgroup:Condition_Cat_AbsSecondOrderIn"] <- -2

但是，无论我如何为主要效果和交互指定 powerSim 函数（参见下面的一些示例），在运行 lastResult()$errors 'object is not a matrix' 时，我都会得到 0% 的功率和以下错误。我知道错误应该是什么意思，但即使将原始数据框和固定效果表转换为矩阵后，错误仍然存​​在，我不确定它指的是什么以及如何获得实际输出。任何帮助将非常感激！

powerSim 函数示例：

powerSim(ModelPower, test=fixed("GroupAbstract", "anova"), nsim=10, seed=1)
powerSim(ModelPower, test=fixed("GroupAbstractOutgroup:Condition_Cat_AbsSecondOrderIn", "anova"), nsim=10, seed=1)
4

0 回答 0