Robust Parameter Estimation for Rational Ordinary Differential Equations

Oren Bassik,Yosef Berman,Soo Go,Hoon Hong,Ilia Ilmer,Alexey Ovchinnikov,Chris Rackauckas,Pedro Soto,Chee Yap
2023-12-17
Abstract:We present a new approach for estimating parameters in rational ODE models from given (measured) time series data. In typical existing approaches, an initial guess for the parameter values is made from a given search interval. Then, in a loop, the corresponding outputs are computed by solving the ODE numerically, followed by computing the error from the given time series data. If the error is small, the loop terminates and the parameter values are returned. Otherwise, heuristics/theories are used to possibly improve the guess and continue the loop. These approaches tend to be non-robust in the sense that their accuracy depend on the search interval and the true parameter values; furthermore, they cannot handle the case where the parameters are locally identifiable. In this paper, we propose a new approach, which does not suffer from the above non-robustness. In particular, it does not require making good initial guesses for the parameter values or specifying search intervals. Instead, it uses differential algebra, interpolation of the data using rational functions, and multivariate polynomial system solving. We also compare the performance of the resulting software with several other estimation software packages.
Mathematical Software,Symbolic Computation,Dynamical Systems,Quantitative Methods
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to perform robust estimation of parameters in the Rational Ordinary Differential Equations (Rational ODEs) model in the given time - series data. Traditional parameter estimation methods usually require an initial guess value and iteratively optimize these parameter values within a certain search interval until the error is minimized. However, the accuracy of this method highly depends on the setting of the search interval and the location of the true parameter values; in addition, they cannot handle the situation of locally identifiable parameters. To overcome the above non - robustness problems, this paper proposes a new method that does not rely on a good initial guess value or a specified search interval. Specifically, the new method utilizes differential algebra, rational function interpolation of data, and the solution techniques of multivariate polynomial systems. Through these means, this method can estimate parameters more accurately and can also provide good estimation results even in the case of locally identifiable parameters. The main steps of the new method include: 1. **Determine the number of parameter values**: Determine the number of parameter values to be estimated according to the model. 2. **Sufficient differentiation of the model**: Differentiate the model a sufficient number of times to generate a system of equations containing sufficient information. 3. **Use data**: Convert the time - series data into rational functions through the interpolation method and substitute these functions into the differentiated system of equations to generate an over - determined polynomial system. 4. **Solve the polynomial system**: Find all possible solutions by squaring the system and filtering the solution set. 5. **Calculate the error**: For each solution, use a numerical method to solve the ODE model and calculate the error between the predicted value and the actual measured value. 6. **Select the optimal solution**: Select the solution with the smallest error as the final parameter estimate. Through this method, the authors of the paper hope to provide a more robust and accurate parameter estimation method, especially suitable for rational ODE models.