# nainstalujte si program R: https://cran.r-project.org/
# nainstalujte si RStudio: https://rstudio.com/products/rstudio/download/#download

# Packages: nainstalujte balicek HardyWeiberg

library(HardyWeinberg)      # načtěte balíček HardyWeinberg
x <- c(SS=66,SF=20,FF=14)   # vytvoření dat - genotypy a jejich počty ve vektoru
HW.test <- HWChisq(x, cc=0, verbose=TRUE) # chí kvadrát test HWE
HW.test <- HWExact(x, verbose=TRUE)       # exaktní test HWE
HWAlltests(x,verbose=TRUE)                # spustí všechny nejčastěji používané testy pro HWE
HWLratio(x,verbose=TRUE, x.linked = FALSE)  # 	Likelihood ratio test pro HWE
HWPerm(x, nperm = 20000, verbose = TRUE, x.linked = FALSE) # permutační test pro HWE