spatstat.explore (version 3.8-2)

circdensity: Density Estimation for Circular Data

Description

Computes a kernel smoothed estimate of the probability density for angular data.

Usage

circdensity(x, sigma = "nrd0", ...,
               bw = NULL,
               weights=NULL,
               unit = c("degree", "radian", "hour", "minute", "other"),
               fullcircle=NULL)

Value

An object of class "density" (produced by

density.default) which can be plotted by plot or by rose.density.

The result also contains information about the arguments

unit and fullcircle which can be used by

rose.density.

Arguments

x

Numeric vector, containing angular data.

sigma

Smoothing bandwidth, or bandwidth selection rule, passed to density.default.

bw

Alternative to sigma for consistency with other functions.

...

Additional arguments passed to density.default, such as kernel and weights.

weights

Optional numeric vector of weights for the data in x.

unit

The unit of angle in which x is expressed. Character string (partially matched).

fullcircle

The period (the value of x that corresponds to one full circle) if unit="other".

Author

Adrian Baddeley , Rolf Turner and Ege Rubak .

Details

The angular values x are smoothed using (by default) the wrapped Gaussian kernel with standard deviation sigma.

See Also

density.default), rose.

Examples

  ang <- runif(1000, max=360)
  rose(circdensity(ang, 12))