Extracting standardized residuals from a gamma(log) glmmTMB

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]]

library(glmmTMB)
dd <- data.frame(g=rgamma(100,shape=4,scale=1))
g1 <- glmmTMB(g~1,family=Gamma(link="log"),data=dd)
resid(g1,"pearson")

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,

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/f9388401f2173b774b16d63d31fff9a1dbc3c48f/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]]

	[[alternative HTML version deleted]]