diptest (version 0.25-2)

dip: Compute Hartigan's Dip Test Statistic for Unimodality

Description

Computes Hartigan's dip test statistic for testing unimodality, and additionally the modal interval.

Usage

dip(x, full.result = FALSE, debug = FALSE)

Arguments

x
numeric; the data.
full.result
logical; if TRUE returns the full result list, see below.
debug
logical; if true, some tracing information is printed (from the C routine).

Value

  • depending on full.result either a number, the dip statistic, or a list with components
  • xthe sorted unname()d data.
  • nlength(x).
  • dipthe dip statistic
  • lo.hiindices into x for lower and higher end of modal interval
  • xl, xulower and upper end of modal interval
  • gcm, lcm(last used) indices for greatest convex minorant and the least concave majorant.
  • mn, mjindex vectors of length n for the GC minorant and the LC majorant respectively.

References

P. M. Hartigan (1985) Computation of the Dip Statistic to Test for Unimodality; Applied Statistics (JRSS C) 34, 320--325. Corresponding (buggy!) Fortran code of AS 217 available from Statlib, http://lib.stat.cmu.edu/apstat/217

J. A. Hartigan and P. M. Hartigan (1985) The Dip Test of Unimodality; Annals of Statistics 13, 70--84.

See Also

isoreg for isotonic regression. qDiptab on how to get approximate P-values.

Examples

Run this code
data(statfaculty)
plot(density(statfaculty))
dip(statfaculty)

str(dip(statfaculty, full = TRUE, debug = TRUE))

data(faithful)
fE <- faithful$eruptions
plot(density(fE))
str(dip(fE, full = TRUE, debug = TRUE))

data(precip)
plot(density(precip))
str(dip(precip, full = TRUE, debug = TRUE))

Run the code above in your browser using DataCamp Workspace