R help - Jul 2003 - unimodality test

Dear R users,

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



Thanks a lot in advance,

J. Sackur
Inserm U562
Orsay, France

>>>>> "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			<><