Description Usage Arguments Details Value Caution Note Author(s) References See Also Examples
When designing a CIR experiment, effort needs to be allocated not only to measuring the magnitude of the harvest, but also to estimating the population sex ratio prior to, and following the harvest. This function examines combinations of different sampling intensities pre- and post-harvest.
1 | make.2d.cir(pop, top.samp = 0.5, samp.incr = 0.1, male.harv = 0.5, fem.harv = 0.05, reps = 10)
|
pop |
Population object being harvested and sampled (likely created by |
top.samp |
Greatest proportion of the population sampled pre- and post-harvest for estimation of sex ratio |
samp.incr |
Step size to be used to iterate through the sampling proportions (both pre- and post-haravest) |
male.harv |
Proportion of males in the population removed during the harvest |
fem.harv |
Proportion of females in the population removed during the harvest |
reps |
Number of simulations to perform at each combination of pre- and post-harvest sampling values. |
This function is a wrapper, that repeatedly calls sim.cir.2.summary
. Exercise some caution
in the number of replicates requested; this can considerably increase the computational burden. Similarly, having
a samp.incr
vanishingly small, will also lead to long waits for results. Reasonable results can be derived
from reps=500
and samp.incr=0.02
.
The resulting data frame consists of
pre-harvest sampling fraction |
Sampling fraction for pre-harvest |
post-harvest sampling fraction |
Sampling fraction for post-harvest |
mean point estimate for sampling combination |
Mean of the |
mean standard error for sampling combination |
Mean of the |
mean coefficient of variation for sampling combination |
Mean of the |
The resulting data frame is comprised of means, calculated after missing values (inadmissable values) have been eliminated. Because distributions of estimates from CIR estimators can be quite skewed, it is worth considering the possibility of examining medians.
To produce contour or heat plots of the resulting data frame, the column of interest must be converted back
to a matrix, e.g., result.mat <- matrix(result$point, nrow=5, ncol=5)
.
Eric Rexstad, RUWPA ericr@mcs.st-and.ac.uk
Borchers, Buckland, and Zucchini (2002), Estimating animal abundance: closed populations. Chapter 5 http://www.ruwpa.st-and.ac.uk/estimating.abundance
sim.cir.2.summary
, two.samp.cir
, and make.twosex.pop
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | example.pop <- make.twosex.pop(abund=300, prop.male=0.4)
example.result <- make.2d.cir(pop=example.pop, top.samp=0.5, samp.incr=0.1, male.harv=0.5,
fem.harv=0.05, reps=100)
# Note, interest here is actually in precision of estimate, not in bias, hence
# matrix is created of CV rather than of point
example.result.mat <- matrix(example.result$cv, nrow=5, ncol=5) # note 5 unique values for pre- and post-
contour(x=seq(0.1, 0.5, by=0.1), y=seq(0.1, 0.5, by=0.1), z=example.result.mat,
main="Example output from function make.2d.cir")
#
# Dont run
# unless you have installed the rgl library
#
library(rgl)
mat2 <- pmin(example.result.mat, 2) # Truncate maximum value to a CV of 2
collut <- terrain.colors(200) # make the plot more colorful
col2 <- collut[mat2*100]
rgl.surface(x=10*(seq(.1,.5, by=.1)),z=10*(seq(.1,.5,by=.1)),mat2,color=col2)
rgl.bbox()
rgl.texts(x=2.5,z=2.5,y=1,"CV for values of p-pre and p-post", justify="center") # label the surface
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.