[R-sig-ME] Profile problem in lme4a

[lancelot]

profile() has not been implemented yet for generalized linear 
mixed-effect models.

R

Le 25/04/2011 06:42, Murray Jorgensen a écrit :
> 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
>
>  > write.table(birds)
> "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
>  > pr2 <- profile(b2.mer)
> Error: is(fm at resp, "lmerResp") is not TRUE
>  > confint(b2.mer)
> Error in object$coefficients : $ operator not defined for this S4 class
>

-- 
Renaud Lancelot
EDENext Project, coordinator: http://www.edenext.eu/

CIRAD, UMR15, Campus International de Baillarguet TA A-DIR / B
F34398 Montpellier

Tel.  +33 4 67 59 37 17  -  Fax  +33 4 67 59 37 95
Secr. +33 4 67 59 37 37  - Cell. +33 6 77 52 08 69