The code for simulation results in Supplementary Fig. 1. 