multimode (version 1.0)

modetree: Mode tree

Description

This function provides the mode tree.

Usage

modetree(data,bws=NULL,gridsize=NULL,cbw1=NULL,cbw2=NULL,display=TRUE,
logbw=FALSE,addplot=FALSE,xlab=NULL,ylab=NULL,col.lines="black")

Arguments

data

Sample in which the mode tree is computed.

bws

Vector or range of bandwidths. If it is a vector of size two, then it is used a grid of bandwidths between the given values. Default lower bandwidth is twice the distance between the grid points used for estimating the density and upper bandwidth equal to the range of the data. Unless it is specified a vector of size greater than two, the number of bandwidths employed is equal to the second element of gridsize.

gridsize

Number of equally spaced points at which the density is to be estimated (first element) and bandwidths used to compute the mode tree (second element). Default is gridsize=c(512,151).

cbw1

Number of modes for which the first critical bandwidth is calculated. This is the first bandwidth used to compute the mode tree when bws is not specified.

cbw2

Number of modes for which the second critical bandwidth is calculated. This is the last bandwidth used to compute the mode tree when bws is not specified.

display

Logical, if TRUE the mode tree plot is plotted. Default TRUE.

logbw

Logical, if TRUE the plot displays and returns the log10 bandwidths. Default logbw=FALSE.

addplot

Logical, if TRUE the plot is added to the current one. Default is FALSE.

xlab

A title for the x axis. See title.

ylab

A title for the y axis. See title.

col.lines

Color employed in the mode tree. If the length is two, the first color is employed for representing the location of the modes and the second one for the splitting of the modes. Default is col.lines="black".

Value

locations

Estimated mode locations for the bandwidths given in the row names.

bandwidths

Employed bandwidths to compute the mode tree.

Details

The mode tree for the sample given in data is computed. For this calculation, a kernel density estimator with Gaussian kernel and bandwidths bws is used. This plot shows with the continuous lines the estimated mode locations at each bandwidth. The horizontal dashed lines indicate the splitting of a mode in more modes.

The NAs will be automatically removed.

References

Minnotte, M. C. and Scott, D. W. (1993). The mode tree: A tool for visualization of nonparametric density features, Journal of Computational and Graphical Statistics, 2, 51--68.

Examples

# NOT RUN {
#Mode tree using a grid of bandwidths between 0.2 and 0.5
set.seed(2016)
data=rnorm(50)
modetree(data,bws=c(0.2,0.5))
# }
# NOT RUN {
#Estimated locations of the five modes in the claw of Marron and Wand (1992)
library(nor1mix)
set.seed(2016)
n<-200
data<-nor1mix::rnorMix(n,MW.nm10)
#Mode tree between the critical bandwidths for 1 and 8 modes
modetree(data,cbw1=1,cbw2=8)
abline(v=1.5);abline(v=-1.5)
# }