使用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进行新水平的预测
我正在尝试拟合一个混合效应模型，然后使用该模型在可能具有不同水平的新数据集上生成估计值。我原以为在新数据集上的估计值将使用估计参数的平均值，但事实并非如此。以下是最小工作示例：library(lme4) d = data.frame(x = rep(1:10, times = 3), ...
Python Statsmodels Mixedlm（混合线性模型）随机效应
我对Statsmodels Mixedlm的输出有点困惑，希望有人能解释一下。 我有一个大型的单户住宅数据集，包括每个属性的前两个销售价格/销售日期。 我已经对整个数据集进行了地理编码，并获取了每个属性的海拔高度。 我正在尝试了解海拔高度与房产价格升值之间的关系在不同城市之间的变化方式。 ...
lmer（来自R包lme4）如何计算对数似然？
我正在尝试理解函数lmer。我已经找到了有关如何使用该命令的大量信息，但几乎没有关于它实际在做什么的信息（除了一些神秘的评论在这里：http://www.bioconductor.org/help/course-materials/2008/PHSIntro/lme4Intro-handout...
使用交叉随机效应的nlme
我正在尝试使用nlme包拟合一个交叉非线性随机效应模型，就像这个问题和这篇邮件列表帖子中提到的那样，将其作为线性随机效应模型。但无论我尝试什么，都会出现错误。以下是一个示例： library(nlme) ##### # simulate data set.seed(18112003) na...
如何在具有多个预测变量的混合模型中绘制随机截距和斜率？
当混合模型拥有多个预测变量时，是否可以绘制随机截距或斜率? 如果只有一个预测变量，我会这样做: #generate one response, two predictors and one factor (random effect) resp<-runif(100,1, 100) ...
R geepack：使用GEE得到过大的估计
我正在使用geepack进行R语言中的逻辑边际模型估计，使用geeglm()函数。但是我得到了垃圾估计值，它们大约比期望值高16个数量级。然而，p值似乎与预期相似。这意味着响应基本上变成了一个阶梯函数。请参见附图 以下是生成该图的代码： require(geepack) data = re...