使用先前估计的值重新启动混合效应模型估计
我正在使用lme4包中的lmer()来估计混合效应模型。这很有效，但现在我想运行一定数量的迭代估计过程，然后通过指定起始值来恢复该过程，就像上次估计过程计算的那样。 根据?lmer的帮助文档，可以通过设置以下参数实现： start - 这些是新的起始值，根据帮助文档，可以从拟合模型中提取...
在ggplot中绘制混合效应模型。
我刚开始学习混合效应模型，希望得到你的帮助。 我已经用 ggplot 绘制了以下图表： ggplot(tempEf,aes(TRTYEAR,CO2effect,group=Myc,col=Myc)) + facet_grid(~N) + geom_smooth(method="lm...
增加新版本lmer的迭代次数？
我刚刚更新了lme4到1.0-4版本，当我运行我的混合效应模型lmer()时，它以前收敛但现在打印出这个警告：Warning message: In (function (fn, par, lower = rep.int(-Inf, n), upper = rep.int(Inf, : ...
在使用anova()函数测试lmer()模型的随机效应时，是否需要将refit=FALSE设置为假？
我目前在测试是否应该在我的lmer模型中包含某些随机效应。为此，我使用anova函数。到目前为止，我的步骤是使用lmer()函数调用拟合模型，并使用REML=TRUE(默认选项)。然后，在其中一个模型上调用anova()，该模型包括要测试的随机效应，而另一个模型则不包括。然而，众所周知，ano...
如何从使用缩放响应拟合的lmer()模型中取消缩放系数
我使用lme4包中的lmer()函数在R中拟合了一个模型。我对因变量进行了缩放： mod <- lmer(scale(Y) ~ X + (X | Z), data = df, ...
lmer（来自R包lme4）如何计算对数似然？
我正在尝试理解函数lmer。我已经找到了有关如何使用该命令的大量信息，但几乎没有关于它实际在做什么的信息（除了一些神秘的评论在这里：http://www.bioconductor.org/help/course-materials/2008/PHSIntro/lme4Intro-handout...
如何在R中使用tryCatch
我希望使用try()或tryCatch()或类似的函数来检测我的模型"fit1"是否存在错误。如果模型没有问题，我想使用"fit1"，否则我想使用"fit2"。 fit1<-glmer(stat ~ dataint + DBH + DBH2 + (1|site_plot), famil...
Lmer模型的标准化系数
我曾经使用下面的代码来计算lmer模型的标准化系数。然而，随着新版lme4的发布，返回对象的结构已经发生了变化。 如何调整stdCoef.lmer函数，使其与新版lme4兼容？ # Install old version of lme 4 install.packages("lme4.0"...
如何在具有多个预测变量的混合模型中绘制随机截距和斜率？
当混合模型拥有多个预测变量时，是否可以绘制随机截距或斜率? 如果只有一个预测变量，我会这样做: #generate one response, two predictors and one factor (random effect) resp<-runif(100,1, 100) ...