g.stats {hierfstat}

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}))

Author(s)

Jerome Goudet, DEE, UNIL, CH-1015 Lausanne Switzerland

[e-mail]

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))