###################################################
### chunk number 1: 
###################################################
options(continue=" ")


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


###################################################
### chunk number 3: 
###################################################
getOption("SweaveHooks")[["fig"]]()
x<-rnorm(300,mean=0,sd=1)
plot(x,type="l",xlab="t")


###################################################
### chunk number 4: 
###################################################
getOption("SweaveHooks")[["fig"]]()
acfWN<-acf(x,main="Vyberova ACF",xlab="k")


###################################################
### chunk number 5: 
###################################################
getOption("SweaveHooks")[["fig"]]()
pacfWN<-pacf(x,main="Vyberova PACF",xlab="k",ylim=c(-0.2,1))


###################################################
### chunk number 6: 
###################################################
Box.test(x)
Box.test(x,type="Ljung-Box")


###################################################
### chunk number 7: 
###################################################
ARphi<-c(0.9,-0.25)
ARroots<-polyroot(c(1,-ARphi))
Mod(ARroots)


###################################################
### chunk number 8: 
###################################################
getOption("SweaveHooks")[["fig"]]()
x<-seq(-1,1,length=1000)
y<-sqrt(1-x^2)
plot(Re(ARroots),Im(ARroots),type="p",pch=19,xlab="Realna cast",ylab="Imaginarni cast",main="Jednotkova kruznice",ylim=c(-3,3),xlim=c(-3,3))
lines(c(x,x),c(y,-y))
abline(h=0)
abline(v=0)
legend("topleft",legend="Koreny procesu AR(2)",pch=19)


###################################################
### chunk number 9: 
###################################################
set.seed(1234)


###################################################
### chunk number 10: 
###################################################
getOption("SweaveHooks")[["fig"]]()
ar.sim<-arima.sim(model=list(ar=ARphi),n=300)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(ar.sim)


###################################################
### chunk number 11: 
###################################################
getOption("SweaveHooks")[["fig"]]()
tARacf<-ARMAacf(ar=ARphi,lag.max=20)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(0:(length(tARacf)-1),tARacf,type="h",main="Teoreticka ACF")
abline(h=0)


###################################################
### chunk number 12: 
###################################################
getOption("SweaveHooks")[["fig"]]()
vARacf<-acf(ar.sim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vARacf,main="Vyberova PACF")


###################################################
### chunk number 13: 
###################################################
getOption("SweaveHooks")[["fig"]]()
tARpacf<-ARMAacf(ar=ARphi,lag.max=20,pacf=TRUE)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(1:length(tARpacf),tARpacf,type="h",main="Teoreticka PACF")
abline(h=0)


###################################################
### chunk number 14: 
###################################################
getOption("SweaveHooks")[["fig"]]()
vARpacf<-pacf(ar.sim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vARpacf,main="Vyberova ACF")


###################################################
### chunk number 15: 
###################################################
ar.fit<-arima(ar.sim,order=c(2,0,0))
print(ar.fit)


###################################################
### chunk number 16: 
###################################################
getOption("SweaveHooks")[["fig"]]()
tsdiag(ar.fit)


###################################################
### chunk number 17: 
###################################################
getOption("SweaveHooks")[["fig"]]()
fileSkript<-paste(data.library,"FunkceM5201.R",sep="")
source(fileSkript)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
PlotPredictARIMA(ar.sim,ar.fit,n.ahead=10)


###################################################
### chunk number 18: 
###################################################
library(forecast)
ar.fitF<-auto.arima(ar.sim,trace=TRUE)
print(ar.fitF)


###################################################
### chunk number 19: 
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(forecast(ar.fitF),n=10)


###################################################
### chunk number 20: 
###################################################
MAtheta<-c(0.75,0.2)
print(MAroots<-polyroot(c(1,MAtheta)))
Mod(MAroots)


###################################################
### chunk number 21: 
###################################################
getOption("SweaveHooks")[["fig"]]()
x<-seq(-1,1,length=1000)
y<-sqrt(1-x^2)
plot(Re(MAroots),Im(MAroots),type="p",pch=19,xlab="Realna cast",ylab="Imaginarni cast",main="Jednotkova kruznice",ylim=c(-3,3),xlim=c(-3,3))
lines(c(x,x),c(y,-y))
abline(h=0)
abline(v=0)
legend("topleft",legend="Koreny procesu MA(2)",pch=19)


###################################################
### chunk number 22: 
###################################################
set.seed(1234)


###################################################
### chunk number 23: 
###################################################
getOption("SweaveHooks")[["fig"]]()
ma.sim<-arima.sim(model=list(ma=MAtheta),n=300)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(ma.sim)


###################################################
### chunk number 24: 
###################################################
getOption("SweaveHooks")[["fig"]]()
tMAacf<-ARMAacf(ma=MAtheta,lag.max=20)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(0:(length(tMAacf)-1),tMAacf,type="h",main="Teoreticka ACF")
abline(h=0)


###################################################
### chunk number 25: 
###################################################
getOption("SweaveHooks")[["fig"]]()
vMAacf<-acf(ma.sim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vMAacf,main="Vyberova ACF")


###################################################
### chunk number 26: 
###################################################
getOption("SweaveHooks")[["fig"]]()
tMApacf<-ARMAacf(ma=MAtheta,lag.max=20,pacf=TRUE)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(1:length(tMApacf),tMApacf,type="h",main="Teoreticka PACF")
abline(h=0)


###################################################
### chunk number 27: 
###################################################
getOption("SweaveHooks")[["fig"]]()
vMApacf<-pacf(ma.sim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vMApacf,main="Vyberova PACF")


###################################################
### chunk number 28: 
###################################################
ma.fit<-arima(ar.sim,order=c(0,0,2))
print(ma.fit)


###################################################
### chunk number 29: 
###################################################
getOption("SweaveHooks")[["fig"]]()
tsdiag(ma.fit)


###################################################
### chunk number 30: 
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
PlotPredictARIMA(ma.sim,ma.fit,n.ahead=10)


###################################################
### chunk number 31: 
###################################################
ma.fitF<-auto.arima(ma.sim,trace=TRUE)
print(ma.fitF)


###################################################
### chunk number 32: 
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(forecast(ma.fitF),n=10)


###################################################
### chunk number 33: 
###################################################
ARphi<-c(0.9,-0.5)
ARroots<-polyroot(c(1,-ARphi))
Mod(ARroots)


###################################################
### chunk number 34: 
###################################################
MAtheta<-c(0.75,0.25)
MAroots<-polyroot(c(1,MAtheta))
Mod(MAroots)


###################################################
### chunk number 35: 
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mfrow=c(1,2))
plot(Re(ARroots),Im(ARroots),type="p",pch=19,xlab="Realna cast",ylab="Imaginarni cast",main="Koreny AR casti",ylim=c(-5,5),xlim=c(-1,8))
lines(c(x,x),c(y,-y))
abline(h=0)
abline(v=0)
plot(Re(MAroots),Im(MAroots),type="p",pch=19,xlab="Realna cast",ylab="Imaginarni cast",main="Koreny MA casti",ylim=c(-5,5),xlim=c(-4.5,4.5))
lines(c(x,x),c(y,-y))
abline(h=0)                                                                                                                         
abline(v=0)


###################################################
### chunk number 36: 
###################################################
set.seed(1234)


###################################################
### chunk number 37: 
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mfrow=c(1,1))
arma.sim<-arima.sim(model=list(ar=ARphi,ma=MAtheta),n=300)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(arma.sim)


###################################################
### chunk number 38: 
###################################################
getOption("SweaveHooks")[["fig"]]()
tARMAacf<-ARMAacf(ar=ARphi,ma=MAtheta,lag.max=20)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(0:(length(tARMAacf)-1),tARMAacf,type="h",main="Teoreticka ACF")
abline(h=0)


###################################################
### chunk number 39: 
###################################################
getOption("SweaveHooks")[["fig"]]()
vARMAacf<-acf(arma.sim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vARMAacf,main="Vyberova ACF")


###################################################
### chunk number 40: 
###################################################
getOption("SweaveHooks")[["fig"]]()
tARMApacf<-ARMAacf(ar=ARphi,ma=MAtheta,lag.max=20,pacf=TRUE)
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(1:length(tARMApacf),tARMApacf,type="h",main="Teoreticka PACF")
abline(h=0)


###################################################
### chunk number 41: 
###################################################
getOption("SweaveHooks")[["fig"]]()
vARMApacf<-pacf(arma.sim,lag.max=20,plot=FALSE)
par(mfrow=c(1,1),mar=c(2,2,3,0)+0.05)
plot(vARMApacf,main="Vyberova PACF")


###################################################
### chunk number 42: 
###################################################
arma.fit<-arima(arma.sim,order=c(2,0,2))
print(arma.fit)


###################################################
### chunk number 43: 
###################################################
getOption("SweaveHooks")[["fig"]]()
tsdiag(arma.fit)


###################################################
### chunk number 44: 
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
PlotPredictARIMA(arma.sim,arma.fit,n.ahead=10)


###################################################
### chunk number 45: 
###################################################
arma.fitF<-auto.arima(arma.sim,trace=TRUE)
print(arma.fitF)


###################################################
### chunk number 46: 
###################################################
getOption("SweaveHooks")[["fig"]]()
par(mfrow=c(1,1),mar=c(2,2,1,0)+0.05)
plot(forecast(arma.fitF),n=10)


