Update `fit_model()` documentation to recommend using 4000 iters

[stephaniereinders]

I just noticed that the example for fit_model() doesn't say that all of our tests used 4000 iterations and dropped the first 1000 iterations for burn-in. Update the documentation accordingly.