 In addition to numerically approximating this integral(e.g.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.