Je veux déterminer s'il y a une différence dans les valeurs p moyennes entre deux groupes. Pour ce faire, j'effectue un test de somme de rang de Wilcoxon (les données ne sont pas normalement distribuées). Jusqu'ici tout va bien. Enfin, je veux calculer la taille d'effet correspondante. Malheureusement, R ne fournit pas cela. Il ne fournit pas non plus de valeur az avec laquelle la taille de l'effet peut être facilement calculée en utilisant: taille de l'effet = z / sqrt (N)
voici un exemple de code R:
a=rep(0:1,each=20) #grouping variable
b=c(rnorm(20, .03,.01), rnorm(20, .02, .009)) #vector of p-values
d=cbind(a,b)
test = wilcox.test(b ~ a, data = d) #perform Wilcoxon rank-sum test
test
Quelqu'un sait-il comment obtenir la taille de l'effet?
Réponses:
L'estimateur qui correspond au test de Wilcoxon est l'estimateur de Hodges-Lehmann; il est retourné en
wilcox.test
utilisant l'conf.int=TRUE
option, sous "différence de lieu".Pour votre exemple:
Pour en savoir plus sur le Wilcoxon et les hypothèses sous-jacentes, et ce qu'il teste réellement, et d'autres estimateurs non paramétriques, ce document est (peut-être) utile: www.stat.umn.edu/geyer/old03/5102/notes/rank.pdf
la source
wilcox.test(b~a,data=d, conf.int=TRUE)$estimate / sqrt(20)
correcte?Obtenez le z de votre formule en
et calculez la taille de l'effet avec votre formule, en définissant N sur 40
la source
coin::wilcoxsign_test
fonction R. Parlez-vous également de la formule OP, ?