[R-sig-ME] Model specification (fixed and random effects)?

ONKELINX, Thierry Thierry.ONKELINX at inbo.be
Tue Apr 22 09:56:03 CEST 2014

Dear Julia,

Your model specification is OK if you assume that there is a linair trend along number of surgeries and tests.

The correlation should be  corAR1(form=~test_time|id/surgery_number)

Best regards,

ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht
+ 32 2 525 02 51
+ 32 54 43 61 85
Thierry.Onkelinx op inbo.be

To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of.
~ Sir Ronald Aylmer Fisher

The plural of anecdote is not data.
~ Roger Brinner

The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

-----Oorspronkelijk bericht-----
Van: r-sig-mixed-models-bounces op r-project.org [mailto:r-sig-mixed-models-bounces op r-project.org] Namens Yla Savh
Verzonden: dinsdag 22 april 2014 9:07
Aan: r-sig-mixed-models op r-project.org
Onderwerp: [R-sig-ME] Model specification (fixed and random effects)?


I would appreciate any help with my model design. Each subject had a few surgeries (from 1 to 5) ("surgery number"). These subjects had a test taken a few times after each surgery (from 1 to 6 times) ("test time"). There was supposed to be a decrease in the dependent variable (test_results) immediately after each surgery.  I coded both variables as continuous (surgery_number and test_time).

I wonder if I specified fixed and random effects correctly. I also think there might be some correlation between the surgery numbers but I am not sure about how to specify it in the model.

x<-lme(test_results~ surgery_number*test_time, random=~1|id/surgery_number) s12<-update(x, correlation=corAR1(form=~1|id/surgery_number))


        [[alternative HTML version deleted]]

* * * * * * * * * * * * * D I S C L A I M E R * * * * * * * * * * * * *
Dit bericht en eventuele bijlagen geven enkel de visie van de schrijver weer en binden het INBO onder geen enkel beding, zolang dit bericht niet bevestigd is door een geldig ondertekend document.
The views expressed in this message and any annex are purely those of the writer and may not be regarded as stating an official position of INBO, as long as the message is not confirmed by a duly signed document.

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