A comprehensive probabilistic solution of random SIS-type epidemiological models using the Random Variable Transformation technique

[EN] This paper provides a complete probabilistic description of SIS-type epidemiological models where all the input parameters (contagion rate, recovery rate and initial conditions) are assumed to be random variables. By applying the Random Variable Transformation technique, the first probability density function, the mean and the variance functions as well as confidence intervals associated to the solution of SIS-type epidemiological models are determined under the general hypothesis that the random inputs have any joint probability density function. The distributions to describe the time until a given proportion of the population remains susceptible and infected are also determined. Lastly, a probabilistic description of the so-called basic reproductive number is included. The theoretical results are applied to an illustrative example showing good fitting.This work has been partially supported by the Ministerio de Economia y Competitividad grants MTM2013-41765-P and TRA2012-36932. Ana Navarro Quiles acknowledges the doctorate scholarship granted by Programa de Ayudas de Investigacion y Desarrollo (PAID), Universitat Politecnica de Valencia.Casabán, M.; Cortés, J.; Navarro-Quiles, A.; Romero, J.; Roselló, M.; Villanueva Micó, RJ. (2016). A comprehensive probabilistic solution of random SIS-type epidemiological models using the Random Variable Transformation technique. Communications in Nonlinear Science and Numerical Simulation. 32:199-210. https://doi.org/10.1016/j.cnsns.2015.08.009S1992103

A comprehensive probabilistic analysis of SIR-type epidemiological models based on full randomized Discrete-Time Markov Chain formulation with applications

[EN] This paper provides a comprehensive probabilistic analysis of a full randomization of approximate SIR-type epidemiological models based on discrete-time Markov chain formulation. The randomization is performed by assuming that all input data (initial conditions, the contagion, and recovering rates involved in the transition matrix) are random variables instead of deterministic constants. In the first part of the paper, we determine explicit expressions for the so called first probability density function of each subpopulation identified as the corresponding states of the Markov chain (susceptible, infected, and recovered) in terms of the probability density function of each input random variable. Afterwards, we obtain the probability density functions of the times until a given proportion of the population remains susceptible, infected, and recovered, respectively. The theoretical analysis is completed by computing explicit expressions of important randomized epidemiological quantities, namely, the basic reproduction number, the effective reproduction number, and the herd immunity threshold. The study is conducted under very general assumptions and taking extensive advantage of the random variable transformation technique. The second part of the paper is devoted to apply our theoretical findings to describe the dynamics of the pandemic influenza in Egypt using simulated data excerpted from the literature. The simulations are complemented with valuable information, which is seldom displayed in epidemiological models. In spite of the nonlinear mathematical nature of SIR epidemiological model, our results show a strong agreement with the approximation via an appropriate randomized Markov chain. A justification in this regard is discussed.Spanish Ministerio de Economia y Competitividad, Grant/Award Number: MTM2017-89664-P; Generalitat Valenciana, Grant/Award Number: APOSTD/2019/128; Ministerio de Economia y Competitividad, Grant/Award Number: MTM2017-89664-PCortés, J.; El-Labany, S.; Navarro-Quiles, A.; Selim, MM.; Slama, H. (2020). A comprehensive probabilistic analysis of SIR-type epidemiological models based on full randomized Discrete-Time Markov Chain formulation with applications. Mathematical Methods in the Applied Sciences. 43(14):8204-8222. https://doi.org/10.1002/mma.6482S820482224314Hamra, G., MacLehose, R., & Richardson, D. (2013). Markov Chain Monte Carlo: an introduction for epidemiologists. International Journal of Epidemiology, 42(2), 627-634. doi:10.1093/ije/dyt043Becker, N. (1981). A General Chain Binomial Model for Infectious Diseases. Biometrics, 37(2), 251. doi:10.2307/2530415Allen, L. J. S. (2010). An Introduction to Stochastic Processes with Applications to Biology. doi:10.1201/b12537Hethcote, H. W. (2000). The Mathematics of Infectious Diseases. SIAM Review, 42(4), 599-653. doi:10.1137/s0036144500371907Brauer, F., & Castillo-Chávez, C. (2001). Mathematical Models in Population Biology and Epidemiology. Texts in Applied Mathematics. doi:10.1007/978-1-4757-3516-1Cortés, J.-C., Navarro-Quiles, A., Romero, J.-V., & Roselló, M.-D. (2018). Some results about randomized binary Markov chains: theory, computing and applications. International Journal of Computer Mathematics, 97(1-2), 141-156. doi:10.1080/00207160.2018.1440290Cortés, J.-C., Navarro-Quiles, A., Romero, J.-V., & Roselló, M.-D. (2017). Randomizing the parameters of a Markov chain to model the stroke disease: A technical generalization of established computational methodologies towards improving real applications. Journal of Computational and Applied Mathematics, 324, 225-240. doi:10.1016/j.cam.2017.04.040Casabán, M.-C., Cortés, J.-C., Romero, J.-V., & Roselló, M.-D. (2015). Probabilistic solution of random SI-type epidemiological models using the Random Variable Transformation technique. Communications in Nonlinear Science and Numerical Simulation, 24(1-3), 86-97. doi:10.1016/j.cnsns.2014.12.016Casabán, M.-C., Cortés, J.-C., Navarro-Quiles, A., Romero, J.-V., Roselló, M.-D., & Villanueva, R.-J. (2016). A comprehensive probabilistic solution of random SIS-type epidemiological models using the random variable transformation technique. Communications in Nonlinear Science and Numerical Simulation, 32, 199-210. doi:10.1016/j.cnsns.2015.08.009Slama, H., Hussein, A., El-Bedwhey, N. A., & Selim, M. M. (2019). An approximate probabilistic solution of a random SIR-type epidemiological model using RVT technique. Applied Mathematics and Computation, 361, 144-156. doi:10.1016/j.amc.2019.05.019Slama, H., El-Bedwhey, N. A., El-Depsy, A., & Selim, M. M. (2017). Solution of the finite Milne problem in stochastic media with RVT Technique. The European Physical Journal Plus, 132(12). doi:10.1140/epjp/i2017-11763-6Kegan, B., & West, R. W. (2005). Modeling the simple epidemic with deterministic differential equations and random initial conditions. Mathematical Biosciences, 195(2), 179-193. doi:10.1016/j.mbs.2005.02.004Dorini, F. A., Cecconello, M. S., & Dorini, L. B. (2016). On the logistic equation subject to uncertainties in the environmental carrying capacity and initial population density. Communications in Nonlinear Science and Numerical Simulation, 33, 160-173. doi:10.1016/j.cnsns.2015.09.009Van den Driessche, P. (2017). Reproduction numbers of infectious disease models. Infectious Disease Modelling, 2(3), 288-303. doi:10.1016/j.idm.2017.06.002Heffernan, J. ., Smith, R. ., & Wahl, L. . (2005). Perspectives on the basic reproductive ratio. Journal of The Royal Society Interface, 2(4), 281-293. doi:10.1098/rsif.2005.0042Khalil, K. M., Abdel-Aziz, M., Nazmy, T. T., & Salem, A.-B. M. (2012). An Agent-Based Modeling for Pandemic Influenza in Egypt. Intelligent Systems Reference Library, 205-218. doi:10.1007/978-3-642-25755-1_1

Combining Polynomial Chaos Expansions and the Random Variable Transformation Technique to Approximate the Density Function of Stochastic Problems, Including Some Epidemiological Models

[EN] In this paper, we deal with computational uncertainty quantification for stochastic models with one random input parameter. The goal of the paper is twofold: First, to approximate the set of probability density functions of the solution stochastic process, and second, to show the capability of our theoretical findings to deal with some important epidemiological models. The approximations are constructed in terms of a polynomial evaluated at the random input parameter, by means of generalized polynomial chaos expansions and the stochastic Galerkin projection technique. The probability density function of the aforementioned univariate polynomial is computed via the random variable transformation method, by taking into account the domains where the polynomial is strictly monotone. The algebraic/exponential convergence of the Galerkin projections gives rapid convergence of these density functions. The examples are based on fundamental epidemiological models formulated via linear and nonlinear differential and difference equations, where one of the input parameters is assumed to be a random variable.This work has been supported by the Spanish Ministerio de Economia y Competitividad grant MTM2017-89664-P. The author Marc Jornet acknowledges the doctorate scholarship granted by Programa de Ayudas de Investigacion y Desarrollo (PAID), Universitat Politecnica de Valencia.Calatayud-Gregori, J.; Chen-Charpentier, BM.; Cortés, J.; Jornet-Sanz, M. (2019). Combining Polynomial Chaos Expansions and the Random Variable Transformation Technique to Approximate the Density Function of Stochastic Problems, Including Some Epidemiological Models. Symmetry (Basel). 11(1):1-28. https://doi.org/10.3390/sym11010043S128111Strand, J. . (1970). Random ordinary differential equations. Journal of Differential Equations, 7(3), 538-553. doi:10.1016/0022-0396(70)90100-2Bharucha-Reid, A. T. (1964). On the theory of random equations. Proceedings of Symposia in Applied Mathematics, 40-69. doi:10.1090/psapm/016/0189071Xiu, D., & Karniadakis, G. E. (2002). The Wiener--Askey Polynomial Chaos for Stochastic Differential Equations. SIAM Journal on Scientific Computing, 24(2), 619-644. doi:10.1137/s1064827501387826Chen-Charpentier, B.-M., Cortés, J.-C., Licea, J.-A., Romero, J.-V., Roselló, M.-D., Santonja, F.-J., & Villanueva, R.-J. (2015). Constructing adaptive generalized polynomial chaos method to measure the uncertainty in continuous models: A computational approach. Mathematics and Computers in Simulation, 109, 113-129. doi:10.1016/j.matcom.2014.09.002Cortés, J.-C., Romero, J.-V., Roselló, M.-D., & Villanueva, R.-J. (2017). Improving adaptive generalized polynomial chaos method to solve nonlinear random differential equations by the random variable transformation technique. Communications in Nonlinear Science and Numerical Simulation, 50, 1-15. doi:10.1016/j.cnsns.2017.02.011Chen-Charpentier, B. M., & Stanescu, D. (2010). Epidemic models with random coefficients. Mathematical and Computer Modelling, 52(7-8), 1004-1010. doi:10.1016/j.mcm.2010.01.014Lucor, D., Su, C.-H., & Karniadakis, G. E. (2004). Generalized polynomial chaos and random oscillators. International Journal for Numerical Methods in Engineering, 60(3), 571-596. doi:10.1002/nme.976Santonja, F., & Chen-Charpentier, B. (2012). Uncertainty Quantification in Simulations of Epidemics Using Polynomial Chaos. Computational and Mathematical Methods in Medicine, 2012, 1-8. doi:10.1155/2012/742086Stanescu, D., & Chen-Charpentier, B. M. (2009). Random coefficient differential equation models for bacterial growth. Mathematical and Computer Modelling, 50(5-6), 885-895. doi:10.1016/j.mcm.2009.05.017Calatayud, J., Cortés, J. C., Jornet, M., & Villanueva, R. J. (2018). Computational uncertainty quantification for random time-discrete epidemiological models using adaptive gPC. Mathematical Methods in the Applied Sciences, 41(18), 9618-9627. doi:10.1002/mma.5315Villegas, M., Augustin, F., Gilg, A., Hmaidi, A., & Wever, U. (2012). Application of the Polynomial Chaos Expansion to the simulation of chemical reactors with uncertainties. Mathematics and Computers in Simulation, 82(5), 805-817. doi:10.1016/j.matcom.2011.12.001Xiu, D., & Em Karniadakis, G. (2002). Modeling uncertainty in steady state diffusion problems via generalized polynomial chaos. Computer Methods in Applied Mechanics and Engineering, 191(43), 4927-4948. doi:10.1016/s0045-7825(02)00421-8Shi, W., & Zhang, C. (2012). Error analysis of generalized polynomial chaos for nonlinear random ordinary differential equations. Applied Numerical Mathematics, 62(12), 1954-1964. doi:10.1016/j.apnum.2012.08.007Calatayud, J., Cortés, J.-C., & Jornet, M. (2018). On the convergence of adaptive gPC for non-linear random difference equations: Theoretical analysis and some practical recommendations. Journal of Nonlinear Sciences and Applications, 11(09), 1077-1084. doi:10.22436/jnsa.011.09.06Casabán, M.-C., Cortés, J.-C., Romero, J.-V., & Roselló, M.-D. (2015). Probabilistic solution of random SI-type epidemiological models using the Random Variable Transformation technique. Communications in Nonlinear Science and Numerical Simulation, 24(1-3), 86-97. doi:10.1016/j.cnsns.2014.12.016Dorini, F. A., Cecconello, M. S., & Dorini, L. B. (2016). On the logistic equation subject to uncertainties in the environmental carrying capacity and initial population density. Communications in Nonlinear Science and Numerical Simulation, 33, 160-173. doi:10.1016/j.cnsns.2015.09.009Dorini, F. A., & Cunha, M. C. C. (2008). Statistical moments of the random linear transport equation. Journal of Computational Physics, 227(19), 8541-8550. doi:10.1016/j.jcp.2008.06.002Hussein, A., & Selim, M. M. (2012). Solution of the stochastic radiative transfer equation with Rayleigh scattering using RVT technique. Applied Mathematics and Computation, 218(13), 7193-7203. doi:10.1016/j.amc.2011.12.088Hussein, A., & Selim, M. M. (2015). Solution of the stochastic generalized shallow-water wave equation using RVT technique. The European Physical Journal Plus, 130(12). doi:10.1140/epjp/i2015-15249-3Hussein, A., & Selim, M. M. (2013). A general analytical solution for the stochastic Milne problem using Karhunen–Loeve (K–L) expansion. Journal of Quantitative Spectroscopy and Radiative Transfer, 125, 84-92. doi:10.1016/j.jqsrt.2013.03.018Xu, Z., Tipireddy, R., & Lin, G. (2016). Analytical approximation and numerical studies of one-dimensional elliptic equation with random coefficients. Applied Mathematical Modelling, 40(9-10), 5542-5559. doi:10.1016/j.apm.2015.12.041Cortés, J.-C., Navarro-Quiles, A., Romero, J.-V., & Roselló, M.-D. (2017). Full solution of random autonomous first-order linear systems of difference equations. Application to construct random phase portrait for planar systems. Applied Mathematics Letters, 68, 150-156. doi:10.1016/j.aml.2016.12.015El-Tawil, M. A. (2005). The approximate solutions of some stochastic differential equations using transformations. Applied Mathematics and Computation, 164(1), 167-178. doi:10.1016/j.amc.2004.04.062Calatayud, J., Cortés, J.-C., & Jornet, M. (2018). The damped pendulum random differential equation: A comprehensive stochastic analysis via the computation of the probability density function. Physica A: Statistical Mechanics and its Applications, 512, 261-279. doi:10.1016/j.physa.2018.08.024Calatayud, J., Cortés, J. C., & Jornet, M. (2018). Uncertainty quantification for random parabolic equations with nonhomogeneous boundary conditions on a bounded domain via the approximation of the probability density function. Mathematical Methods in the Applied Sciences, 42(17), 5649-5667. doi:10.1002/mma.5333Cortés, J.-C., Navarro-Quiles, A., Romero, J.-V., & Roselló, M.-D. (2018). Solving second-order linear differential equations with random analytic coefficients about ordinary points: A full probabilistic solution by the first probability density function. Applied Mathematics and Computation, 331, 33-45. doi:10.1016/j.amc.2018.02.051Casabán, M.-C., Cortés, J.-C., Navarro-Quiles, A., Romero, J.-V., Roselló, M.-D., & Villanueva, R.-J. (2016). A comprehensive probabilistic solution of random SIS-type epidemiological models using the random variable transformation technique. Communications in Nonlinear Science and Numerical Simulation, 32, 199-210. doi:10.1016/j.cnsns.2015.08.009Kegan, B., & West, R. W. (2005). Modeling the simple epidemic with deterministic differential equations and random initial conditions. Mathematical Biosciences, 194(2), 217-231. doi:10.1016/j.mbs.2005.02.002Crestaux, T., Le Maıˆtre, O., & Martinez, J.-M. (2009). Polynomial chaos expansion for sensitivity analysis. Reliability Engineering & System Safety, 94(7), 1161-1172. doi:10.1016/j.ress.2008.10.008Sudret, B. (2008). Global sensitivity analysis using polynomial chaos expansions. Reliability Engineering & System Safety, 93(7), 964-979. doi:10.1016/j.ress.2007.04.002Chen-Charpentier, B. M., Cortés, J.-C., Romero, J.-V., & Roselló, M.-D. (2013). Some recommendations for applying gPC (generalized polynomial chaos) to modeling: An analysis through the Airy random differential equation. Applied Mathematics and Computation, 219(9), 4208-4218. doi:10.1016/j.amc.2012.11.007Ernst, O. G., Mugler, A., Starkloff, H.-J., & Ullmann, E. (2011). On the convergence of generalized polynomial chaos expansions. ESAIM: Mathematical Modelling and Numerical Analysis, 46(2), 317-339. doi:10.1051/m2an/2011045Giraud, L., Langou, J., & Rozloznik, M. (2005). The loss of orthogonality in the Gram-Schmidt orthogonalization process. Computers & Mathematics with Applications, 50(7), 1069-1075. doi:10.1016/j.camwa.2005.08.009Marzouk, Y. M., Najm, H. N., & Rahn, L. A. (2007). Stochastic spectral methods for efficient Bayesian solution of inverse problems. Journal of Computational Physics, 224(2), 560-586. doi:10.1016/j.jcp.2006.10.010Marzouk, Y., & Xiu, D. (2009). A Stochastic Collocation Approach to Bayesian Inference in Inverse Problems. Communications in Computational Physics, 6(4), 826-847. doi:10.4208/cicp.2009.v6.p826SCOTT, D. W. (1979). On optimal and data-based histograms. Biometrika, 66(3), 605-610. doi:10.1093/biomet/66.3.605National Spanish Health Survey (Encuesta Nacional de Salud de España, ENSE)http://pestadistico.inteligenciadegestion.msssi.es/publicoSNS/comun/ArbolNodos.asp

Probabilistic solution of random SI-type epidemiological models using the Random Variable Transformation technique

[EN] This paper presents a full probabilistic description of the solution of random SI-type epidemiological models which are based on nonlinear differential equations. This description consists of determining: the first probability density function of the solution in terms of the density functions of the diffusion coefficient and the initial condition, which are assumed to be independent random variables; the expectation and variance functions of the solution as well as confidence intervals and, finally, the distribution of time until a given proportion of susceptibles remains in the population. The obtained formulas are general since they are valid regardless the probability distributions assigned to the random inputs. We also present a pair of illustrative examples including in one of them the application of the theoretical results to model the diffusion of a technology using real data.This work has been partially supported by the Ministerio de Economia y Competitividad Grants MTM2013-41765-P and TRA2012-36932.Casabán Bartual, MC.; Cortés López, JC.; Romero Bauset, JV.; Roselló Ferragud, MD. (2015). Probabilistic solution of random SI-type epidemiological models using the Random Variable Transformation technique. Communications in Nonlinear Science and Numerical Simulation. 24(1):86-97. https://doi.org/10.1016/j.cnsns.2014.12.016S869724

Computational methods for random differential equations: probability density function and estimation of the parameters

[EN] Mathematical models based on deterministic differential equations do not take into account the inherent uncertainty of the physical phenomenon (in a wide sense) under study. In addition, inaccuracies in the collected data often arise due to errors in the measurements. It thus becomes necessary to treat the input parameters of the model as random quantities, in the form of random variables or stochastic processes. This gives rise to the study of random ordinary and partial differential equations. The computation of the probability density function of the stochastic solution is important for uncertainty quantification of the model output. Although such computation is a difficult objective in general, certain stochastic expansions for the model coefficients allow faithful representations for the stochastic solution, which permits approximating its density function. In this regard, Karhunen-Loève and generalized polynomial chaos expansions become powerful tools for the density approximation. Also, methods based on discretizations from finite difference numerical schemes permit approximating the stochastic solution, therefore its probability density function. The main part of this dissertation aims at approximating the probability density function of important mathematical models with uncertainties in their formulation. Specifically, in this thesis we study, in the stochastic sense, the following models that arise in different scientific areas: in Physics, the model for the damped pendulum; in Biology and Epidemiology, the models for logistic growth and Bertalanffy, as well as epidemiological models; and in Thermodynamics, the heat partial differential equation. We rely on Karhunen-Loève and generalized polynomial chaos expansions and on finite difference schemes for the density approximation of the solution. These techniques are only applicable when we have a forward model in which the input parameters have certain probability distributions already set. When the model coefficients are estimated from collected data, we have an inverse problem. The Bayesian inference approach allows estimating the probability distribution of the model parameters from their prior probability distribution and the likelihood of the data. Uncertainty quantification for the model output is then carried out using the posterior predictive distribution. In this regard, the last part of the thesis shows the estimation of the distributions of the model parameters from experimental data on bacteria growth. To do so, a hybrid method that combines Bayesian parameter estimation and generalized polynomial chaos expansions is used.[ES] Los modelos matemáticos basados en ecuaciones diferenciales deterministas no tienen en cuenta la incertidumbre inherente del fenómeno físico (en un sentido amplio) bajo estudio. Además, a menudo se producen inexactitudes en los datos recopilados debido a errores en las mediciones. Por lo tanto, es necesario tratar los parámetros de entrada del modelo como cantidades aleatorias, en forma de variables aleatorias o procesos estocásticos. Esto da lugar al estudio de las ecuaciones diferenciales aleatorias. El cálculo de la función de densidad de probabilidad de la solución estocástica es importante en la cuantificación de la incertidumbre de la respuesta del modelo. Aunque dicho cálculo es un objetivo difícil en general, ciertas expansiones estocásticas para los coeficientes del modelo dan lugar a representaciones fieles de la solución estocástica, lo que permite aproximar su función de densidad. En este sentido, las expansiones de Karhunen-Loève y de caos polinomial generalizado constituyen herramientas para dicha aproximación de la densidad. Además, los métodos basados en discretizaciones de esquemas numéricos de diferencias finitas permiten aproximar la solución estocástica, por lo tanto, su función de densidad de probabilidad. La parte principal de esta disertación tiene como objetivo aproximar la función de densidad de probabilidad de modelos matemáticos importantes con incertidumbre en su formulación. Concretamente, en esta memoria se estudian, en un sentido estocástico, los siguientes modelos que aparecen en diferentes áreas científicas: en Física, el modelo del péndulo amortiguado; en Biología y Epidemiología, los modelos de crecimiento logístico y de Bertalanffy, así como modelos de tipo epidemiológico; y en Termodinámica, la ecuación en derivadas parciales del calor. Utilizamos expansiones de Karhunen-Loève y de caos polinomial generalizado y esquemas de diferencias finitas para la aproximación de la densidad de la solución. Estas técnicas solo son aplicables cuando tenemos un modelo directo en el que los parámetros de entrada ya tienen determinadas distribuciones de probabilidad establecidas. Cuando los coeficientes del modelo se estiman a partir de los datos recopilados, tenemos un problema inverso. El enfoque de inferencia Bayesiana permite estimar la distribución de probabilidad de los parámetros del modelo a partir de su distribución de probabilidad previa y la verosimilitud de los datos. La cuantificación de la incertidumbre para la respuesta del modelo se lleva a cabo utilizando la distribución predictiva a posteriori. En este sentido, la última parte de la tesis muestra la estimación de las distribuciones de los parámetros del modelo a partir de datos experimentales sobre el crecimiento de bacterias. Para hacerlo, se utiliza un método híbrido que combina la estimación de parámetros Bayesianos y los desarrollos de caos polinomial generalizado.[CA] Els models matemàtics basats en equacions diferencials deterministes no tenen en compte la incertesa inherent al fenomen físic (en un sentit ampli) sota estudi. A més a més, sovint es produeixen inexactituds en les dades recollides a causa d'errors de mesurament. Es fa així necessari tractar els paràmetres d'entrada del model com a quantitats aleatòries, en forma de variables aleatòries o processos estocàstics. Açò dóna lloc a l'estudi de les equacions diferencials aleatòries. El càlcul de la funció de densitat de probabilitat de la solució estocàstica és important per a quantificar la incertesa de la sortida del model. Tot i que, en general, aquest càlcul és un objectiu difícil d'assolir, certes expansions estocàstiques dels coeficients del model donen lloc a representacions fidels de la solució estocàstica, el que permet aproximar la seua funció de densitat. En aquest sentit, les expansions de Karhunen-Loève i de caos polinomial generalitzat esdevenen eines per a l'esmentada aproximació de la densitat. A més a més, els mètodes basats en discretitzacions mitjançant esquemes numèrics de diferències finites permeten aproximar la solució estocàstica, per tant la seua funció de densitat de probabilitat. La part principal d'aquesta dissertació té com a objectiu aproximar la funció de densitat de probabilitat d'importants models matemàtics amb incerteses en la seua formulació. Concretament, en aquesta memòria s'estudien, en un sentit estocàstic, els següents models que apareixen en diferents àrees científiques: en Física, el model del pèndol amortit; en Biologia i Epidemiologia, els models de creixement logístic i de Bertalanffy, així com models de tipus epidemiològic; i en Termodinàmica, l'equació en derivades parcials de la calor. Per a l'aproximació de la densitat de la solució, ens basem en expansions de Karhunen-Loève i de caos polinomial generalitzat i en esquemes de diferències finites. Aquestes tècniques només són aplicables quan tenim un model cap avant en què els paràmetres d'entrada tenen ja determinades distribucions de probabilitat. Quan els coeficients del model s'estimen a partir de les dades recollides, tenim un problema invers. L'enfocament de la inferència Bayesiana permet estimar la distribució de probabilitat dels paràmetres del model a partir de la seua distribució de probabilitat prèvia i la versemblança de les dades. La quantificació de la incertesa per a la resposta del model es fa mitjançant la distribució predictiva a posteriori. En aquest sentit, l'última part de la tesi mostra l'estimació de les distribucions dels paràmetres del model a partir de dades experimentals sobre el creixement de bacteris. Per a fer-ho, s'utilitza un mètode híbrid que combina l'estimació de paràmetres Bayesiana i els desenvolupaments de caos polinomial generalitzat.This work has been supported by the Spanish Ministerio de Economía y Competitividad grant MTM2017–89664–P.Calatayud Gregori, J. (2020). Computational methods for random differential equations: probability density function and estimation of the parameters [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/138396TESISPremios Extraordinarios de tesis doctorale

Lp-solution to the random linear delay differential equation with stochastic forcing term

[EN] This paper aims at extending a previous contribution dealing with the random autonomous-homogeneous linear differential equation with discrete delay tau > 0, by adding a random forcing term f(t) that varies with time: x'(t) = ax(t) + bx(t-tau) + f(t), t >= 0, with initial condition x(t) = g(t), -tau <= t <= 0. The coefficients a and b are assumed to be random variables, while the forcing term f(t) and the initial condition g(t) are stochastic processes on their respective time domains. The equation is regarded in the Lebesgue space L-p of random variables with finite p-th moment. The deterministic solution constructed with the method of steps and the method of variation of constants, which involves the delayed exponential function, is proved to be an L-p-solution, under certain assumptions on the random data. This proof requires the extension of the deterministic Leibniz's integral rule for differentiation to the random scenario. Finally, we also prove that, when the delay tau tends to 0, the random delay equation tends in L-p to a random equation with no delay. Numerical experiments illustrate how our methodology permits determining the main statistics of the solution process, thereby allowing for uncertainty quantification.This work has been supported by the Spanish Ministerio de Economia, Industria y Competitividad (MINECO), the Agencia Estatal de Investigacion (AEI) and Fondo Europeo de Desarrollo Regional (FEDER UE) grant MTM2017-89664-P.Cortés, J.; Jornet, M. (2020). Lp-solution to the random linear delay differential equation with stochastic forcing term. Mathematics. 8(6):1-16. https://doi.org/10.3390/math8061013S11686Xiu, D., & Karniadakis, G. E. (2004). Supersensitivity due to uncertain boundary conditions. International Journal for Numerical Methods in Engineering, 61(12), 2114-2138. doi:10.1002/nme.1152Casabán, M.-C., Cortés, J.-C., Navarro-Quiles, A., Romero, J.-V., Roselló, M.-D., & Villanueva, R.-J. (2016). A comprehensive probabilistic solution of random SIS-type epidemiological models using the random variable transformation technique. Communications in Nonlinear Science and Numerical Simulation, 32, 199-210. doi:10.1016/j.cnsns.2015.08.009Strand, J. . (1970). Random ordinary differential equations. Journal of Differential Equations, 7(3), 538-553. doi:10.1016/0022-0396(70)90100-2Villafuerte, L., Braumann, C. A., Cortés, J.-C., & Jódar, L. (2010). Random differential operational calculus: Theory and applications. Computers & Mathematics with Applications, 59(1), 115-125. doi:10.1016/j.camwa.2009.08.061Cortés, J.-C., Jódar, L., Roselló, M.-D., & Villafuerte, L. (2012). Solving initial and two-point boundary value linear random differential equations: A mean square approach. Applied Mathematics and Computation, 219(4), 2204-2211. doi:10.1016/j.amc.2012.08.066Calatayud, J., Cortés, J.-C., Jornet, M., & Villafuerte, L. (2018). Random non-autonomous second order linear differential equations: mean square analytic solutions and their statistical properties. Advances in Difference Equations, 2018(1). doi:10.1186/s13662-018-1848-8Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). Improving the Approximation of the First- and Second-Order Statistics of the Response Stochastic Process to the Random Legendre Differential Equation. Mediterranean Journal of Mathematics, 16(3). doi:10.1007/s00009-019-1338-6Licea, J. A., Villafuerte, L., & Chen-Charpentier, B. M. (2013). Analytic and numerical solutions of a Riccati differential equation with random coefficients. Journal of Computational and Applied Mathematics, 239, 208-219. doi:10.1016/j.cam.2012.09.040Burgos, C., Calatayud, J., Cortés, J.-C., & Villafuerte, L. (2018). Solving a class of random non-autonomous linear fractional differential equations by means of a generalized mean square convergent power series. Applied Mathematics Letters, 78, 95-104. doi:10.1016/j.aml.2017.11.009Nouri, K., & Ranjbar, H. (2014). Mean Square Convergence of the Numerical Solution of Random Differential Equations. Mediterranean Journal of Mathematics, 12(3), 1123-1140. doi:10.1007/s00009-014-0452-8Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). Random differential equations with discrete delay. Stochastic Analysis and Applications, 37(5), 699-707. doi:10.1080/07362994.2019.1608833Calatayud, J., Cortés, J.-C., & Jornet, M. (2019). $$\mathrm {L}^p$$-calculus Approach to the Random Autonomous Linear Differential Equation with Discrete Delay. Mediterranean Journal of Mathematics, 16(4). doi:10.1007/s00009-019-1370-6Caraballo, T., Cortés, J.-C., & Navarro-Quiles, A. (2019). Applying the random variable transformation method to solve a class of random linear differential equation with discrete delay. Applied Mathematics and Computation, 356, 198-218. doi:10.1016/j.amc.2019.03.048Zhou, T. (2014). A Stochastic Collocation Method for Delay Differential Equations with Random Input. Advances in Applied Mathematics and Mechanics, 6(4), 403-418. doi:10.4208/aamm.2012.m38Shi, W., & Zhang, C. (2017). Generalized polynomial chaos for nonlinear random delay differential equations. Applied Numerical Mathematics, 115, 16-31. doi:10.1016/j.apnum.2016.12.004Khusainov, D. Y., Ivanov, A. F., & Kovarzh, I. V. (2009). Solution of one heat equation with delay. Nonlinear Oscillations, 12(2), 260-282. doi:10.1007/s11072-009-0075-3Shaikhet, L. (2016). Stability of equilibrium states of a nonlinear delay differential equation with stochastic perturbations. International Journal of Robust and Nonlinear Control, 27(6), 915-924. doi:10.1002/rnc.3605Benhadri, M., & Zeghdoudi, H. (2018). Mean square asymptotic stability in nonlinear stochastic neutral Volterra-Levin equations with Poisson jumps and variable delays. Functiones et Approximatio Commentarii Mathematici, 58(2). doi:10.7169/facm/1657Santonja, F.-J., & Shaikhet, L. (2012). Analysing Social Epidemics by Delayed Stochastic Models. Discrete Dynamics in Nature and Society, 2012, 1-13. doi:10.1155/2012/530472Liu, L., & Caraballo, T. (2018). Analysis of a Stochastic 2D-Navier–Stokes Model with Infinite Delay. Journal of Dynamics and Differential Equations, 31(4), 2249-2274. doi:10.1007/s10884-018-9703-xLupulescu, V., & Abbas, U. (2011). Fuzzy delay differential equations. Fuzzy Optimization and Decision Making, 11(1), 99-111. doi:10.1007/s10700-011-9112-7Krapivsky, P. L., Luck, J. M., & Mallick, K. (2011). On stochastic differential equations with random delay. Journal of Statistical Mechanics: Theory and Experiment, 2011(10), P10008. doi:10.1088/1742-5468/2011/10/p10008GARRIDO-ATIENZA, M. J., OGROWSKY, A., & SCHMALFUSS, B. (2011). RANDOM DIFFERENTIAL EQUATIONS WITH RANDOM DELAYS. Stochastics and Dynamics, 11(02n03), 369-388. doi:10.1142/s0219493711003358Cortés, J.-C., Villafuerte, L., & Burgos, C. (2017). A Mean Square Chain Rule and its Application in Solving the Random Chebyshev Differential Equation. Mediterranean Journal of Mathematics, 14(1). doi:10.1007/s00009-017-0853-6Cortés, J. C., Jódar, L., & Villafuerte, L. (2007). Numerical solution of random differential equations: A mean square approach. Mathematical and Computer Modelling, 45(7-8), 757-765. doi:10.1016/j.mcm.2006.07.017Braumann, C. A., Cortés, J.-C., Jódar, L., & Villafuerte, L. (2018). On the random gamma function: Theory and computing. Journal of Computational and Applied Mathematics, 335, 142-155. doi:10.1016/j.cam.2017.11.045Khusainov, D. Y., & Pokojovy, M. (2015). Solving the Linear 1D Thermoelasticity Equations with Pure Delay. International Journal of Mathematics and Mathematical Sciences, 2015, 1-11. doi:10.1155/2015/47926

Big data clustering: Data preprocessing, variable selection, and dimension reduction

[no abstract available