*********************************************************** * STT-4100/STT-7230 * Planification des expériences * * Exemple du chapitre 6: Les plans hiérarchisés * Données tirées de Milliken & Johnson, 2009, p. 627. ***********************************************************; *********************** * Lecture des données * ***********************; data insect; input cie $ prod Nmoust @@; label cie="Compagnie"; label prod="Produit"; label Nmoust="Nombre d'insectes"; datalines; A 1 151 B 1 140 D 1 79 A 1 135 B 1 152 D 1 74 A 1 137 B 1 133 D 1 73 A 2 118 B 2 151 D 2 67 A 2 132 B 2 132 D 2 78 A 2 135 B 2 139 D 2 63 A 3 131 C 1 96 D 3 90 A 3 137 C 1 108 D 3 81 A 3 121 C 1 94 D 3 96 C 2 84 D 4 83 C 2 87 D 4 89 C 2 82 D 4 94 ; ***************************** * Diagramme de dispersion * *****************************; proc sgplot data=insect; scatter x=cie y=Nmoust/group=prod markerattrs=(symbol=Circlefilled); xaxis offsetmin=0.1 offsetmax=0.1; run; ***************************** * Analyse avec Proc GLM : * équivalente à Proc MIXED ici, * car aucun facteur aléatoire. *****************************; proc glm data=insect; class cie prod; model Nmoust= cie prod(cie); run; ************************************* * Analyse avec Proc MIXED * privilégiée dans tout le chapitre 6 *************************************; proc mixed data=insect; class cie prod; lsmeans cie*prod/pdiff; run; proc mixed data=insect covtest cl; class cie prod; *model Nmoust = cie*prod; * anova à un facteur pour comparer les 11 produits; model Nmoust = cie prod(cie)/e3 outpm=out vciry; * résidus marginaux; *model Nmoust = cie prod(cie)/ e3 solution outp=out2; * résidus conditionnels; lsmeans cie/pdiff; lsmeans prod(cie)/ pdiff; lsmeans prod(cie)/slice=cie pdiff; run; ************************************* * Analyse des résidus: * outpm et outp donneront la même chose ici * car il n'y a pas de facteur aléatoire, * donc pas de corrélation à enlever; *************************************; ods select TestsForNormality; proc univariate data=out normal; var ScaledResid; run;