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

###################################################
### code chunk number 1: M5201cv05.Rnw:61-62
###################################################
options(continue=" ")


###################################################
### code chunk number 2: M5201cv05.Rnw:65-67
###################################################
lattice.options(default.theme = "col.whitebg")
options(SweaveHooks = list(fig = function() par(mfrow=c(1,1))),width=80)


###################################################
### code chunk number 3: M5201cv05.Rnw:400-405
###################################################
fileDat<-paste(data.library,"UKgas.txt",sep="")
con<-file(fileDat)
(popis<-readLines(con,n=1))
close(con)
UKgas <- scan(fileDat,skip=1)


###################################################
### code chunk number 4: M5201cv05.Rnw:414-417
###################################################
UKgasTS<-ts(UKgas,start=1963,frequency=4)
par(mar=c(2,2,1,0)+0.5)
plot(UKgasTS)


###################################################
### code chunk number 5: M5201cv05.Rnw:489-493
###################################################
fileSkript<-paste(data.library,"FunkceM5201.R",sep="")
source(fileSkript)
outp<-powtr(UKgas,seglen=8,figure=TRUE,location="median",variability="iqr")
str(outp)


###################################################
### code chunk number 6: M5201cv05.Rnw:525-526
###################################################
outp<-powtr(UKgas,seglen=8,figure2=TRUE,location="median",variability="iqr")


###################################################
### code chunk number 7: M5201cv05.Rnw:547-548
###################################################
outp<-powtr(UKgas,seglen=8,figure3=TRUE,location="median",variability="iqr")


###################################################
### code chunk number 8: M5201cv05.Rnw:569-570
###################################################
outp<-powtr(UKgas,seglen=4,figure3=TRUE,location="median",variability="iqr")


###################################################
### code chunk number 9: M5201cv05.Rnw:587-588
###################################################
outp<-powtr(UKgas,seglen=6,figure3=TRUE,location="median",variability="iqr")


###################################################
### code chunk number 10: M5201cv05.Rnw:605-606
###################################################
outp<-powtr(UKgas,seglen=10,figure3=TRUE,location="median",variability="iqr")


###################################################
### code chunk number 11: M5201cv05.Rnw:623-624
###################################################
outp<-powtr(UKgas,seglen=12,figure3=TRUE,location="median",variability="iqr")


###################################################
### code chunk number 12: M5201cv05.Rnw:672-674
###################################################
HistFit(UKgas)
mtext("original values",side=3,line=-0.5,cex=0.95)


###################################################
### code chunk number 13: M5201cv05.Rnw:688-689
###################################################
boxplotSegments(UKgas^(-0.5),seglen=8)


###################################################
### code chunk number 14: M5201cv05.Rnw:699-700
###################################################
HistFit(UKgas^(-0.5),xlab=expression(y == x^-0.5))


###################################################
### code chunk number 15: M5201cv05.Rnw:714-715
###################################################
boxplotSegments(UKgas^(-1/3),seglen=8)


###################################################
### code chunk number 16: M5201cv05.Rnw:725-726
###################################################
HistFit(UKgas^(-1/3),xlab=expression(y == x^-0.33))


###################################################
### code chunk number 17: M5201cv05.Rnw:740-741
###################################################
boxplotSegments(UKgas^(-0.25),seglen=8)


###################################################
### code chunk number 18: M5201cv05.Rnw:751-752
###################################################
HistFit(UKgas^(-0.25),xlab=expression(y == x^-0.25))


###################################################
### code chunk number 19: M5201cv05.Rnw:780-794
###################################################
y <- as.vector(UKgasTS)
x <- as.vector(time(UKgasTS))
n<-length(y)
nn<-300
data<-data.frame(x,y)
LinTrend <- lm(y ~ x+I(x^2),data=data)
print(summary(LinTrend))
new <- data.frame(x = seq(x[1],x[n],length.out=nn))
pred.w.plim <- predict(LinTrend, new, interval="prediction")
pred.w.clim <- predict(LinTrend, new, interval="confidence")
par(mar=c(2,2,1,0)+0.5)
matplot(new$x,cbind(pred.w.clim, pred.w.plim[,-1]),col=c(2,3,3,4,4),
        lty=c(1,2,2,3,3), type="l", ylab="predicted y")
lines(x,y)


###################################################
### code chunk number 20: M5201cv05.Rnw:811-813
###################################################
x<-resid(LinTrend)
boxplotSegments(x,seglen=8)


###################################################
### code chunk number 21: M5201cv05.Rnw:831-833
###################################################
x<-resid(LinTrend)+coef(LinTrend)[1]
outp<-powtr(x,seglen=8,figure=TRUE,location="median",variability="iqr")


###################################################
### code chunk number 22: M5201cv05.Rnw:852-853
###################################################
HistFit(outp$transfx,xlab=x)


###################################################
### code chunk number 23: M5201cv05.Rnw:876-882
###################################################
library(car)
TT<-time(UKgasTS)
data<-data.frame(TIME=TT-mean(TT),X=UKgas)
transf1<-powerTransform(X ~ TIME+I(TIME^2), data=data)
summary(transf1)
str(transf1)


###################################################
### code chunk number 24: M5201cv05.Rnw:891-892
###################################################
print(transf1$roundlam)


###################################################
### code chunk number 25: M5201cv05.Rnw:904-911
###################################################
m1 <- lm(X ~ TIME+I(TIME^2), data=data)
summary(m1)
outT <- powerTransform(m1)
summary(outT)
print(outT$roundlam)
m2<-update(m1,basicPower(outT$y,outT$roundlam)~.)
summary(m2)


###################################################
### code chunk number 26: M5201cv05.Rnw:920-922
###################################################
res<-resid(m2)
boxplotSegments(res,seglen=8,xlab="Residuals of Box-Cox model")


###################################################
### code chunk number 27: M5201cv05.Rnw:933-934
###################################################
HistFit(res,xlab="Residuals of Box-Cox model")


###################################################
### code chunk number 28: M5201cv05.Rnw:955-962
###################################################
m3 <- lm(X ~ TIME, data=data)
summary(m3)
outT <- powerTransform(m3)
summary(outT)
print(outT$roundlam)
m4<-update(m3,basicPower(outT$y,outT$roundlam)~.)
summary(m4)


###################################################
### code chunk number 29: M5201cv05.Rnw:976-978
###################################################
res<-resid(m4)
boxplotSegments(res,seglen=8,xlab="Residuals of Box-Cox model")


###################################################
### code chunk number 30: M5201cv05.Rnw:989-990
###################################################
HistFit(res,xlab="Residuals of Box-Cox model")


