1 pièce(s) jointe(s)
Anova successives et boucle for
Bonjour,
Si je poste ici, c'est que je suis confronté au problème suivant :
je travaille sur des communautés de champignons du sol, et dans le design expérimental que j'ai réalisé, j'ai regardé l'effet de différents traitements sur des parcelles de sol. Cela a abouti à l'obtention d'un tableau de ce style.
Pièce jointe 209045
pour voir si l'effet traitement est significatif sur l'otu2 je réalise l'opération suivante :
Code:
1 2
| otu2 <- lm(traitement~otu2, data=OTU_R_Ready)
anova(otu2) |
et apparait la sortie telle que :
Code:
1 2 3 4 5 6
| Analysis of Variance Table
Response: traitement
Df Sum Sq Mean Sq F value Pr(>F)
otu2 1 3.8 3.797 0.3119 0.5783
Residuals 70 852.2 12.174 |
N'ayant pas vraiment l'envie de retaper à la main 82 fois l’opération précédente pour avoir le résultat pour les autres otu, je me suis dit qu'il était sans doute possible de réaliser une boucle for pour automatiser le processus.
Cependant je manque de connaissances en la matière. Auriez-vous une idée de ce que cela pourrait donner ?
En vous remerciant par avance de toute l'aide que vous pourrez m'apporter.