/* TWO SAMPLE PROFILE EXAMPLE */

DM 'LOG;CLEAR;OUT;CLEAR;';
Options ls=79 Compress = No sasautos='~/gauss/sas_macros/';

/* Data from Morrison, Page 195 */

DATA Schizo_raw;
  INFILE 'schizo.dat';
  INPUT grp a1-a4;
  acid1=sqrt(a1); acid2=sqrt(a2); acid3=sqrt(a3); acid4=sqrt(a4);
  psi=(acid2+acid3)/2-acid4;

%transform(data = Schizo_raw, response=a1-a4, class = grp, explanatory = grp,
    levels=4, polydeg=3); 


proc glm data = Schizo_raw;
  class grp;
  model acid1-acid4 = grp/nouni;
  lsmeans grp;
  repeated Time 4 (0 15 30 45) polynomial/printe printm;
  manova H=grp M=(1 -.5 -.5 0);
run;

proc glm data = Schizo_raw;
    class grp;
    model psi=grp;
    estimate 'psi diff' grp 1 -1;
run;