Modal regions
Evaluate the modal regions for a data set. Only the version for circular data is implemented.
modal.region(x, ...)
## Default S3 method:
modal.region(x, ...)
## S3 method for class 'circular'
modal.region(x, z=NULL, q=0.95, bw, adjust = 1,
type = c("K", "L"), kernel = c("vonmises", "wrappednormal"),
na.rm = FALSE, step=0.01, eps.lower=10^(-4), eps.upper=10^(-4), ...)x |
numeric or an object of class |
z |
numeric or object of class |
q |
numeric in the interval [0,1]. The quantile of the modal region. |
bw |
the smoothing bandwidth to be used. When the |
adjust |
the bandwidth used is actually |
type |
Not Yet Used. |
kernel |
a character string giving the smoothing kernel to be
used. This must be one of |
na.rm |
logical; if |
step |
numeric. Used in the construction of the regular grid |
eps.lower,eps.upper |
the cut point in the density is searched in the interval [min(density)*(1+eps.lower),max(density)*(1-eps.upper)]. |
... |
further arguments passed to the next methods. |
Only the version for circular data is actually implemented.
A list of class modal.region.circular with the following elements
zeros |
extremes of modal regions, possible as a matrix |
areas |
a list with two components: |
density |
the object from function |
q |
the modal region order as in input. |
level |
the cut point at the density scale. |
Claudio Agostinelli
L.G.R. Oliveira-Santos, C.A. Zucco and C. Agostinelli (2013) Using conditional circular kernel density functions to test hypotheses on animal circadian activity. Animal Behaviour, 85(1) 269-280.
x <- rvonmises(100, circular(pi), 10) res <- modal.region(x, bw=50) plot(res)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.