いくつかの lmer モデルからいくつかの ggplot ファセットに適合する係数に基づいて、いくつかの線をプロットしようとしています。ファセットごとに異なる線をプロットしたいと思います。これを達成する最良の方法は何ですか？
フィットを取得するために、次のアプローチを使用しています。これは二項 lmer モデルであるためです。
model.coefs <- fixef(model)
curve( invlogit( cbind(1, x) %*% model.coefs ), add=TRUE )
これは、この質問から取られています。変数に保存すると、グラフにプロットしたい一連の x 座標と y 座標が曲線関数によって得られます。
ベースRグラフ関数を使用して適切なプロットを作成することはできますが、代わりにggplotを使用したいと思います。したがって、これを実行して 2 つの異なるモデルの 2 つの異なるプロットを取得するのではなく、ここでは 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を含む) をまとめようとしていることにあると思いますが、これは簡単に解決できるはずです。ここで間違っている、または行方不明になっているのは何ですか?