HWf: Computation of inbreeding coefficient
In HardyWeinberg: Statistical Tests and Graphics for Hardy-Weinberg Equilibrium
HWf R Documentation
Computation of inbreeding coefficient
Description
HWf computes the inbreeding coefficient for sample
of genotype counts, or a matrix of genotype counts.
Usage
HWf(X)
Arguments
X
a vector or matrix of genotype counts (AA, AB, BB)
Details
For monomorphic markers a warning is issued, and the estimate for the
inbreeding coefficient is NaN.
Value
Returns a single inbreeding coefficient (intraclass correlation coefficient), if X is a single sample, or a vector of inbreeding coefficients, if X is a matrix with genotype counts.
Author(s)
Jan Graffelman jan.graffelman@upc.edu
References
Crow, J. F. and Kimura, M. (1970) An introduction to population genetics theory.
Harper & Row, publishers, New York
See Also
HWChisq
Examples
#
# A single sample
#
x <- c(MM=298,MN=489,NN=213)
fhat <- HWf(x)
fhat
#
# Multiple samples
#
set.seed(123)
X <- HWData(nm=100,n=1000)
fhat <- HWf(X)
HardyWeinberg documentation built on May 29, 2024, 6:17 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
| HWf | R Documentation |
Computation of inbreeding coefficient
Description
HWf computes the inbreeding coefficient for sample
of genotype counts, or a matrix of genotype counts.
Usage
HWf(X)
Arguments
X |
a vector or matrix of genotype counts (AA, AB, BB) |
Details
For monomorphic markers a warning is issued, and the estimate for the inbreeding coefficient is NaN.
Value
Returns a single inbreeding coefficient (intraclass correlation coefficient), if X is a single sample, or a vector of inbreeding coefficients, if X is a matrix with genotype counts.
Author(s)
Jan Graffelman jan.graffelman@upc.edu
References
Crow, J. F. and Kimura, M. (1970) An introduction to population genetics theory. Harper & Row, publishers, New York
See Also
HWChisq
Examples
#
# A single sample
#
x <- c(MM=298,MN=489,NN=213)
fhat <- HWf(x)
fhat
#
# Multiple samples
#
set.seed(123)
X <- HWData(nm=100,n=1000)
fhat <- HWf(X)
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.