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