A quick scan of the source tells me this should be straightforward, it will require one extra line to dump any Parameters from initial into the list of parameters of the model, and an extra line or two in __call__ to correctly feed the provided guesses to the ode-solver.
In relation to this SO question, it might be worth re-investigating how to allow Parameters to be allowed as initial values, e.g.
A quick scan of the source tells me this should be straightforward, it will require one extra line to dump any Parameters from initial into the list of parameters of the model, and an extra line or two in
__call__
to correctly feed the provided guesses to the ode-solver.