data a1;
input gen env y;
cards;
1 1 8
1 2 9
1 3 10
1 4 11
1 5 12
2 1 12
2 2 11
2 3 10
2 4 9
2 5 8
*given overall sums of squares for the interaction;
proc glm;classes gen env;
model y=env gen env*gen/ss1;

proc sort data=a1;by env;
proc means noprint;by env;var y;
output out=m1 mean=my css=sy;
data m2;set m1; sy=sqrt(sy);
proc means noprint data=m2;var sy;
output css=scalee;
proc print;
run;

proc sort data=a1; by gen;
proc means noprint;by gen;var y;
output out=m1 mean=my css=sy;
data m2;set m1; sy=sqrt(sy);
proc means noprint data=m2;var sy;
output css=scaleg;
proc print;
run;
quit;