模型中的线性混合模型的平均系数?修复不再有效

时间:2017-10-26 05:13:32

标签: r lme4 mixed-models glmulti

我正在使用glmulti包对lme4中混合模型的固定效果进行变量选择。我在检索由this thread中的包的作者解决的系数和置信区间时遇到了同样的问题。即使用coefcoef.multi会出现check.names错误,并且在调用NULL方法时系数会列为predict。所以我尝试了上面链接的线程中列出的解决方案,使用:

setMethod('getfit', 'merMod', function(object, ...) {
summ=summary(object)$coef
summ1=summ[,1:2]
if (length(dimnames(summ)[[1]])==1) {
    summ1=matrix(summ1, nr=1, dimnames=list(c("(Intercept)"),c("Estimate","Std.    Error")))
}
cbind(summ1, df=rep(10000,length(fixef(object))))
})

我修复了原始帖子中遗漏的"并运行了代码。但是,现在而不是

  

data.frame(...,check.names = FALSE)中的错误:参数暗示   不同的行数:1,0

我为每个模型都会收到此错误...

  

计算Satterthwaite近似值时出错。输出   lme4包的返回摘要从lme4返回一些   lmerTest中发生了计算错误

我正在使用lmerTest,如果glmulti无法从模型中提取正确的信息,它就会失败并不会让我感到惊讶。所以真的是错误的前两行可能是应该关注的内容。

原始修补程序的说明位于开发人员网站here上。显然这个软件包暂时没有更新,是的,我应该学习一个新的软件包......但在那之前我希望能解决这个问题。我会直接通过他的网站联系开发人员。但是,与此同时,有没有人试过这个并找到了解决办法?

lme4 glmulti rJava和其他相关套餐已全部更新至最新版本。

0 个答案:

没有答案