我正在尝试根据某些 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?]**)
我怀疑问题的一部分是我试图一次将几个示例组合在一起(包括这个、这个和这个),但这感觉应该很容易解决。我在这里做错了什么或错过了什么?