Displaying 2 results from an estimated 2 matches for "nplants".
Did you mean:
plants
2005 Feb 01
3
polynomials REML and ML in nlme
...copy
of Bates and Pinheiro. It is probably quite obvious but the answer might
be of general interest.
If I fit a fixed effect with an added quadratic term and then do it as
an orthogonal polynomial using maximum likelihood I get the expected
result- they have the same logLik.
mod2a<-lme(wthole~nplants+I(nplants^2),data=d3,random=~1|field/subplot,m
ethod="ML")
mod2b<-lme(wthole~poly(nplants,2),data=d3,random=~1|field/subplot,method
="ML")
> anova(mod2a,mod2b)
Model df AIC BIC logLik
mod2a 1 6 6698.231 6723.869 -3343.116
mod2b 2 6 6698.231 6...
2004 Feb 17
3
parse error in GLMM function
Hi R-Helpers:
I?m trying to use the function GLMM from lme4 package, (R-1.8.1, Windows
98),and I get the following error:
> pd5 = GLMM(nplant~sitio+
+ fert+
+ remo+
+ sitio:fert+
+ remo:sitio+
+ remo:fert+
+ remo:fert:sitio
+ data=datos,
+ family=binomial,
+