[R-sig-ME] profile and/or bootstrapped CIs for GAMM random effects

John Morrongiello john.morrongiello at unimelb.edu.au
Wed Sep 2 00:51:54 CEST 2015

Hi all
I thought I'd re-post this message from last week as it may have got lost in the weekend traffic:

I'm having trouble estimating CIs for random effects (i.e. CODE) from a gamm4 model fit with a gamma distribution,
My model is:

##this returns
Error in mkMerMod(rho = environment(devfun), opt = opt, reTrms = b$reTrms,  :
  unused argument (devFunOnly = TRUE)

##this returns
Error in mkNewReTrms(object, newdata, compReForm, na.action = na.action,  :
  random effects specified in re.form that were not present in original model
Are these errors arising because I'm using a gamma distribution? Is there another way to get the random effect CIs?

I can get random effect CIs for the following model:

confint(c2$mer,method="boot")##same error message as above

Thanks for your time

Dr. John R. Morrongiello
School of BioSciences
University of Melbourne
Victoria 3010, Australia
T: +61 3 8344 8929
M: +61 403 338 554
E: john.morrongiello at unimelb.edu.au<mailto:%20jmorrongiell at unimelb.edu.au>
W: morrongiellolab.com<http://morrongiellolab.com/>

	[[alternative HTML version deleted]]

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