R中多重比较的成对检验(Nemenyi检验)

Question

我正在使用PMCMR包执行博士后Kruskal-Nemenyi测试。当我使用默认设置运行测试时：

posthoc.kruskal.nemenyi.test(x=coastal$HIGH_MORTGAGE, g=coastal$SIZECLASS, method="Tukey")

我得到以下结果和警告：

Pairwise comparisons using Tukey and Kramer (Nemenyi) test  
                   with Tukey-Dist approximation for independent samples

data:  coastal$HIGH_MORTGAGE and coastal$SIZECLASS 

       Large Medium
Medium 0.931 -     
Small  0.746 0.078 

P value adjustment method: none

Warning message:
In posthoc.kruskal.nemenyi.test.default(x = coastal$HIGH_MORTGAGE,  :
  Ties are present, p-values are not corrected.*

当我运行测试将分布更改为Chisq以对关系进行校正时，我仍然得到相同的结果，并且没有使用卡方分布。

posthoc.kruskal.nemenyi.test(x=coastal$HIGH_MORTGAGE, g=coastal$SIZECLASS, method="Chisq")

Pairwise comparisons using Tukey and Kramer (Nemenyi) test  
                   with Tukey-Dist approximation for independent samples

data:  coastal$HIGH_MORTGAGE and coastal$SIZECLASS

       Large Medium
Medium 0.931 -     
Small  0.746 0.078

P value adjustment method: none 

Warning message:
In posthoc.kruskal.nemenyi.test.default(x = coastal$HIGH_MORTGAGE,  :
  Ties are present, p-values are not corrected.

我想知道包中是否有bug，或者是否有我不知道的方法来解决这个问题。

Answer 1

从PMCMR1.0版本到PMCMR1.1 (和> 1.1)，语法略有变化，所以它是dist，而不是method，它是：

posthoc.kruskal.nemenyi.test( x, g, dist = c("Tukey", "Chisquare"), ...)

或

posthoc.kruskal.nemenyi.test(formula, data, subset, na.action, dist =
 c("Tukey", "Chisquare"), ...)

包含示例的vignette已在PMCMR1.3版本中进行了相应更新。