Thierry,
The first suggestion worked. Thank you very much.
*Ben Caldwell*

University of California, Berkeley
137 Mulford Hall #3114
Berkeley, CA 94720
Office 223 Mulford Hall
(510)859-3358



On Fri, Apr 29, 2011 at 1:52 AM, ONKELINX, Thierry <Thierry.ONKELINX@inbo.be
> wrote:

> Dear Ben,
>
> Are site, transect and plot factors? And do they have unique id's?
>
> You could try this
>
> rws30.UL$site <- factor(rws30.UL$site)
> rws30.UL$transect <- interaction(rws30.UL$site, rws30.UL$transect, drop =
> TRUE)
> rws30.UL$plot <- interaction(rws30.UL$site, rws30.UL$transect,
> rws30.UL$plot, drop = TRUE)
> modelincrBS<-glmer(l.ru.ba.incr~shigo.av+pre.f.crwn.length+bak.thick.bh+Date+slope.pos.num+dens.T+dbh+leaf.area+can.pos.num
> +(1|site/transect/plot),
>  data=rws30.UL, family=gaussian, na.action=na.omit)
>
> Or
>
> modelincrBS<-glmer(l.ru.ba.incr~shigo.av+pre.f.crwn.length+bak.thick.bh+Date+slope.pos.num+dens.T+dbh+leaf.area+can.pos.num
> +(1|site) + (1|transect) + (1|plot),
>  data=rws30.UL, family=gaussian, na.action=na.omit)
>
> Best regards,
>
> Thierry
>
>
> ----------------------------------------------------------------------------
> ir. Thierry Onkelinx
> Instituut voor natuur- en bosonderzoek
> team Biometrie & Kwaliteitszorg
> Gaverstraat 4
> 9500 Geraardsbergen
> Belgium
>
> Research Institute for Nature and Forest
> team Biometrics & Quality Assurance
> Gaverstraat 4
> 9500 Geraardsbergen
> Belgium
>
> tel. + 32 54/436 185
> Thierry.Onkelinx@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-help-bounces@r-project.org
> > [mailto:r-help-bounces@r-project.org] Namens Benjamin Caldwell
> > Verzonden: vrijdag 29 april 2011 0:37
> > Aan: r-help
> > Onderwerp: [R] using lme4 with three nested random effects
> >
> > Hi all,
> > I'm trying to fit models for data with three levels of nested random
> > effects: site/transect/plot. For example,
> >
> > modelincrBS<-glmer(l.ru.ba.incr~shigo.av+pre.f.crwn.length+bar
> > k.thick.bh+Date+slope.pos.num+dens.T+dbh+leaf.area+can.pos.num
> > +(1|site/transect/plot),
> > data=rws30.UL, family=gaussian, na.action=na.omit)
> >
> > but I get the following error:
> >
> > Error: length(f1) == length(f2) is not TRUE In addition:
> > Warning messages:
> > 1: In plot:(transect:site) :
> >   numerical expression has 92 elements: only the first used
> > 2: In plot:(transect:site) :
> >   numerical expression has 92 elements: only the first used
> >
> > The formulation works for two nested effects (e.g. 1|site/transect)
> >
> > I can get it to run in lme
> > modelincrBS<-lme(l.ru.ba.incr~shigo.av+pre.f.crwn.length+bark.
> > thick.bh+Date+slope.pos.num+dens.T+dbh+leaf.area+can.pos.num,
> > data=rws30.UL, random=(~1| site/transect/plot),na.action=na.omit)
> >
> > but I can't specify a distribution family in that package.
> >
> > Any help much appreciated.
> >
> > Ben Caldwell
> >
> > *
> > *
> >
> >       [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help@r-project.org 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.
> >
>

	[[alternative HTML version deleted]]

