Estimate the log-likelihood using Gaussian Quadrature (multidimensional grid)

calc.2LL(fit, nnodes.gq = 8, nsd.gq = 4)

## Arguments

fit saemFit fit number of nodes to use for the Gaussian quadrature when computing the likelihood with this method (defaults to 1, equivalent to the Laplaclian likelihood) span (in SD) over which to integrate when computing the likelihood by Gaussian quadrature. Defaults to 3 (eg 3 times the SD)