[R-sig-ME] Profile problem in lme4a

[Murray Jorgensen]

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

-- 
Dr Murray Jorgensen      http://www.stats.waikato.ac.nz/Staff/maj.html
Department of Statistics, University of Waikato, Hamilton, New Zealand
Email: maj at waikato.ac.nz    majorgensen at ihug.co.nz      Fax 7 838 4155
Phone  +64 7 838 4773 wk    Home +64 7 825 0441   Mobile 021 0200 8350