Programming the statistical procedures from SAS

Multiple repeated measures glimmix

Accepted Solution Solved
Reply
New Contributor
Posts: 3
Accepted Solution

Multiple repeated measures glimmix

I'm trying to run an analysis in glimmix of some ecological data thats throwing me for a loop, but I'm probabaly overthinking it. The data were collected over 4 years, at 4 blocks where each block contained 1 of 4 different vegetative types. Because the same plots were visited each year, I know I can account for potential covaraince at this level using a random effect of study plot or as a repeated measure using one of the following

 

random int/ subject=plot;

OR

random year/ subject=plot;

 

However, these plots were also visited twice within a given year. How do I account for this second level of repeated-ness? Is it just as simple as :

 

random visit/ subject=plot*yr?

 

and can this kind of statement be combined with one of the above?


Accepted Solutions
Solution
‎01-05-2017 02:05 PM
Respected Advisor
Posts: 2,655

Re: Multiple repeated measures glimmix

Here's an idea.  Throw out year.  If you have two visits to a plot per year, then fit visit 1 through 8 as a single repeated measure.  You can then look at 'year' effects by combining adjacent visits in an LSMESTIMATE statement.

 

And if the visits are not equally spaced in time, consider fitting a spline to both the fixed effect and as a random effect (see example 45.15 Comparing Multiple B Splines and Example 45.6 Radial Smoothing of Repeated Measures Data.

 

Steve Denham

View solution in original post


All Replies
Respected Advisor
Posts: 4,754

Re: Multiple repeated measures glimmix

Shouldn't visits be called season instead and considered as a fixed effect?

PG
New Contributor
Posts: 3

Re: Multiple repeated measures glimmix

I don't think so, they're repeated measurements of the same sample unit. But I'm open to it if you can explain why.
Respected Advisor
Posts: 4,754

Re: Multiple repeated measures glimmix

Well, I am missing some context here.

 

But ecological systems typically got through a cycle during the year. Observations from the Spring may be systematically different from observations from the Fall. Random parameter distributions available in GLIMMIX are unimodal, they won't fit too well. If your visits occurred at different moments of a systematic cycle, then it would be preferable to model that cycle instead of considering them as independent random realizations of a single distribution.  

PG
New Contributor
Posts: 3

Re: Multiple repeated measures glimmix

I think in this case I can assume that they're not systematically different. The data are sort of a maximum number of unique small mammals trapped during each survey. In this case, we know individuals trapped early in the year show up in surveys later in the year, and that fall abundances ought to be larger than spring abundances because there is the potential for multiple litters. For this reason, I know there is covariance among the two samples within a year and and I am just searching for a way to model that out, in addition to the covariance associated with returning to the exact same plots year after year.

Solution
‎01-05-2017 02:05 PM
Respected Advisor
Posts: 2,655

Re: Multiple repeated measures glimmix

Here's an idea.  Throw out year.  If you have two visits to a plot per year, then fit visit 1 through 8 as a single repeated measure.  You can then look at 'year' effects by combining adjacent visits in an LSMESTIMATE statement.

 

And if the visits are not equally spaced in time, consider fitting a spline to both the fixed effect and as a random effect (see example 45.15 Comparing Multiple B Splines and Example 45.6 Radial Smoothing of Repeated Measures Data.

 

Steve Denham

☑ This topic is solved.

Need further help from the community? Please ask a new question.

Discussion stats
  • 5 replies
  • 249 views
  • 0 likes
  • 3 in conversation