Re: [R] crossed random fx nlme lme4

From: Emilio A. Laca <ealaca_at_ucdavis.edu>
Date: Mon 18 Jul 2005 - 07:32:38 EST

random = pdBlocked( list( pdIdent( ~ y), pdIdent( ~ observer - 1), pdIdent( ~ set - 1) )

gave me output only when the data was a groupedData. The results are different depending on whether I specify observer or set as grouping factor.
I am missing something ...

Simon, Prof Bates, thanks for taking the time to reply.

eal

On Jul 15, 2005, at 6:35 AM, Douglas Bates wrote:

> On 7/13/05, Simon Blomberg <blomsp@ozemail.com.au> wrote:
>
>> At 09:35 AM 14/07/2005, Emilio A. Laca wrote:
>>
>>> I need to specify a model similar to this
>>>
>>> lme.formula(fixed = sqrt(lbPerAc) ~ y + season + y:season, data =
>>> cy,
>>> random = ~y | observer/set, correlation = corARMA(q = 6))
>>>
>>> except that observer and set are actually crossed instead of nested.
>>>
>>
>> Does this work for you? (following P&B pp 162-3 and an R-help archive
>> search on "crossed random effects")...
>>
>> fit <- lme(sqrt(lbPerAc) ~ y * season, random=list(pdBlocked
>> (pdIdent(~y),
>> pdIdent(observer-1), pdIdent(set-1))), correlation=corARMA(q = 6),
>> data=cy)
>>
>> lme isn't very well set up for crossed random effects. It's easier
>> in lmer.
>> I don't think lmer can handle alternative correlation structures yet,
>> though. (Prof. Bates?)
>>
>
> Exactly. Thanks for the summary.
>
>
>>
>> HTH,
>>
>> Simon.
>>
>>
>>
>>> observer and set are factors
>>> y and lbPerAc are numeric
>>>
>>> If you know how to do it or have suggestions for reading I will be
>>> grateful.
>>>
>>>
>>> eal
>>>
>>> ps I have already read Pinheiro & Bates, the jan 05 newsletter, and
>>> several postings.
>>>
>>> ______________________________________________
>>> R-help@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
>>>
>>
>> Simon Blomberg, B.Sc.(Hons.), Ph.D, M.App.Stat.
>> Centre for Resource and Environmental Studies
>> The Australian National University
>> Canberra ACT 0200
>> Australia
>> T: +61 2 6125 7800 email: Simon.Blomberg_at_anu.edu.au
>> F: +61 2 6125 0757
>> CRICOS Provider # 00120C
>>
>> ______________________________________________
>> R-help@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
>>
>



R-help@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 Received on Mon Jul 18 07:38:37 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:33:46 EST