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.2 Index]