从 lme4 mer 模型对象中提取随机效应方差
我有一个具有固定和随机效果的 mer 对象。如何提取随机效应的方差估计？这是我的问题的简化版本。 study <- lmer(Reaction ~ Days + (1|Subject), dat…
为什么 log() 语句会导致 lmList 失败？
在解决这个问题时：lme4和nlme包之间的干扰OP和我确定将响应向量包装在对 lme4::lmList 的调用内的 log() 中会导致 lmList 失败。 例如， my.lmList.m…
lme4 和 nlme 包之间的干扰
我有这个脚本： mod<-lmList(log(indwgt)~log(lngtclass)| sex, data=mbsp) 两周前就有效了。今天我感觉 Error in eval(expr, envir, enclos) : objec…
获取 nlme 或 lme4 中固定效应的广义最小二乘均值
可以使用 model.tables 函数获得 aov 对象的最小二乘均值及其标准误差： npk.aov <- aov(yield ~ block + N*P*K, npk) model.tables(npk.aov, "means"…
lme4 中的多元线性混合模型
我想知道如何使用 lme4 拟合多元线性混合模型。我使用以下代码拟合单变量线性混合模型： library(lme4) lmer.m1 <- lmer(Y1~A*B+(1|Block)+(1|Block:A…
有什么方法可以从 lme4 mer 模型拟合对象生成 LaTeX 表吗？
有谁知道如何从 lme4 mer 对象生成一个好的出版质量的 LaTeX 表？ xtable 方法（包 xtable）和 latex 方法（包 Hmisc）都不知道如何处理 mer 对象。 …
从 lmList 对象中提取 RMSE
我有以下代码从 lmList 对象中提取系数。 library(lme4) library(plyr) lm.model <- lmList(Y ~ X | eventID, df) param <- ldply(lm.model, coef) 这…
将 sim() 与 lmer() 一起使用
我使用相同的预测变量运行了两个多级逻辑回归，但有两个不同的响应： fruitMLM <- lmer(InsuffFruit ~ Income + HDI + Income:HDI + (1 + Income | Co…
LanguageR/ lme4：对具有随机相关参数的模型的 p 值有任何帮助吗？
我使用包 languageR 来实现混合效果模型，其语法在本文末尾。我可以使用 pvals.fnc 获取模型 1 和 3（hd_lmer1 和 hd_lmer2）的 p 值。将其与模型二一…
使用 Lattice （或其他东西）在 R 中绘制 lme4 的回归结果
感谢 之前的答案，我已经使用 lme4 进行了回归。现在我已经对每个州进行了回归拟合，我想使用点阵来绘制每个州的 QQ 图。我还想以网格格式绘制每个状…
