Hi All, I have a dataset with 51 genotypes and 2 treatments. I fit a model with proc mixed as proc mixed data=data ; class treatment genotype ; model trait = treatment genotype treatment*genotype / htype=1 ddfm=satterth ; lsmestimate treatment*genotype '(geno1-geno2)trt1 vs (geno1-geno2)trt2' 1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0; lsmestimate treatment*genotype '(geno1-geno3)trt1 vs (geno1-geno3)trt2' 1 0 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0; run; I would like to compare the difference between one genotype with the rest of the genotypes in treatment1 and the difference between one genotype with the rest of the genotypes in treatment2 Ho: (geno1 - geno2)treatment1 = (geno1 - geno2)treatment2 Is there any way to do this without having to write hundreds of individual constrasts? Thanks.
... View more