๐Ÿ“‰ ๐€๐๐Ž๐•๐€ ๐š̀ ๐๐ž๐ฎ๐ฑ ๐Ÿ๐š๐œ๐ญ๐ž๐ฎ๐ซ๐ฌ (๐“๐ฐ๐จ-๐–๐š๐ฒ ๐€๐๐Ž๐•๐€) ๐ฌ๐จ๐ฎ๐ฌ ๐‘ : ๐†๐ฎ๐ข๐๐ž ๐œ๐ฅ๐š๐ข๐ซ ๐ž๐ญ ๐š๐ฉ๐ฉ๐ฅ๐ข๐ช๐ฎรฉ

Dans cet exemple, nous analysons :

·         Variable dรฉpendante : ๐‘๐‘œ๐‘š๐‘๐‘Ÿ๐‘’ ๐‘‘๐‘’ ๐‘Ž๐‘๐‘Ž๐‘Ÿ๐‘–๐‘’๐‘›๐‘  (Mites)

·         Facteur 1 : ๐‘‡๐‘ฆ๐‘๐‘’ ๐‘‘๐‘’ ๐‘ ๐‘–๐‘ก๐‘’ (Brรปlรฉ / Non brรปlรฉ)

·         Facteur 2 : ๐‘†๐‘’๐‘ฅ๐‘’ (Mรขle / Femelle)

Voici un guide pas ร  pas, rรฉdigรฉ dans un style clair et accessible aux chercheurs dรฉbutants.

๐™๐™ฃ ๐™ก๐™ž๐™ ๐™š ๐™š๐™ฉ ๐™ช๐™ฃ ๐™‹๐™–๐™ง๐™ฉ๐™–๐™œ๐™š de ce post avant d'aller plus loin nous feront plaisir.


๐Ÿ”ด ๐Ÿ. ๐๐ซ๐ž́๐ฉ๐š๐ซ๐š๐ญ๐ข๐จ๐ง ๐๐ž๐ฌ ๐๐จ๐ง๐ง๐ž́๐ž๐ฌ (๐ž́๐ญ๐š๐ฉ๐ž ๐ฅ๐š ๐ฉ๐ฅ๐ฎ๐ฌ ๐ข๐ฆ๐ฉ๐จ๐ซ๐ญ๐š๐ง๐ญ๐ž)

Les donnรฉes doivent รชtre en format long :

Site

Sexe

Mites

Burned

Male

95

Burned

Male

100

Burned

Female

30

Burned

Female

32

Unburned

Male

75

Unburned

Male

78

Unburned

Female

25

Unburned

Female

27

Une ligne = une observation
Chaque facteur dans une colonne sรฉparรฉe
Variable rรฉponse dans une seule colonne


๐Ÿ”ด ๐Ÿ. ๐ˆ๐ฆ๐ฉ๐จ๐ซ๐ญ๐ž๐ซ ๐ฅ๐ž๐ฌ ๐๐จ๐ง๐ง๐ž́๐ž๐ฌ ๐๐š๐ง๐ฌ ๐‘

data <- read.csv("mites_data.csv")

Ou crรฉation manuelle :

data <- data.frame(
  Site = factor(c("Burned","Burned","Burned","Burned",
                  "Unburned","Unburned","Unburned","Unburned")),
  Sex = factor(c("Male","Male","Female","Female",
                 "Male","Male","Female","Female")),
  Mites = c(95,100,30,32,75,78,25,27)
)

๐Ÿ”ด ๐Ÿ‘. ๐•๐ž́๐ซ๐ข๐Ÿ๐ข๐ž๐ซ ๐ฅ๐ž๐ฌ ๐ก๐ฒ๐ฉ๐จ๐ญ๐ก๐ž̀๐ฌ๐ž๐ฌ ๐๐ž ๐ฅ๐€๐๐Ž๐•๐€

(a) Normalitรฉ des rรฉsidus

model <- aov(Mites ~ Site * Sex, data = data)
shapiro.test(residuals(model))

p > 0,05 → normalitรฉ respectรฉe

(b) Homogรฉnรฉitรฉ des variances

install.packages("car")
library(car)
leveneTest(Mites ~ Site * Sex, data = data)

p > 0,05 → variances homogรจnes


๐Ÿ”ด ๐Ÿ’. ๐‹๐š๐ง๐œ๐ž๐ซ ๐ฅ๐€๐๐Ž๐•๐€ ๐š̀ ๐๐ž๐ฎ๐ฑ ๐Ÿ๐š๐œ๐ญ๐ž๐ฎ๐ซ๐ฌ (๐ž́๐ญ๐š๐ฉ๐ž ๐œ๐ž๐ง๐ญ๐ซ๐š๐ฅ๐ž)

anova_result <- aov(Mites ~ Site * Sex, data = data)
summary(anova_result)

Le rรฉsultat fournit trois effets principaux :

1 Effet Site
2
Effet Sexe
3
Interaction Site × Sexe


๐Ÿ”ด ๐Ÿ“. ๐ˆ๐ง๐ญ๐ž๐ซ๐ฉ๐ซ๐ž́๐ญ๐š๐ญ๐ข๐จ๐ง ๐๐ž๐ฌ ๐ซ๐ž́๐ฌ๐ฎ๐ฅ๐ญ๐š๐ญ๐ฌ

Si :

·         Effet Site significatif (p < 0,05)
→ Le nombre d’acariens diffรจre entre les sites brรปlรฉs et non brรปlรฉs

·         Effet Sexe significatif (p < 0,05)
→ Les mรขles et les femelles prรฉsentent des diffรฉrences d’abondance

·         Interaction significative (p < 0,05)
→ L’effet du site dรฉpend du sexe (rรฉsultat biologiquement trรจs important)

Si l’interaction est significative, on l’interprรจte toujours en prioritรฉ.


๐Ÿ”ด ๐Ÿ”. ๐“๐ž๐ฌ๐ญ ๐ฉ๐จ๐ฌ๐ญ-๐ก๐จ๐œ (๐“๐ฎ๐ค๐ž๐ฒ ๐‡๐’๐ƒ)

Si l’ANOVA est significative :

TukeyHSD(anova_result)

Ce test permet d’expliquer les lettres (a, b, c) affichรฉes dans les figures scientifiques.


๐Ÿ”ด ๐Ÿ•. ๐†๐ซ๐š๐ฉ๐ก๐ข๐ช๐ฎ๐ž ๐๐ข๐ง๐ญ๐ž๐ซ๐š๐œ๐ญ๐ข๐จ๐ง

interaction.plot(
  x.factor = data$Site,
  trace.factor = data$Sex,
  response = data$Mites,
  col = c("purple", "gold"),
  lwd = 2,
  ylab = "Nombre d'acariens",
  xlab = "Type de site"
)

Ce graphique permet de visualiser l’existence (ou non) d’une interaction entre les deux facteurs.


๐ŸŽฏ Conclusion

L’ANOVA ร  deux facteurs est un outil puissant pour analyser simultanรฉment deux variables explicatives et leur interaction. Elle est largement utilisรฉe en biologie, agronomie, sciences sociales et recherche expรฉrimentale.


Si vous avez trouvรฉ cette publication utile, n'hรฉsitez pas ร  ๐’๐’‚ ๐’๐’Š๐’Œ๐’†๐’“ ๐’†๐’• ร  ๐’๐’‚ ๐’‘๐’‚๐’“๐’•๐’‚๐’ˆ๐’†๐’“ avec vos amis et collรจgues ! Pour mieux apprendre l’utilisation des logiciel et modรจles statistiques, nous vous invitons ร  prendre part ร  la prochaine session de notre formation en ๐™€๐™˜๐™ค๐™ฃ๐™ค๐™ขรฉ๐™ฉ๐™ง๐™ž๐™š ๐™š๐™ฉ ๐™๐™š๐™˜๐™๐™ฃ๐™ž๐™ฆ๐™ช๐™š๐™จ ๐™Œ๐™ช๐™–๐™ฃ๐™ฉ๐™ž๐™ฉ๐™–๐™ฉ๐™ž๐™ซ๐™š๐™จ https://forms.gle/mopoGRWKwyTtbEJ16



#ANOVA
#StatistiquesAvecR
#AnalyseDeDonnรฉes
#MรฉthodesQuantitatives
#DataScience
#RechercheScientifique
#FormationStatistique
#RProgramming

 

Commentaires

Posts les plus consultรฉs de ce blog

ร‰conomรฉtrie des donnรฉes de panel: de la thรฉorie ร  la pratique