我试图运行 LMM 允许时间成为随机效果:

lmer(efficacy ~ time + (1+time|id), data=long, REML=TRUE)

但是,我收到以下错误消息:

错误:观察次数(= 3948)<=术语随机效应数(= 4496)(1次| id);随机效应参数和残差方差(或尺度参数)可能无法识别

我有三个时间点,但并非所有受试者都拥有所有3个时间点(磨损)的数据,但我对结果变量(效力)的遗漏计数只有13个缺失值,而且我的整个数据集中的总数为 3961 participants . 如果有人能够对此有所了解或指出我正确的方向,我会非常感激 .

谢谢 .