## (a)
# nacteni dat
data <- read.delim('./Data/test-koncetiny-trup.txt')
str(data) 
sum(is.na(data)) # zjisteni poctu NA hodnot - 0 znamena, ze zadne NA v souboru nemame

data$sex <- factor(data$sex)
data.f <- data[data$sex == "f", ] # vyber pouze zeny

dim(data.f) # dimenze datove tabulky: 100 3 znamena, ze mame 3 promenne (znaky) pro 100 objektu (tj. mame 100 zen)

## (b)
lowex <- data.f$lowex.L
tru <- data.f$tru.L
# vypocet kvantilu (dolni, median, horni) pro koncetinu
kv <- quantile(lowex, probs = c(0.25, 0.5, 0.75), type = 2)
# sikmost a spicatost
sikmost <- e1071::skewness(lowex, type = 3)    
spicatost <- e1071::kurtosis(lowex, type = 3)  

# sestaveni tabulky koncetina
tab <- data.frame(m = mean(lowex), s = sd(lowex), min = min(lowex),
                  dolni.kv = kv[1], median = kv[2], horni.kv = kv[3], 
                  max = max(lowex), sikmost, spicatost,
                  row.names = c('koncetina'))
(tab <- round(tab, digits = 3)) # zaokrouhleni

# charakteristiky pro trup
kvt <- quantile(tru, probs = c(0.25, 0.5, 0.75), type = 2)
sikmost <- e1071::skewness(tru, type = 3)    # 
spicatost <- e1071::kurtosis(tru, type = 3)  # 

# sestaveni tabulky trup
tabt <- data.frame(m = mean(tru), s = sd(tru), min = min(tru),
                  dolni.kv = kvt[1], median = kvt[2], horni.kv = kvt[3], 
                  max = max(tru), sikmost, spicatost,
                  row.names = c('trup'))
(tabt <- round(tabt, digits = 3))
rbind(tab,tabt) # spojeni tabulky do jedne, vysledek:

#                m      s min dolni.kv median horni.kv  max sikmost spicatost
# koncetina 940.50 45.471 836    909.0    943      969 1076   0.196    -0.017
# trup      424.47 31.502 358    401.5    419      450  492   0.234    -0.738

## interpretace: prumerna delka prave dolni koncetiny u zen je 940.5 mm se smerodatnou odchylkou 45.47
## namerene delky dolni koncetiny se pohybovaly v rozsahu 836 - 1076 mm, s dolnim kvartilem 909, 
## tj. 25% zen melo dolni koncetinu o delce 909 mm nebo kratsi, median 943 - polovina zen v souboru 
## mela delku dolni koncetiny mensi nebo rovnu 943 mm, horni kvartil 969 - 75% zen melo dolni koncetinu
## delky 969 mm nebo kratsi. Sikmost 0.196 ukazuje na mirne kladne zesikmeni dat (prodlouzeny pravy konec)
## zaporna spicatost (-0.0170) ukazuje na zplosteni dat, ovsem vzhledem k male hodnote pouze velmi mirne

## analogicky pro delku trupu, tam ovsem spicatost vychazi vyrazne vyssi, takze jde o znacne zplostela data
## a take vice (kladne) zesikmena

## (c)
par(mfrow = c(1,2), mar = c(4,5,1,1)) # nastaveni 2 obrazku vedle sebe
boxplot(data.f$lowex.L, border = "forestgreen",
        col = "darkseagreen2", xlab = "delka dolni koncetiny (mm)")
boxplot(data.f$tru.L, border = "purple", col = "plum",
        xlab = "delka trupu (mm)")
# odlehlou (extremni) hodnotu obsahuje pouze diagram pro delku koncetiny
