### R code from vignette source 'M5201cv01new.Rnw'
### Encoding: CP1250

###################################################
### code chunk number 1: M5201cv01new.Rnw:38-39
###################################################
options(continue=" ")


###################################################
### code chunk number 2: M5201cv01new.Rnw:42-44
###################################################
lattice.options(default.theme = "col.whitebg")
options(SweaveHooks = list(fig = function() par(mfrow=c(1,1))),width=90)


###################################################
### code chunk number 3: M5201cv01new.Rnw:78-79 (eval = FALSE)
###################################################
## pracka <- scan(file="clipboard", sep=";", dec=",")


###################################################
### code chunk number 4: M5201cv01new.Rnw:89-91 (eval = FALSE)
###################################################
## pocitac <- scan(file="clipboard", sep=";", dec=",")
## 


###################################################
### code chunk number 5: M5201cv01new.Rnw:92-94
###################################################
pracka<-c(63.7, 65.6, 69.7, 72.5, 78.1, 80.1, 81.0, 84.1, 86.4, 88.7, 89.9, 91.2, 91.7, 91.9, 93.4)
pocitac<-c(5.2, 5.3, 7.2, 8.5, 11.5, 13.6, 17.5, 21.4, 25.2, 28.4, 34.2, 39.7, 41.9, 48.6, 55.6)


###################################################
### code chunk number 6: M5201cv01new.Rnw:102-103
###################################################
vybaveni <- data.frame(rok = 1993:2007, pracka, pocitac)


###################################################
### code chunk number 7: M5201cv01new.Rnw:108-110
###################################################
str(vybaveni)



###################################################
### code chunk number 8: M5201cv01new.Rnw:116-118
###################################################
head(vybaveni)
tail(vybaveni)


###################################################
### code chunk number 9: M5201cv01new.Rnw:132-133
###################################################
vybaveniTS <- ts(vybaveni[,-1], start = 1993, frequency = 1)


###################################################
### code chunk number 10: M5201cv01new.Rnw:138-139
###################################################
str(vybaveniTS)


###################################################
### code chunk number 11: M5201cv01new.Rnw:149-150
###################################################
data.library<-"C:/Users/Levakova/Documents/DSP/cviceni/M5201/data/"


###################################################
### code chunk number 12: M5201cv01new.Rnw:162-163
###################################################
fileDat <- paste(data.library, "airtransport.txt" ,sep = "")


###################################################
### code chunk number 13: M5201cv01new.Rnw:171-172
###################################################
airtrans <- read.table(fileDat, header = TRUE)


###################################################
### code chunk number 14: M5201cv01new.Rnw:176-179
###################################################
str(airtrans)
head(airtrans)
tail(airtrans)


###################################################
### code chunk number 15: M5201cv01new.Rnw:186-188
###################################################
airtransTS<-ts(airtrans[,-1], start = 2006, frequency = 12) 
str(airtransTS)


###################################################
### code chunk number 16: M5201cv01new.Rnw:210-214
###################################################
fileTxt <- paste(data.library, "tourism_israel.inf", sep = "")
con <- file(fileTxt)
(popis <- readLines(con))
close(con)


###################################################
### code chunk number 17: M5201cv01new.Rnw:221-224
###################################################
fileDat <- paste(data.library, "tourism_israel.txt", sep = "")
turiste <- read.table(fileDat, header = FALSE)
str(turiste)


###################################################
### code chunk number 18: M5201cv01new.Rnw:229-232
###################################################
names(turiste) <- c("time", "America", "Europe", "Other countries")
head(turiste)
tail(turiste)


###################################################
### code chunk number 19: M5201cv01new.Rnw:237-239
###################################################
turisteTS<-ts(turiste[,-1], start = 1991, frequency = 12) 
str(airtransTS)


###################################################
### code chunk number 20: M5201cv01new.Rnw:253-254
###################################################
getOption("SweaveHooks")[["fig"]]()
plot(vybaveniTS, type = "o")


###################################################
### code chunk number 21: M5201cv01new.Rnw:266-267
###################################################
getOption("SweaveHooks")[["fig"]]()
plot(vybaveniTS,type="o",yax.flip=TRUE)


###################################################
### code chunk number 22: M5201cv01new.Rnw:281-283
###################################################
library(Hmisc)
describe(vybaveni[,2:3])


###################################################
### code chunk number 23: M5201cv01new.Rnw:298-299
###################################################
getOption("SweaveHooks")[["fig"]]()
plot(airtransTS, type = "o", pch = 20, cex = 3)


###################################################
### code chunk number 24: M5201cv01new.Rnw:314-319
###################################################
getOption("SweaveHooks")[["fig"]]()
k <- dim(airtransTS)[2]
plot(airtransTS,type = "o", plot.type = "single", lty = 1:k,pch=1:k)
legend("topleft",
legend=c("Hungary","Austria","Slovakia","Czech Republic"),
lty = 1:k,pch=1:k, bty = "n")


###################################################
### code chunk number 25: M5201cv01new.Rnw:333-341
###################################################
getOption("SweaveHooks")[["fig"]]()
layout(matrix(c(1,2), nrow = 1), widths = c(0.8, 0.22))
par(mar = c(5, 4, 4, 2) + 0.1)
plot(airtransTS,type = "o", plot.type = "single", lty = 1:k,pch=1:k)
par(mar = c(5, 0, 4, 2) + 0.1)
plot(1:3, rnorm(3), pch = 1, lty = 1, ylim=c(-2,2), type = "n", axes = FALSE, ann = FALSE)
legend("left",
legend=c("Hungary","Austria","Slovakia","Czech Republic"),
lty = 1:k,pch=1:k, bty = "n")


###################################################
### code chunk number 26: M5201cv01new.Rnw:355-356
###################################################
describe(airtrans$Czech.Republic)


###################################################
### code chunk number 27: M5201cv01new.Rnw:363-364
###################################################
options(SweaveHooks = list(fig = function() par(mfrow=c(1,1))),width=70)


###################################################
### code chunk number 28: M5201cv01new.Rnw:369-382
###################################################
getOption("SweaveHooks")[["fig"]]()
x <- airtransTS[,4]
par(mar = c(4,2,1,0) + 0.75)
h <- hist(x, probability = TRUE, breaks = "FD",
        xlab = "Number of passengers",
        main = "Histogram, Kernel Density Estimate, Normal Curve")
xfit <- seq(min(x), max(x), length=512)
yfit <- dnorm(xfit, mean = mean(x), sd = sd(x))
lines(xfit, yfit, col = "dodgerblue", lwd = 2)
lines(density(x,n=512), lwd = 2, col = "red", lty = 2)
rug(x, side = 1, ticksize = 0.02, col = "grey20")
legend("topright",
legend=c("Histogram","Kernel Density Estimate","Normal Curve"),
col = c("black","red","dodgerblue"), lty = c(1,2,1), bty = "n")


###################################################
### code chunk number 29: M5201cv01new.Rnw:401-402
###################################################
getOption("SweaveHooks")[["fig"]]()
plot(turisteTS, type = "o", pch = 20, cex = 3)


###################################################
### code chunk number 30: M5201cv01new.Rnw:412-416
###################################################
getOption("SweaveHooks")[["fig"]]()
k <- dim(turisteTS)[2]
plot(turisteTS,type = "o", plot.type = "single", col = 1:k, pch = 20, cex = 3)
legend("top",
legend=c("America","Europe","Other countries"), lty = c(1,1,1), col = 1:k, bty = "n")


###################################################
### code chunk number 31: M5201cv01new.Rnw:429-430
###################################################
describe(turiste$Europe)


###################################################
### code chunk number 32: M5201cv01new.Rnw:436-437
###################################################
options(SweaveHooks = list(fig = function() par(mfrow=c(1,1))),width=70)


###################################################
### code chunk number 33: M5201cv01new.Rnw:442-455
###################################################
getOption("SweaveHooks")[["fig"]]()
x <- turisteTS[,2]
par(mar = c(4,2,1,0) + 0.75)
h <- hist(x, probability = TRUE, breaks = "FD",
        xlab = "Number of tourists from Europe",
        main = "Histogram, Kernel Density Estimate, Normal Curve")
xfit <- seq(min(x), max(x), length=512)
yfit <- dnorm(xfit, mean = mean(x), sd = sd(x))
lines(xfit, yfit, col = "dodgerblue", lwd = 2)
lines(density(x,n=512), lwd = 2, col = "red", lty = 2)
rug(x, side = 1, ticksize = 0.02, col = "grey20")
legend("topright",
legend=c("Histogram","Kernel Density Estimate","Normal Curve"),
col = c("black","red","dodgerblue"), lty = c(1,2,1), bty = "n")


