Data analysis code for aggregated data as reported in the manuscript