Valor crítico de la prueba de rango con signo de una muestra de Wilcoxon en R

Question

Intento encontrar el valor crítico para la prueba de rango con signo de una muestra de Wilcoxon. Actualmente, puedo encontrar el valor utilizando tablas. He buscado en qwilcox() en R, pero parece que esto da los valores críticos para la prueba de dos muestras de Wilcoxon (o prueba de Mann-Whitney). ¿Existe una función en R que puedo utilizar para calcular este valor crítico?

Answer 1

Puede utilizar la función qsignrank() función. Ejemplo:

> qsignrank(.025, 10, lower.tail=FALSE)
46

Esto significa que para un tamaño de muestra de 10 y un dos caras con un nivel de significación del 5%, el estadístico de la prueba debe ser superior a 46 (es decir, 47 o más) para ser estadísticamente significativo. Ejemplo de datos:

> set.seed(1)
> x = rnorm(10, .5)
> wilcox.test(x)

    Wilcoxon signed rank test

data:  x
V = 47, p-value = 0.04883
alternative hypothesis: true location is not equal to 0

En este caso, la estadística de la prueba es 47, y significativa al nivel del 5%.

Tenga en cuenta que para una prueba de dos caras, el estadístico de prueba devuelto por qsignrank() es el más grande de los dos estadísticos de prueba posibles. Por ejemplo, wilcox.test(-x) da un estadístico de prueba de 8, que puede transformarse en 47 mediante $\frac{10\cdot 11}{2}-8$ .