如何解决MEEM(object, conLin, control$niterEM) 中的错误:在级别 0,块 1 处反向求解中的奇异性...线性混合模型,nlme,R
m1 =lme(fixed=Hour~Age*Ethnicity*GHUsedFlag*Disability,random=~Age|Id,data=mydf)
在我的模型中,我正在检查多级交互。年龄-数字重复测量、种族-特征因子、GHUsedFlag-逻辑、残疾-特征-因子、小时-数字、结果变量。
在模型中引入 Disability
时出现此错误。无法真正理解为什么会发生错误以及如何应对。任何帮助表示赞赏。
我在 R 中使用 nlme 包。
MEEM(object,conLin,control$niterEM) 中的错误: 在级别 0,块 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)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。