A nonlinear Gauss-Seidel type algorithm is proposed for computing the maximum posterior estimates of the random effects in a generalized linear mixed model. We show that the algorithm converges in virtually all typical situations of generalized linear mixed models. A numerical example shows the...