capture log close
log using blood, replace

use blood, clear
list, clean

* estimate model
gllamm y A O B, i(id) offset(lnN) fam(poiss) comp(clus ind cA cAB cO cB) init nocons eform

* obtain predicted counts
gllapred mu, mu
list y mu, clean

* Perform goodness-of-fit test by estimating saturated model
 * store current log-likelihood
local ll=e(ll)

 * dummies for the four counts:
gen n=_n if y>0
quietly tab n, gen(dummy)

 * estimate saturated model using poisson
quietly poisson y dummy*, nocons

 * compare log-likelihoods
local deviance = 2*(e(ll)-`ll')
local pval = chiprob(`deviance',1)
disp "deviance = " `deviance' ", p-value = " `pval'

log close
exit