RMultBinary function

Simulating a multivariate Bernoulli distribution