Generation of multiple binary and continuous non-normal variables simultaneously
given the marginal characteristics and association structure based on the methodology
proposed by Demirtas et al. (2012) <doi:10.1002/sim.5362>.