Dear R-users,
Sorry if someone came out with a similar question but after one day of
searching I am giving up:
Does anyone know how to plot the original points used in a lm model and the
set of resulting regression lines generated by the model?
This is how I do it using the "plot" and "lines" functions but I would like
to do it using lattice and I cannot find a way to plot in the same panel
points and lines that come from different datasets.
Note: I would like to use this sort of "predict" approach rather than using
the regression coefficients of the model as in complex models I get confused
when I have to combine the coefficients to build up each regression line.
Many thanks in advance and sorry for bothering you.
Javier
fm<-lm(c4~SEX+AREA+c3+SEX:AREA+SEX:c3+AREA:c3,data=mydata)
mygrid<-expand.grid(SEX=c("male","female"),AREA=c("area1","area2","area3"),c3=seq(-2.5,2.5,length.out=20))
pred_fm<-predict(fm,newdata=mygrid)
plot(mydata$c3, mydata$c4, col = ifelse(mydata$SEX== "male", "blue",
"red"), pch = ifelse(mydata$SEX=="male",2,1),xlab="c3",ylab="c4")
lines(seq(-2.5, 2.5, length.out = 20), pred_fm[mygrid$SEX == "female" &
mygrid$AREA == "area3"] , col = "red", lwd = 2, lty = 3)
lines(seq(-2.5, 2.5, length.out = 20), pred_fm[mygrid$SEX == "female" &
mygrid$AREA == "area1"] , col = "red", lwd = 2, lty = 1)
lines(seq(-2.5, 2.5, length.out = 20), pred_fm[mygrid$SEX == "female" &
mygrid$AREA == "area2"] , col = "red", lwd = 2, lty = 2)
lines(seq(-2.5, 2.5, length.out = 20), pred_fm[mygrid$SEX == "male" &
mygrid$AREA == "area3"] , col = "blue", lwd = 2, lty = 3)
lines(seq(-2.5, 2.5, length.out = 20), pred_fm[mygrid$SEX == "male" &
mygrid$AREA == "area1"] , col = "blue", lwd = 2, lty = 1)
lines(seq(-2.5, 2.5, length.out = 20), pred_fm[mygrid$SEX == "male" &
mygrid$AREA == "area2"] , col = "blue", lwd = 2, lty = 2)
--
View this message in context:
http://www.nabble.com/Plotting-regression-lines-and-points-from-lm-using-lattice-tp21052486p21052486.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
