OPTIONS LS=72;
TITLE 'ANCOVA ANALYSIS OF SBP DATA';
DATA SBP;
   INPUT AGE SBP SEX;
CARDS;
39 144 1
45 138 1
47 145 1
65 162 1
46 142 1
67 170 1
42 124 1
67 158 1
56 154 1
64 162 1
56 150 1
59 140 1
34 110 1
42 128 1
48 130 1
45 135 1
17 114 1
20 116 1
19 124 1
36 136 1
50 142 1
39 120 1
21 120 1
44 160 1
53 158 1
63 144 1
29 130 1
25 125 1
69 175 1
41 158 0
60 185 0
41 152 0
47 159 0
66 176 0
47 156 0
68 184 0
43 138 0
68 172 0
57 168 0
65 176 0
57 164 0
61 154 0
36 124 0
44 142 0
50 144 0
47 149 0
19 128 0
22 130 0
21 138 0
38 150 0
52 156 0
41 134 0
18 134 0
51 174 0
55 174 0
65 158 0
33 144 0
23 139 0
70 180 0
56 165 0
62 172 0
51 160 0
48 157 0
59 170 0
40 153 0
35 148 0
33 140 0
26 132 0
61 169 0
;
PROC PRINT;

PROC PLOT DATA=SBP;
     PLOT SBP*AGE=SEX;
RUN;

PROC GLM DATA=SBP;
     MODEL SBP=AGE;
RUN;

DATA SBP1;
SET SBP;
    IF SEX=1;
PROC PRINT;
PROC MEANS;
PROC GLM DATA=SBP1;
     MODEL SBP=AGE;
RUN;

DATA SBP0;
SET SBP;
    IF SEX=0;
PROC PRINT;
PROC MEANS;
PROC GLM DATA=SBP0;
     MODEL SBP=AGE;
RUN;

PROC GLM DATA=SBP;
   CLASS SEX;
   MODEL SBP=SEX AGE*SEX /SS1;
   TITLE 'SEPARATE INTERCEPTS AND SEPARATE SLOPES, TESTING ZERO SLOPE';
RUN;

PROC GLM DATA=SBP;
   CLASS SEX;
   MODEL SBP=AGE SEX AGE*SEX / SS1 SS3;
   TITLE 'SEPARATE INTERCEPTS AND SEPARATE SLOPES, TESTING EQUAL SLOPE';
RUN;


PROC GLM DATA=SBP;
     CLASS SEX;
     MODEL SBP=AGE SEX / SS3;
     TITLE 'FITTING MODEL WITH SEPARATE INTERCEPTS AND EQUAL SLOPES';
RUN;

PROC GLM DATA=SBP;
     CLASS SEX;
     MODEL SBP=AGE SEX / SOLUTION;
     LSMEANS SEX / STDERR PDIFF COV OUT=ADJMEANS;
     TITLE 'FITTING MODEL WITH SEPARATE INTERCEPTS AND EQUAL SLOPES';
RUN;

PROC PRINT DATA=ADJMEANS;
RUN;

PROC GLM DATA=SBP;
     CLASS SEX;
     MODEL SBP=AGE SEX / P CLI CLM;
     TITLE 'FITTING MODEL WITH SEPARATE INTERCEPTS AND EQUAL SLOPES';
RUN;