Hello Rick, Thanks for your reply. I've attached the data. I tried both glimmix and proc gee, but I still got stuck. @Rick_SAS frt=fertilizer I collected data for 27 weeks for the same plants sprayed with different fertilizers, and I am looking to see if there is any significant difference between fertilizer treatments on plant height and biomass. Proc gee Data =Waterford; class frt week; model ph=frt; repeated subject=frt/type=ar(1) covb ; lsmeans frt*week/ cl diff plot=diff Adjust=Tukey lines; run; Error: Effects used in the LSMEANS statement must have appeared previously in the MODEL statement. Proc glimmix Data =Waterford; class frt week; model ph = frt week frt*week/link=log dist=multinomial; random week/subject=frt type=cs residual; lsmeans frt*week / Adjust = Tukey lines slicediff= week; output out=second predicted=pred residual=resid residual(noblup)=mresid student=studentresid student(noblup)=smresid; run; For the multinomial distribution only the following links are possible: CUMLOGIT, CUMPROBIT, CUMLOGLOG, CUMCLL, GLOGIT. Thank you again.
... View more