Discussion :

[bouda53]

Bonjour,

J'aimerais automatiser le calcul de test de wilcox.test

Exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

wilcox.test(fichier$parametre1[fichier$condition=="pd" ],fichier$parametre1[fichier$condition=="avc"])

Sachant que j'ai 8 paramètres et 6 conditions a tester et que je voudrais extraire uniquement les p.value dans une matrice qui devrait être croisée.

Merci à vous.

[dev_ggy]

Bonjour Bouda,

Par manque de clarté de ta question, j'essaie de répondre selon ce que j'en comprends. J'ai pris deux conditions et teste toutes les combinaisons avec chacune des variables et conditions. Concernant la matrice, comme j'ai qu'une double entrée et non une triple avec deux conditions et une variable, j'ai mis toutes les conditions en colonne. De plus, dans le fichier j'ai supposé que tes variables 1 à 6 étaient placées dans un tableau aux positions 1 à 6. Le code n'est peut-être pas optimisé et ne correspond peut-être pas tout à fait à ton problème mais j'espère qu'il t'aidera dans ta recherche.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
condition<-c("pd", "avc")
Mat <- matrix(length(condition)*length(condition)*6, ncol = (length(condition)+length(condition))
for (i in 1:length(condition)){
  for (j in 1:length(condition)){
   for (k in 1:6){
    Mat(k, i+(j-1)*length(condition))<-wilcox.test(fichier[k,fichier$condition==i ],fichier[k,fichier$condition==j])$p.value
    }
  } 
}

[bouda53]

Bonjours et merci pour votre réponse;

J'aimerais savoir dans le script que vous m'avez posté à quel moment les paramètres sont pris en compte?

Je viens de faire un test en prenant en considération mes 6 conditions, et j'aimerais savoir à quoi correspond for (k in 1:6), sachant que vous ne considérez que 2 conditions?