1. Add MNL_InitData.R to generate a common dataset about choice outcomes
2. print acceptance rate of discreteMH function