As per the code samples in DOI:10.1111/j.1096-0031.2009.00255.x
As per the code samples in DOI:10.1111/j.1096-0031.2009.00255.x