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

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


###################################################
### code chunk number 2: M5201cv08.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: M5201cv08.Rnw:191-194
###################################################
ARphi<-c(0.9)
ARroots<-polyroot(c(1,-ARphi))
Mod(ARroots)


###################################################
### code chunk number 4: M5201cv08.Rnw:212-213
###################################################
set.seed(123456)


###################################################
### code chunk number 5: M5201cv08.Rnw:217-221
###################################################
getOption("SweaveHooks")[["fig"]]()
ARphi<-c(0.9)
ari.sim<-arima.sim(model=list(order=c(1,1,0),ar=ARphi),n=300)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(ari.sim)


###################################################
### code chunk number 6: M5201cv08.Rnw:232-235
###################################################
getOption("SweaveHooks")[["fig"]]()
vARIacf<-acf(ari.sim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vARIacf,main="Vyberova ACF")


###################################################
### code chunk number 7: M5201cv08.Rnw:242-245
###################################################
getOption("SweaveHooks")[["fig"]]()
vARIpacf<-pacf(ari.sim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vARIpacf,main="Vyberova PACF")


###################################################
### code chunk number 8: M5201cv08.Rnw:258-260
###################################################
ari.fit<-arima(ari.sim,order=c(1,1,0))
print(ari.fit)


###################################################
### code chunk number 9: M5201cv08.Rnw:268-272
###################################################
getOption("SweaveHooks")[["fig"]]()
fileFun<-paste(data.library,"FunkceM5201.R",sep="")
source(fileFun)
par(mfrow = c(1, 1), mar = c(2, 2, 1, 0) + 0.05)
PlotPredictARIMA(ari.sim, ari.fit, n.ahead = 10)


###################################################
### code chunk number 10: M5201cv08.Rnw:286-289
###################################################
getOption("SweaveHooks")[["fig"]]()
dif.ari<-diff(ari.sim)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(dif.ari)


###################################################
### code chunk number 11: M5201cv08.Rnw:300-303
###################################################
getOption("SweaveHooks")[["fig"]]()
difARIacf<-acf(dif.ari,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(difARIacf,main="Vyberova ACF")


###################################################
### code chunk number 12: M5201cv08.Rnw:312-315
###################################################
getOption("SweaveHooks")[["fig"]]()
vARIpacf<-pacf(dif.ari,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vARIpacf,main="Vyberova PACF")


###################################################
### code chunk number 13: M5201cv08.Rnw:327-329
###################################################
dif.ari.fit<-arima(dif.ari,order=c(1,0,0))
print(dif.ari.fit)


###################################################
### code chunk number 14: M5201cv08.Rnw:377-379
###################################################
MAtheta<-c(0.75,0.2)
Mod(polyroot(c(1,MAtheta)))


###################################################
### code chunk number 15: M5201cv08.Rnw:386-387
###################################################
set.seed(987654)


###################################################
### code chunk number 16: M5201cv08.Rnw:391-394
###################################################
getOption("SweaveHooks")[["fig"]]()
ima.sim<-arima.sim(model=list(order=c(0,1,2),ma=MAtheta),n=300)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(ima.sim)


###################################################
### code chunk number 17: M5201cv08.Rnw:405-408
###################################################
getOption("SweaveHooks")[["fig"]]()
vIMAacf<-acf(ima.sim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vIMAacf,main="Vyberova ACF")


###################################################
### code chunk number 18: M5201cv08.Rnw:417-420
###################################################
getOption("SweaveHooks")[["fig"]]()
vIMApacf<-pacf(ima.sim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vIMApacf,main="Vyberova PACF")


###################################################
### code chunk number 19: M5201cv08.Rnw:434-439
###################################################
getOption("SweaveHooks")[["fig"]]()
ima.fit<-arima(ima.sim,order=c(0,1,2))
print(ima.fit)

par(mfrow = c(1, 1), mar = c(2, 2, 1, 0) + 0.05)
PlotPredictARIMA(ima.sim, ima.fit, n.ahead = 10)


###################################################
### code chunk number 20: M5201cv08.Rnw:450-453
###################################################
getOption("SweaveHooks")[["fig"]]()
dif.ima<-diff(ima.sim)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(dif.ima)


###################################################
### code chunk number 21: M5201cv08.Rnw:464-467
###################################################
getOption("SweaveHooks")[["fig"]]()
difIMAacf<-acf(dif.ima,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(difIMAacf,main="Vyberova ACF")


###################################################
### code chunk number 22: M5201cv08.Rnw:475-478
###################################################
getOption("SweaveHooks")[["fig"]]()
vIMApacf<-pacf(dif.ima,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vIMApacf,main="Vyberova PACF")


###################################################
### code chunk number 23: M5201cv08.Rnw:495-497
###################################################
dif.ima.fit<-arima(dif.ima,order=c(0,0,2))
print(dif.ima.fit)


###################################################
### code chunk number 24: M5201cv08.Rnw:521-525
###################################################
ARphi<-c(0.9,-0.25)
Mod(polyroot(c(1,-ARphi)))
MAtheta<-c(0.75,0.2)
Mod(polyroot(c(1,MAtheta)))


###################################################
### code chunk number 25: M5201cv08.Rnw:532-533
###################################################
set.seed(123)


###################################################
### code chunk number 26: M5201cv08.Rnw:537-540
###################################################
getOption("SweaveHooks")[["fig"]]()
arimaSim<-arima.sim(model=list(order=c(2,1,2),ar=ARphi,ma=MAtheta),n=300)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(arimaSim)


###################################################
### code chunk number 27: M5201cv08.Rnw:551-554
###################################################
getOption("SweaveHooks")[["fig"]]()
vARIMAacf<-acf(arimaSim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vARIMAacf,main="Vyberova ACF")


###################################################
### code chunk number 28: M5201cv08.Rnw:561-564
###################################################
getOption("SweaveHooks")[["fig"]]()
vARIMApacf<-pacf(arimaSim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vARIMApacf,main="Vyberova PACF")


###################################################
### code chunk number 29: M5201cv08.Rnw:581-586
###################################################
getOption("SweaveHooks")[["fig"]]()
arima.fit<-arima(arimaSim,order=c(2,1,2))
print(arima.fit)

par(mfrow = c(1, 1), mar = c(2, 2, 1, 0) + 0.05)
PlotPredictARIMA(arimaSim, arima.fit, n.ahead = 10)


###################################################
### code chunk number 30: M5201cv08.Rnw:598-601
###################################################
getOption("SweaveHooks")[["fig"]]()
dif.arima<-diff(arimaSim)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(dif.arima)


###################################################
### code chunk number 31: M5201cv08.Rnw:614-617
###################################################
getOption("SweaveHooks")[["fig"]]()
difARIMAacf<-acf(dif.arima,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(difARIMAacf,main="Vyberova ACF")


###################################################
### code chunk number 32: M5201cv08.Rnw:625-628
###################################################
getOption("SweaveHooks")[["fig"]]()
vARIMApacf<-pacf(dif.arima,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vARIMApacf,main="Vyberova PACF")


###################################################
### code chunk number 33: M5201cv08.Rnw:641-643
###################################################
dif.arima.fit<-arima(dif.arima,order=c(2,0,2))
print(dif.arima.fit)


###################################################
### code chunk number 34: M5201cv08.Rnw:658-663
###################################################
getOption("SweaveHooks")[["fig"]]()
fileDat<-paste(data.library,"beer.txt",sep="")
beer<-read.table(fileDat,header=TRUE)
beerTS <- ts(beer, start = 1958, frequency = 12)
par(mar=c(2,2,3,0)+0.05)
plot(beerTS,main="Produkce piva v Australii")


###################################################
### code chunk number 35: M5201cv08.Rnw:679-681
###################################################
ima.beer<-arima(beerTS,order=c(0,1,1))
ima.beer


###################################################
### code chunk number 36: M5201cv08.Rnw:692-694
###################################################
beer1992<-predict(ima.beer,n.ahead=12)
sum(beer1992$pred)


###################################################
### code chunk number 37: M5201cv08.Rnw:703-704
###################################################
getOption("SweaveHooks")[["fig"]]()
acf(resid(ima.beer))


###################################################
### code chunk number 38: M5201cv08.Rnw:732-736
###################################################
getOption("SweaveHooks")[["fig"]]()
fileDat<-paste(data.library,"gnp96.dat",sep="")
gnp96<-read.table(fileDat)
gnpTS<-ts(gnp96[,2],start=1947,frequency=4)
plot(gnpTS)


###################################################
### code chunk number 39: M5201cv08.Rnw:750-751
###################################################
getOption("SweaveHooks")[["fig"]]()
acf(gnpTS,lag.max=50)


###################################################
### code chunk number 40: M5201cv08.Rnw:761-763
###################################################
getOption("SweaveHooks")[["fig"]]()
gnp.dif<-diff(gnpTS)
plot(gnp.dif)


###################################################
### code chunk number 41: M5201cv08.Rnw:784-786
###################################################
getOption("SweaveHooks")[["fig"]]()
gnpgr<-diff(log(gnpTS))
plot(gnpgr)


###################################################
### code chunk number 42: M5201cv08.Rnw:798-801
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mfrow=c(2,1),mar=c(2,2,3,0)+0.05)
acf(gnpgr,lag.max=24)
pacf(gnpgr,lag.max=24,main="")


###################################################
### code chunk number 43: M5201cv08.Rnw:820-824
###################################################
ma.gnpgr<-arima(gnpgr,order=c(0,0,2))
ma.gnpgr
ar.gnpgr<-arima(gnpgr,order=c(1,0,0))
ar.gnpgr


###################################################
### code chunk number 44: M5201cv08.Rnw:855-856
###################################################
ARMAtoMA(ar=0.35,ma=0,10)


###################################################
### code chunk number 45: M5201cv08.Rnw:869-870
###################################################
getOption("SweaveHooks")[["fig"]]()
PlotPredictARIMA(gnpgr,ma.gnpgr,n.ahead=30)


