## pdf("F:\\Rmisc\\graphs\\winepop.pdf") winepop <- read.table("http://www.stat.ufl.edu/~winner/data/winepop.dat",header=F, col.names=c("year","tpop","pop5","pop14","pop24","pop34","pop44","pop54","pop64","pop65p","wine")) attach(winepop) adultpop <- (tpop-pop5-pop14)/1000 wine.mod1 <- lm(wine ~ adultpop) summary(wine.mod1) library(lmtest) dwtest(wine ~ adultpop) library(nlme) wine.mod1a <- gls(wine ~ adultpop, method='ML') summary(wine.mod1a) wine.mod2 <- gls(wine ~ adultpop, correlation=corAR1(), method='ML') summary(wine.mod2) wine.mod3 <- gls(wine ~ adultpop, correlation=corARMA(p=1,q=0), method='ML') summary(wine.mod3) plot(ACF(wine.mod1a,maxLag=30), alpha=0.05) ## dev.off()