๐ ๐๐๐๐๐ ๐̀ ๐๐๐ฎ๐ฑ ๐๐๐๐ญ๐๐ฎ๐ซ๐ฌ (๐๐ฐ๐จ-๐๐๐ฒ ๐๐๐๐๐) ๐ฌ๐จ๐ฎ๐ฌ ๐ : ๐๐ฎ๐ข๐๐ ๐๐ฅ๐๐ข๐ซ ๐๐ญ ๐๐ฉ๐ฉ๐ฅ๐ข๐ช๐ฎรฉ
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
.png)
Commentaires
Enregistrer un commentaire