Bonjour
J'ai un tableau avec plusieurs colonnes (échantillons) et ces valeurs j'aimerais comparer les échantillons entre eux avec un test de wilcoxon, comment automatiser cette test pour comparer tous les échantillons en 1 seul pas? ou je dois le faire deux par deux?
MERCI
Bonjour.
C'est pas clair ! les tests de Wilcoxon permettent de comparer 2 échantillons (appariés, ou pas, suivant le test). Si tu veux "comparer" (quoique ça signifie) plusieurs échantillons, ce n'est plus un Wilcoxon qu'il faut. Et si tu veux les comparer 2 à deux, tu fais le nombre de tests nécessaires.
Tout ça dépend de ce que tu veux savoir ou mettre en évidence.
Cordialement.
Bonjour
oui en effet je veux les comparer deux a deux et je cherche une façon de faire ces comparaisons en 1 seule pas avec un script,
par example la function tapply permets de faire le même calcul (par exemple: moyenne, indice de diversite, variance etc) pour toutes les échantillons d'un tableau
moi je veux une function pour le dire a R de comparer ces colonnes entre elles deux a deux sans pour autant devoir taper le code wilcox.test (a,b,....) a chaque fois pour chaque couple. J'essaie avec la function tapply et ça ne marche pas...
MERCI d'avance.
tu peux toujours faire avec des boucles. Supposons que seule la P-value t'intéresse. supposons aussi que ton tableau s'appelle x.
n<-ncol(x)
pvalues<-matrix(1,ncol=n,nrow=n)
for(i in 1:n)for(j in (i+1):n) {
pvalues[i,j]<-wilcox.test[x[,i],x[,j])$p.value
p.values[j,i]<-p.values[i,j])
}
bonjour,
j'ai des messages d’erreur avec votre code et je ne sais pas très bien utiliser les boucles,
> for(i in 1:n)for(j in (i+1):n) {
+ pvalues[i,j]<-wilcox.test[x[,i],x[,j])$p.value
Erreur : ')' inattendu(e) in:
"for(i in 1:n)for(j in (i+1):n) {
pvalues[i,j]<-wilcox.test[x[,i],x[,j])"
> p.values[j,i]<-p.values[i,j])
Erreur : ')' inattendu(e) in "p.values[j,i]<-p.values[i,j])"
> }
MERCI
ah oui j'ai oublié la parenthèse ouvrante dans l'appel à wilcox.test.
merci !
si non j'ai trouve une maniéré de réaliser le test sans avoir recours aux boucles:
with(data,pairwise.wilcox.test (y,x, p.adj="none"))
MERCI
