只适用于1组的线性混合模型

Question

是否可以强制具有随机效果的lmer模型只适用于只有一个级别的数据？我们希望在极少数情况下保持相同的模型结构，我们的数据只包含1个分组级别。以下说明错误。

library(lme4)
#> Loading required package: Matrix

sleepstudy$Subject <- as.character(sleepstudy$Subject)
ss <- sleepstudy[sleepstudy$Subject == "308", ]

m1 <- lmer(Reaction~Days+(1|Subject), ss)
#> Error: grouping factors must have > 1 sampled level

需要注意的是，我们正在修正方差（参见上一个问题：Fixing variance values in lme4）。因此，我们不需要估计方差。