散布図とallEffects plot. 以下にデータの例をいくつか示します。
library("effects")
x<-c(1,2,6,7,4,3,5)
y<-c(5,6,3,6,9,4,4)
a<-as.factor(c(1,1,1,2,2,2,2))
xyplot(y ~ x|a)
lm.y<-lm(y~x*a)
plot(allEffects(lm.y))
明確にするために、プロットにxyplotレイヤーを追加したいと思います。plot(allEfects...)アドバイスをいただければ幸いです。
このような ggplot2 パッケージをお勧めします。
ggplot(data, aes(y=y, x=x)) +
geom_point(shape = 1) +
geom_smooth(method="lm", fullrange=T) +
facet_grid(~a) +
theme(aspect.ratio = 1)
プロットを重ねることができます。
そうは言っても、 allEffects() 情報を組み込む方法がわからないため、これは完全な回答ではありません。他の誰かがそれについてより良いフィードバックを提供できるかもしれません.