theta.h

0th

Percentile

Population Parameter THETA using Homozygosity

This function computes the population parameter THETA using the homozygosity (or mean heterozygosity) from gene frequencies.

Keywords
manip, univar
Usage
theta.h(x, standard.error = FALSE)
Arguments
x
a vector or a factor.
standard.error
a logical indicating whether the standard error of the estimated theta should be returned (TRUE), the default being FALSE.
Details

The argument x can be either a factor or a vector. If it is a factor, then it is taken to give the individual alleles in the population. If it is a numeric vector, then its values are taken to be the numbers of each allele in the population. If it is a non-numeric vector, it is a coerced as a factor.

The standard error is computed with an approximation due to Chakraborty and Weiss (1991).

Value

  • a numeric vector of length one with the estimated theta (the default), or of length two if the standard error is returned (standard.error = TRUE).

References

Zouros, E. (1979) Mutation rates, population sizes and amounts of electrophoretic variation at enzyme loci in natural populations. Genetics, 92, 623--646. Chakraborty, R. and Weiss, K. M. (1991) Genetic variation of the mitochondrial DNA genome in American Indians is at mutation-drift equilibrium. American Journal of Human Genetics, 86, 497--506.

See Also

heterozygosity, theta.s, theta.k

Aliases
  • theta.h
Documentation reproduced from package ape, version 2.1-1, License: GPL (>= 2)

Community examples

Looks like there are no examples yet.