# Nastaveni pracovniho adresare (working directory)
setwd("C:/Jirka/Environment/Výuka/VLBS0621 Biostatistika pro P-PooL/Jaro 2022")


# Nacteni souboru zelvy do datoveho ramce
zelvy<-read.csv(file="02_zelvy.csv",header=TRUE,dec=",",sep=";",skip=0,encoding="Windows-1250")


# Zakladni popisne statistiky promennych v souboru
summary(zelvy)
summary(zelvy)[,5]
summary(zelvy)[,6]


# Frekvencni tabulka
table(zelvy[,5])


# Vypocet konkretnich statistik na jednotlivych sloupcich
mean(zelvy[,1])
sd(zelvy[,1])

mean(zelvy[,"Šířka"])
sd(zelvy[,"Výška"])

mean(zelvy$Výška)
sd(zelvy$Výška)

vek<-zelvy$Věk
mean(vek)
sd(vek)


# Histogram
hist(vek,main="Histogram věku želv",xlab="Věk",ylab="Počet jedinců",breaks=15)


# Hezci histogram
install.packages("genasis")
library(genasis)
genhistogram(vek,col="black",emboss=3,main="Histogram věku želv",xlab="Věk",ylab="Počet jedinců",distr="")