7.3 Exponentielles Glätten (exponential smoothing)

# 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()