turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Proc Mixed - Covary a baseline value

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

04-27-2017 11:37 PM

I'm hoping someone might be able to help me adjust the lsmeans for a baseline value. That is, DV was measured 8 times: 0 min, 15 min, 30 min, 45 min, 60 min, 120 min, and 90 min, and we'd like to adjust the lsmeans for the baseline value (0 min). Is this possible?

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to eileenwright

04-28-2017 12:07 AM

With option / AT time=0 ?

PG

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to eileenwright

05-03-2017 01:44 AM

An analysis of covariance model might work, if what you want is to estimate the means at 15, 30, 45, 60, 90, and 120 minutes for a common value of baseline. (The baseline time=0 value would be the covariate.)

Because (presumably) of the repeated measures on a subject at 0, ..., 120 min, the model would be more complicated than your standard ANCOVA (which would have data for time=0 and time= a second value). In particular, the relationship between data at time 15 and time 0 might be stronger than the relationship between data at time 120 and time 0 because noise intrudes as time passes--in other words, the slope of the regression of the response on baseline (time=0) might decrease with later times.

This is speculative and untested, but I would consider

proc glimmix data=have;

class time subjectID;

model response = time baseline time*baseline;

random time / subject=subjectID type=<some covariance structure, maybe ar(1)> residual;

lsmeans time / at mean; /* or some other value of baseline */

run;

Issues to consider are

-- the nature of the relationship between response and baseline at each time (e.g., linear)

-- an appropriate covariance structure for the repeated measures within a subject

-- normality and homogeneity of variance (assuming normal distribution)

-- what a sensible "common" value for baseline might be, in context

Alternatively, you could compute a variable that represents deviance from baseline, either absolute (i.e., subtract the baseline value) or relative (e.g., divide by the baseline value). Or maybe a random coefficients model.