使用先前估计的值重新启动混合效应模型估计
我正在使用lme4包中的lmer()来估计混合效应模型。这很有效，但现在我想运行一定数量的迭代估计过程，然后通过指定起始值来恢复该过程，就像上次估计过程计算的那样。 根据?lmer的帮助文档，可以通过设置以下参数实现： start - 这些是新的起始值，根据帮助文档，可以从拟合模型中提取...
从lme4 mer模型对象中提取随机效应方差
我有一个mer对象，其中包含固定效应和随机效应。如何提取随机效应的方差估计值？以下是我问题的简化版本。study <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy) study 这会产生很长的输出 - 在这种情况下不算太长。...
R中的交叉重复效应和AR1协方差结构的线性混合模型
我拥有来自参与者（part）的受试者内生理数据，他们都在三轮（round）中查看刺激物（阅读报纸），每个轮次有五份报纸（paper），每个报纸中都有不同数量的访问（visit）。我有两个固定因素（CONDhier和CONDSbund）加上交互作用来预测生理状态（例如EDA），通常是自回归的。我...
使用qqmath或dotplot从lmer（lme4包）绘制随机效应：如何使其看起来漂亮？
qqmath函数使用lmer软件包的输出来生成随机效应的优秀毛虫图。也就是说，qqmath非常适合绘制分层模型中拦截的估计值及其误差的图表。下面是使用lme4软件包中名为Dyestuff的内置数据的lmer和qqmath函数示例。该代码将生成分层模型以及使用ggmath函数生成漂亮的图表。li...
如何在混合效应模型中获取系数及其置信区间？
在 lm 和 glm 模型中，我使用函数 coef 和 confint 来达到目标:m = lm(resp ~ 0 + var1 + var1:var2) # var1 categorical, var2 continuous coef(m) confint(m) 我现在将随机效应添加到模型中...
lme4::lmer报告“固定效应模型矩阵秩不足”，我需要修复吗？如何修复？
我正尝试使用混合效应模型，将其余的列作为预测变量来预测F2_difference，但是我收到了一个错误信息： fixed-effect model matrix is rank deficient so dropping 7 columns / coefficients. 根据这个...
lme4开发版本的收敛误差
我正在尝试使用lme4的开发版本进行混合效应模型的功效分析，并参考此教程。我注意到在教程中，lme4会出现收敛错误：## Warning: Model failed to converge with max|grad| = 0.00187101 (tol = ## 0.001) 当我使用我的数...
如何使用lme4比较没有随机效应的模型和有随机效应的模型？
我可以使用nlme包中的gls()函数来建立没有随机效应的mod1模型。 然后，我可以使用lme()函数来建立包含随机效应的mod2模型，并通过AIC对比mod1和mod2。 mod1 = gls(response ~ fixed1 + fixed2, method="REML", data...
有没有一种方法可以从`glmer`对象中获得“边际效应”？
我正在使用glmer估计随机效应logit模型，并且希望报告自变量的边际效应。对于glm模型，包mfx可以帮助计算边际效应。是否有任何用于glmer对象的包或函数？ 谢谢您的帮助。 下面给出一个可重现的示例## mydata <- read.csv("http://www.at...
如何从lme4 mer模型拟合对象生成LaTeX表格？
有没有人知道如何从lme4的mer对象中生成漂亮的出版质量LaTeX表格？无论是xtable方法（包xtable）还是latex方法（包Hmisc），都不知道如何处理mer对象。例如，给定以下拟合：library(lme4) fm1 <- lmer(Reaction ~ Days...