simmixstock1: Simulate marker frequencies and distributions in a mixed...

In mixstock: functions for mixed stock analysis

Description

Functions for simulating marker frequency distributions and samples in source and mixed populations

Usage

simmixstock1(sampsize = NULL, true.freq = matrix(c(0.65, 0.33, 0.01, 
    0.01, 0.33, 0.65, 0.01, 0.01), ncol = 2), true.contrib = c(0.9, 
    0.1), boot = FALSE, param = FALSE, data = NULL, rseed = 1004, 
    nboot = 1000, chainlen = NULL, ests = c("cmlboot.nonpar", 
        "cmlboot.par", "umlboot.nonpar", "umlboot.par", "mcmc"), 
    verbose = FALSE, contrib.only = FALSE)
sim.mark.freq(H,R,g.mark,g.source)

Arguments

sampsize	total sampsize: half from mixed population, (1/(2R)) from each source
true.freq	matrix of marker frequencies in sources
true.contrib	contributions from each source to source population
boot	bootstrap existing data?
param	parametric bootstrap?
data	original data set to bootstrap
nboot	number of bootstrap samples
chainlen	chain length for MCMC
ests	list of estimates to produce (parametric or nonparametric bootstrap for CML or UML estimation, MCMC)
H	number of markers
R	number of sources
g.mark	geometric distribution parameter for marker frequency
g.source	geometric distribution parameter for source contribution
rseed	random number seed
contrib.only	save only source contributions in MCMC chain results?
verbose	verbose output?

Value

sim.mark.freq just returns an HxR matrix of marker simmixstock1 returns a list with a genboot result for each type of estimate requested;

Author(s)

Ben Bolker

See Also

genboot, mysumvec

mixstock documentation built on May 2, 2019, 6:48 p.m.