Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- matplot(matrix(log(AirPassengers), ncol=12), type='l', axes=F, ylab='log(Passengers)', xlab='Month')
- axis(1, at=1:12, labels=month.abb)
- axis(2)
- box()
- title('Monthly air passenger data 1949:1960')
- ap <- data.frame('lflights'= log(c(AirPassengers)), month=month.abb, year=rep(1949:1960, each=12))
- ap$month.n <- match(ap$month, month.abb)
- ap$monthly.diff <- lm(lflights ~ factor(year), data=ap)$residuals
- matplot(matrix(ap$monthly.diff, ncol=12), type='l', axes=F, ylab='log(Passengers)', xlab='Month')
- axis(1, at=1:12, labels=month.abb)
- axis(2)
- box()
- title('Monthly residuals for air passenger data 1949:1960')
- library(splines)
- ap$monthly.pred <- lm(monthly.diff~bs(month.n, degree=2, knots = c(5)), data=ap)$fitted
- lines(1:12, ap$monthly.pred[1:12], lwd=2)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement