[R-sig-ME] fixing the value of some parameters in an lme or lmer model

Asher Strauss asher.strauss at gmail.com
Wed Mar 29 10:39:46 CEST 2017

Dear list,

I am trying to fit a two level hierarchical model using the lme/lmer
I have two levels in my data: 3 time points (codded: week 0, week 4, week
8) nested within subjects (ID). I would like to set the intercept of each
patient to the value of the outcome at week 0 (the first observation), or
in other words force the model fit for each subject to pass through the
subject's baseline score on the outcome variable.

Here are the formulas (in informal notation):

Level 1: (time)

Level 2: (subjects)

I underlined and made bold the values I would like to force.

Can this be done in nlme or lme4 packages?

Thank you very much!!!

Asher Strauss

	[[alternative HTML version deleted]]

More information about the R-sig-mixed-models mailing list