activity (version 1.0)

compareCkern: Compare circular distributions.

Description

Test probability that two sets of circular observations come from the same distribution.

Usage

compareCkern(y1, y2, reps = 1000)

Arguments

y1,y2

Numeric vectors of radian data.

reps

Number of bootstrap iterations.

Value

A named 2-element vector: Overlap = observed overlap index; p = probability observed index arose by chance.

Details

Bootstrap test calculates overlap index 2 (see references) for the observed data samples, then generates a null distribution of overlap indices using data sampled randomly with replacement from the combined data. This randomised distribution is then used to estimate the probability that the observed overlap arose by chance.

References

Ridout, M.S. & Linkie, M. (2009) Estimating overlap of daily activity patterns from camera trap data. Journal of Agricultural Biological and Environmental Statistics, 14, 322-337.

See Also

overlapEst

Examples

#Example with bootstrap reps limited to speed up
data(BCItime)
tPaca <- 2*pi*BCItime$time[BCItime$species=="paca"]
tRat <- 2*pi*BCItime$time[BCItime$species=="rat"]
compareCkern(tPaca,tRat,reps=10)