# Nastavení pracovního adresáře
setwd("C:/Jirka/Environment/Výuka/Bi7541 Analýza dat na PC/2017 jaro")

# Načtení dat ve formátu csv oddělených středníkem
data<-read.csv2(file="04_vzorce.csv",header=FALSE,dec=",",sep=";",skip=5)

# Počet radku a sloupcu souboru
nrow(data)
ncol(data)

# Pojmenovani sloupcu
colnames(data)
colnames(data)<-c("kategorie","cislo","nazev","","mnozstvi")

# Vyber pouze vhodnych sloupcu
data2<-data[,c(1,2,3,5)]

# Datovy typ sloupce
class(data2[,1])
class(data2[,"mnozstvi"])

# Agregacni funkce
mean(data2[,"mnozstvi"])
mean(data2$mnozstvi)

median(data2$mnozstvi)

quantile(data2$mnozstvi,0.25)

min(data2$mnozstvi)
max(data2$mnozstvi)

# Vypocet rozsahu
rozsah=max(data2$mnozstvi)-min(data2$mnozstvi)

# Distribuce do 8 skupin
distribuce<-c(NA,NA,NA,NA,NA,NA,NA,NA)
hranice<-c(0:8)*rozsah/8+min(data2$mnozstvi)

for (i in 1:8) {
  distribuce[i]<-length(data2$mnozstvi[which(data2$mnozstvi> hranice[i]&
                                             data2$mnozstvi<=hranice[i+1])])
}

barplot(distribuce,xlab="Odpad",ylab="Počet")


t<-wilcox.test(c(1,2,3,4,5),c(45,2,1,4,3),method="paired")
t$estimate