Matthias Schmidt
2007-Aug-15 13:21 UTC
[R] problems with specifying an offset term in lmer and lmer2
Dear all,
I just want to integrate an offset in a linear mixed model applying
functions lmer and lmer2 from library(lme4). It seems that lmer just
ignores the offset term and lmer2 does the fit but replies an error if
calling the fitted object:
#lmer with offset
lmeTNr_offset<-lmer(res_id~1+offset(fit.gam)+(1|TNr),Bu_WP)
lmeTNr_offset
Linear mixed-effects model fit by REML
Formula: res_id ~ 1 + offset(fit.gam) + (1 | TNr)
Data: Bu_WP
AIC BIC logLik MLdeviance REMLdeviance
-8869 -8852 4436 -8883 -8873
Random effects:
Groups Name Variance Std.Dev.
TNr 0.016274 0.12757
Residual 0.033355 0.18263
number of obs: 25208, groups: TNr, 4812
Fixed effects:
Estimate Std. Error t value
(Intercept) -0.008815 0.002342 -3.764
#results in the same output as the specification without an offset:
#lmer without offset
lmeTNr<-lmer(res_id~1+(1|TNr),Bu_WP)
lmeTNr
Linear mixed-effects model fit by REML
Formula: res_id ~ 1 + (1 | TNr)
Data: Bu_WP
AIC BIC logLik MLdeviance REMLdeviance
-8869 -8852 4436 -8883 -8873
Random effects:
Groups Name Variance Std.Dev.
TNr 0.016274 0.12757
Residual 0.033355 0.18263
number of obs: 25208, groups: TNr, 4812
Fixed effects:
Estimate Std. Error t value
(Intercept) -0.008815 0.002342 -3.764
#lmer2 with offset
lme2TNr_offset<-offset<-lmer2(res_id~1+offset(fit.gam)+(1|TNr),Bu_WP)
#calling lme2TNr_offset results in an error:
error in checkSlotAssignment(object, name, value) :
assignment of an object of class "array" is not valid for
slot
"offset" in an object of class "summary.lmer2"; is(value,
"numeric") is
not TRUE
thanks for any help
best regards
Matthias
-------------------------------------------
Dr. Matthias Schmidt
Nordwestdeutsche Forstliche Versuchsanstalt
Abteilung Waldwachstum
Gr?tzelstr. 2
37079 G?ttingen
TEL.: 0551/69401-110
Fax: 0551/69401-160
The Northwest German Forest Research Institute
Department of Forest Growth
Gr?tzelstr. 2
37079 G?ttingen
TEL.: 0551/69401-110
Fax: 0551/69401-160