1. update code for simulating dataset (now the dataset better mimics the real one)
2. organize code in Step 0---now there are two parallel ones, Step0a for preparing real data and Step0b for generating simulated data
3. Step 1 and 2 are updated accordingly to be seamlessly integrated with output of 0a or 0b
4. Complete Step 1 and 2 with 8 models in both