# -------------------------
# P R I K L A D   2.1
# -------------------------

dir()
data <- read.delim('17-anova-newb.txt', sep = ',')
head(data, 5)
sum(is.na(data))
data <- na.omit(data)
dim(data)

# -------------------------
# P R I K L A D   2.2
# -------------------------

data$edu.M <- factor(data$edu.M, labels = c('ZS', 'SS', 'SSm', 'VS'))
head(data)
data$weight.K <- factor(data$weight.K, labels = c('nizka', 'norma', 'vysoka'))
head(data)

# -------------------------
# P R I K L A D   2.3
# -------------------------

edu <- data$edu.M
nj <- as.numeric(table(edu))
sum(edu == 'ZS')
pj <- nj / 1276
Nj <- cumsum(nj)
Fj  <- cumsum(pj)

edu.name <- c('ZS', 'SS', 'SSm', 'VS')
edu.rada <- data.frame(nj, pj, Nj, Fj, row.names = edu.name)
round(edu.rada, 4)

# -------------------------
# P R I K L A D   2.4
# -------------------------

source('Sbirka-AS-I-2018-funkce.R')
round(variacni.rada(edu, row.names = edu.name), 4)

barplot(nj, col = 'yellow', border = 'red',
        density = 30, ylim = c(0, 500),
        xlab = 'nejvyssi vzdelani',
        ylab = 'absolutni cetnost',
        names = edu.name, las = 1)

abline(h = seq(0, 500, by = 100), col = 'grey80', lty = 2)
box(bty = 'o', col = 'red', lty = 2)

rel.barplot(nj)
rel.barplot(nj, names = edu.name, col = c('red', 'blue', 'green', 'black'))
# -------------------------
# P R I K L A D   2.5
# -------------------------

wei <- data$weight.K

KT.abs <- table(edu, wei)
KT.abs
nj. <- apply(KT.abs, 1, sum)


KT.abs <- cbind(KT.abs, nj.)
n.k <- apply(KT.abs, 2, sum)
KT.abs <- rbind(KT.abs, n.k)

KT.rel <- KT.abs / 1276
round(KT.rel, digits = 4)

KT.abs <- table(edu, wei)

# -------------------------
# P R I K L A D   2.6
# -------------------------

radek.p <- prop.table(KT.abs, margin = 1)
round(radek.p, 4)

sloupec.p <- prop.table(KT.abs, margin = 2)
round(sloupec.p, 4)

# -------------------------
# P R I K L A D   2.7
# -------------------------
data <- read.delim('01-one-sample-mean-skull-mf.txt', sep = '\t')
head(data, n = 4)
sum(is.na(data))
data <- na.omit(data)
dim(data)

# -------------------------
# P R I K L A D   2.8
# -------------------------

skull.BM <- data[data$sex == 'm', 'skull.B' ] 
skull.BM
(n <- length(skull.BM))
range(skull.BM)

(r <- round(1 + 3.3*log10(n)))

(b <- seq(from = 123, to = 150, by = 3))

(centr <- seq(from = 124.5, to = 148.5, by = 3))

par(mfrow = c(1,2))
hist(skull.BM, breaks = b, ylim = c(0, 52), col = 'yellow',
     border = 'red', density = 20, 
     xlab = 'nejvetsi sirka mozkovny (mm)',
     ylab = 'absolutni cetnost', main = '', las = 1,
     axes = F)

axis(side = 1, centr)
axis(side = 2, las = 1)
box(bty = 'o')

boxplot(skull.BM, type = 2, horizontal = T, col = 'khaki',
        border = 'darkred', medcol = 'red',
        xlab = 'nejvetsi sirka mozkovny (mm)')