Hi

I fitted the linear mixed-effects individual tree model with three levels of
random parameter (Stand, Plot and Growth-Period).  I would like to make
prediction at population averaged and localize by stand level only (at this
point) using n number of randomly selected trees.

fm <-     lme(log(di) ~ log(dbh)+dbh^2+BA+REL,
                    random=~1|Stand/Plot/GP,
                    data=tree,
                    method="REML",
                    weights=varPower())

Here is what I would like to do:
- randomly pick at least n samples from each level from the newdata with
same hierarchical data structure i.e. Stand/Plot/GP
- generate random effects coefficients for each level (at least at the stand
level for now)
- make prediction using newly generated random effects coefficient for the
entire newdata
- also predict population level average

Is there any library I could call.  Any help would be highly appreciated.
Thank you so much

regards,
Prasid

	[[alternative HTML version deleted]]


