7.2.1 Lineares Modell

		# Zeitreihenanalyse: # beer Daten von # http://134.76.173.220/beer.zip beer <- read.csv2( " beer.csv " ) beer <- ts(beer[,1],start=1956,freq=12) # Lineare Regression library(rrcov) pdf("ts4.pdf",width=5,height=4) par(mar=c(4,4,2,1)) logbeer<-log(beer) t <- (1:length(beer))/12+1956 plot(logbeer,ylab="log(Produzierte Biermenge)",xlab="Zeit",col=gray(0.5)) res <- ltsReg(logbeer~t) lines(t,res$fit,col=4,lwd=1) dev.off()
		# Zeitreihenanalyse: # beer Daten von # http://134.76.173.220/beer.zip beer <- read.csv2( " beer.csv " ) beer <- ts(beer[,1],start=1956,freq=12) # Quadratische Regression logbeer<-log(beer) t <- (1:length(beer))/12+1956 t2<-t^2 pdf("ts5.pdf",width=5,height=4) par(mar=c(4,4,2,1)) plot(logbeer,ylab="log(Produzierte Biermenge)",xlab="Zeit",col=gray(0.5)) res <- ltsReg(logbeer~t+t2) lines(t,res$fit,col=4,lwd=1) dev.off()