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--tp20399460p20399460.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.