glmer and glmmPQL

Dear Raguel,

Make sure that both site and transect are factors.

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

-----Oorspronkelijk bericht-----
Van: r-sig-mixed-models-bounces at r-project.org [mailto:r-sig-mixed-models-bounces at r-project.org] Namens Raquel Benavides
Verzonden: woensdag 23 november 2011 16:59
Aan: r-sig-mixed-models at r-project.org
Onderwerp: [R-sig-ME] glmer and glmmPQL

Hi, 

I am a beginner in mixed models in R. In my design I have 3 nested random factors. Running glmer everything is fine when I just include the first level, but when I try to include the second level I have this error message:

model_seed2<-glmer(seedling~altitude+slope+herbs+shrubs+rocks+fec+AB+light+(
1|site/transect),data=datos,family=poisson)

Error: length(f1) == length(f2) is not TRUE

In addition: Warning messages:

1: In transect:site :

  numerical expression has 1319 elements: only the first used

2: In transect:site :

  numerical expression has 1319 elements: only the first used

 

What does it means? I try this model with glmmPQL and everything was fine, but I have no idea why I have problems with glmer and no with glmPQL.Any clue?

Thanks a lot!

 



_______________________________________________
R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models