R中线性混合模型的功率仿真

Question

我需要进行线性混合模型模拟以获得不同样本量的功效。

我的模型是：

评分 = y

固定效果，x = Ring

随机效应 = 参与者

我试过的代码如下。它只返回“基于 100 次模拟，（0 个警告，100 个错误）” alpha = 0.05, nrow = 2000' ....

谢谢！！

#create a dataframe
library(lmerTest)
library(simr)
library(tidyverse)
Ring = c('Ring', 'NoRing')
#from 1 to 10 (11 is not included).
Ring = rep(Ring, times = 1000)
attractiveness = floor(runif(10, min=1, max=11)) #this creates random numbers
#from 1 to 10 (11 is not included).

participants<-rep(factor(1:100),each=20)
targetID = rep(c(1,2,3,4,5,6,7,8,9,10), each= 2)
targetImage= rep(targetID, times= 100)
Ratings = rep(attractiveness, times = 200)
data<-data.frame(participants, Ring, targetImage, Ratings)
#parameters for the model:

## Intercept and slopes for ring
fixed <- c(3, 0.5)

## Random intercepts for participants 
rand <- 0.5

## residual variance
res <- 2
model <- makeLmer(Ratings ~ Ring + (1|participants), fixef=fixed, 
                  VarCorr=rand, sigma=res, data=data)
sim_treat <- powerSim(model, nsim=100, test = fcompare(Ratings~Ring))
sim_treat