Dear community,
I would really like to know if the way I wrote my lsmestimates is correct. I'm combining the data of 2 experiments in a weighted analysis. I have 6 treatments applied at rate 1, 2, 3 and 4. So I want to compare trt 2 vs trt 3 at rate 1, and so on. In this case my dependent variable is the fresh shoot weight of watermelons.
Proc glimmix data=one; class Exp Blk Trt Rate; model Fshoot_weight= Trt|Rate/dist=lognormal ddfm=kr; random intercept/subject=Blk; weight Fshoot_weight_wt; lsmestimate Trt*Rate 'Trt 2 vs 3 @ rate=1' 1 11 -11 1 1 1, 'Trt 2 vs 4 @ rate=1' 1 11 1 -11 1 1, 'Trt 2 vs 5 @ rate=1' 1 11 1 1 -11 1/adjust=t elsm ; run;
I would greatly appreciate if you could check if these estimates statements make sense.
Thanks a lot!
Caroline
... View more