Estimating random coefficient logit models with full covariance matrix: comparing performance of mixed logit and laplace approximation methods