Bifurcation analysis of models with uncertain function specification: how should we proceed?

When we investigate the bifurcation structure of models of natural phenomena, we usually assume that all model functions are mathematically specified and that the only existing uncertainty is with respect to the parameters of these functions. In this case, we can split the parameter space into domains corresponding to qualitatively similar dynamics, separated by bifurcation hypersurfaces. On the other hand, in the biological sciences, the exact shape of the model functions is often unknown, and only some qualitative properties of the functions can be specified: mathematically, we can consider that the unknown functions belong to a specific class of functions. However, the use of two different functions belonging to the same class can result in qualitatively different dynamical behaviour in the model and different types of bifurcation. In the literature, the conventional way to avoid such ambiguity is to narrow the class of unknown functions, which allows us to keep patterns of dynamical behaviour consistent for varying functions. The main shortcoming of this approach is that the restrictions on the model functions are often given by cumbersome expressions and are strictly model-dependent: biologically, they are meaningless. In this paper, we suggest a new framework (based on the ODE paradigm) which allows us to investigate deterministic biological models in which the mathematical formulation of some functions is unspecified except for some generic qualitative properties. We demonstrate that in such models, the conventional idea of revealing a concrete bifurcation structure becomes irrelevant: we can only describe bifurcations with a certain probability. We then propose a method to define the probability of a bifurcation taking place when there is uncertainty in the parameterisation in our model. As an illustrative example, we consider a generic predator-prey model where the use of different parameterisations of the logistic-type prey growth function can result in different dynamics in terms of the type of the Hopf bifurcation through which the coexistence equilibrium loses stability. Using this system, we demonstrate a framework for evaluating the probability of having a supercritical or subcritical Hopf bifurcation.