data glmxiyin5_12;
	input s $ d $ n $ y @@;
	label s="品種" d="密度" n="氮肥" y="小區產量";
cards;
s1 2500 n1 18.4 s1 2500 n1 20.3 s1 2500 n1 19.7
s1 2500 n2 29.1 s1 2500 n2 31.6 s1 2500 n2 30.7
s1 3500 n1 23.7 s1 3500 n1 27.1 s1 3500 n1 25.4
s1 3500 n2 38.9 s1 3500 n2 42.3 s1 3500 n2 40.6
s2 2500 n1 24.7 s2 2500 n1 26.2 s2 2500 n1 25.3
s2 2500 n2 38.9 s2 2500 n2 37.6 s2 2500 n2 39.1
s2 3500 n1 31.3 s2 3500 n1 34.5 s2 3500 n1 33.8
s2 3500 n2 50.4 s2 3500 n2 52.3 s2 3500 n2 48.7
;
run;
proc anova;
	class s d n;
	model y=s d n s*d s*n d*n s*d*n;
	means s d n s*d s*n d*n s*d*n/bon;
run;
proc glm;
	class s d n;
	model y=s d n s*d s*n d*n;
    contrast '（2500-3500）/s1 '  d  1 -1  s*d   1  -1 ;
    contrast '（2500-3500）/s2 '  d  1 -1  s*d   0   0  1 -1 ;
    contrast '（n1-n2）/s1 '      n  1 -1  s*n   1  -1 ;
    contrast '（n1-n2）/s2 '      n  1 -1  s*n   0   0  1 -1 ;
    contrast '（n1-n2）/2500 '    n  1 -1  d*n   1  -1 ;
    contrast '（n1-n2）/3500 '    n  1 -1  d*n   0   0  1 -1 ;
run;
quit;