Has anyone seen an R function or code to perform the Two-Dimenstional Kolmogorov-Smirnov (2DKS) test described in this reference ...

James E. Garvey, Elizabeth A. Marschall, Russell A. Wright (1998) FROM STAR CHARTS TO STONEFLIES: DETECTING RELATIONSHIPS IN CONTINUOUS BIVARIATE DATA. Ecology: Vol. 79, No. 2, pp. 442-447.

There have been (at least) three previous questions (basically one thread) on the R-help mailing list but either without answers or to links that are no longer live ..

http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg22028.html
http://tolstoy.newcastle.edu.au/R/help/04/06/1060.html
http://r.789695.n4.nabble.com/2D-kolmogorov-td849793.html#a849793

Apparently there is a routine in numerical recipes (http://www.fizyka.umk.pl/nrbook/c14-7.pdf) and a compiled version on the lead authors website but I have been unable to find an implementation in R.

Does anyone know of such an implementation?  Thank you in advance for any help you can offer.

	[[alternative HTML version deleted]]

