filename hdcp 'h:\HAL\Courses\Stat210\handicap.txt';
data discrim;
    infile hdcp firstobs=2;
    input score code;
proc sort;
    by code;
proc gplot;
    plot score*code;
proc means mean stddev min max n;
    var score;
    by code;
proc glm order=data;
    class code;
    model score = code;
    means code / lsd bon tukey scheffe;
    contrast 'all' code 4 -1 -1 -1 -1;
run;