*********************************************************** * STT-4100/STT-7230 * Planification des expériences * * Exemple 3 du chapitre 6: Les plans hiérarchisés * Données tirées de Milliken & Johnson, 2009, p. 147. ***********************************************************; *********************** * Lecture des données * ***********************; data confort; input temp chambre sexe2 confort @@; label temp ="Température"; label confort ="Niveau de confort"; if sexe2=1 then sexe="Homme"; if sexe2=2 then sexe="Femme"; drop sexe2; datalines; 15 1 1 5 20 1 1 8 25 1 1 12 15 1 1 4 20 1 1 8 25 1 1 8 15 1 2 1 20 1 2 10 25 1 2 11 15 1 2 2 20 1 2 7 25 1 2 13 15 2 1 5 20 2 1 6 25 2 1 8 15 2 1 4 20 2 1 3 25 2 1 7 15 2 2 5 20 2 2 8 25 2 2 8 15 2 2 5 20 2 2 8 25 2 2 8 15 3 1 4 20 3 1 5 25 3 1 6 15 3 1 2 20 3 1 7 25 3 1 6 15 3 2 1 20 3 2 8 25 3 2 6 15 3 2 3 20 3 2 8 25 3 2 7 ; *************************** * Diagramme de dispersion * ***************************; proc sgplot data=confort; scatter x=temp y=confort/group=sexe markerattrs=(symbol=Circlefilled); xaxis offsetmin=0.1 offsetmax=0.1; run; ************************************************* * Courbes d'interaction pour les facteurs fixes *************************************************; proc sort data=confort; by temp sexe; proc means data=confort noprint; by temp sexe; var confort; output out=out mean=conf_moy; run; title "Courbes d'interaction entre les facteurs sexe et température"; proc sgplot data=out; series x=temp y=conf_moy/group=sexe markers markerattrs=(symbol=Circlefilled); xaxis offsetmin=0.1 offsetmax=0.1; run; title ; ************************************* * Analyse avec Proc MIXED *************************************; proc sort data=confort; by chambre; proc mixed data=confort covtest cl; class temp chambre sexe; model confort= temp sexe temp*sexe/e3 outpm=out vciry; random chambre(temp)/v ; lsmeans temp*sexe/ slice=temp pdiff; run; * pour le test du rapport de vraisemblances sur composante de variance; proc mixed data=confort ; class temp sexe; model confort= temp sexe temp*sexe; run; ************************************* * Analyse des residus; *************************************; ods select TestsForNormality; proc univariate data=out normal; var ScaledResid; run;