title: kruskal-wallis contrast test in r
body: <p>All,</p>
  
  <p>I would like to perform the equivalent of TukeyHSD on the rank ordering median shift test that such as kruskal wallis</p>
  
  <pre><code>X=matrix(c(1,1,1,1,2,2,2,4,4,4,4,4,1,3,6,9,4,6,8,10,1,2,1,3),ncol=2)
  anova=aov(X[,2]~factor(X[,1]))
  TukeyHSD(anova)
  ## Tukey multiple comparisons of means
  ## 95% family-wise confidence level
  ##
  ## Fit: aov(formula = X[, 2] ~ factor(X[, 1]))
  ##
  ## $`factor(X[, 1])`
  ## diff lwr upr p adj
  ## 2-1 1.25 -5.927068 8.427068 0.8794664
  ## 4-1 -1.35 -7.653691 4.953691 0.8246844
  ## 4-2 -2.60 -9.462589 4.262589 0.5617125
  kruskal.test(X[,2]~factor(X[,1]))
  ##
  ## Kruskal-Wallis rank sum test
  ##
  ## data: X[, 2] by factor(X[, 1])
  ## Kruskal-Wallis chi-squared = 1.7325, df = 2, p-value = 0.4205
  </code></pre>
  
  <p>I would like now to analyze the contrasts. Please help. Thanks.
  Rik</p>