Ho i risultati dello stesso test applicato a due campioni indipendenti:
x <- c(17, 12, 13, 16, 9, 19, 21, 12, 18, 17)
y <- c(10, 6, 15, 9, 8, 11, 8, 16, 13, 7, 5, 14)
E voglio calcolare un test di somma dei ranghi di Wilcoxon.
Quando calcolo la statistica a mano, ottengo: T W = ∑ rango ( X i ) = 156,5
Quando lascio eseguire a R un wilcox.test(x, y, correct = F)
, ottengo:
W = 101.5
Perché? Non dovrebbe la statistica essere restituito solo quando eseguo un firmato rank test con ? O fraintendere il test di somma dei ranghi?paired = T
Come posso dire a R di emettere
Come parte dei risultati del test, non attraverso qualcosa come:
dat <- data.frame(v = c(x, y), s = factor(rep(c("x", "y"), c(10, 12))))
dat$r <- rank(dat$v)
T.W <- sum(dat$r[dat$s == "x"])
Ho fatto una domanda di follow-up sul significato dei diversi modi per calcolare la statistica del test per il test di somma dei ranghi di Wilcoxon