### R code from vignette source 'C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv09.Rnw'
### Encoding: CP1250

###################################################
### code chunk number 1: M5201cv09.Rnw:83-84
###################################################
options(continue=" ")


###################################################
### code chunk number 2: M5201cv09.Rnw:87-89
###################################################
lattice.options(default.theme = "col.whitebg")
options(SweaveHooks = list(fig = function() par(mfrow=c(1,1))),width=70)


###################################################
### code chunk number 3: M5201cv09.Rnw:485-487
###################################################
fileDat<-paste(data.library,"wool.txt",sep="")
wool<-scan(fileDat)


###################################################
### code chunk number 4: M5201cv09.Rnw:494-498
###################################################
getOption("SweaveHooks")[["fig"]]()
woolTS<-ts(wool,start=1965,frequency=4)
str(woolTS)
par(mar=c(2,2,0,0)+0.5)
plot(woolTS)


###################################################
### code chunk number 5: M5201cv09.Rnw:526-531
###################################################
x<-woolTS
model1 <- HoltWinters(x,beta=FALSE,gamma=FALSE)
summary(model1)
model1$alpha
coefficients(model1)


###################################################
### code chunk number 6: M5201cv09.Rnw:552-556
###################################################
getOption("SweaveHooks")[["fig"]]()
pred1<-predict(model1, n.ahead=12, prediction.interval=TRUE)
summary(pred1)
par(mar=c(2,2,1,0)+0.5)
plot(model1, predicted.values=pred1)


###################################################
### code chunk number 7: M5201cv09.Rnw:590-595
###################################################
model2 <- HoltWinters(x,gamma=FALSE)
summary(model2)
model2$alpha
model2$beta
coefficients(model2)


###################################################
### code chunk number 8: M5201cv09.Rnw:608-612
###################################################
getOption("SweaveHooks")[["fig"]]()
pred2 <- predict(model2, n.ahead=12, prediction.interval=TRUE)
summary(pred2)
par(mar=c(2,2,1,0)+0.5)
plot(model2, predicted.values=pred2)


###################################################
### code chunk number 9: M5201cv09.Rnw:650-656
###################################################
model3 <- HoltWinters(x)
summary(model3)
model3$alpha
model3$beta
model3$gamma
coefficients(model3)


###################################################
### code chunk number 10: M5201cv09.Rnw:663-666
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mar=c(2,2,0,0)+0.5)
plot(1:4,coefficients(model3)[3:6],type="o")
for(k in 1:4) abline(v=k,col="gray",lty=1)


###################################################
### code chunk number 11: M5201cv09.Rnw:686-690
###################################################
getOption("SweaveHooks")[["fig"]]()
pred3 <- predict(model3, n.ahead=12, prediction.interval=T)
summary(pred3)
par(mar=c(2,2,1,0)+0.5)
plot(model3, predicted.values=pred3)


###################################################
### code chunk number 12: M5201cv09.Rnw:704-705
###################################################
library(forecast)


###################################################
### code chunk number 13: M5201cv09.Rnw:718-721
###################################################
x<-woolTS
modelETS1<-ets(x,damped=FALSE,additive.only=TRUE)
summary(modelETS1)


###################################################
### code chunk number 14: M5201cv09.Rnw:725-726
###################################################
modelETS1$par


###################################################
### code chunk number 15: M5201cv09.Rnw:743-745
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mar=c(2,2,1,0)+0.5)
plot(modelETS1)


###################################################
### code chunk number 16: M5201cv09.Rnw:756-758
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mar=c(2,2,1,0)+0.5)
plot(modelETS1,plot.type="single",col=1:3,ylab="")


###################################################
### code chunk number 17: M5201cv09.Rnw:773-777
###################################################
getOption("SweaveHooks")[["fig"]]()
predETS1 <- forecast(modelETS1)
summary(predETS1)
par(mar=c(2,2,1,0)+0.5)
plot(predETS1)


###################################################
### code chunk number 18: M5201cv09.Rnw:791-794
###################################################
x<-woolTS
modelETS2 <- ets(x)
summary(modelETS2)


###################################################
### code chunk number 19: M5201cv09.Rnw:801-803
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mar=c(2,2,1,0)+0.5)
plot(modelETS2)


###################################################
### code chunk number 20: M5201cv09.Rnw:813-815
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mar=c(2,2,1,0)+0.5)
plot(modelETS2,plot.type="single",col=1:3,ylab="")


###################################################
### code chunk number 21: M5201cv09.Rnw:826-830
###################################################
getOption("SweaveHooks")[["fig"]]()
predETS2 <- forecast(modelETS2)
summary(predETS2)
par(mar=c(2,2,1,0)+0.5)
plot(predETS2)


