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

John Morrongiello jrmorrongiello at gmail.com
Wed Aug 26 15:23:18 CEST 2015


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

c1<-gamm4((longdist)~s(time,k=4),random=~(1|CODE),data=catfish1,family=Gamma(link=log))

confint(c1$mer,method="profile")
##this returns
Error in mkMerMod(rho = environment(devfun), opt = opt, reTrms = b$reTrms,
:
  unused argument (devFunOnly = TRUE)

confint(c1$mer,method="boot")
##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:
c2<-gamm4((longdist)~s(time,k=4),random=~(1|CODE),data=catfish1,family=gaussian)

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

Thanks for your time
John

	[[alternative HTML version deleted]]



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