网友案例
options ls=72; goptions reset=all;
title "Fabrics";
Data Youden;
do row=1 to 3; do column = 1 to 7; input trt $ Y @@;
output; end; end;
cards;
A 57 B 60 C 39 D 42 E 36 F 27 G 69
B 72 C 39 D 33 E 45 F 45 G 81 A 54
D 21 E 24 F 27 G 63 A 51 B 60 C 27
;
ods output means =means lsmeans=lsmeans_FIX;
proc glm data=Youden; class row column trt;
model Y = row column trt/solution;
means trt; lsmeans trt/e pdiff stderr; run; quit;
proc gplot data=Youden;
plot Y*trt=column/legend; plot Y*column=trt/legend;
symbol1 v=dot i=none;
label row = "Machine"; label column = "Speed"; run;quit;
** suppose row=day, column = technician (random) **;
ods output lsmeans=lsmeans_RAN;
proc mixed data=Youden; class row column trt;
model Y = trt/solution;
random row column;
lsmeans trt/e pdiff;
run; quit;
|