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

Earthworm management in tropical agroecosystems

In agrosystem field experiments, earthworm inoculation did not impede depletion of soil organic stocks in most cases, in spite of increased carbon inputs through enhanced primary production. Slight evidence of soil organic matter (SOM) protection was found in poorly structured soil, such as a yam plot in Ivory Coast (soil sieved before experimentation), and a pasture plot on Martinique. Aggregation inherited from past earthworm activities probably maintains SOM protection after earthworms have disappeared ; longer term experiments are necessary to observe C dynamics after the disappearance of inherited earthworm structures. In two experiments with maize in Ivory Coast and Peru, the activity of earthworms led to a small increase in the incorporation of organic matter from surface mulch in the SOM. Most of the C incorporated into the SOM originated from root material, and earthworm activities only slightly modified this pattern. Earthworm activity had significant effects on the distribution of C among particle size fractions. The general trend was a depletion of large (greater than 50 micrometers) particles and an accumulation of small (less than 2 micrometers) particles. Nutrient depletion in low-input cropping systems was not impeded by earthworm activities ; at Yurimaguas, some signs of a better conservation of K were noted after 3 years in the traditional rotation. (Résumé d'auteur

Control of membrane barrier during bacterial type-III protein secretion

Type-III secretion systems (T3SSs) of the bacterial flagellum and the evolutionarily related injectisome are capable of translocating proteins with a remarkable speed of several thousand amino acids per second. Here, we investigate how T3SSs are able to transport proteins at such a high rate while preventing the leakage of small molecules. Our mutational and evolutionary analyses demonstrate that an ensemble of conserved methionine residues at the cytoplasmic side of the T3SS channel create a deformable gasket (M-gasket) around fast-moving substrates undergoing export. The unique physicochemical features of the M-gasket are crucial to preserve the membrane barrier, to accommodate local conformational changes during active secretion, and to maintain stability of the secretion pore in cooperation with a plug domain (R-plug) and a network of salt-bridges. The conservation of the M-gasket, R-plug, and salt-bridge network suggests a universal mechanism by which the membrane integrity is maintained during high-speed protein translocation in all T3SSs.Peer Reviewe

Study of Tau-pair Production in Photon-Photon Collisions at LEP and Limits on the Anomalous Electromagnetic Moments of the Tau Lepton

Tau-pair production in the process e+e- -> e+e-tau+tau- was studied using data collected by the DELPHI experiment at LEP2 during the years 1997 - 2000. The corresponding integrated luminosity is 650 pb^{-1}. The values of the cross-section obtained are found to be in agreement with QED predictions. Limits on the anomalous magnetic and electric dipole moments of the tau lepton are deduced.Comment: 20 pages, 9 figures, Accepted by Eur. Phys. J.

CP asymmetry in B→ϕKSB \to \phi K_S in a general two-Higgs-doublet model with fourth-generation quarks

We discuss the time-dependent CP asymmetry of decay $B \to \phi K_S$ in an extension of the Standard Model with both two Higgs doublets and additional fourth-generation quarks. We show that although the Standard Model with two-Higgs-doublet and the Standard model with fourth generation quarks alone are not likely to largely change the effective $\sin 2 \beta$ from the decay of $B \to \phi K_S$, the model with both additional Higgs doublet and fourth-generation quarks can easily account for the possible large negative value of $\sin 2 \beta$ without conflicting with other experimental constraints. In this model, additional large CP violating effects may arise from the flavor changing Yukawa interactions between neutral Higgs bosons and the heavy fourth generation down type quark, which can modify the QCD penguin contributions. With the constraints obtained from $b \to s \bar{s} s$ processes such as $B \to X_s \gamma$ and $\Delta m_{B_s^0}$, this model can lead to the effective $\sin 2 \beta$ to be as large as $- 0.4$ in the CP asymmetry of $B \to \phi K_S$.Comment: 13 pages, 5 figures, references added, to appear in Eur.Phys.J.

Energy dependence of Cronin momentum in saturation model for p+Ap+A and A+AA+A collisions

We calculate $\sqrt{s}$ dependence of Cronin momentum for $p+A$ and $A+A$ collisions in saturation model. We show that this dependence is consistent with expectation from formula which was obtained using simple dimentional consideration. This can be used to test validity of saturation model (and distinguish among its variants) and measure $x$ dependence of saturation momentum from experimental data.Comment: LaTeX2e, 12 pages, 8 figure

Search for composite and exotic fermions at LEP 2

A search for unstable heavy fermions with the DELPHI detector at LEP is reported. Sequential and non-canonical leptons, as well as excited leptons and quarks, are considered. The data analysed correspond to an integrated luminosity of about 48 pb^{-1} at an e^+e^- centre-of-mass energy of 183 GeV and about 20 pb^{-1} equally shared between the centre-of-mass energies of 172 GeV and 161 GeV. The search for pair-produced new leptons establishes 95% confidence level mass limits in the region between 70 GeV/c^2 and 90 GeV/c^2, depending on the channel. The search for singly produced excited leptons and quarks establishes upper limits on the ratio of the coupling of the excited fermio

Impact of lopinavir/ritonavir use on antiretroviral resistance in recent clinical practice

Objectives This observational study was requested by French health authorities to determine the impact of lopinavir/ritonavir (Kaletra®) on antiretroviral resistance in clinical practice. Virological failures of lopinavir/ritonavir and their effects on the resistance to protease inhibitors and reverse transcriptase inhibitors were evaluated in protease inhibitor-experienced patients.Patients and methods Virological failure was defined as an HIV-1 plasma viral load >50 copies/mL after at least 3 months of lopinavir/ritonavir-containing antiretroviral therapy. For all patients, a resistance genotypic test was available at failure and before lopinavir/ritonavir treatment. Data from 72 patients with inclusion criteria were studied. Results The mean viral load at baseline was 4 log10 copies/mL (1.6–6.5). Mutations in the protease gene significantly selected between baseline and failure were L10V, K20R, L33F, M36I, I47V, I54V, A71V and I85V (P < 0.05). Patients who had more than seven protease inhibitor mutations at baseline showed a significantly increased risk of occurrence of protease inhibitor mutations. The proportion of viruses susceptible to atazanavir, fosamprenavir and darunavir decreased significantly between baseline and failure (P < 0.05). Among patients with a virus susceptible to atazanavir at day 0, 26% (n = 14) exhibited a virus resistant or possibly resistant at the time of failure. This proportion was 32% (n = 16) for fosamprenavir and 16% (n = 7) for darunavir. Conclusions A darunavir-based regimen appears to be a sequential option in the case of lopinavir/ritonavir failure. To compare and determine the best treatment sequencing, similar studies should be performed for darunavir/ritonavir and atazanavir/ritonavir

Search for lightest neutralino and stau pair production in light gravitino scenarios with stau NLSP

Promptly decaying lightest neutralinos and long-lived staus are searched for in the context of light gravitino scenarios. It is assumed that the stau is the next to lightest supersymmetric particle (NLSP) and that the lightest neutralino is the next to NLSP (NNLSP). Data collected with the Delphi detector at centre-of-mass energies from 161 to 183 \GeV are analysed. No evidence of the production of these particles is found. Hence, lower mass limits for both kinds of particles are set at 95% C.L.. The mass of gaugino-like neutralinos is found to be greater than 71.5 GeV/c^2. In the search for long-lived stau, masses less than 70.0 to 77.5 \GeVcc are excluded for gravitino masses from 10 to 150 \eVcc . Combining this search with the searches for stable heavy leptons and Minimal Supersymmetric Standard Model staus a lower limit of 68.5 \GeVcc may be set for the stau mas