[R] unimodality test

Martin Maechler maechler at stat.math.ethz.ch
Mon Jul 14 12:17:43 CEST 2003


replying to my own posting of last Friday:
>>>>> "MM" == Martin Maechler <maechler at stat.math.ethz.ch> writes:

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

    MM> I had got a version with Fortran and S-plus from Dario
    MM> Ringach (@ NYU.edu) in 1994 (from what I see) and had
    MM> worked on it in 2000, made it into an R package back then.
(yes)

    MM> The reason it hasn't made its way to CRAN was
    MM> that the Fortran code (which I f2c'ed to C) still has
    MM> bugs (leading to segmentation faults) that I've not yet
    MM> found time to debug.

not quite true.

    MM> Let me have a look at it before making it available  <....>

I've put a new source package diptest (0.9-1) on CRAN's incoming.
For the interested ones, there's a version of it also available
as
	ftp://stat.ethz.ch/U/maechler/R/diptest_0.9-1.tar.gz

	(yes, only for those who can install packages from source!)

Note that because of vacation and conferences, it may take a
bit longer before this will be available from CRAN and its mirrors 
(and even a bit longer before "pre-compiled" availability for Windows).

Regards,

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