mssumstats: Calculates several summary statistics from microsats

In mastoffel/sealABC: useful functions for the abc pinniped analysis

Description

Calculates several summary statistics from microsats

Usage

mssumstats(data, by_pop = NULL, start_geno = NULL, mratio = c("strict",
  "loose"), rarefaction = FALSE, nresamp = 1000, nind = NULL,
  nloc = NULL)

Arguments

data	microsatellites, whereby every locus is represented by two adjacent columns. Or a gtypes object from strataG
by_pop	name of population variable. If specified, all summary statistics will be calculated within populations
start_geno	integer, specifying the first column with genotypes. If now specified, all column are expected to be genotypes
mratio	defaults to "strict". if "loose", the mratio will be calculated differently, see ?m_ratio
rarefaction	if TRUE, calculates mean and sd number of alleles as mean of n_samp individuals and n_loc loci over nresamp bootstraps
nloc	number of loci to subsample
nsamp	number of samples to subssample
nboot	number of bootstraps

Examples

data(fur_seal)
out <- mssumstats(fur_seal, by_pop = "pop", start_geno = 4, mratio = "loose", rarefaction = TRUE,
nresamp = 10, nind = 20, nloc = 5)

mastoffel/sealABC documentation built on May 21, 2019, 12:43 p.m.