hierfstat (version 0.04-14)

g.stats: Calculates likelihood-ratio G-statistic on contingency table

Description

Calculates the likelihood ratio G-statistic on a contingency table of alleles at one locus X sampling unit. The sampling unit could be any hierarchical level

Usage

g.stats(data,diploid=TRUE)

Arguments

data

a two-column data frame. The first column contains the sampling unit, the second the genotypes

diploid

Whether the data are from diploid (default) organisms

Value

• obs
  Observed contingency table
• exp
  Expected number of allelic observations
• X.squared
  The chi-squared statistics, $\sum\frac{(O-E)^2}{E}$
• g.stats
  The likelihood ratio statistics, $2 \sum(O \log(\frac{O}{E}))$

References

Goudet J., Raymond, M., DeMeeus, T. and Rousset F. (1996) Testing differentiation in diploid populations. Genetics. 144: 1933-1940 Goudet J. (2005). Hierfstat, a package for R to compute and test variance components and F-statistics. Molecular Ecology Notes. 5:184-186 Petit E., Balloux F. and Goudet J.(2001) Sex-biased dispersal in a migratory bat: A characterization using sex-specific demographic parameters. Evolution 55: 635-640.

See Also

g.stats.glob.

Examples

data(gtrunchier)
attach(gtrunchier)
g.stats(data.frame(Patch,L21.V))