Simulate correlated binary variables

Hi,

I would like to simulate correlated binary variables.

Theoretically, I should have been able to do the simulation from the following paper

https://journals.sagepub.com/doi/pdf...867X1501500118

and use the command proposed in the aforementioned paper

rbinary

Unfortunately, this command does not work. I get an error message

file "rbinary_simudata.dta" not found
r(601);



Do you have any suggestions? How can I overcome the problem of this command? Or, do you have any recommendation of an alternative way to simulate correlated binary variables?

Many Thanks,
Michalis