[Rd] kruskalmc(pgirmess) (PR#10664)

From: <sommerfeld_at_mpil-ploen.mpg.de>
Date: Wed, 30 Jan 2008 11:20:22 +0100 (CET)


Full_Name: Ralf Sommerfeld
Version: 2.3.1
OS: Windows XP
Submission from: (NULL) (193.175.214.49)

The kruskalmc function runs a multiple comparison test after a Kruskal-Wallis test. I used the function for treatment comparison (default) as compared to a comparison against a control. The alignment of names and obs.dif is however in the wrong order!
++++++++++++++
The correct ordering (lines 13-19) is as follows:

         for (i in 2:(length(name))) {
             for (j in 1:(i-1)) {
                 vname <- c(vname, paste(name[j], "-", name[i], 
                   sep = ""))
		 indices <- rbind(indices, c( j, i))
             }
         }

instead of:

        for (i in 1:(length(name) - 1)) {
            for (j in (i + 1):length(name)) {
                vname <- c(vname, paste(name[i], "-", name[j], 
                  sep = ""))
                indices <- rbind(indices, c(i, j))
            }
        }

+++++++++++++++

I did not check the algorithm for ordering a two-sided or one-sided test against a control.

All the best,
Ralf



R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Wed 30 Jan 2008 - 17:59:59 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 30 Jan 2008 - 18:30:14 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive