/*************  bread.sas ************************
 Subsampling Example from NKNW Ch. 28 (Table 28.9)
**************************************************/

options linesize = 79;

proc format;  /* value labels used in data step below */
     value tempfmt 1 = 'Low'  2 = 'Medium' 3 = 'High';
data school;
     infile 'bread.dat';
     input temp batch crusty;
     label crusty = 'Crustiness of bread';
     format temp tempfmt.;

proc glm;
     class temp batch;
     model crusty = temp batch(temp);
     random batch(temp) / test;
     means temp batch(temp);

proc glm; /* Pretending both effects random */
     class temp batch;
     model crusty = temp batch(temp);
     random temp batch(temp) / test;