https://github.com/cran/HardyWeinberg
Raw File
Tip revision: 59bc105de5b1eda230f8d1ce6ac1c7e09cfe1ad3 authored by Jan Graffelman on 06 April 2024, 09:53:07 UTC
version 1.7.8
Tip revision: 59bc105
HWChisqMat.R
HWChisqMat <- function(X,...) {
# applys chi-square test for HWE to all rows of X.
pvalvec <- NULL
chisqvec <- NULL
Dvec <- NULL
for (i in 1:nrow(X)) {
   out <- HWChisq(X[i,],...)
   pvalvec <- c(pvalvec,out$pval)
   chisqvec <- c(chisqvec,out$chisq)
   Dvec <- c(Dvec,out$D)   
}
return(list(pvalvec=pvalvec,chisqvec=chisqvec,Dvec=Dvec))
}

back to top