Example 4 for PROC ANOVA
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: ANOVAEX4 */
/* TITLE: Example 4 for PROC ANOVA */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: analysis of variance, balanced data, design */
/* PROCS: ANOVA */
/* DATA: */
/* */
/* SUPPORT: sasrdt */
/* REF: PROC ANOVA, EXAMPLE 4. */
/* MISC: */
/****************************************************************/
/* Latin Square Split Plot Design ------------------------------*/
title1 'Sugar Beet Varieties';
title3 'Latin Square Split-Plot Design';
data Beets;
do Harvest=1 to 2;
do Rep=1 to 6;
do Column=1 to 6;
input Variety Y @;
output;
end;
end;
end;
datalines;
3 19.1 6 18.3 5 19.6 1 18.6 2 18.2 4 18.5
6 18.1 2 19.5 4 17.6 3 18.7 1 18.7 5 19.9
1 18.1 5 20.2 6 18.5 4 20.1 3 18.6 2 19.2
2 19.1 3 18.8 1 18.7 5 20.2 4 18.6 6 18.5
4 17.5 1 18.1 2 18.7 6 18.2 5 20.4 3 18.5
5 17.7 4 17.8 3 17.4 2 17.0 6 17.6 1 17.6
3 16.2 6 17.0 5 18.1 1 16.6 2 17.7 4 16.3
6 16.0 2 15.3 4 16.0 3 17.1 1 16.5 5 17.6
1 16.5 5 18.1 6 16.7 4 16.2 3 16.7 2 17.3
2 17.5 3 16.0 1 16.4 5 18.0 4 16.6 6 16.1
4 15.7 1 16.1 2 16.7 6 16.3 5 17.8 3 16.2
5 18.3 4 16.6 3 16.4 2 17.6 6 17.1 1 16.5
;
/* Harvest: Split Plot on Original Latin Square for Whole Plots */
proc anova data=Beets;
class Column Rep Variety Harvest;
model Y=Rep Column Variety Rep*Column*Variety
Harvest Harvest*Rep
Harvest*Variety;
test h=Rep Column Variety e=Rep*Column*Variety;
test h=Harvest e=Harvest*Rep;
run;