This seems much easier to do by using blackbox optimization (via Optim) with the actual functional form you want that provides a closure over your data. It’s true that you can make this into a LsqFit problem, but the more general approach of blackbox optimization is easier once you’re used to it.
1 Like