https://github.com/ehm-lab/2010_gasparrini_StatMed_Rcode