 via Gauss–Hermite quadrature), methods motivated by laplace approximation have been proposed.For example, the penalized quasi-likelihood method, which essentially involves repeatedly fitting (i.e.doubly iterative) a weighted normal mixed model with a working variate, is implemented by various commercial and open source statistical programs.