https://github.com/ehm-lab/2011_gasparrini_JSS_Rcode