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

In hierfstat: Estimation and Tests of Hierarchical F-Statistics

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 (O-E)^2/E
g.stats	The likelihood ratio statistics, 2 Sum(O Log(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))

hierfstat documentation built on May 6, 2022, 1:05 a.m.