1

我正在尝试根据某些 lmer 模型拟合某些 ggplot 方面的系数绘制一些线。我想为每个方面绘制不同的线。实现这一目标的最佳方法是什么?

为了获得我的拟合,我使用以下方法,因为它是二项式 lmer 模型:

model.coefs <- fixef(model)
curve( invlogit( cbind(1, x) %*% model.coefs ), add=TRUE )

这是取自这个问题。保存到变量时,曲线函数会给出一堆我想在图表上绘制的 x 和 y 坐标。

我可以使用基本的 R 绘图函数让它生成一个不错的绘图,但我宁愿选择 ggplot。因此,与其为两个不同的模型获取两个不同的图,而不是在此处标记为模型A 和模型B:

plot(x, y)
modelA.coefs <- fixef(modelA)
curve( invlogit( cbind(1, x) %*% modelA.coefs ), add=TRUE )

plot(x, y)
modelB.coefs <- fixef(modelB)
curve( invlogit( cbind(1, x) %*% modelB.coefs ), add=TRUE )

我正试图让它在 ggplot 中工作。

ggplot(dataset)+
aes(x=x, y=y)+
geom_point()+
facet_wrap(~groupingFactor)+
geom_line(subset(dataset, groupFactor=="factorA"), **[would the solution be to feed this line something here?]**)

我怀疑问题的一部分是我试图一次将几个示例组合在一起(包括这个这个这个),但这感觉应该很容易解决。我在这里做错了什么或错过了什么?

4

0 回答 0