simmixstock1 {mixstock}R Documentation

Simulate marker frequencies and distributions in a mixed stock

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


[Package mixstock version 0.9 Index]