sim.cir.2.summary: Simulate 2-occasion CIR multiple times and return a summary

In dill/wisp: Wildlife Simulation Package (wisp)

Description

Perform 2-occasion change-in-ratio (CIR) according to users' specifications, and return mean metrics from the simulation

Usage

sim.cir.2.summary(popn, frac.prehunt, frac.posthunt, frac.harv.male, frac.harv.fem, B)

Arguments

popn	A WiSP population object, perhaps created by make.twosex.pop
frac.prehunt	Proportion [0,1] of the population sampled prior to harvest season to estimate sex composition
frac.posthunt	Proportion [0,1] of the population sampled following the harvest season to estimate sex composition
frac.harv.male	Proportion [0,1] of the males harvested during the hunting season
frac.harv.fem	Proportion [0,1] of the females harvested during the hunting season
B	Number of simulated replications to perform

Details

The population (popn) remains unchanged for all simulations; i.e., the same population is used for all replicates. Hence none of the stochasticity in the results is attributable to changing populations.

Value

A data frame is returned consisting of the following summary statistics:

$point.mean	Mean of the B point estimates generated (after ignoring missing values)
$se.mean	Standard error of the B point estimates generated (after ignoring missing values)
$cv.mean	Coefficient of variation (not percentage) of the B point estimates generated (after ignoring missing values)

Note

Consult sim.cir.2 for the manner in which missing values are generated. The number of replicates requested may not be equal to the number of replicates realized because of the presence of missing values. The utility of this function is to concatenate results from multiple sets of simulations into a form for plotting.

Author(s)

Eric Rexstad, RUWPA [e-mail]

References

Borchers, Buckland, and Zucchini (2002), Estimating animal abundance: closed populations. Chapter 5 http://www.ruwpa.st-and.ac.uk/estimating.abundance

See Also

two.samp.cir, sim.cir.2, make.twosex.pop

Examples

library(wisp)
example <- make.twosex.pop(abund=200, prop.male=0.4)
example.result <- sim.cir.2.summary(popn=example, frac.prehunt=0.2, frac.posthunt=0.4,
                  frac.harv.male=0.6, frac.harv.fem=0.02, B=200)
example.result

dill/wisp documentation built on May 15, 2019, 8:31 a.m.