Posterior a observar en varios repositorios y sitios webs al concluir encontramos la respuesta que te mostramos ahora.
Solución:
Otro enfoque además de kruskal::agricolae mencionado por Marek, es la prueba Nemenyi-Damico-Wolfe-Dunn implementada en la página de ayuda para oneway_test en el paquete de moneda que usa multcomp. Usando la configuración de Hadley y reduciendo el valor B= para la función aproximada () para que finalice en un tiempo finito:
#updated translation of help page implementation of NDWD
NDWD <-
independence_test(dv ~ iv, data = sum_codings1, distribution = approximate(B = 10000),
ytrafo = function(data) trafo(data, numeric_trafo = rank_trafo),
xtrafo = mcp_trafo(iv = "Tukey"))
### global p-value
print(pvalue(NDWD))
### sites (I = II) != (III = IV) at alpha = 0.01 (page 244)
print(pvalue(NDWD, method = "single-step"))
Los resultados más estables en ese conjunto de datos más grande pueden requerir aumentar el valor B y aumentar la paciencia del usuario.
Ene: 2012: Recientemente hubo una publicación en R-help que afirmaba resultados inesperados de este método, así que reenvié ese correo electrónico al mantenedor. Mark Difford dijo que había confirmado los problemas y ofreció pruebas alternativas con el paquete nparcomp: https://stat.ethz.ch/pipermail/r-help/2012-January/300100.html
También hubo en la misma semana un par de otras sugerencias sobre ayuda para contrastes post-hoc con pruebas KW: kruskalmc sugerido por Mario Garrido Escudero y
rms::polr seguido por rms::contrasts sugerido por Frank Harrell https://stat.ethz.ch/pipermail/r-help/2012-January/300329.html
Nov 2015: acuerdo con toto_tico que el código de la página de ayuda del paquete de monedas se ha cambiado en los años intermedios. El ?independence_test La página de ayuda ahora ofrece una prueba KW multivariable y la ?oneway_test La página de ayuda ha reemplazado su implementación anterior con el código anterior usando el independence_test función.
Puede usar el paquete PMCMR. Hay más información al respecto.
Spelling_Grades <- c(90,87,89,90,75,88,97,99,78,85,72,76,77,79,70)
Methods <- c("A","A","A","A","B","B","B","B","B","B","C","C","C","C","C")
kruskalmc(Spelling_Grades~Methods)
#This method doesn't accept characters that's why I've changed the methods to integer
Methods <- c(1,1,1,1,2,2,2,2,2,2,3,3,3,3,3)
posthoc.kruskal.nemenyi.test(Spelling_Grades~Methods)
Los dos métodos anteriores dan los mismos resultados.
Acuérdate de que tienes la opción de añadir una tasación objetiva .