[R] Linear Mixed Models with nlme, more than one random effect

[Giulio Di Giovanni]

Thanks a lot !!!

I was about to reply to my previous email that I found a solution via lme4 !
But really thanks.




>From: joris.dewolf at cropdesign.com
>To: "Giulio Di Giovanni" <perimessaggini at hotmail.com>
>CC: r-help at stat.math.ethz.ch, r-help-bounces at stat.math.ethz.ch
>Subject: Re: [R] Linear Mixed Models with nlme, more than one random effect
>Date: Wed, 20 Jun 2007 15:17:42 +0200
>
>
>
>Guilio,
>
>Have a look at Rnew volume 5/1 (http://cran.r-project.org/doc/Rnews/) where
>Doug Bates explains this nicely. Condider using lme4 for your purpose.
>
>Joris
>
>
>
>
>
>
>
>
>
>              "Giulio Di
>              Giovanni"
>              <perimessaggini at h                                          To
>              otmail.com>               r-help at stat.math.ethz.ch
>              Sent by:                                                   cc
>              r-help-bounces at st
>              at.math.ethz.ch                                       Subject
>                                        [R] Linear Mixed Models with nlme,
>                                        more than one random effect
>              20/06/2007 14:09
>
>
>
>
>
>
>
>
>
>
>Hi, I' trying to learn how to use lme for Linear Mixed Models  and I have a
>
>problem when I have to include more than one random effect in my model. I
>know that this could be a stupid question to ask, but I'm not able to solve
>
>it by myself... One example: if my model is
>
>response = operator + block + day
>with operator and block as fixed effects and day as random effect, I use
>
>res.lme <- lme(resp ~ oper + block , random=~1|day)
>
>If I want to include also another random effect, as "experiment", what I
>should do ?
>This effect doesn't have to be nested, at the and I would like to have the
>COV matrix using (if I'm not wrong) getVarCov function.
>
>Thanks in advance for any help or suggestions, I'm a beginner on this
>field...
>
>Davide
>
>_________________________________________________________________
>Cinema, Tv, Gossip e Orsoscopo…Tutto su MSN intrattenimento!
>
>______________________________________________
>R-help at stat.math.ethz.ch mailing list
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide
>http://www.R-project.org/posting-guide.html
>and provide commented, minimal, self-contained, reproducible code.

_________________________________________________________________
Calcio, Quiz, Sudoku, Scacchi
 Inizia la sfida su Messenger, GRATIS!