sim_joint_dist function

Simulate category joint distribution