[R-sig-ME] Problem with bootMer and models with offset

Martin Hecht martin.hecht at iqb.hu-berlin.de
Mon Oct 13 19:11:23 CEST 2014


Hi,

I am using the latest version of R (3.1.1) and lme4 (1.1-7) .

I'm getting some weird implausible bootstrap results from bootMer for a
model with an offset (models without offset work fine).

I've attached an Rdata file with the returned object from lmer (lmerObj)
and the returned object from bootMer (booted)

Any help would be appreciated very much.


the syntax is:

booted <- bootMer ( x = lmerObj , FUN =
function(fit){return(c(fixef(fit),getME(fit,'theta')))} , nsim = 100 ,
seed = 65835 , verbose = TRUE )

> booted$t0[1]
(Intercept)
    2.606621

> mean(booted$t[,1])
[1] -78.21035

> boot::boot.ci ( booted , conf = 0.95 , type = "basic" , index = 1 )
BOOTSTRAP CONFIDENCE INTERVAL CALCULATIONS
Based on 100 bootstrap replicates

CALL :
boot::boot.ci(boot.out = booted, conf = 0.95, type = "basic",
    index = 1)

Intervals :
Level      Basic        
95%   (69.864, 95.425 ) 
Calculations and Intervals on Original Scale
Some basic intervals may be unstable



More information about the R-sig-mixed-models mailing list