## 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;