网友案例
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;