###################################################
### chunk number 1: 
###################################################
#line 46 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
options(continue=" ")


###################################################
### chunk number 2: 
###################################################
#line 50 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
lattice.options(default.theme = "col.whitebg")
options(SweaveHooks = list(fig = function() par(mfrow=c(1,1))),width=80)


###################################################
### chunk number 3: 
###################################################
#line 84 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
fileDat<-paste(data.library, "eriel.dat", sep = "")
(TXT<-paste(scan(fileDat, what=character(), nlines=1), collapse = " "))


###################################################
### chunk number 4: 
###################################################
#line 94 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
eriel <- read.table(fileDat, header=F, skip=2)
str(eriel)
head(eriel)


###################################################
### chunk number 5: 
###################################################
getOption("SweaveHooks")[["fig"]]()
#line 104 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
erielTS<-ts(eriel, start = 1921, frequency = 12)
par(mar = c(2,2,1,0) + 0.5)
plot(erielTS, main = TXT, cex.main = 0.85)


###################################################
### chunk number 6: 
###################################################
#line 120 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
Time <- time(erielTS)
model1 <- lm(erielTS ~ Time)
summary(model1)


###################################################
### chunk number 7: 
###################################################
#line 145 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
anova(model1)


###################################################
### chunk number 8: 
###################################################
#line 168 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
n<-length(Time)
gridt <- seq(Time[1],Time[n], length.out=300)
pred1 <- predict(model1, data.frame(Time=gridt), se.fit = T)
str(pred1)


###################################################
### chunk number 9: 
###################################################
#line 179 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
alpha <- 0.05
CV <- qt(1-alpha/2,pred1$df)


###################################################
### chunk number 10: 
###################################################
getOption("SweaveHooks")[["fig"]]()
#line 188 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
par(mar=c(2,2,1,0)+0.5)
plot(erielTS, main = TXT, cex.main = 0.85)
matlines(gridt, cbind(pred1$fit,pred1$fit-CV*pred1$se.fit,pred1$fit+CV*pred1$se.fit),lty = c(1,2,2), col = c("red","gray35","gray35"))


###################################################
### chunk number 11: 
###################################################
#line 203 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
Time<-time(erielTS)
model5<-lm(erielTS~Time+I(Time^2)+I(Time^3)+I(Time^4)+I(Time^5))
summary(model5)


###################################################
### chunk number 12: 
###################################################
#line 215 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
delta <- 0.5*(start(erielTS)[1] + end(erielTS)[1])
TimeC <- Time - delta
model5<-lm(erielTS ~ TimeC + I(TimeC^2) + I(TimeC^3)+ I(TimeC^4) + I(TimeC^5))
summary(model5)


###################################################
### chunk number 13: 
###################################################
#line 224 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
anova(model5)


###################################################
### chunk number 14: 
###################################################
#line 236 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
anova(model5,model1)


###################################################
### chunk number 15: 
###################################################
getOption("SweaveHooks")[["fig"]]()
#line 259 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
gridtC <- gridt-delta
pred5<-predict(model5,data.frame(TimeC=gridtC), se.fit= T)
alpha <- 0.05
CV <- qt(1-alpha/2, pred5$df)
par(mar=c(2,2,1,0)+0.5)
plot(erielTS, main = TXT, cex.main = 0.85)
matlines(gridt, cbind(pred5$fit,pred5$fit-CV*pred5$se.fit, pred5$fit+CV*pred5$se.fit),
        lty = c(1,2,2), col = c("red","gray35","gray35"))


###################################################
### chunk number 16: 
###################################################
#line 280 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
model9<-lm(erielTS~TimeC+I(TimeC^2)+I(TimeC^3)+I(TimeC^4)+I(TimeC^5)+I(TimeC^6)+I(TimeC^7)+I(TimeC^8)+I(TimeC^9))
summary(model9)
anova(model9)


###################################################
### chunk number 17: 
###################################################
#line 289 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
anova(model9,model5)


###################################################
### chunk number 18: 
###################################################
getOption("SweaveHooks")[["fig"]]()
#line 305 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
gridtC <- gridt-delta
pred9 <- predict(model9,data.frame(TimeC=gridtC),se.fit=T)
alpha <- 0.05
CV <- qt(1-alpha/2,pred9$df)
par(mar=c(2,2,1,0)+0.5)
plot(erielTS, main = TXT, cex.main = 0.85)
matlines(gridt, cbind(pred9$fit,pred9$fit-CV*pred9$se.fit,pred9$fit+CV*pred9$se.fit),
        lty = c(1,2,2), col = c("red","gray35","gray35"))


###################################################
### chunk number 19: 
###################################################
#line 432 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
fileSkript<-paste(script.library,"fndptr.R",sep="")
source(fileSkript)


###################################################
### chunk number 20: 
###################################################
getOption("SweaveHooks")[["fig"]]()
#line 447 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
fndptr(erielTS)


###################################################
### chunk number 21: 
###################################################
getOption("SweaveHooks")[["fig"]]()
#line 465 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
fndptr(erielTS,maxdgr=20)


###################################################
### chunk number 22: 
###################################################
#line 479 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
model18<-lm(erielTS~TimeC+I(TimeC^2)+I(TimeC^3)+I(TimeC^4)+I(TimeC^5)+I(TimeC^6)+
   I(TimeC^7)+I(TimeC^8)+I(TimeC^9)+I(TimeC^10)+I(TimeC^11)+I(TimeC^12)+I(TimeC^13)+
   I(TimeC^14)+I(TimeC^15)+I(TimeC^16)+I(TimeC^17)+I(TimeC^18))
summary(model18)
anova(model18)


###################################################
### chunk number 23: 
###################################################
#line 491 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
anova(model18,model9)


###################################################
### chunk number 24: 
###################################################
#line 506 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
model17<-lm(erielTS~TimeC+I(TimeC^2)+I(TimeC^3)+I(TimeC^4)+I(TimeC^5)+I(TimeC^6)+
   I(TimeC^7)+I(TimeC^8)+I(TimeC^9)+I(TimeC^10)+I(TimeC^11)+I(TimeC^12)+I(TimeC^13)+
   I(TimeC^14)+I(TimeC^15)+I(TimeC^16)+I(TimeC^17))
anova(model18,model17)


###################################################
### chunk number 25: 
###################################################
getOption("SweaveHooks")[["fig"]]()
#line 524 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
pred18<-predict(model18,data.frame(TimeC=gridtC),se.fit=T)
alpha<-0.05
CV<-qt(1-alpha/2,pred18$df)
par(mar=c(2,2,1,0)+0.5)
plot(erielTS,main=TXT,cex.main=0.85)
matlines(gridt,cbind(pred18$fit,pred18$fit-CV*pred18$se.fit,pred18$fit+CV*pred18$se.fit),
        lty=c(1,2,2),col=c("red","gray35","gray35"))


###################################################
### chunk number 26: 
###################################################
getOption("SweaveHooks")[["fig"]]()
#line 553 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
par(mfrow=c(1,2),mar=c(2,2,1,0)+0.5)
plot(model9,which=c(1,2))


###################################################
### chunk number 27: 
###################################################
getOption("SweaveHooks")[["fig"]]()
#line 577 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
par(mfrow=c(1,2),mar=c(2,2,1,0)+0.5)
plot(model18,which=c(1,2))


###################################################
### chunk number 28: 
###################################################
getOption("SweaveHooks")[["fig"]]()
#line 598 "C:/Users/Marie/Documents/DSP/cviceni/M5201/M5201cv03.Rnw"
par(mfrow=c(1,1))
par(mar=c(2,2,1,0)+0.5)
plot(erielTS,main=TXT,cex.main=0.85)
matlines(gridt,cbind(pred1$fit,pred5$fit,pred9$fit,pred18$fit),
        lty=c(1,1,1,1),col=c("magenta","dodgerblue","darkgreen","red"))
legend("bottomright",legend=paste("dgr =",c(1,5,9,18)),bty="n",
       lty=c(1,1,1,1),col=c("magenta","dodgerblue","darkgreen","red"))


