错误边界（单数）拟合：参见 lme4 Model 中的 isSingular

Question

我有以下结构：

'data.frame':   1041 obs. of  4 variables:
 $ RST_AST: num  0.45 0.576 0.45 0.45 0.675 ...
 $ subj   : chr  "4" "13" "24" "25" ...
 $ domain : chr  "perceptual" "perceptual" "perceptual" "perceptual" ...
 $ ce     : chr  "compromise" "compromise" "compromise" "compromise" ...

现在我想使用 lme4 package 运行多级 Model ，如下所示：

library(lme4) #load package

ce_model <- lmer(RST_AST ~ domain + ce+ (1|subj), data = df_RST_AST)

然后我收到错误消息：

boundary (singular) fit: see ?isSingular

我不知道出了什么问题？

有人能帮我吗？

Answer 1

您的变量之一很可能没有被评估为对方差有贡献。 最有可能的是，“罪魁祸首”是随机变量（在您的情况下是 subj）。 如果您查看 model 的 output（只需打印 ce_model 并按回车键或摘要（ce_model）），应该有一个变量的方差和标准偏差为零或几乎为零。 有几种方法可以解决这个问题（请参阅上面的任何链接）。 我建议分析错误变量，它有什么样的分布，可以转换，...等等。 也许您需要对数据进行一个或多个通用模型。 查看您的 model 或型号也可能会有所帮助； 但是，如果没有进一步的信息，很难说。