Thanks for your fast reply, PG It was supposed to test whether there were significant changes in behaviour in each group between the baseline and intervention period. I did not receive explanations why the GLIMMIX procedure was most fitted for the analysis, but I guess it is because the data were non-parametric, it involved repeated measures, there were differences in sample sizes (4, 6 and 7 animals) and because not every subject was observed equally often, and the length in baseline and intervention periods differed. This is the code I actually run: data dataset; set dataset; run; Proc univariate plot normal; var (behaviour x) run; Proc glimmix data=dataset; class chimp cgroup intervention period; model behaviour x = intervention*period intervention period/ link=logit dist=binomial ddfm=kr; random _residual_; random intercept / subject=chimp(cgroup) type=un; lsmeans intervention*period intervention period/pdiff adjust=bon; lsmeans intervention*period intervention period/cl ilink; run;
... View more