Code to replicate (partially) the NMMAPS dataset and fit the relevant model.