data example5_3;
input lit dos med y;
cards;
1	1	3	2.33
1	2	2	4.33
1	3	1	2.94
1	4	4	5.11
2	1	2	2.78
2	2	1	2.72
2	3	4	4.33
2	4	3	3.06
3	1	1	1.88
3	2	4	3.56
3	3	3	2.86
3	4	2	3.89
4	1	4	4.54
4	2	3	2.28
4	3	2	4.39
4	4	1	2.72
;
/*正态性检验*/
proc univariate normal;
class med;
var y;
run;
/*方差齐性检验*/
proc glm;
class lit dos med ;
model y=med;
means med/hovtest;
run;
/*方差分析*/
proc glm;
class lit dos med ;
model y=med dos lit;
means med/hovtest;
lsmeans med/tdiff adjust=bon;
run;
