MEEM(object, conLin, control$niterEM) 中的错误：在第 0 级，块 1 处反向求解中的奇点...线性混合模型，nlme，R答案

【问题标题】：Error in MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1... Linear mixed model, nlme, RMEEM(object, conLin, control$niterEM) 中的错误：在第 0 级，块 1 处反向求解中的奇点...线性混合模型，nlme，R
【发布时间】：2021-09-02 19:50:01
【问题描述】：

m1 =lme(fixed=Hour~Age*Ethnicity*GHUsedFlag*Disability,
   random=~Age|Id, data=mydf)

在我的模型中，我正在检查多级交互。年龄数字重复测量、种族特征因素、GHUsedFlaglogical、残疾特征因素、小时数字、结果变量。在模型中引入 Disability 时出现此错误。无法真正理解为什么会发生错误以及如何应对它。任何帮助表示赞赏。我在 R 中使用 nlme 包。

MEEM(object, conLin, control$niterEM) 中的错误：在第 0 级，第 1 块的反向求解中的奇点

【问题讨论】：

标签： r model nlme mixed

【解决方案1】：

这几乎总是由模型中的共线性引起的。特别是，如果数据中不存在您的因子变量的组合，就会发生这种情况（我认为）。 (1) 试试lme4::lmer() 的模型，看看你是否收到关于秩不足的警告：

lme4::lmer(Hour~Age*Ethnicity*GHUsedFlag*Disability + (Age|Id), data=mydf)

（2）建立固定效应模型矩阵，使用caret::findLinearCombos

X <- model.matrix(~Age*Ethnicity*GHUsedFlag*Disability, data=mydf)
caret::findLinearCombos(X)

【讨论】：

谢谢本。我会试一试的。