skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Direct Estimation of Parameters in ODE Models Using WENDy: Weak-Form Estimation of Nonlinear Dynamics
Abstract We introduce the Weak-form Estimation of Nonlinear Dynamics (WENDy) method for estimating model parameters for non-linear systems of ODEs. Without relying on any numerical differential equation solvers, WENDy computes accurate estimates and is robust to large (biologically relevant) levels of measurement noise. For low dimensional systems with modest amounts of data, WENDy is competitive with conventional forward solver-based nonlinear least squares methods in terms of speed and accuracy. For both higher dimensional systems and stiff systems, WENDy is typically both faster (often by orders of magnitude) and more accurate than forward solver-based approaches. The core mathematical idea involves an efficient conversion of the strong form representation of a model to its weak form, and then solving a regression problem to perform parameter inference. The core statistical idea rests on the Errors-In-Variables framework, which necessitates the use of the iteratively reweighted least squares algorithm. Further improvements are obtained by using orthonormal test functions, created from a set of$$C^{\infty }$$ C bump functions of varying support sizes.We demonstrate the high robustness and computational efficiency by applying WENDy to estimate parameters in some common models from population biology, neuroscience, and biochemistry, including logistic growth, Lotka-Volterra, FitzHugh-Nagumo, Hindmarsh-Rose, and a Protein Transduction Benchmark model. Software and code for reproducing the examples is available athttps://github.com/MathBioCU/WENDy.  more » « less
Award ID(s):
2054085
PAR ID:
10523729
Author(s) / Creator(s):
; ;
Publisher / Repository:
springer
Date Published:
Journal Name:
Bulletin of Mathematical Biology
Volume:
85
Issue:
11
ISSN:
0092-8240
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract By means of a unifying measure-theoretic approach, we establish lower bounds on the Hausdorff dimension of the space-time set which can support anomalous dissipation for weak solutions of fluid equations, both in the presence or absence of a physical boundary. Boundary dissipation, which can occur at both the time and the spatial boundary, is analyzed by suitably modifying the Duchon & Robert interior distributional approach. One implication of our results is that any bounded Euler solution (compressible or incompressible) arising as a zero viscosity limit of Navier–Stokes solutions cannot have anomalous dissipation supported on a set of dimension smaller than that of the space. This result is sharp, as demonstrated by entropy-producing shock solutions of compressible Euler (Drivas and Eyink in Commun Math Phys 359(2):733–763, 2018.https://doi.org/10.1007/s00220-017-3078-4; Majda in Am Math Soc 43(281):93, 1983.https://doi.org/10.1090/memo/0281) and by recent constructions of dissipative incompressible Euler solutions (Brue and De Lellis in Commun Math Phys 400(3):1507–1533, 2023.https://doi.org/10.1007/s00220-022-04626-0 624; Brue et al. in Commun Pure App Anal, 2023), as well as passive scalars (Colombo et al. in Ann PDE 9(2):21–48, 2023.https://doi.org/10.1007/s40818-023-00162-9; Drivas et al. in Arch Ration Mech Anal 243(3):1151–1180, 2022.https://doi.org/10.1007/s00205-021-01736-2). For$$L^q_tL^r_x$$ L t q L x r suitable Leray–Hopf solutions of the$$d-$$ d - dimensional Navier–Stokes equation we prove a bound of the dissipation in terms of the Parabolic Hausdorff measure$$\mathcal {P}^{s}$$ P s , which gives$$s=d-2$$ s = d - 2 as soon as the solution lies in the Prodi–Serrin class. In the three-dimensional case, this matches with the Caffarelli–Kohn–Nirenberg partial regularity. 
    more » « less
  2. Abstract Deep Neural Networks (DNNs) are increasingly deployed in critical applications, where ensuring their safety and robustness is paramount. We present$$_\text {CAV25}$$ CAV 25 , a high-performance DNN verification tool that uses the DPLL(T) framework and supports a wide-range of network architectures and activation functions. Since its debut in VNN-COMP’23, in which it achieved the New Participant Award and ranked 4th overall,$$_\text {CAV25}$$ CAV 25 has advanced significantly, achieving second place in VNN-COMP’24. This paper presents and evaluates the latest development of$$_\text {CAV25}$$ CAV 25 , focusing on the versatility, ease of use, and competitive performance of the tool.$$_\text {CAV25}$$ CAV 25 is available at:https://github.com/dynaroars/neuralsat. 
    more » « less
  3. Abstract This article revisits the problem of global well-posedness for the generalized parabolic Anderson model on$$\mathbb {R}^+\times \mathbb {T}^2$$ R + × T 2 within the framework of paracontrolled calculus (Gubinelli et al. in Forum Math, 2015). The model is given by the equation:$$\begin{aligned} (\partial _t-\Delta ) u=F(u)\eta \end{aligned}$$ ( t - Δ ) u = F ( u ) η where$$\eta \in C^{-1-\kappa }$$ η C - 1 - κ with$$1/6>\kappa >0$$ 1 / 6 > κ > 0 , and$$F\in C_b^2(\mathbb {R})$$ F C b 2 ( R ) . Assume that$$\eta \in C^{-1-\kappa }$$ η C - 1 - κ and can be lifted to enhanced noise, we derive new a priori bounds. The key idea follows from the recent work by Chandra et al. (A priori bounds for 2-d generalised Parabolic Anderson Model,,2024), to represent the leading error term as a transport type term, and our techniques encompass the paracontrolled calculus, the maximum principle, and the localization approach (i.e. high-low frequency argument). 
    more » « less
  4. Abstract Let$$p_{1},\ldots ,p_{n}$$ p 1 , , p n be a set of points in the unit square and let$$T_{1},\ldots ,T_{n}$$ T 1 , , T n be a set of$$\delta $$ δ -tubes such that$$T_{j}$$ T j passes through$$p_{j}$$ p j . We prove a lower bound for the number of incidences between the points and tubes under a natural regularity condition (similar to Frostman regularity). As a consequence, we show that in any configuration of points$$p_{1},\ldots , p_{n} \in [0,1]^{2}$$ p 1 , , p n [ 0 , 1 ] 2 along with a line$$\ell _{j}$$ j through each point$$p_{j}$$ p j , there exist$$j\neq k$$ j k for which$$d(p_{j}, \ell _{k}) \lesssim n^{-2/3+o(1)}$$ d ( p j , k ) n 2 / 3 + o ( 1 ) . It follows from the latter result that any set of$$n$$ n points in the unit square contains three points forming a triangle of area at most$$n^{-7/6+o(1)}$$ n 7 / 6 + o ( 1 ) . This new upper bound for Heilbronn’s triangle problem attains the high-low limit established in our previous work arXiv:2305.18253. 
    more » « less
  5. Abstract We study certain one-parameter families of exponential sums of Airy–Laurent type. Their general theory was developed in Katz and Tiep (Airy sheaves of Laurent type: an introduction,https://web.math.princeton.edu/~nmk/kt31_11sept.pdf). In the present paper, we make use of that general theory to compute monodromy groups in some particularly simple families (in the sense of “simple to remember), realizing Weyl groups of type$$E_6$$ E 6 and$$E_8$$ E 8
    more » « less