Uncertainty quantification for hybrid random logistic models with harvesting via density functions

The so-called logistic model with harvesting, p′(t)=rp(t)(1−p(t)K)−c(t)p(t), p(t0)=p0, is a classical ecological model that has been extensively studied and applied in the deterministic setting. It has also been studied, to some extent, in the stochastic framework using the Itô Calculus by formulating a Stochastic Differential Equation whose uncertainty is driven by the Gaussian white noise. In this paper, we present a new approach, based on the so-called theory of Random Differential Equations, that permits treating all model parameters as a random vector with an arbitrary join probability distribution (so, not just Gaussian). We take extensive advantage of the Random Variable Transformation method to probabilistically solve the full randomized version of the above logistic model with harvesting. It is done by exactly computing the first probability density function of the solution assuming that all model parameters are continuous random variables with an arbitrary join probability density function. The probabilistic solution is obtained in three relevant scenarios where the harvesting or influence function is mathematically described by discontinuous parametric stochastic processes having a biological meaning. The probabilistic analysis also includes the computation of the probability density function of the nontrivial equilibrium state, as well as the probability that stability is reached. All these results are new and extend their deterministic counterpart under very general assumptions. The theoretical findings are illustrated via two numerical examples. Finally, we show a detailed example where results are applied to describe the dynamics of stock of fishes over time using real data.