search for: lgaplag

Displaying 1 result from an estimated 1 matches for "lgaplag".

2011 Dec 13
2
Problem with ploting fitted values
...y I can only see the original values. Below is the code with data library: / library(forecast) data(AirPassengers) AP <- AirPassengers class(AP) start(AP) end(AP) frequency(AP) lgAP <- log(AP) t<-2:length(AP) APsL<-log(ts(lgAP, start=1949, frequency=12)) final<-(ts.intersect(lgAP,lgAPlag=lag(lgAP,-1), t=t, sin12 = sin(2*pi*t/12), cos12 = cos(2*pi*t/12),dframe=TRUE)) summary(AIRlm <- lm(lgAP ~ lgAPlag + t + sin12 + cos12, data=final)) plot(lgAP, main="Log of Air Passengers",type="l", col=4, lty=2, lwd=2) lines(t, AIRlm$fitted, col=2, lwd=2) legend("...