options linesize = 80 ;
footnote "~john-c/5421/snedecor.cochran.sas &sysdate &systime" ;

data alfalfa ;
infile 'snedecor-cochran.table12.12.1' ;
     length variety $7 date $1 ;
     input variety date block yield ;

Ladak     A   1   2.17
Ladak     A   2   1.88
Ladak     A   3   1.62
Ladak     A   4   2.34
Ladak     A   5   1.58
Ladak     A   6   1.66
Ladak     B   1   1.58
Ladak     B   2   1.26
Ladak     B   3   1.22
Ladak     B   4   1.59
Ladak     B   5   1.25
Ladak     B   6   0.94
Ladak     C   1   2.29
Ladak     C   2   1.60
Ladak     C   3   1.67
Ladak     C   4   1.91
Ladak     C   5   1.39
Ladak     C   6   1.12
Ladak     D   1   2.23
Ladak     D   2   2.01
Ladak     D   3   1.82
Ladak     D   4   2.10
Ladak     D   5   1.66
Ladak     D   6   1.10
Cossack   A   1   2.33
Cossack   A   2   2.01
Cossack   A   3   1.70
Cossack   A   4   1.78
Cossack   A   5   1.42
Cossack   A   6   1.35
Cossack   B   1   1.38
Cossack   B   2   1.30
Cossack   B   3   1.85
Cossack   B   4   1.09
Cossack   B   5   1.13
Cossack   B   6   1.06
Cossack   C   1   1.86
Cossack   C   2   1.70
Cossack   C   3   1.81
Cossack   C   4   1.54
Cossack   C   5   1.67
Cossack   C   6   0.88
Cossack   D   1   2.27
Cossack   D   2   1.81
Cossack   D   3   2.01
Cossack   D   4   1.40
Cossack   D   5   1.31
Cossack   D   6   1.06
Ranger    A   1   1.75
Ranger    A   2   1.95
Ranger    A   3   2.13
Ranger    A   4   1.78
Ranger    A   5   1.31
Ranger    A   6   1.30
Ranger    B   1   1.52
Ranger    B   2   1.47
Ranger    B   3   1.80
Ranger    B   4   1.37
Ranger    B   5   1.01
Ranger    B   6   1.31
Ranger    C   1   1.55
Ranger    C   2   1.61
Ranger    C   3   1.82
Ranger    C   4   1.56
Ranger    C   5   1.23
Ranger    C   6   1.13
Ranger    D   1   1.56
Ranger    D   2   1.72
Ranger    D   3   1.99
Ranger    D   4   1.55
Ranger    D   5   1.51
Ranger    D   6   1.33

run ;

proc glm data = alfalfa ;
     class variety date block ;
     model yield = variety block variety*block date variety*date
                   / solution ;
title1 'PROC GLM for Data From Snedecor-Cochran, 6th Edition' ;
title2 'Table 12.12.1 on Alfalfa Yield: Split Plot' ;
run ;