[R] - dunn.test gives strange output

Joachim Audenaert Joachim.Audenaert at pcsierteelt.be
Wed Apr 22 13:14:06 CEST 2015


Dear all,

I have a question concerning my output of the dunn.test function in R. I 
like to compare different datasets, which are not distributed normally, so 
I use the Dunn.test to do pairwise comparison. I have 2 questions 
concerning the output:

- why are my groupnames changed, the output gives 4 times swirskii, 
although my groupnames are longer (see my dataset)
- secondly when I check the p-values, I see something very odd: I see that 
p values where the same groups are compared sometimes indicate that they 
are different, why are not all those p-values = 1 ?

This is my dataset:

structure(list(controle = c(111, 88, 216, 169), chemie = c(47, 31, 35, 
30), IPM = c(0, 0, 0, 1), gallicus = c(102, 152, 102, 75), swirskii3 = 
c(1, 0, 0, 0), swirskiiA = c(0, 0, 1, 3), swirskiiP = c(0, 0, 1, 6), 
swirskii1x = c(12, 2, 75, 46)), .Names = c("controle", "chemie", "IPM", 
"gallicus", "swirskii3", "swirskiiA", "swirskiiP", "swirskii1x"), 
row.names = c(NA, 4L), class = "data.frame")

I get the following output:

 DUNN <- dunn.test(value, variable, method = "none",kw=TRUE)
 
  Kruskal-Wallis rank sum test

data: value and variable
Kruskal-Wallis chi-squared = 26.8894, df = 7, p-value = 0


                        Comparison of value by variable  
                                (No adjustment)  
Col Mean-|
Row Mean |   controle     chemie        IPM   gallicus   swirskii swirskii
---------+------------------------------------------------------------------
  chemie |  -1.341044  -3.410083  -2.069039  -0.325682   1.015361 3.084401
         |     0.0900     0.0003     0.0193     0.3723     0.1550 0.0010
         |
     IPM |  -3.410083  -2.069039  -0.325682   1.015361   3.084401 
-3.410083
         |     0.0003     0.0193     0.3723     0.1550     0.0010 0.0003
         |
gallicus |  -0.325682   1.015361   3.084401  -3.410083  -2.069039 0.000000
         |     0.3723     0.1550     0.0010     0.0003     0.0193 0.5000
         |
swirskii |  -3.410083  -2.069039   0.000000  -3.084401  -3.007770 
-1.666726
         |     0.0003     0.0193     0.5000     0.0010     0.0013 0.0478
         |
swirskii |  -3.007770  -1.666726   0.402313  -2.682088   0.402313 
-2.969454
         |     0.0013     0.0478     0.3437     0.0037     0.3437     
0.0015
         |
swirskii |  -2.969454  -1.628410   0.440628  -2.643772   0.440628 0.038315
         |     0.0015     0.0517     0.3297     0.0041     0.3297 0.4847
         |
swirskii |  -1.475148  -0.134104   1.934935  -1.149466   1.934935 1.532621
         |     0.0701     0.4467     0.0265     0.1252     0.0265 0.0627
Col Mean-|
Row Mean |   swirskii
---------+-----------
  chemie |  -3.410083
         |     0.0003
         |
     IPM |  -2.069039
         |     0.0193
         |
gallicus |  -3.084401
         |     0.0010
         |
swirskii |   0.402313
         |     0.3437
         |
swirskii |  -1.628410
         |     0.0517
         |
swirskii |  -1.475148
         |     0.0701
         |
swirskii |   1.494306
         |     0.0675 




Met vriendelijke groeten - With kind regards,

Joachim Audenaert 
onderzoeker gewasbescherming - crop protection researcher

PCS | proefcentrum voor sierteelt - ornamental plant research

Schaessestraat 18, 9070 Destelbergen, België
T: +32 (0)9 353 94 71 | F: +32 (0)9 353 94 95
E: joachim.audenaert at pcsierteelt.be | W: www.pcsierteelt.be   

Heb je je individuele begeleiding bemesting (CVBB) al aangevraagd? | Het 
PCS op LinkedIn
Disclaimer | Please consider the environment before printing. Think green, 
keep it on the screen!
	[[alternative HTML version deleted]]



More information about the R-help mailing list