Universal Differential Equations while also estimating non NN parameters

I have been using the the tutorial on Automatically Discover Missing Physics by Embedding Machine Learning into Differential Equations as a textbook to build my ODEs with embedded NN’s (It’s awesome). I also want to optimize the parameters for the rest of my ODE (the ones not in the NN). I am thinking I will iterate between updating the NN parameters while keeping the remaining ODE parameters fixed and then the ODE parameters while keeping the NN parameters fixed.

Is there any precedent for this? Perhaps some good examples or tips that one might offer?

Thanks,

DS

The bottom of the tutorial does this. After fitting the UDE and doing the symbolic regression, it then does a refit on the full set of parameters.

The thing you need to be careful about is that the NN can capture “anything”, so your parameters will be biased downwards as a compensation of that. There isn’t good literature on this, but I’ve been looking to do some things with a student to showcase different regularizations as a fix for this behavior.