data;
infile 'a:\stat706\ancova.dat';
input displace force brand $;
proc glm;
class brand;
model displace=force brand;
lsmeans brand/pdiff stderr;
run;