Advertisement
Guest User

Untitled

a guest
Jul 26th, 2016
61
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.74 KB | None | 0 0
  1. matplot(matrix(log(AirPassengers), ncol=12), type='l', axes=F, ylab='log(Passengers)', xlab='Month')
  2. axis(1, at=1:12, labels=month.abb)
  3. axis(2)
  4. box()
  5. title('Monthly air passenger data 1949:1960')
  6.  
  7. ap <- data.frame('lflights'= log(c(AirPassengers)), month=month.abb, year=rep(1949:1960, each=12))
  8. ap$month.n <- match(ap$month, month.abb)
  9. ap$monthly.diff <- lm(lflights ~ factor(year), data=ap)$residuals
  10. matplot(matrix(ap$monthly.diff, ncol=12), type='l', axes=F, ylab='log(Passengers)', xlab='Month')
  11. axis(1, at=1:12, labels=month.abb)
  12. axis(2)
  13. box()
  14. title('Monthly residuals for air passenger data 1949:1960')
  15.  
  16. library(splines)
  17. ap$monthly.pred <- lm(monthly.diff~bs(month.n, degree=2, knots = c(5)), data=ap)$fitted
  18. lines(1:12, ap$monthly.pred[1:12], lwd=2)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement