simmixstock0: Simulate a sample from sources and mixed stock

In bbolker/mixstock: Functions for Mixed Stock Analysis

Simulate a sample from sources and mixed stock

Description

Given information on marker frequencies in sources, the contributions of each source to a mixed stock, and sample sizes, generates multinomial samples of markers from sources and the mixed stock.

Usage

simmixstock0(sourcefreq, input.freq, sourcesampsize, mixsampsize, rseed=NULL)

Arguments

sourcefreq	Matrix (markers=rows, source=cols) of true marker frequencies in sources
input.freq	True contributions from each source to the mixed stock
sourcesampsize	Sample size from (each) source
mixsampsize	Sample size from mixed stock
rseed	Random-number seed

Value

Multinomial sample.

mixsamp	Sample from mixed stock (vector)
sourcesamp	Sample from sources (matrix, markers=rows, source=cols)

...

Author(s)

Ben Bolker

Examples

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)
x <- simmixstock0(true.freq,true.contrib,sourcesampsize=50,mixsampsize=100,rseed=1004)

bbolker/mixstock documentation built on July 23, 2024, 12:18 p.m.