Hi, I have the following question. I've carried out a study in which 96 animals were involved. They were spread over four rounds, so in each round 24 animals were followed up. In each round, animals were spread over 4 treatments, so in each round I had 6 animals per treatment. Treatments in the four rounds were the same. Each round lasted for two weeks and animals were observed on day 2, 4, 6, 8, 10 and 12. I want to use proc glimmix to verify differences between the four treatments. However, I have repeated measures for each animals (6 different days) and I'm not sure how to include this. At first I used the following procedure: Proc glimmix data=dataset; class animal treatment day; model y = treatment day/dist=binary; random animal; lsmeans treatment/pdiff; run; However, in this way I'm not taking into account the repeated measures per animal. I feel like I need to change the random statement but I'm not sure how. Do you have any suggestions?
... View more