simData function

Simulate data from a multivariate joint model