It is a standard approach to consider the maximum likelihood estimation procedure for the estimation of parameters in statistical modelling. The sample likelihood function has a closed form representation only if the two densities in the integrand are conjugate to each other. In case of any non- conjugate pair, no closed form representation exists. In such situations, we need to approximate the integral by making use of some numerical techniques. A first or second order Laplace approximation or the (adaptive) Gauss-Hermite quadrature method can be applied in order to get an ...