scripts for a simulation study in (bielejec et al 2014)