Hi all,
I am getting an error message with profile(). Perhaps someone can
explain why?
I am running R 2.12.1 on Vista Home Basic. Package lme4a has been loaded
from lme4a_0.999375-63.zip.
Regards, Murray Jorgensen
"N" "AR" "EL" "DEc" "NNI" "Island"
"1" 36 0.33 1.26 36 14 "Chiles"
"2" 30 0.5 1.17 234 13 "Las Papas-Coconuco"
"3" 37 2.03 1.06 543 83 "Sumapaz"
"4" 35 0.99 1.9 551 23 "Tolima-Quindio"
"5" 11 0.03 0.46 773 45 "Paramilo"
"6" 21 2.17 2 801 14 "Cocuy"
"7" 11 0.22 0.7 950 14 "Pamploma"
"8" 13 0.14 0.74 958 5 "Cachira"
"9" 17 0.05 0.61 995 29 "Tama"
"10" 13 0.07 0.66 1065 55 "Batallon"
"11" 29 1.8 1.5 1167 35 "Merida"
"12" 4 0.17 0.75 1182 75 "Perija"
"13" 18 0.61 2.28 1238 75 "Santa Marta"
"14" 15 0.07 0.55 1380 35 "Cende"
> library(lme4a)
> obs <- as.factor(1:14)
> b2.mer <- lmer(N ~ AR + DEc + (1|obs), family = poisson, data = birds)
> print(b2.mer, corr = TRUE)
Generalized linear mixed model fit by maximum likelihood ['merMod']
Family: poisson
Formula: N ~ AR + DEc + (1 | obs)
Data: birds
AIC BIC logLik deviance
32.1190 34.6753 -12.0595 24.1190
Random effects:
Groups Name Variance Std.Dev.
obs (Intercept) 0.00478 0.06913
Number of obs: 14, groups: obs, 14
Fixed effects:
Estimate Std. Error z value
(Intercept) 3.3979212 0.1493763 22.747
AR 0.3286582 0.0765854 4.291
DEc -0.0007958 0.0001594 -4.994
Correlation of Fixed Effects:
(Intr) AR
AR -0.468
DEc -0.796 0.035
Error: is(fm at resp, "lmerResp") is not TRUE
Error in object$coefficients : $ operator not defined for this S4 class