Libname quirk 'k:\sboker\temp\GSA2001\';

title1 '***********************************************';
title2 ' Univariate statistics for GSAMULT             ';
title3 '***********************************************';

PROC MEANS DATA=quirk.gsamult N MEAN STD MIN MAX;
RUN;

title1 '***********************************************';
title2 ' Univariate statistics for GSAMULT3            ';
title3 '***********************************************';

PROC MEANS DATA=quirk.gsamult3 N MEAN STD MIN MAX;
RUN;

title1 '***********************************************';
title2 ' Univariate statistics for GSADATA2            ';
title3 '***********************************************';

PROC MEANS DATA=quirk.gsadata2 N MEAN STD MIN MAX;
RUN;

DATA quirk.iqincome;
    SET quirk.gsamult;
    ciq = y1-110;
    agecat = -1;
    IF age > 75 THEN agecat = 1;
    income = 50000 + (2000 * ciq) + (2000 * eit);
    IF agecat = 1 THEN income = 50000 - (2000 * ciq) + (2000 * eit);
    KEEP  
       id ciq income age agecat eit;
RUN;

PROC SORT DATA=quirk.iqincome;
    BY agecat;
RUN;

PROC MEANS DATA=quirk.iqincome N MEAN STD MIN MAX;
    BY agecat;
RUN;

PROC GLM DATA=quirk.iqincome;
     CLASS agecat;
     MODEL income = ciq agecat;
RUN;

PROC MIXED DATA=quirk.iqincome;
  MODEL income=ciq / S;
  RANDOM int ciq / SUB=age TYPE=un gcorr;
RUN;


PROC GLM DATA=quirk.iqincome;
     CLASS agecat;
     MODEL income = ciq agecat ciq*agecat;
RUN;

PROC MIXED DATA=quirk.iqincome;
  MODEL income=ciq agecat ciq*agecat / S;
  RANDOM int ciq / SUB=age TYPE=un gcorr;
RUN;