rbifactor function

simulate from bi-factor copula model