Thank you Daniel, I did try to change xlim's and ylim's, but the fitted line still doesn't appear on the plot. Below is a small expample I've created (I'm sorry for using names, but I het confused with x's and y's):
>mat<-matrix(c(1,2,3,4,5,6,7,8,1,2,3,4,1,2,3,4,1,5,2,7,3,7,5,8),nrow=8,ncol=3) >mats<-ts(mat[,3],start=c(1,1),freq=4) >quart<-mat[,2] >time<-seq(1,length(mats),1) >nmat<-data.frame(time=c(9),quart=c(1)) >quar<-factor(nmat$quart) >t1<-lm(mats~time) >plot(mats,xlim=c(1,3.25),ylim=c(0,15)) >lines(fitted(t1)~time,col="blue") >points(c(3),predict(t1,newdata=nmat),col="red") Unfortunately for me, this example works and does plot fitted line and predicted points. I really don't understand why I can't make the same commands work when trying with my real life data. Anyway, thanks for getting back to me earlier on. Regards Ula uep wrote: > > I'm trying to plot a time series (quarterly observations, seasonal effect, > T=56),regression line and predicted point on the same plot. > I'm using following commands: >> plot(qdts,xlim=c(1982,1997)) >> lines(fitted(m2)~time,data=qd,col="red") >> points(predict(m2,newdata=nqd),col="blue") > > where: > qdts<-ts(qd$durables,start=c(1982,1),freq=4) > m2<-lm(qdts~time+quar) > time<-seq(1,length(qdts),1) > nqd<-data.frame(time=c(57),quar=c(1)) > > Unfortunately once I enter all three commands the fitted line and > predicted point don't appear on the plot. What surprises me is that I am > not getting any error message. > > What I am doing wrong? > Many thanks in advance > > Ula > > -- View this message in context: http://www.nabble.com/Plotting-Standard-Regression-Fit-Am-I-wrong-or-a-bug--tp20399460p20401091.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.
