Bonjour,
Alors voila mon problème, après une anova j'ai des résidus non-homoscedastiques...Et j'essaye tout pour transformer mes données (Boxcox (lambda=1 sniff...), racines, log, arcsin...) mais rien ne marche...Que dois je faire?
Si ca peut aider voila mon jeux de données:(rendement de pois en fonction de traitement ou pas (0/1) en N, K et P sur 6 blocs
block N P K yield
1 1b 0N 1P 1K 49.5
2 1b 1N 1P 0K 62.8
3 1b 0N 0P 0K 46.8
4 1b 1N 0P 1K 57.0
5 2b 1N 0P 0K 59.8
6 2b 1N 1P 1K 58.5
7 2b 0N 0P 1K 55.5
8 2b 0N 1P 0K 56.0
9 3b 0N 1P 0K 62.8
10 3b 1N 1P 1K 55.8
11 3b 1N 0P 0K 69.5
12 3b 0N 0P 1K 55.0
13 4b 1N 0P 0K 62.0
14 4b 1N 1P 1K 48.8
15 4b 0N 0P 1K 45.5
16 4b 0N 1P 0K 44.2
17 5b 1N 1P 0K 52.0
18 5b 0N 0P 0K 51.5
19 5b 1N 0P 1K 49.8
20 5b 0N 1P 1K 48.8
21 6b 1N 0P 1K 57.2
22 6b 1N 1P 0K 59.0
23 6b 0N 1P 1K 53.2
24 6b 0N 0P 0K 56.0
Les résultats de l'anova :
> aovrdt=aov(yield~block+N*P*K)
> summary(aovrdt) ##Effet azote, effet bloc et effet K mais pas d'effet de P et des interactions###
Df Sum Sq Mean Sq F value Pr(>F)
block 5 343.29 68.66 4.4467 0.015939 *
N 1 189.28 189.28 12.2587 0.004372 **
P 1 8.40 8.40 0.5441 0.474904
K 1 95.20 95.20 6.1657 0.028795 *
N:P 1 21.28 21.28 1.3783 0.263165
N:K 1 33.14 33.14 2.1460 0.168648
P:K 1 0.48 0.48 0.0312 0.862752
Residuals 12 185.29 15.44
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>
ET les resultats du test de levene:
leveneaov<-aov((abs(aovrdt$residuals))~bl ock+N*P*K, data = levenetable)
> summary(leveneaov
+ )
Df Sum Sq Mean Sq F value Pr(>F)
block 5 11.4347 2.2869 1.0548 0.43096
N 1 0.8067 0.8067 0.3721 0.55325
P 1 0.2334 0.2334 0.1076 0.74850
K 1 17.3400 17.3400 7.9980 0.01523 *
N:P 1 0.7350 0.7350 0.3390 0.57118
N:K 1 0.1612 0.1612 0.0743 0.78976
P:K 1 0.8313 0.8313 0.3834 0.54735
Residuals 12 26.0167 2.1681
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>
-----
