Attach a data-generating function to an existing model_spec. This enables
the simulate() method for parametric bootstrap and power analysis.
Examples
spec <- model_spec_normal()
spec <- add_simulator(spec, function(n, mu, sigma) rnorm(n, mu, sigma))
simulate(spec, n = 10, par = c(mu = 5, sigma = 2))