...
- Read the original data, create a subset table to include only the columns with variables that you want to simulate and remove any rows that include missing values
- Calculate the mean value and the variance-covariance matrix (using the cov() function) of the independent (explanatory) variables
- Generate random independent variables that follow the multivariate normal distribution (using the mvrnorm() function), having the same mean and the same pairwise correlations as the mean and the pairwise correlations of the original independent variables
- Apply logistic regression models (using the glm() function with family='binomial') to estimate the coefficients that represent the relationship between each binary dependent variable and the independent variables of the original dataset
- Use the estimated coefficients to calculate the log odds of the predicted variables and use them to generate random binary variables (using the rbinom() function) with probability equal to each log odd
- Combine the simulated dependent and independent variables in one new table
...