我试图根据某些lmer模型拟合到某些ggplot方面的系数绘制一些线条。我想为每个方面绘制不同的线条。实现这一目标的最佳方式是什么?
为了获得我的契合,我使用以下方法,因为它是二项式lmer模型:
model.coefs <- fixef(model)
curve( invlogit( cbind(1, x) %*% model.coefs ), add=TRUE )
这取自this question。当保存到变量时,曲线函数给出了一堆x和y坐标,我想在我的图上绘制。
我能够使用基本R图形函数来生成一个像样的图,但我宁愿选择ggplot。因此,不要这样做两个不同模型的两个不同的图,在这里标记为modelA和modelB:
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?]**)
我怀疑部分问题是我试图同时将几个例子拼凑在一起(包括this one,this one和this one),但这感觉应该很简单锻炼。我在这里做错了什么或错过了什么?