Learn R Programming

mrds (version 2.3.0)

varn: Compute empirical variance of encounter rate

Description

Computes one of a series of possible variance estimates for the observed encounter rate for a set of sample measurements (e.g., line lengths) and number of observations per sample.

Usage

varn(lvec,nvec,type)

covn(lvec, groups1, groups2, type)

Value

Variance of encounter rate as defined by arguments

Arguments

lvec

vector of sample measurements (e.g., line lengths)

nvec

vector of number observed

type

choice of variance estimator to use for encounter rate

groups1

vector of number of groups observed

groups2

vector of number of individuals observed

Author

Jeff Laake, David L Miller

Details

The choice of type follows the notation of Fewster et al. (2009) in that there are 8 choices of encounter rate variance that can be computed for lines and one for points:

R2

random line placement with unequal line lengths (design-assisted estimator)

R3

random line placement, model-assisted estimator, based on true contagion process

R4

random line placement, model-assisted estimator, based on apparent contagion process

S1

systematic line placement, post-stratification with no strata overlap

S2

systematic line placement, post-stratification with no strata overlap, variances weighted by line length per stratum

O1

systematic line placement, post-stratification with overlapping strata (akin to S1)

O2

systematic line placement, post-stratification with overlapping strata (weighted by line length per stratum, akin to S2)

O3

systematic line placement, post-stratification with overlapping strata, model-assisted estimator with trend in encounter rate with line length

P2

random point placement, potentially unequal number of visits per point, design-based estimator

P3

random point placement, potentially unequal number of visits per point, model-based estimator

Default value is "R2", shown in Fewster et al. (2009) to have good performance for completely random designs for lines. For systematic parallel line transect designs, Fewster et al. recommend "O2". For point transects the default is "P2" (but "P3" is also available).

For the systematic estimators, pairs are assigned in the order they are given in the lengths and groups vectors.

References

Fewster, R.M., S.T. Buckland, K.P. Burnham, D.L. Borchers, P.E. Jupp, J.L. Laake and L. Thomas. 2009. Estimating the encounter rate variance in distance sampling. Biometrics 65: 225-236.