https://github.com/gasparrini/2010_gasparrini_StatMed_Rcode