08-20-2013 06:00 PM
I want to set up doubly repeated measures using GLIMMIX but I haven't used GLIMMIX for repeated measures. The model so far I have now is
Proc glimmix data = wheat namelen=200 plots=residualpanel (conditional marginal);
class Year Block Plot Cultivar Nrate Sampling;
model DM = Nrate|Cultivar|Part|Sampling/ddfm=kr ;
Random Year Block(Year);
It's split-plot design with Nrate as whole plot (3 levels) and cultivar as split-plot (4 levels). Number of years are 3 and blocks (replications) are four. So, each year there were 48 plots. From each plot (48), plant samples were taken and divided into two "Part" (Head and stem). Their dry matter (DM) was measured. These samples were taken four times (Sampling) during the season.
The two parts (head and stem) have residuals effect and so is from sampling. Therefore Part and Sampling should be set as repeated measures but I don't know how to do it.
08-21-2013 03:27 PM
Check out Example 41.5 in the PROC GLIMMIX documentation. It gives a way to jointly model variables of different distributions. In your case, you can set up a joint model, but specify the same distribution. This is as close to the Kronecker doubly repeated measures method in MIXED as you are going to be able to do in GLIMMIX.
If the doubly repeated measures were nested (say week within month, or month and season), I'd say create a compound measure and sort out the comparisons using LSMESTIMATE statements.
08-22-2013 04:10 PM
Thanks for the direction. It seems the possible solution is to create multiple G-side effects but only one R-side random effect.