pdf("I:\\data\\tamparev_r.pdf") tamparev <- read.fwf("I:\\data\\tamparev.txt", width=c(8,8,8,8), col.names=c("pct_occ", "ave_rr", "city", "month")) tamparev <- tamparev[1:72,] attach(tamparev) revpar <- pct_occ*ave_rr/100 x1<-cos(2*3.14159*1*month/12); x2<-sin(2*3.14159*1*month/12); x3<-cos(2*3.14159*2*month/12); x4<-sin(2*3.14159*2*month/12); x5<-cos(2*3.14159*3*month/12); x6<-sin(2*3.14159*3*month/12); x7<-cos(2*3.14159*4*month/12); x8<-sin(2*3.14159*4*month/12); x9<-cos(2*3.14159*5*month/12); x10<-sin(2*3.14159*5*month/12); x11<-cos(3.14159*month); x12<-month; model1 <- lm(revpar~x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11+x12) summary(model1) anova(model1) model2 <- lm(revpar~x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x12) summary(model2) anova(model2) anova(model1,model2) plot(month,revpar,pch=16,xlim=c(0,73),ylim=c(20,100),main="Tampa - Order=6 Trig Reg - Predicted") lines(month,predict(model1)) abline(v=seq(3,63,12)) plot(month,residuals(model1),pch=16,main="Tampa - Order=6 Trig Reg - Residuals") abline(h=0) abline(v=seq(3,63,12)) model3 <- lm(revpar~x1+x2+x3+x4+x5+x6+x7+x8+x12) summary(model3) anova(model3) anova(model2,model3) plot(month,revpar,pch=16,xlim=c(0,73),ylim=c(20,100),main="Tampa - Order=4 Trig Reg - Predicted") lines(month,predict(model3)) abline(v=seq(3,63,12)) plot(month,residuals(model3),pch=16,main="Tampa - Order=4 Trig Reg - Residuals") abline(h=0) abline(v=seq(3,63,12)) dev.off()