########## Funkce pro konstrukci tabulky souhrnných charakteristik (cz)################################
souhrn.stat=function(values,groups,sign.digits,digits.in) {
# values ... datový vektor (vektor typu Numeric)
# groups ... skupiny (vektor typu Factor)
# sign.digits ... počet des. míst, na něž se zaokrouhlují vypočtené souhrnné statistiky
# digits.in ... počet des. míst složek datového vektoru

# Výpočet charakteristik
rozsah=tapply(values,groups,function(x) sum(!is.na(x)))
minimum=tapply(values,groups,min,na.rm=TRUE)
dolni.kv=tapply(values,groups,quantile,probs=0.25,na.rm=TRUE)
prumer=tapply(values,groups,mean,na.rm=TRUE)
median=tapply(values,groups,quantile,probs=0.5,na.rm=TRUE)
horni.kv=tapply(values,groups,quantile,probs=0.75,na.rm=TRUE)
maximum=tapply(values,groups,max,na.rm=TRUE)
smer.odch=tapply(values,groups,sd,na.rm=TRUE)
var.koef=100*smer.odch/prumer
sikmost=tapply(values,groups,skewness,na.rm=TRUE)
spicatost=tapply(values,groups,kurtosis,na.rm=TRUE)
hradby.min=dolni.kv-1.5*(horni.kv-dolni.kv)
hradby.max=horni.kv+1.5*(horni.kv-dolni.kv)

# Zaokrouhlení charakteristik
dolni.kv=formatC(dolni.kv,digits=sign.digits, format="f")
prumer=formatC(prumer,digits=sign.digits, format="f")
median=formatC(median,digits=sign.digits, format="f")
horni.kv=formatC(horni.kv,digits=sign.digits, format="f")
smer.odch= formatC(smer.odch,digits=sign.digits, format="f")
var.koef= formatC(var.koef,digits=1, format="f")
sikmost= formatC(sikmost,digits=3, format="f")
spicatost= formatC(spicatost,digits=3, format="f")
hradby.min= formatC(hradby.min,digits=digits.in+1, format="f")
hradby.max= formatC(hradby.max,digits=digits.in+1, format="f")

# Sloučení výsledků do tabulky (rbind ... slučuje řádky)
souhrn=rbind(rozsah,minimum,dolni.kv,prumer,median,horni.kv,maximum,smer.odch,var.koef,sikmost,spicatost,hradby.min,hradby.max)
# Definování popisků řádků tabulky
souhrn.popis=c("rozsah","minimum","dolní kvartil","průměr","medián","horní kvartil","maximum","směrodatná odchylka","variační koeficient (%)","šikmost","špičatost","dolní mez (vnitřní hr.)","horní mez (vnitřní hr.)")
rownames(souhrn)=souhrn.popis

# Definování výstupu funkce
return(souhrn)}

##########################################################################################################

# Generování tabulky souhrnných statistik s využitím funkce souhrn.stat
souhrn.po_zak=souhrn.stat(values=data_po_zak$po_zak,groups=data_po_zak$vyrobce,sign.digits=3,digits.in=0)