lines.density.circular {circular}R Documentation

Add a Plot for Kernel Density Estimation for Circular Data

Description

The lines add a plot for density.circular objects.

Usage

## S3 method for class 'density.circular':
lines(x, type = "l", zero.line = TRUE, points.plot = FALSE, points.col = 1, points.pch = 1, points.cex = 1, plot.type = c("circle", "line"), bins = NULL, shrink = 1, tcl = 0.025, sep = 0.025, join = TRUE, nosort = FALSE, plot.info = NULL, zero = NULL, rotation = NULL, ...)

Arguments

x an object of class density.circular.
type plotting parameter with useful default.
zero.line logical; if TRUE, add a base line at y = 0. Used when plot.type is "line".
points.plot logical; if TRUE original data are added to the plot.
points.col, points.pch, points.cex parameters used to draw the points.
plot.type type of the plot.
bins number of ticks to plot.
shrink parameter that controls the size of the plotted circle. Default is 1. Larger values shrink the circle, while smaller values enlarge the circle.
tcl length of the ticks.
sep constant used to specify the distance between stacked points. Default is 0.025; smaller values will create smaller spaces.
join logical: should the first and the last point joined.
nosort logical: should the data sort before plotting. Defaults is to sort.
plot.info an object from plot.circular that contains information on the zero, the rotation and next.points.
zero the zero of the plot. Ignored if plot.info is provided.
rotation the rotation of the plot. Ignored if plot.info is provided.
... futher parameters passed to lines.default.

Value

A list with information on the plot: zero, rotation and next.points.

Author(s)

Claudio Agostinelli

See Also

density.circular and plot.density.circular

Examples

set.seed(1234)
x <- rvonmises(n=100, mu=circular(pi), kappa=2)
y <- rvonmises(n=100, mu=circular(pi/2), kappa=2)
resx <- density(x, bw=25)
res <- plot(resx, points.plot=TRUE, xlim=c(-1.5,1), ylim=c(-1.1, 1.5))
resy <- density(y, bw=25)
lines(resy, points.plot=TRUE, col=2, points.col=2, plot.info=res)

[Package circular version 0.3-8 Index]