A very big thank you both @SteveDenham and @lvm for your advices. I have just revised my codes, and could then you want to help to correct my codes below. However, how to compare each treatment to others ? meanwhile I can't use Duncan in Proc Mixed, can you ? If so, I propose to use "adjust = DUNNET" to compare with control, am I wrong ? but I didn't understand how to interpret it in my paper with a,b,c, notation, it's looked very different results of "means A B A*B / DUNCAN e=A(Block)" in Proc Anova. Please help me... data sp; input A $ B $ Block Y; datalines; A1 B1 1 72 A1 B1 2 63 A1 B2 1 45 A1 B2 2 30 A1 B3 1 4 A1 B3 2 7 A2 B1 1 23 A2 B1 2 27 A2 B2 1 85 A2 B2 2 65 A2 B3 1 55 A2 B3 2 58 A3 B1 1 5 A3 B1 2 8 A3 B2 1 23 A3 B2 2 21 A3 B3 1 15 A3 B3 2 19 A4 B1 1 30 A4 B1 2 39 A4 B2 1 42 A4 B2 2 44 A4 B3 1 12 A4 B3 2 18 ; proc mixed; class A B Block; model Y = A B A*B; random Block A*Block; lsmeans A|B / adjust=DUNNETT; run;
... View more