powered by
HWf computes the inbreeding coefficient for a sample of genotypes.
HWf
HWf(X)
a vector of genotype counts (AA, AB, BB)
Returns the inbreeding coefficient (intraclass correlation coefficient)
For monomorphic markers a warning is issued, and the estimate for the inbreeding coefficient is set to zero.
Crow, J. F. and Kimura, M. (1970) An introduction to population genetics theory. Harper \& Row, publishers, New York
HWChisq
# NOT RUN { x <- c(MM=298,MN=489,NN=213) f <- HWf(x) cat("Inbreeding coefficient: ",f,"\n") # }
Run the code above in your browser using DataLab