Source: http://www2.sas.com/proceedings/sugi25/25/aa/25p003.pdf