Re: [R] model specification using lme

2016-05-30 Thread Thierry Onkelinx
Dear Hanna, None of the models are correct is you want the same random intercept for the different methods but different random slope per method. You can random = ~ 1 + time:method | individual The easiest way to get alpha_0 and tau_i is to apply post-hoc contrasts. That is fairly easy to do

Re: [R] model specification using lme

2016-05-29 Thread Ben Bolker
li li gmail.com> writes: > > Hi all, > For the following data, I consider the following random intercept and > random slope model. Denote as y_ijk the response value from *j*th > individual within *i*th method at time point *k*. Assume the following > model for y_ijk: > > y_ijk=

[R] model specification using lme

2016-05-29 Thread li li
Hi all, For the following data, I consider the following random intercept and random slope model. Denote as y_ijk the response value from *j*th individual within *i*th method at time point *k*. Assume the following model for y_ijk: y_ijk= (alpha_0+ tau_i +a_j(i))+(beta_i+b_j(i)) T_k +