locus_table: Create a table of summary statistics per locus.

Description

Create a table of summary statistics per locus.

Usage

locus_table(x, index = "simpson", lev = "allele", population = "ALL",
  information = TRUE)

Arguments

a genind or genclone object.

index

Which diversity index to use. Choices are

"simpson"(Default) to give Simpson's index
"shannon"to give the Shannon-Wiener index
"invsimpson"to give the Inverse Simpson's index aka the Stoddar

lev

At what level do you want to analyze diversity? Choices are "allele" (Default) or "genotype".

population

Select the populations to be analyzed. This is the parameter sublist passed on to the function popsub. Defaults to "ALL".

information

When TRUE (Default), this will print out a header of information to the R console.

Value

a table with 4 columns indicating the Number of alleles/genotypes observed, Diversity index chosen, Nei's 1978 expected heterozygosity, and Evenness.

References

Jari Oksanen, F. Guillaume Blanchet, Roeland Kindt, Pierre Legendre, Peter R. Minchin, R. B. O'Hara, Gavin L. Simpson, Peter Solymos, M. Henry H. Stevens, and Helene Wagner. vegan: Community Ecology Package, 2012. R package version 2.0-5.

Niklaus J. Gr"unwald, Stephen B. Goodwin, Michael G. Milgroom, and William E. Fry. Analysis of genotypic diversity data for populations of microorganisms. Phytopathology, 93(6):738-46, 2003

J.A. Ludwig and J.F. Reynolds. Statistical Ecology. A Primer on Methods and Computing. New York USA: John Wiley and Sons, 1988.

E.C. Pielou. Ecological Diversity. Wiley, 1975.

J.A. Stoddart and J.F. Taylor. Genotypic diversity: estimation and prediction in samples. Genetics, 118(4):705-11, 1988.

Masatoshi Nei. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics, 89(3):583-590, 1978.

Claude Elwood Shannon. A mathematical theory of communication. Bell Systems Technical Journal, 27:379-423,623-656, 1948

Examples

Run this code

# Analyze locus statistics for the North American population of P. infestans.
data(Pinf)
locus_table(Pinf, population = "North America")