[R-sig-ME] Extracting standardized residuals from a gamma(log) glmmTMB
Mollie Brooks
mollieebrooks at gmail.com
Tue Oct 17 17:33:48 CEST 2017
Hi Aoibheann,
That seems to be a bug in glmmTMB since the Gamma family does include a
variance function.
I think we need to change this line of code https://github.com/
glmmTMB/glmmTMB/blob/f9388401f2173b774b16d63d31fff9
a1dbc3c48f/glmmTMB/R/methods.R#L483
Can you add an issue on GitHub (https://github.com/glmmTMB/glmmTMB/issues)?
If not, I could do it.
cheers,
Mollie
On Tue, Oct 17, 2017 at 3:54 PM, Aoibheann Gaughran <gaughra at tcd.ie> wrote:
> Hi all,
>
> I have run a glmmTMB with (family="Gamma",link="log") and now want to
> validate the model, but am unable to extract the standardised residuals
> with sresid <- resid(mod1tmb, "pearson"). It produces the following error
> message:
>
> Error in residuals.glmmTMB(mod1tmb, "pearson") :
> variance function undefined for family ‘Gamma’; cannot compute Pearson
> residuals
>
>
> Would anyone be able to assist me with this?
>
> Many thanks,
>
> --
> Aoibheann Gaughran
>
> Behavioural and Evolutionary Ecology Research Group
> Zoology Building
> School of Natural Sciences
> Trinity College Dublin
> Dublin 2
> Ireland
> Phone: +353 (86) 3812615
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models
mailing list