data a1; input grp monwalk; cards; 1 9.00 1 9.50 1 9.75 1 10.00 1 13.00 1 9.50 2 11.00 2 10.00 2 10.00 2 11.75 2 10.50 2 15.00 3 11.50 3 12.00 3 9.00 3 11.50 3 13.25 3 13.00 4 13.25 4 11.50 4 12.00 4 13.50 4 11.50 ; run; data a2; set a1; if grp=4 then group='4: control'; if grp=1 then group='1: active'; if grp=2 then group='2: passive'; if grp=3 then group='3: none'; if grp=1 then treat=1; else treat=0; run; proc gplot data=a2; plot monwalk*group / haxis=axis1 vaxis=axis2 frame; axis1 offset=(3,3) label=(H=3) value=(H=2) minor=NONE; axis2 offset=(1,1) value=(H=2) minor=NONE label=(A=90 R=0 H=3); symbol1 v=circle c=black i=none l=1 h=2 r=1; run; proc glm data=a2; class group; model monwalk=group; *means group / hovtest=bartlett; *means group / hovtest=bartlett t bon; *lsmeans group / adjust=bon pdiff cl; means group / bon; run; proc ttest data=a2; class treat; var monwalk; run; proc glm data=a2; class treat; model monwalk=treat / solution; run; proc npar1way wilcoxon data=a2; class group; var monwalk; run; proc ttest data=a2; where grp=1 or grp=4; class group; var monwalk; run;