options ls=80 nodate;

data teacher;
input y x1 x2 x3 x4 gender $;
z=(gender='f');
cards;
489 81 151 45.50 43.61 m
423 68 156 46.45 44.69 m
507 80 165 76.50 54.57 m
467 107 149 55.50 43.27 m
340 43 134 49.40 49.21 m
524 129 163 72.00 49.96 m
488 139 159 86.20 53.05 m
445 88 135 64.00 49.51 m
388 99 141 44.15 39.57 m
410 69 125 59.00 55.66 f
569 57 131 31.75 63.97 f
425 77 141 80.50 45.32 f
344 81 122 75.50 46.67 f
324 19 141 49.00 41.21 f
505 53 152 49.35 43.83 f
235 77 141 60.75 41.61 f
501 76 132 41.25 64.57 f
600 65 157 50.75 42.41 f
;

proc print;
title 'teacher effectiveness data';
run;

proc glm;
model y = x1 x2 x3 x4 z/xpx i;
run;