setwd("C:/Jirka/Environment/Výuka/E7541 Analýza dat na PC/2022 podzim")

#install.packages("xlsx")
library("xlsx")

# Kontingenční tabulka - kongenery PCB v regionech UNEP
data1.0<-read.xlsx("12_kontingencni_tabulky_zadani.xlsx",2,startRow=12)
data1.1<-data1.0[1:1190,3:11]
colnames(data1.1)<-c("region","stat","lokalita","matrice","submatrice","skupina","latka","koncentrace","jednotka" )
data1<-data1.1

tabulka1<-table(data1$region,data1$latka)
chisq.test(tabulka1)

# alternativni použití Fisherova exaktního testu
fisher.test(tabulka1) # vrací chybu
fisher.test(tabulka1,workspace=1e6)



## X2 test závislosti nemocnosti na pohlaví
data2.0<-read.xlsx("12_kontingencni_tabulky_zadani.xlsx",3,startRow=10)
data2.1<-data2.0[1:200,3:4]
colnames(data2.1)<-c("nemoc","pohlavi")
data2<-data2.1

tabulka2<-table(data2$nemoc,data2$pohlavi)
chisq.test(tabulka2)

# alternativni použití Fisherova exaktního testu
fisher.test(tabulka2)