Scaling Bayesian inference of mixed multinomial logit models to very large datasets