mrgsim function

Simulate from a model object