options nocenter ls=78 ps=120;

data brouwer;
infile 'brouwer.dat' missover;
input loghet gen s1 $ s2 $ i2 s3 $ i3 s4 $ i4;
  if gen=4 then delete;
  
proc glm;
  classes s1 s2 s3;
  model loghet = gen s1 i2*s2 i3*s3/solution;
  random s1/test;