setwd("C:/Jirka/Environment/Výuka/E7541 Analýza dat na PC/2023 podzim")
library("xlsx")
data0<-read.xlsx("10_procviceni.xlsx",1)

data1<-data0[46:35255,2:4]

class(data1[,2])
data1[,2]<-as.numeric(data1[,2])

ks.test(data1[which(data1[,3]=="CHRY"),2], 
        "pnorm",
        mean(data1[which(data1[,3]=="CHRY"),2]),
        sd(data1[which(data1[,3]=="CHRY"),2]))

wilcox.test(data1[which(data1[,3]=="CHRY"),2],
            data1[which(data1[,3]=="BAP"),2],
            paired = FALSE)

par(mar=c(2,1,1,1))

boxplot(log10(data1[which(data1[,3]=="CHRY"),2]),
        log10(data1[which(data1[,3]=="BAP"),2]))