1.  Ought to work.  I'm more used to the usual notation of A*B, rather than A(B) for the model and lsmeans statements, but they should work.  I am curious why you would use type3 analysis rather than the default REML with this design, though. 2.  A BY statement would give separate analyses by year.  Be sure the data are sorted appropriately. 3.  But phenology is a random effect, so if you want to do this, you will have to construct ESTIMATE statements to compare the BLUPs at each phenology.  This could become tedious. Steve Denham ... View more