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

Question

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

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

Error in MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1

Answer 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)

Error in MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1... Linear mixed model, nlme, R