Unable to run post hoc test after macOS Ventura upgrade!

I need to run non-parametric test and posthoc test on my data, I was using PMCMRPlus package, but after upgrading MacOS to Ventura RStudio doesn't run anymore. Does anyone have any suggestions for another package for this posthoc test?

Error: 'posthoc.kruskal.nemenyi.test.default' is defunct.
Use 'PMCMRplus::kwAllPairsNemenyiTest' instead.
See help("Defunct") and help("PMCMR-defunct").

library(PMCMRplus)
ans <- posthoc.kruskal.nemenyi.test.default(count ~ spray, data = InsectSprays)
#> Error in posthoc.kruskal.nemenyi.test.default(count ~ spray, data = InsectSprays): could not find function "posthoc.kruskal.nemenyi.test.default"
ans <- kwAllPairsNemenyiTest(count ~ spray, data = InsectSprays)
#> Warning in kwAllPairsNemenyiTest.default(c(10, 7, 20, 14, 14, 12, 10, 23, : Ties
#> are present, p-values are not corrected.
summary(ans)
#> 
#>  Pairwise comparisons using Tukey-Kramer-Nemenyi all-pairs test with Tukey-Dist approximation
#> data: count by spray
#> alternative hypothesis: two.sided
#> P value adjustment method: single-step
#> H0
#>            q value   Pr(>|q|)    
#> B - A == 0   0.441 0.99960707    
#> C - A == 0   6.738 2.7914e-05 ***
#> D - A == 0   4.400 0.02292762   *
#> E - A == 0   5.435 0.00169497  **
#> F - A == 0   0.572 0.99861029    
#> C - B == 0   7.179 5.7061e-06 ***
#> D - B == 0   4.841 0.00812908  **
#> E - B == 0   5.876 0.00046606 ***
#> F - B == 0   0.131 0.99999905    
#> D - C == 0   2.338 0.56300179    
#> E - C == 0   1.303 0.94108652    
#> F - C == 0   7.311 3.4955e-06 ***
#> E - D == 0   1.035 0.97809294    
#> F - D == 0   4.973 0.00584731  **
#> F - E == 0   6.007 0.00031142 ***
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Created on 2023-01-21 with reprex v2.0.2

You should

update.packages("PMCMRplus")

update R to v.4.1 or 4.2 and run from the R.App.

Thank you so much!
but now I have another problem, kruskal.test is significant and the posthoc test doesn't shows difference. Could it be due to the message that "p-values are not corrected" using kwAllPairsNemenyiTest?

kruskal.test(DESEMP$TCE~trat)

Kruskal-Wallis rank sum test

data: DESEMP$TCE by trat
Kruskal-Wallis chi-squared = 26.446, df = 14, p-value = 0.02271

ansTCE <- kwAllPairsNemenyiTest(DESEMP$TCE~trat, dist = c("Tukey", "Chisquare"))
Warning message:
In kwAllPairsNemenyiTest.default(c(71.07142857, 32.88095238, -33.19047619, :
Ties are present, p-values are not corrected.
summary(ansTCE)

Pairwise comparisons using Tukey-Kramer-Nemenyi all-pairs test with Tukey-Dist approximation

data: DESEMP$TCE by trat
alternative hypothesis: two.sided
P value adjustment method: single-step
H0
q value Pr(>|q|)
ALIMENTA15E - ALIMENTA15 == 0 1.539 0.999075
ALIMENTA22 - ALIMENTA15 == 0 0.308 1.000000
ALIMENTA30 - ALIMENTA15 == 0 0.440 1.000000
ALIMENTA7 - ALIMENTA15 == 0 0.088 1.000000
REALIMENTA15 - ALIMENTA15 == 0 0.659 1.000000
REALIMENTA15E - ALIMENTA15 == 0 0.747 1.000000
REALIMENTA22 - ALIMENTA15 == 0 0.176 1.000000
REALIMENTA30 - ALIMENTA15 == 0 1.846 0.993658
REALIMENTA7 - ALIMENTA15 == 0 0.528 1.000000
RESTRI15 - ALIMENTA15 == 0 2.813 0.805134
RESTRI15E - ALIMENTA15 == 0 2.769 0.822891
RESTRI22 - ALIMENTA15 == 0 1.934 0.989997
RESTRI30 - ALIMENTA15 == 0 1.583 0.998739
RESTRI7 - ALIMENTA15 == 0 3.165 0.637022
ALIMENTA22 - ALIMENTA15E == 0 1.231 0.999931
ALIMENTA30 - ALIMENTA15E == 0 1.099 0.999983
ALIMENTA7 - ALIMENTA15E == 0 1.451 0.999523
REALIMENTA15 - ALIMENTA15E == 0 2.198 0.967975
REALIMENTA15E - ALIMENTA15E == 0 2.286 0.955532
REALIMENTA22 - ALIMENTA15E == 0 1.714 0.997027
REALIMENTA30 - ALIMENTA15E == 0 3.385 0.519656
REALIMENTA7 - ALIMENTA15E == 0 2.066 0.981451
RESTRI15 - ALIMENTA15E == 0 4.352 0.125190
RESTRI15E - ALIMENTA15E == 0 4.308 0.135941
RESTRI22 - ALIMENTA15E == 0 3.473 0.472960
RESTRI30 - ALIMENTA15E == 0 3.121 0.659930
RESTRI7 - ALIMENTA15E == 0 4.704 0.061275 .
ALIMENTA30 - ALIMENTA22 == 0 0.132 1.000000
ALIMENTA7 - ALIMENTA22 == 0 0.220 1.000000
REALIMENTA15 - ALIMENTA22 == 0 0.967 0.999997
REALIMENTA15E - ALIMENTA22 == 0 1.055 0.999990
REALIMENTA22 - ALIMENTA22 == 0 0.484 1.000000
REALIMENTA30 - ALIMENTA22 == 0 2.154 0.973106
REALIMENTA7 - ALIMENTA22 == 0 0.835 0.999999
RESTRI15 - ALIMENTA22 == 0 3.121 0.659930
RESTRI15E - ALIMENTA22 == 0 3.077 0.682452
RESTRI22 - ALIMENTA22 == 0 2.242 0.962134
RESTRI30 - ALIMENTA22 == 0 1.890 0.991998
RESTRI7 - ALIMENTA22 == 0 3.473 0.472960
ALIMENTA7 - ALIMENTA30 == 0 0.352 1.000000
REALIMENTA15 - ALIMENTA30 == 0 1.099 0.999983
REALIMENTA15E - ALIMENTA30 == 0 1.187 0.999956
REALIMENTA22 - ALIMENTA30 == 0 0.615 1.000000
REALIMENTA30 - ALIMENTA30 == 0 2.286 0.955532
REALIMENTA7 - ALIMENTA30 == 0 0.967 0.999997
RESTRI15 - ALIMENTA30 == 0 3.253 0.590383
RESTRI15E - ALIMENTA30 == 0 3.209 0.613811
RESTRI22 - ALIMENTA30 == 0 2.374 0.939854
RESTRI30 - ALIMENTA30 == 0 2.022 0.984775
RESTRI7 - ALIMENTA30 == 0 3.605 0.405116
REALIMENTA15 - ALIMENTA7 == 0 0.747 1.000000
REALIMENTA15E - ALIMENTA7 == 0 0.835 0.999999
REALIMENTA22 - ALIMENTA7 == 0 0.264 1.000000
REALIMENTA30 - ALIMENTA7 == 0 1.934 0.989997
REALIMENTA7 - ALIMENTA7 == 0 0.615 1.000000
RESTRI15 - ALIMENTA7 == 0 2.901 0.767082
RESTRI15E - ALIMENTA7 == 0 2.857 0.786515
RESTRI22 - ALIMENTA7 == 0 2.022 0.984775
RESTRI30 - ALIMENTA7 == 0 1.670 0.997739
RESTRI7 - ALIMENTA7 == 0 3.253 0.590383
REALIMENTA15E - REALIMENTA15 == 0 0.088 1.000000
REALIMENTA22 - REALIMENTA15 == 0 0.484 1.000000
REALIMENTA30 - REALIMENTA15 == 0 1.187 0.999956
REALIMENTA7 - REALIMENTA15 == 0 0.132 1.000000
RESTRI15 - REALIMENTA15 == 0 2.154 0.973106
RESTRI15E - REALIMENTA15 == 0 2.110 0.977580
RESTRI22 - REALIMENTA15 == 0 1.275 0.999895
RESTRI30 - REALIMENTA15 == 0 0.923 0.999998
RESTRI7 - REALIMENTA15 == 0 2.506 0.909575
REALIMENTA22 - REALIMENTA15E == 0 0.571 1.000000
REALIMENTA30 - REALIMENTA15E == 0 1.099 0.999983
REALIMENTA7 - REALIMENTA15E == 0 0.220 1.000000
RESTRI15 - REALIMENTA15E == 0 2.066 0.981451
RESTRI15E - REALIMENTA15E == 0 2.022 0.984775
RESTRI22 - REALIMENTA15E == 0 1.187 0.999956
RESTRI30 - REALIMENTA15E == 0 0.835 0.999999
RESTRI7 - REALIMENTA15E == 0 2.418 0.930696
REALIMENTA30 - REALIMENTA22 == 0 1.670 0.997739
REALIMENTA7 - REALIMENTA22 == 0 0.352 1.000000
RESTRI15 - REALIMENTA22 == 0 2.638 0.870606
RESTRI15E - REALIMENTA22 == 0 2.594 0.884579
RESTRI22 - REALIMENTA22 == 0 1.758 0.996133
RESTRI30 - REALIMENTA22 == 0 1.407 0.999665
RESTRI7 - REALIMENTA22 == 0 2.989 0.726010
REALIMENTA7 - REALIMENTA30 == 0 1.319 0.999842
RESTRI15 - REALIMENTA30 == 0 0.967 0.999997
RESTRI15E - REALIMENTA30 == 0 0.923 0.999998
RESTRI22 - REALIMENTA30 == 0 0.088 1.000000
RESTRI30 - REALIMENTA30 == 0 0.264 1.000000
RESTRI7 - REALIMENTA30 == 0 1.319 0.999842
RESTRI15 - REALIMENTA7 == 0 2.286 0.955532
RESTRI15E - REALIMENTA7 == 0 2.242 0.962134
RESTRI22 - REALIMENTA7 == 0 1.407 0.999665
RESTRI30 - REALIMENTA7 == 0 1.055 0.999990
RESTRI7 - REALIMENTA7 == 0 2.638 0.870606
RESTRI15E - RESTRI15 == 0 0.044 1.000000
RESTRI22 - RESTRI15 == 0 0.879 0.999999
RESTRI30 - RESTRI15 == 0 1.231 0.999931
RESTRI7 - RESTRI15 == 0 0.352 1.000000
RESTRI22 - RESTRI15E == 0 0.835 0.999999
RESTRI30 - RESTRI15E == 0 1.187 0.999956
RESTRI7 - RESTRI15E == 0 0.396 1.000000
RESTRI30 - RESTRI22 == 0 0.352 1.000000
RESTRI7 - RESTRI22 == 0 1.231 0.999931
RESTRI7 - RESTRI30 == 0 1.583 0.998739

Signif. codes: 0 ‘’ 0.001 ‘’ 0.01 ‘’ 0.05 ‘.’ 0.1 ‘ ’ 1

signifies that there were tied ranks. For an explanation of how this affects the way in which test statistics are calculated, see this paper by the package author.

This topic was automatically closed 7 days after the last reply. New replies are no longer allowed.

If you have a query related to it or one of the replies, start a new topic and refer back with a link.