options ps=60 ls=80 pageno=1 formdlim='_';

data birds;
input calcium hormone $ gender $;
cards;
32.0 yes male
23.8 yes male
28.8 yes male
25.0 yes male
29.3 yes male
39.1 yes female
26.2 yes female
21.3 yes female
35.8 yes female
40.2 yes female
14.5 no male
11.0 no male
10.8 no male
14.3 no male
10.0 no male
16.5 no female
18.4 no female
12.7 no female
14.0 no female
12.8 no female
;

proc glm;
class hormone gender;
model calcium = hormone gender hormone*gender;
means hormone;
lsmeans hormone/pdiff cl;
run;