[R-sig-ME] taking in account results of a gmml in, despite, of error warning about memory?

glenda mendieta glendamendieta at gmail.com
Thu Dec 1 22:16:37 CET 2011


> Message: 1
> Date: Wed, 30 Nov 2011 11:06:25 +0000
> From: "ONKELINX, Thierry"<Thierry.ONKELINX at inbo.be>
> To:"r-sig-mixed-models at r-project.org"
> 	<r-sig-mixed-models at r-project.org>
> Subject: Re: [R-sig-ME] taking in account results of a gmml in
> 	despite, of error warning about memory?
> Message-ID:
> 	<AA818EAD2576BC488B4F623941DA742756FD665A at inbomail.inbo.be>
> Content-Type: text/plain; charset="us-ascii"
>
> Dear Glenda,
>
> Your (0+spp|tree) random effect requires a 89 x 89 matrix , or 3916 (co)variances to be estimated. You don't have enough data to support that kind of model. Furthermore it will probably require large design matrices and thus lead to the out of memory errors you got. Maybe (1|tree) + (1|spp) is a doable random effects model that is still sensible. This is a cross effect of tree and spp. You even could consider (1|tree) + (1|spp) + (1|tree:spp). (1|tree) = effect of tree regardsless of spp. (1|spp) = effect of spp regardless of tree. (1|tree:spp) = combined effects of tree and spp.

Thank you very much Thierry, I am trying the last option you proposed, 
seemed to me as the most sensitive approach.
Cheers,

Glenda
> Best regards,
>
> Thierry
>
>
> 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
> Belgium
> Thierry.Onkelinx at inbo.be
> www.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




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