options nocenter ps=100 ls=79;

data a;
   infile 'company.dat' firstobs=2;
   input steel site type $ reps score;

proc glm;
   class steel site type reps;
   model score = steel | site | type;
   means steel type steel*type / lsd;
run;