data IleumSucraseActivity;
length treatment $20;
input treatment$ block IleumKm IleumVmax;
Non-Antibiotics 5 46.18 28.06
Non-Antibiotics 3 55.61 24.64
Non-Antibiotics 4 58.71 27.35
Non-Antibiotics 3 47.32 22.15
Non-Antibiotics 4 43.61 12.53
Non-Antibiotics 8 52.00 25.94
Non-Antibiotics 11 54.65 21.34
Non-Antibiotics 5 64.54 25.73
Non-Antibiotics 10 51.27 24.36
Non-Antibiotics 5 57.61 23.04
Non-Antibiotics 5 48.37 32.42
Non-Antibiotics 9 56.04 40.59
Non-Antibiotics 6 53.89 11.33
Antibiotics 8 39.02 27.74
Antibiotics 4 55.56 16.92
Antibiotics 7 56.97 25.41
Antibiotics 9 58.29 31.49
Antibiotics 4 53.18 28.04
Antibiotics 11 41.89 43.48
Antibiotics 5 40.35 38.12
Antibiotics 3 54.84 40.96
Antibiotics 10 47.66 32.51
Antibiotics 6 61.14 35.47
Antibiotics 5 52.42 38.02
;
run;
Proc print data=Activity;
run;
proc glimmix data=Activity nobound;
class treatment block;
model IleumKm = treatment;
random block/subject=treatment*block type=cs residual;
random _residual_/subject=treatment*block group=treatment;
lsmeans treatment/tdiff lines;
covtest homogeneity;
run;
结果:
ERROR: QUANEW Optimization cannot be completed.