Update MNL Scenario 0, 1, 3M, 3A, 3Reg:
1. new discrete MH algorithm (using discretized Normal proposal)
2. use the right rgamma (instead of rpois) for simulating lambda in the Gibbs sampling iterations