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
profile and/or bootstrapped CIs for GAMM random effects
1 message · John Morrongiello