We study multi-patch epidemic models where individuals may migrate from one patch to another in either of the susceptible, exposed/latent, infectious and recovered states. We assume that infections occur both locally with a rate that depends on the patch as well as ``from distance" from all the other patches. The exposed and infectious periods have general distributions, and are not affected by the possible migrations of the individuals. The migration processes in either of the three states are assumed to be Markovian, and independent of the exposed and infectious periods. We establish a functional law of large number (FLLN) and a function central limit theorem (FCLT) for the susceptible, exposed/latent, infectious and recovered processes. In the FLLN, the limit is determined by a set of Volterra integral equations. In the special case of deterministic exposed and infectious periods, the limit becomes a system of ODEs with delays. In the FCLT, the limit is given by a set of stochastic Volterra integral equations driven by a sum of independent Brownian motions and continuous Gaussian processes with an explicit covariance structure.
{"title":"Multi-patch epidemic models with general exposed and infectious periods","authors":"G. Pang, É. Pardoux","doi":"10.1051/ps/2023003","DOIUrl":"https://doi.org/10.1051/ps/2023003","url":null,"abstract":"We study multi-patch epidemic models where individuals may migrate from one patch to another in either of the susceptible, exposed/latent, infectious and recovered states. We assume that infections occur both locally with a rate that depends on the patch as well as ``from distance\" from all the other patches. The exposed and infectious periods have general distributions, and are not affected by the possible migrations of the individuals. The migration processes in either of the three states are assumed to be Markovian, and independent of the exposed and infectious periods. We establish a functional law of large number (FLLN) and a function central limit theorem (FCLT) for the susceptible, exposed/latent, infectious and recovered processes. In the FLLN, the limit is determined by a set of Volterra integral equations. In the special case of deterministic exposed and infectious periods, the limit becomes a system of ODEs with delays. In the FCLT, the limit is given by a set of stochastic Volterra integral equations driven by a sum of independent Brownian motions and continuous Gaussian processes with an explicit covariance structure.","PeriodicalId":51249,"journal":{"name":"Esaim-Probability and Statistics","volume":"50 1","pages":""},"PeriodicalIF":0.4,"publicationDate":"2020-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90650622","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We study a kinetic stochastic model with a non-linear time-inhomogeneous drag force and a Brownian-type random force. More precisely, the Kolmogorov type diffusion [[EQUATION]] is considered : here [[EQUATION]] is the position of the particle and [[EQUATION]] is its velocity and is solution of a stochastic differential equation driven by a one-dimensional Brownian motion, with the drift of the form [[EQUATION]] . The function F satisfies some homogeneity condition and [[EQUATION]] is positive. The behaviour of the process in large time is proved by using stochastic analysis tools.
{"title":"Asymptotic behaviour for a time-inhomogeneous Kolmogorov type diffusion","authors":"M. Gradinaru, Emeline Luirard","doi":"10.1051/ps/2022014","DOIUrl":"https://doi.org/10.1051/ps/2022014","url":null,"abstract":"We study a kinetic stochastic model with a non-linear time-inhomogeneous drag force and a Brownian-type random force. More precisely, the Kolmogorov type diffusion [[EQUATION]] is considered : here [[EQUATION]] is the position of the particle and [[EQUATION]] is its velocity and is solution of a stochastic differential equation driven by a one-dimensional Brownian motion, with the drift of the form [[EQUATION]] . The function F satisfies some homogeneity condition and [[EQUATION]] is positive. The behaviour of the process in large time is proved by using stochastic analysis tools.","PeriodicalId":51249,"journal":{"name":"Esaim-Probability and Statistics","volume":"17 1","pages":""},"PeriodicalIF":0.4,"publicationDate":"2020-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77101056","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This paper is devoted to the estimation of a partial graphical model with a structural Bayesian penalization. Precisely, we are interested in the linear regression setting where the estimation is made through the direct links between potentially high-dimensional predictors and multiple responses, since it is known that Gaussian graphical models enable to exhibit direct links only, whereas coefficients in linear regressions contain both direct and indirect relations (due e.g. to strong correlations among the variables). A smooth penalty reflecting a generalized Gaussian Bayesian prior on the covariates is added, either enforcing patterns (like row structures) in the direct links or regulating the joint influence of predictors. We give a theoretical guarantee for our method, taking the form of an upper bound on the estimation error arising with high probability, provided that the model is suitably regularized. Empirical studies on synthetic data and a real dataset are conducted.
{"title":"A partial graphical model with a structural prior on the direct links between predictors and responses","authors":"Eunice Okome Obiang, Pascal J'ez'equel, F. Proïa","doi":"10.1051/PS/2021010","DOIUrl":"https://doi.org/10.1051/PS/2021010","url":null,"abstract":"This paper is devoted to the estimation of a partial graphical model with a structural Bayesian penalization. Precisely, we are interested in the linear regression setting where the estimation is made through the direct links between potentially high-dimensional predictors and multiple responses, since it is known that Gaussian graphical models enable to exhibit direct links only, whereas coefficients in linear regressions contain both direct and indirect relations (due e.g. to strong correlations among the variables). A smooth penalty reflecting a generalized Gaussian Bayesian prior on the covariates is added, either enforcing patterns (like row structures) in the direct links or regulating the joint influence of predictors. We give a theoretical guarantee for our method, taking the form of an upper bound on the estimation error arising with high probability, provided that the model is suitably regularized. Empirical studies on synthetic data and a real dataset are conducted.","PeriodicalId":51249,"journal":{"name":"Esaim-Probability and Statistics","volume":"11 1","pages":""},"PeriodicalIF":0.4,"publicationDate":"2020-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75899947","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Valentin De Bortoli, A. Desolneux, B. Galerne, Arthur Leclaire
In this paper, we introduce a notion of spatial redundancy in Gaussian random fields. This study is motivated by applications of the a contrario method in image processing. We define similarity functions on local windows in random fields over discrete or continuous domains. We derive explicit Gaussian asymptotics for the distribution of similarity functions when computed on Gaussian random fields. Moreover, for the special case of the squared L2 norm, we give non-asymptotic expressions in both discrete and continuous periodic settings. Finally, we present fast and accurate approximations of these non-asymptotic expressions using moment methods and matrix projections.
{"title":"Redundancy in Gaussian random fields","authors":"Valentin De Bortoli, A. Desolneux, B. Galerne, Arthur Leclaire","doi":"10.1051/PS/2020010","DOIUrl":"https://doi.org/10.1051/PS/2020010","url":null,"abstract":"In this paper, we introduce a notion of spatial redundancy in Gaussian random fields. This study is motivated by applications of the a contrario method in image processing. We define similarity functions on local windows in random fields over discrete or continuous domains. We derive explicit Gaussian asymptotics for the distribution of similarity functions when computed on Gaussian random fields. Moreover, for the special case of the squared L2 norm, we give non-asymptotic expressions in both discrete and continuous periodic settings. Finally, we present fast and accurate approximations of these non-asymptotic expressions using moment methods and matrix projections.","PeriodicalId":51249,"journal":{"name":"Esaim-Probability and Statistics","volume":"95 4 1","pages":""},"PeriodicalIF":0.4,"publicationDate":"2020-03-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80878546","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this paper, we study backward stochastic Volterra integral equations introduced in Lin [Stochastic Anal. Appl. 20 (2002) 165–183] and Yong [Stochastic Process. Appl. 116 (2006) 779–795] and extend the existence, uniqueness or comparison results for general filtration as in Papapantoleon et al. [Electron. J. Probab. 23 (2018) EJP240] (not only Brownian-Poisson setting). We also consider Lp-data and explore the time regularity of the solution in the Itô setting, which is also new in this jump setting.
{"title":"Backward stochastic Volterra integral equations with jumps in a general filtration","authors":"A. Popier","doi":"10.1051/PS/2021006","DOIUrl":"https://doi.org/10.1051/PS/2021006","url":null,"abstract":"In this paper, we study backward stochastic Volterra integral equations introduced in Lin [Stochastic Anal. Appl. 20 (2002) 165–183] and Yong [Stochastic Process. Appl. 116 (2006) 779–795] and extend the existence, uniqueness or comparison results for general filtration as in Papapantoleon et al. [Electron. J. Probab. 23 (2018) EJP240] (not only Brownian-Poisson setting). We also consider Lp-data and explore the time regularity of the solution in the Itô setting, which is also new in this jump setting.","PeriodicalId":51249,"journal":{"name":"Esaim-Probability and Statistics","volume":"89 1","pages":""},"PeriodicalIF":0.4,"publicationDate":"2020-02-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76712206","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We prove a simple criterion of exponential tightness for sequences of Gaussian r.v.’s with values in a separable Banach space from which we deduce a general result of Large Deviations which allows easily to obtain LD estimates in various situations.
{"title":"Tightness and exponential tightness of Gaussian probabilities","authors":"P. Baldi","doi":"10.1051/ps/2020003","DOIUrl":"https://doi.org/10.1051/ps/2020003","url":null,"abstract":"We prove a simple criterion of exponential tightness for sequences of Gaussian r.v.’s with values in a separable Banach space from which we deduce a general result of Large Deviations which allows easily to obtain LD estimates in various situations.","PeriodicalId":51249,"journal":{"name":"Esaim-Probability and Statistics","volume":"28 1","pages":""},"PeriodicalIF":0.4,"publicationDate":"2020-01-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81082360","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
A quenched large deviation principle for Brownian motion in a stationary potential is proved. As the proofs are based on a method developed by Sznitman [Comm. Pure Appl. Math. 47 (1994) 1655–1688] for Brownian motion among obstacles with compact support no regularity conditions on the potential is needed. In particular, the sufficient conditions are verified by potentials with polynomially decaying correlations such as the classical potentials studied by Pastur [Teoret. Mat. Fiz. 32 (1977) 88–95] and Fukushima [J. Stat. Phys. 133 (2008) 639–657] and the potentials recently introduced by Lacoin [Ann. Inst. Henri Poincaré Probab. Stat. 48 (2012) 1010–1028; 1029–1048].
{"title":"Large deviations for Brownian motion in a random potential","authors":"D. Boivin, Thi Thu Hien Lê","doi":"10.1051/ps/2020007","DOIUrl":"https://doi.org/10.1051/ps/2020007","url":null,"abstract":"A quenched large deviation principle for Brownian motion in a stationary potential is proved. As the proofs are based on a method developed by Sznitman [Comm. Pure Appl. Math. 47 (1994) 1655–1688] for Brownian motion among obstacles with compact support no regularity conditions on the potential is needed. In particular, the sufficient conditions are verified by potentials with polynomially decaying correlations such as the classical potentials studied by Pastur [Teoret. Mat. Fiz. 32 (1977) 88–95] and Fukushima [J. Stat. Phys. 133 (2008) 639–657] and the potentials recently introduced by Lacoin [Ann. Inst. Henri Poincaré Probab. Stat. 48 (2012) 1010–1028; 1029–1048].","PeriodicalId":51249,"journal":{"name":"Esaim-Probability and Statistics","volume":"48 1","pages":""},"PeriodicalIF":0.4,"publicationDate":"2020-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78244589","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Consider a continuous time particle system ηt = (ηt(k), k ∈ 𝕃), indexed by a lattice 𝕃 which will be either ℤ, ℤ∕nℤ, a segment {1, ⋯ , n}, or ℤd, and taking its values in the set Eκ𝕃 where Eκ = {0, ⋯ , κ − 1} for some fixed κ ∈{∞, 2, 3, ⋯ }. Assume that the Markovian evolution of the particle system (PS) is driven by some translation invariant local dynamics with bounded range, encoded by a jump rate matrix ⊤. These are standard settings, satisfied by the TASEP, the voter models, the contact processes. The aim of this paper is to provide some sufficient and/or necessary conditions on the matrix ⊤ so that this Markov process admits some simple invariant distribution, as a product measure (if 𝕃 is any of the spaces mentioned above), the law of a Markov process indexed by ℤ or [1, n] ∩ ℤ (if 𝕃 = ℤ or {1, …, n}), or a Gibbs measure if 𝕃 = ℤ/nℤ. Multiple applications follow: efficient ways to find invariant Markov laws for a given jump rate matrix or to prove that none exists. The voter models and the contact processes are shown not to possess any Markov laws as invariant distribution (for any memory m). (As usual, a random process X indexed by ℤ or ℕ is said to be a Markov chain with memory m ∈ {0, 1, 2, ⋯ } if ℙ(Xk ∈ A | Xk−i, i ≥ 1) = ℙ(Xk ∈ A | Xk−i, 1 ≤ i ≤ m), for any k.) We also prove that some models close to these models do. We exhibit PS admitting hidden Markov chains as invariant distribution and design many PS on ℤ2, with jump rates indexed by 2 × 2 squares, admitting product invariant measures.
{"title":"Invariant measures of interacting particle systems: Algebraic aspects","authors":"Luis Fredes, J. Marckert","doi":"10.1051/ps/2020008","DOIUrl":"https://doi.org/10.1051/ps/2020008","url":null,"abstract":"Consider a continuous time particle system ηt = (ηt(k), k ∈ 𝕃), indexed by a lattice 𝕃 which will be either ℤ, ℤ∕nℤ, a segment {1, ⋯ , n}, or ℤd, and taking its values in the set Eκ𝕃 where Eκ = {0, ⋯ , κ − 1} for some fixed κ ∈{∞, 2, 3, ⋯ }. Assume that the Markovian evolution of the particle system (PS) is driven by some translation invariant local dynamics with bounded range, encoded by a jump rate matrix ⊤. These are standard settings, satisfied by the TASEP, the voter models, the contact processes. The aim of this paper is to provide some sufficient and/or necessary conditions on the matrix ⊤ so that this Markov process admits some simple invariant distribution, as a product measure (if 𝕃 is any of the spaces mentioned above), the law of a Markov process indexed by ℤ or [1, n] ∩ ℤ (if 𝕃 = ℤ or {1, …, n}), or a Gibbs measure if 𝕃 = ℤ/nℤ. Multiple applications follow: efficient ways to find invariant Markov laws for a given jump rate matrix or to prove that none exists. The voter models and the contact processes are shown not to possess any Markov laws as invariant distribution (for any memory m). (As usual, a random process X indexed by ℤ or ℕ is said to be a Markov chain with memory m ∈ {0, 1, 2, ⋯ } if ℙ(Xk ∈ A | Xk−i, i ≥ 1) = ℙ(Xk ∈ A | Xk−i, 1 ≤ i ≤ m), for any k.) We also prove that some models close to these models do. We exhibit PS admitting hidden Markov chains as invariant distribution and design many PS on ℤ2, with jump rates indexed by 2 × 2 squares, admitting product invariant measures.","PeriodicalId":51249,"journal":{"name":"Esaim-Probability and Statistics","volume":"41 1","pages":""},"PeriodicalIF":0.4,"publicationDate":"2020-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75833649","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Given p ∈ (1, 2), we study 𝕃p -solutions of a reflected backward stochastic differential equation with jumps (RBSDEJ) whose generator g is Lipschitz continuous in (y , z , u ). Based on a general comparison theorem as well as the optimal stopping theory for uniformly integrable processes under jump filtration, we show that such a RBSDEJ with p -integrable parameters admits a unique 𝕃p solution via a fixed-point argument. The Y -component of the unique 𝕃p solution can be viewed as the Snell envelope of the reflecting obstacle 𝔏 under g -evaluations, and the first time Y meets 𝔏 is an optimal stopping time for maximizing the g -evaluation of reward 𝔏.
{"title":"𝕃p solutions of reflected backward stochastic differential equations with jumps","authors":"Song Yao","doi":"10.1051/ps/2020026","DOIUrl":"https://doi.org/10.1051/ps/2020026","url":null,"abstract":"Given p ∈ (1, 2), we study 𝕃p -solutions of a reflected backward stochastic differential equation with jumps (RBSDEJ) whose generator g is Lipschitz continuous in (y , z , u ). Based on a general comparison theorem as well as the optimal stopping theory for uniformly integrable processes under jump filtration, we show that such a RBSDEJ with p -integrable parameters admits a unique 𝕃p solution via a fixed-point argument. The Y -component of the unique 𝕃p solution can be viewed as the Snell envelope of the reflecting obstacle 𝔏 under g -evaluations, and the first time Y meets 𝔏 is an optimal stopping time for maximizing the g -evaluation of reward 𝔏.","PeriodicalId":51249,"journal":{"name":"Esaim-Probability and Statistics","volume":"75 1","pages":"935-962"},"PeriodicalIF":0.4,"publicationDate":"2020-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91253160","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
J. Azaïs, F. Bachoc, A. Lagnoux, Thi Mong Ngoc Nguyen
We consider the semi-parametric estimation of the scale parameter of the variogram of a one-dimensional Gaussian process with known smoothness. We suggest an estimator based both on quadratic variations and the moment method. We provide asymptotic approximations of the mean and variance of this estimator, together with asymptotic normality results, for a large class of Gaussian processes. We allow for general mean functions, provide minimax upper bounds and study the aggregation of several estimators based on various variation sequences. In extensive simulation studies, we show that the asymptotic results accurately depict the finite-sample situations already for small to moderate sample sizes. We also compare various variation sequences and highlight the efficiency of the aggregation procedure.
{"title":"Semi-parametric estimation of the variogram scale parameter of a Gaussian process with stationary increments","authors":"J. Azaïs, F. Bachoc, A. Lagnoux, Thi Mong Ngoc Nguyen","doi":"10.1051/PS/2020021","DOIUrl":"https://doi.org/10.1051/PS/2020021","url":null,"abstract":"We consider the semi-parametric estimation of the scale parameter of the variogram of a one-dimensional Gaussian process with known smoothness. We suggest an estimator based both on quadratic variations and the moment method. We provide asymptotic approximations of the mean and variance of this estimator, together with asymptotic normality results, for a large class of Gaussian processes. We allow for general mean functions, provide minimax upper bounds and study the aggregation of several estimators based on various variation sequences. In extensive simulation studies, we show that the asymptotic results accurately depict the finite-sample situations already for small to moderate sample sizes. We also compare various variation sequences and highlight the efficiency of the aggregation procedure.","PeriodicalId":51249,"journal":{"name":"Esaim-Probability and Statistics","volume":"36 1","pages":""},"PeriodicalIF":0.4,"publicationDate":"2020-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75221144","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}