![]() |
# Zeitreihenanalyse:
# beer Daten von # http://134.76.173.220/beer.zip beer <- read.csv2( " beer.csv " ) beer <- ts(beer[,1],start=1956,freq=12) pdf("ts7.pdf",width=8,height=6) par(mfrow=c(2,1)) par(mar=c(4,4,2,1)) plot(beer,ylab="Produzierte Biermenge",xlab="Zeit",col=gray(0.5),xlim=c(1956,1999)) beer.hw <- HoltWinters(beer) lines(beer.hw$fitted[,1],col=4) plot(beer,ylab="Produzierte Biermenge",xlab="Zeit",col=gray(0.5),xlim=c(1956,1999)) lines(predict(beer.hw,n.ahead=48),col=4,lwd=2) dev.off() |