[R-sig-ME] Bug in afex, pbkrtest or lme4 in parametric bootstrapping?
Tom Wenseleers
Tom.Wenseleers at bio.kuleuven.be
Wed Aug 21 19:47:19 CEST 2013
Dear all,
I was just trying to use parametric bootstrapping to calculate significance levels in a Poisson mixed model using afex/pbkrtest. Using the data in attachment and the code
data_alates=read.csv("alates_raw_FINAL3.csv", header=TRUE, fill=TRUE)
mixed(ALATES~ANT_TENDED*MELEZITOSE+NR_START+(1|CLONES)+(1|PLANT)+(1|REPLICATE),type=3,family=poisson,method="PB",args.test = list(nsim = 10),data=data_alates)
I get the following error message though:
Fitting 6 lmer() models:
[......]
Obtaining 5 p-values:
[
Error in (function (classes, fdef, mtable) :
unable to find an inherited method for function 'refit' for signature '"glmerMod"'
In addition: Warning message:
In mixed(ALATES ~ ANT_TENDED * MELEZITOSE + NR_START + (1 | CLONES) + :
Numerical variables NOT centered on 0 (i.e., likely bogus results if in interactions): NR_START
Any thoughts perhaps what might be wrong, and what workaround I could perhaps use to make this work?
Cheers,
Tom
PS I was using the following package versions:
sessionInfo()
R version 3.0.0 (2013-04-03)
Platform: x86_64-w64-mingw32/x64 (64-bit)
locale:
[1] LC_COLLATE=English_United Kingdom.1252 LC_CTYPE=English_United Kingdom.1252 LC_MONETARY=English_United Kingdom.1252
[4] LC_NUMERIC=C LC_TIME=English_United Kingdom.1252
attached base packages:
[1] grid stats4 splines parallel stats graphics grDevices utils datasets methods base
other attached packages:
[1] lsmeans_1.10-00 plyr_1.8 effects_2.2-4 colorspace_1.2-2 nlme_3.1-109 afex_0.6-77 stringr_0.6.2 reshape2_1.2.2
[9] car_2.0-16 nnet_7.3-6 coin_1.0-22 modeltools_0.2-19 mvtnorm_0.9-9994 survival_2.37-4 pbkrtest_0.3-5.1 lme4_1.1-0
[17] Matrix_1.0-12 lattice_0.20-15 glmmADMB_0.7.7 R2admb_0.7.5.3 MASS_7.3-26 devtools_1.2
loaded via a namespace (and not attached):
[1] digest_0.6.3 evaluate_0.4.3 httr_0.2 memoise_0.1 minqa_1.2.1 multcomp_1.2-18 RCurl_1.95-4.1 tools_3.0.0 whisker_0.3-2
More information about the R-sig-mixed-models
mailing list