Dear R, The stage 1 (nAGQ=0) use the initial values for theta by default and initial values for beta is obtained from a generalized linear model fit. The stage 2 initial values are theta estimated from stage 1 and initial values for beta the same ones from stage 1, right? Thank's! Lia [[alternative HTML version deleted]]