使用lme4在R中进行功率分析

问题描述

我正在尝试按照Westfall教程来确定参加研究的最低人数(https://jakewestfall.shinyapps.io/crossedpower/)。

采用平衡设计,我试图将他在网站中提供的lme代码应用于我的数据:

blkdmor_fit <- lmer(ia_rt ~ contrast * Clength * Cpredict * Cfrequency + (contrast|participants) + (contrast|Word_id),data=blkdmor)

但是,当我运行模型时(我有两个参与者的试验数据,并且是两个条件的子集),我收到此错误消息:

Errore:观察次数(= 701)

重点是,我肯定知道Word_id中的单词是549。所以问题是,由于某种原因,lme将它们乘以2,因此它们的数量当然大于观察数。有人知道为什么lme会以这种方式表现这种乘法吗?它会乘以Word_id x参与者,还是每个随机效果?无论如何,我不明白原因。

非常感谢任何想要帮助的人!

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)