1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/hrm.R
\name{hrm}
\alias{hrm}
\alias{HRM}
\title{Bickel's half-range mode estimator}
\usage{
hrm(x, bw = NULL, ...)
}
\arguments{
\item{x}{numeric. Vector of observations.}
\item{bw}{numeric. The bandwidth to be used. Should belong to (0, 1].
This gives the fraction of the observations to consider at
each step of the iterative algorithm.}
\item{...}{Additional arguments.}
}
\value{
A numeric value is returned, the mode estimate.
}
\description{
SINCE THIS FUNCTION USED TO DEPEND ON THE BIOCONDUCTOR PACKAGE 'GENEFILTER',
IT IS CURRENTLY DEFUNCT.
This function computes Bickel's half range mode estimator
described in Bickel (2002). It is a wrapper around the function
\code{half.range.mode} from package \pkg{genefilter}.
}
\details{
The mode estimator is computed by iteratively identifying
densest half ranges. A densest half range is an interval
whose width equals half the current range, and which
contains the maximal number of observations.
The subset of observations falling in the selected
densest half range is then used to compute a new range,
and the procedure is iterated.
}
\note{
The user may call \code{hrm} through
\code{mlv(x, method = "hrm", bw, ...)}.
}
\examples{
\dontrun{
# Unimodal distribution
x <- rgamma(1000, shape = 31.9)
## True mode
gammaMode(shape = 31.9)
## Estimate of the mode
hrm(x, bw = 0.4)
mlv(x, method = "hrm", bw = 0.4)
}
}
\references{
\itemize{
\item Bickel D.R. (2002).
Robust estimators of the mode and skewness of continuous data.
\emph{Computational Statistics and Data Analysis}, \bold{39}:153-163.
\item Hedges S.B. and Shah P. (2003).
Comparison of mode estimation methods and application in molecular clock analysis.
\emph{BMC Bioinformatics}, \bold{4}:31-41.
\item Bickel D.R. and Fruehwirth R. (2006).
On a Fast, Robust Estimator of the Mode:
Comparisons to Other Robust Estimators with Applications.
\emph{Computational Statistics and Data Analysis}, \bold{50}(12):3500-3530.
}
}
\seealso{
\code{\link[modeest]{mlv}()} for general mode estimation;
\code{\link[modeest]{hsm}()} for the half sample mode;
\code{\link[modeest]{venter}()} for the Venter mode estimate.
}
\author{
The C and R code are due to Richard Bourgon \email{bourgon@stat.berkeley.edu},
see package \pkg{genefilter}. The algorithm is described in Bickel (2002).
}
|
