I am trying to get bootstrapped predictions from a logistic glmer (lme4) object and compare them to the “observed” independent variable from each bootstrap iteration. I can use bootMer to get the bootstrapped predictions:

bb <- bootMer(glmm.M35, FUN=function(x)predict(x,re.form=NA,newdata=df.fit,type="response"), nsim=10) # test with just 10 draws

but it doesn’t appear that bootMer saves the replicate datasets. Can anyone help me modify the bootMer function to output the bootstrap data frames or at least the independent variable (y) from each? Also any advice on handling the occasions when an iteration fails to converge in bootMer would be appreciated.

Thanks,
Dan
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Daniel Hocking
Department of Environmental Conservation
Northeast Climate Science Center
University of Massachusetts

http://www.danieljhocking.wordpress.com
dhocking@umass.edu<mailto:dhocking@umass.edu>
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

	[[alternative HTML version deleted]]

