LPCM (version 0.41-6)

coverage: Coverage and self-coverage plots.

Description

These functions compute coverages (for any principal object), and self-coverages (only for local principal curves, these may be used for bandwidth selection).

Usage

coverage.raw(X, vec, tau, weights=1, plot.type="p", print=FALSE,
      label=NULL,...)
coverage(X, vec, taumin=0.02, taumax, gridsize=25, weights=1,
      plot.type="o", print=TRUE,...)
lpc.coverage(object, taumin=0.02, taumax, gridsize=25, quick=TRUE,
      plot.type="o", print=TRUE, ...)
lpc.self.coverage(X,  taumin=0.02, taumax=0.5,   gridsize=25, x0, mult=1, 
     way = "two", scaled=TRUE,  weights=1, pen=2, depth=1,
     control=lpc.control(boundary=0, cross=FALSE),   quick=TRUE,
     plot.type="o", print=TRUE, ... )
select.self.coverage(self, sens=0.02, from, plot.type="o", auto=FALSE)

Arguments

X

a $N \times d$ data matrix.

object

An object of type lpc or lpc.spline.

vec

A matrix with $d$ columns. The rows contain the points which make up the fitted object.

tau

tube size.

taumin

Minimal tube size.

taumax

Maximal tube size.

weights

An optional vector of weights. If weights are specified, then the coverage is the weighted mean of the indicator functions for falling within the tube. The function lpc.coverage does not have a weights argument, a

label

Experimental option; don't use.

gridsize

The number of different tube sizes to consider.

quick

If TRUE, an approximate coverage curve is provided by computing distances between data points and the curve through the closest local centers or mass; whereas with FALSE we use the distances of the points when projected orthogonally onto the splin

self

An object of class self, or a matrix with two colums providing a self-coverage curve.

sens

Sensitivity parameter: only bandwidths corresponding to curvature values of the self-coverage curve, which are less than the negative value of sens, are selected.

from

Minimum coverage for bandwidth selection. Default: 1/3 for clustering, 2/3 for principal curves.

auto

Unimplemented. If TRUE, the sensitivity is decreased until a bandwidth is selected.

plot.type

If set to 0, no plotted output is given. Otherwise, an appropriate plot is provided, using the plotting type as specified.

print

If TRUE, coverage values are printed on the screen as soon as computed. This is quite helpful especially if gridsize is large.

x0, way, scaled, pen, mult, depth, control

LPC parameters as outlined in lpc and lpc.control.

...

Optional graphical parameters passed to the corresponding plotting functions.

Value

• A list of items, and a plot (unless plot.type=0).

  For function lpc.self.coverage, the item $select recommends suitable bandwidths for the use in lpc. These correspond to points of strong negative curvature (implemented via second differences) of the self-coverage curve.

Details

See Einbeck, Tutz, and Evers (2005). This paper uses `quick' coverage curves.

Functions coverage.raw and coverage can be used for any object fitted by an unsupervised learning technique (for instance, HS principal curves, or even clustering algorithms), while the functions prefixing with lpc. can only be used for local principal curves. The function lpc.coverage is a wrapper around coverage which takes directly a fitted lpc object, rather than a data matrix.

Function select.self.coverage is called by lpc.self.coverage. It extracts suitable bandwidths from the self-coverage curve, and produces a plot. The function can be called directly by the user if the sensitivity sens or the minimum coverage from are to be modified.

References

Einbeck, J., Tutz, G., & Evers, L. (2005). Local principal curves. Statistics and Computing 15, 301-313.

See Also

lpc

Examples

data(gvessel)
gvessel.lpc <- lpc(gvessel[,c(2,4,5)], h=0.11,  x0=c(35, 1870, 6.3))
lpc.coverage(gvessel.lpc, gridsize=10, print=FALSE)
lpc.self.coverage(gvessel[,c(2,4,5)], x0=c(35, 1870, 6.3), print=FALSE)

data(calspeedflow)
fitms <- ms(calspeedflow[,3:4], h=0.1)
coverage(fitms$data, fitms$cluster.center)