Per eseguire in R il test di Wilcoxon per due campioni è necesario utilizzare la funzione wilcox.test.
Vediamo un esempio di applicazione del test a due code su due campioni da (66 elementi) e db (40 elementi) non appaiati
wilcox.test(da,db,alternative="two.sided",paired=FALSE)
il risultato sarà nel nostro caso
Wilcoxon rank sum test with continuity correction data: da and db W = 1142.5, p-value = 0.2482 alternative hypothesis: true location shift is not equal to 0
Da notare che vengno restituiti sia il P-value che il valore della statistica test W. La funzione wilcox.test prende in input molti parametri, alcuni dei quali opzionali; nell'esempio ne abbiamo utilizzati due:
Se inseriamo il risultato del test in una variabile invece che stamparlo a video otteniamo una struttura come la seguente
structure(list(statistic = structure(1142.5, .Names = "W"), parameter = NULL, p.value = 0.248158967552218, null.value = structure(0, .Names = "location shift"), alternative = "two.sided", method = "Wilcoxon rank sum test with continuity correction", data.name = "da and db"), .Names = c("statistic", "parameter", "p.value", "null.value", "alternative", "method", "data.name" ), class = "htest")