Description Usage Arguments Details Value Examples

View source: R/activity_code_1_1.r

Fits a circular kernel density to radian time-of-day data and estimates activity level from this distribution. Optionally bootstraps the distribution, in which case SEs and confidence limits are also stored for activity level and PDF.

1 2 |

`dat` |
A numeric vector of radian time-of-day data. |

`wt` |
A numeric vector of weights for each |

`reps` |
Number of boostrap iterations to perform. Ignored if sample=="none". |

`bw` |
Numeric value for kernel bandwidth. If NULL, calculated internally. |

`adj` |
Numeric bandwidth adjustment multiplier. |

`sample` |
Character string defining sampling method for bootstrapping errors (see details). |

`show` |
Logical whether or not to show a progress bar while bootstrapping. |

The bandwidth adjustment multiplier `adj`

is provided to allow
exploration of the effect of adjusting the internally calculated bandwidth on
accuracy of activity level estimates. The alternative bootstrapping methods
defined by `sample`

are:
`data`

: sample from the data;
`model`

: sample from the fitted probability density distribution;
`none`

: no bootstrapping.
Confidence interval coverage seems to be better at large sample size
(greater than 100-200) using `"model"`

, but better at small sample size
when using `"data"`

. The reason for this needs further investigation.

An object of type `actmod`

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | ```
#Fit without confidence limits
data(BCItime)
tdat <- 2*pi*BCItime$time[BCItime$species=="ocelot"]
mod1 <- fitact(tdat, sample="none")
plot(mod1)
#Fit with confidence limits (limited reps to speed up)
mod2 <- fitact(tdat, reps=10)
plot(mod2)
#Fit weighted function to correct for detection radius 1.21 times higher
#by day than by night, assuming day between pi/2 (6am) and pi*2/3 (6pm)
weight <- 1/ifelse(tdat>pi/2 & tdat<pi*3/2, 1.2, 1)
mod3 <- fitact(tdat, wt=weight, sample="none")
plot(mod3)
#Overplot unweighted version for comparison
plot(mod1, add=TRUE, lcol=3)
``` |

activity documentation built on May 30, 2017, 8:18 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.