[R] unimodality test

Martin Maechler maechler at stat.math.ethz.ch
Fri Jul 11 16:04:07 CEST 2003


>>>>> "Jerome" == Jerome Sackur <sackur at heraclite.ens.fr>
>>>>>     on Fri, 11 Jul 2003 12:17:33 +0200 (MET DST) writes:

    Jerome> Dear R users, I am interested in uni- bi-
    Jerome> multimodality tests, for analysing reaction times
    Jerome> data. I was lead to Hartigan's dip test
    Jerome> (Ann. Statistics, 13, 1985, pp. 70-84, Applied
    Jerome> Statistics, 34, 1985, 320-325). Not being a
    Jerome> programmer I am unable to translate the Fortran code
    Jerome> given in ref. 2 into a R function. I'd be glad to
    Jerome> learn that someone already did it, or has devised a
    Jerome> better solution for this kind of problem..

I had got a version with Fortran and S-plus from Dario Ringach
(@ NYU.edu) in 1994 (from what I see) and had worked on it in 2000,
made it into an R package back then.
The reason it hasn't made its way to CRAN was that the Fortran code
(which I f2c'ed to C) still has bugs (leading to segmentation
faults)  that I've not yet found time to debug.

Let me have a look at it before making it available (not on
CRAN) but via FTP as a source package.

Martin Maechler <maechler at stat.math.ethz.ch>	http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum  LEO C16	Leonhardstr. 27
ETH (Federal Inst. Technology)	8092 Zurich	SWITZERLAND
phone: x-41-1-632-3408		fax: ...-1228			<><




More information about the R-help mailing list