/* Example of Two-Period Crossover MJ 32.3 */
options linesize=80 pagesize=50;
data cross;
   infile 'cross2.dat';
   input seq person a b;
data b; set cross;
   trt = 'A'; y = a; output;
   trt = 'B'; y = b; output;
   drop a b;
proc glm;
   class trt seq person;
   model y = seq person(seq) trt trt*seq;
   random person(seq) / test;
   test h=seq e=person(seq);
   lsmeans seq / e=person(seq);
   lsmeans trt trt*seq;
proc mixed;
   class trt seq person;
   model y = seq trt trt*seq / ddfm = satterth;
   random person(seq);
   lsmeans seq trt;
   lsmeans trt*seq / pdiff;