我对R中的格子库中的xyplot有疑问。
以下是数据示例:
set.seed(4)
mydata <- data.frame(x.data = rnorm(50),
y.data = rnorm(50),
type = rep(c("A","B"), 50))
head(mydata)
mod <- lm(x.data ~ y.data*type, data= mydata)
summary(mod)
xyplot(y.data + fitted(mod) ~ x.data, groups= type, data= mydata, auto.key=F)
如何仅将回归线添加到A型数据并且仅在图的拟合(mod)部分中添加。是否可以仅使用从最小值到最大值的回归线?
感谢您的帮助。我真的很感激。
答案 0 :(得分:1)
library(lattice)
library(latticeExtra)
set.seed(4)
mydata <- data.frame(x.data = rnorm(50),
y.data = rnorm(50),
type = rep(c("A","B"), 50))
head(mydata)
mod <- lm(x.data ~ y.data*type, data= mydata)
p1 <- predict(mod, newdata = data.frame(y.data = range(mydata$y.data), type = "A"))
xyplot(y.data + fitted(mod) ~ x.data, groups= type, data= mydata, auto.key=F) +
layer(panel.lines(range(mydata$x.data), p1, col = "black", lty = 1), rows = 1, columns = 2)