simulate an ordinal dependent variable that is correlated with multiple existing variables

Question

I\'d like to simulate an ordinal dependent variable that is correlated with multiple existing variables:

set.seed(12345)

pop <- data.fram