multi-start capability?

[bbolker]

Can't remember whether we've discussed this before/whether this is already documented somewhere. It might be useful to provide the capability (or at least to illustrate in a vignette somewhere) to do a series of fits from distributed initial parameter values. (There are of course way too many possibilities for how the range of parameter values could be specified: list, range for uniform variation, Latin hypercube, etc.)