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

[Ben Bolker]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

  Is there any chance you could post a reproducible example ... ?

  (Out of curiosity, do you have strong reasons to use a Gamma/log-link
rather than a log-Normal model?)

  Ben Bolker

On 15-09-01 06:51 PM, John Morrongiello wrote:
> 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:
> 
> 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
> 
> -- 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]]
> 
> _______________________________________________ 
> R-sig-mixed-models at r-project.org mailing list 
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
> 

-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.11 (GNU/Linux)

iQEcBAEBAgAGBQJV5kZ3AAoJEOCV5YRblxUHSRIH/jSfi8qVSxuFx9a+r6jy6zhQ
gq2pLAdbR3/gImwlXWyuJAzJRpPgd/oEYOP7l88iyLMOZNCo5hT5rUpsxzE0yxCH
3fsW1Pdy5+uL3As+jS3WYaECZp8wQt/Rh1tmV0ON/rVY7f+6Oa7ZGhc4Hjbg8+a8
MOgDMObnhVj+W2fd0WMgCSNczXUSKbFkgSEHCLoc2gt3MkvVDgL+mV5BI8XrYLob
u83k39Lu/Wp+BmkJoaVSSME3aQM1bpt+oL6wAQS/L0PFkMPIcLQDJ3F7drUrFxNo
UdCuHZ0OxeS2UEvpMIB8qSMb8FjJmKT64Zh0V3rrHXq4Hn1lh6CsHvBt9cmhcuY=
=fo9G
-----END PGP SIGNATURE-----