# [R-sig-ME] glmmTMB pearson residuals

Stephanie Rivest @rive046 @ending from uott@w@@c@
Tue Sep 18 20:20:59 CEST 2018

```Hello,

I am working with glmmTMB and have a question about calculating the pearson
residuals. Following the math laid out in Zuur et al. (2012) I can
understand how the following code works for models fit with the family
nbinom2:

v <- family(model)\$variance
p <- predict(model,zitype="zprob")  ## z-i probability
mu <- predict(model,zitype="conditional")  ## mean of conditional
distribution
pred <- predict(model,zitype="response") ## (1-p)*mu
k <- sigma(m5) ## dispersion parameter
pvar <- (1-p)*v(mu,k)+mu^2*(p^2+p)
pearson_resid <- (data.frame\$Response-pred)/sqrt(pvar)

What I am unsure about is how this changes if the family being used is
nbinom1?

For the poisson family, does the following code make sense for calculating
pearson resids?

p <- predict(model,zitype="zprob")  ## z-i probability
mu <- predict(model,zitype="conditional")  ## mean of conditional
distribution
pred <- predict(model,zitype="response") ## (1-p)*mu
pvar <- (1-p)*(mu + p*mu^2)
pearson_resid <- (data.frame\$Response - pred) / sqrt(pvar)

Any help is greatly appreciated!

Stephanie Rivest
Ph.D. Candidate | Candidate au Doctorat
Dept. of Biology | Dép. de Biologie
University of Ottawa | Université d'Ottawa

[[alternative HTML version deleted]]

```