je dois enseigner les tests en BTS agricole et je suis en train de m'y initier (et tout seul, c'est un peu dur ;-)).
Je tente de le faire à l'aide de R.
J'ai pris les données (essais de 7 modalités avec 4 répétitions (blocs)) dans le rapport d'un ancien élève et je tente de refaire son analyse d'expérimentation.
Les données sont ici (séparateur de données ; et séparateur décimal ,):
http://dl.dropbox.com/u/1478925/expe.csv
Un extrait du rapport là :
http://dl.dropbox.com/u/1478925/expe.pdf
L'élève a utilisé StatBox. J'aimerai trouver les équivalents R et Rcmdr
J'ai converti les données numériques Modalité et bloc en facteurs.
L'analyse commence par :
INTERACTION TRAITEMENTS*BLOCS
SCE test de TUKEY = 24.676 PROBA = 0.340022
Comment faire ce test ?
Plus loin, on teste la normalité des résidus avec :
INDICES DE NORMALITE (coefficients de K.PEARSON)
SYMETRIE (valeur idéale théorique = 0) : BETA 1 = 0.017 PROB : 0.76596
APLATISSEMENT (valeur idéale théorique = 3) : BETA 2 = 2.726 PROBA : 0.74979
Comment faire ce test ?
Puis, on teste si les résidus par modalité puis par bloc sont égaux avec :
ECARTS-TYPES DES RESIDUS
ECARTS-TYPES FACTEUR 1 = F1
1 (F1n1) 2 (F1n2) 3 (F1n3) 4 (F1n4) 5 (F1n5) 6 (F1n6) 7 (F1n7)
3,45 5,9 1,12 5,1 7,17 3,86 3,45
KHI2 = 7.771 PROB =0.25452
ECARTS-TYPES BLOCS = BLOC
1 (b1) 2 (b2) 3 (b3) 4 (b4)
4,59 3,22 3,97 5,38
KHI2 = 1.547 PROB =0.67572
Comment faire ce test ?
Puis j'ai essayé de faire le test de comparaisons des moyennes.
Lorsque je demande une comparaison de moyenne (ANOVA à deux facteurs), Rcmdr écrit :
Code : Tout sélectionner
AnovaModel.2 <- (lm(TonnageParHa ~ Bloc*Modalité, data=Dataset))
Anova(AnovaModel.2)
La deuxième ligne provoque une erreur :
Code : Tout sélectionner
NA/NaN/Inf dans appel à une fonction externe (argument 1)
Si je mets un a minuscule, j'ai (et cela ressemble à ce que je cherche) :
Code : Tout sélectionner
> anova(AnovaModel.2)
Analysis of Variance Table
Response: TonnageParHa
Df Sum Sq Mean Sq F value Pr(>F)
Bloc 3 74.09 24.70
Modalité 6 128.90 21.48
Bloc:Modalité 18 456.71 25.37
Residuals 0 0.00
mais je n'ai pas les valeurs F value et Pr(>F)...
Pourtant, je les obtiens avec une ANOVA à un facteur, mais pas pour la ligne Bloc:Modalité. Est-ce la valeur du premier test (INTERACTION TRAITEMENTS*BLOCS, SCE test de TUKEY = 24.676 PROBA = 0.340022)
à la fin de cet exemple, on trouve la puissance du test.
comment l'obtenir ?
Il y a encore plein d'essais dans ce rapport, je pense que vos réponses pourraient me mettre sur la voie pour les comprendre.
Merci d'avance pour toutes vos idées, conseils et réponses,
Fabrice Leclercq.