Julian Hofstadler , Krzysztof Łatuszyński , Gareth Roberts , Daniel Rudolf
{"title":"Almost sure convergence rates of adaptive increasingly rare Markov chain Monte Carlo","authors":"Julian Hofstadler , Krzysztof Łatuszyński , Gareth Roberts , Daniel Rudolf","doi":"10.1016/j.spa.2026.104905","DOIUrl":"10.1016/j.spa.2026.104905","url":null,"abstract":"<div><div>We consider adaptive increasingly rare Markov chain Monte Carlo (MCMC) algorithms, which are adaptive MCMC methods, where the adaptation concerning the “past” happens less and less frequently over time. Under a contraction assumption with respect to a Wasserstein-like function we deduce upper bounds of the convergence rate of Monte Carlo sums taking a renormalisation factor into account that is “almost” the one that appears in a law of the iterated logarithm. We demonstrate the applicability of our results by considering different settings, among which are those of simultaneous geometric and uniform ergodicity. All proofs are carried out on an augmented state space, including the classical non-augmented setting as a special case. In contrast to other adaptive MCMC limit theory, some technical assumptions, like diminishing adaptation, are not needed.</div></div>","PeriodicalId":51160,"journal":{"name":"Stochastic Processes and their Applications","volume":"196 ","pages":"Article 104905"},"PeriodicalIF":1.2,"publicationDate":"2026-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146154144","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Conditional delta-method for resampling empirical processes in multiple sample problems","authors":"Merle Munko , Dennis Dobler","doi":"10.1016/j.spa.2026.104885","DOIUrl":"10.1016/j.spa.2026.104885","url":null,"abstract":"<div><div>The functional delta-method has a wide range of applications in statistics. Applications on functionals of empirical processes yield various limit results for classical statistics. To improve the finite sample properties of statistical inference procedures that are based on the limit results, resampling procedures such as random permutation and bootstrap methods are a popular solution. In order to analyze the behaviour of the functionals of the resampling empirical processes, corresponding conditional functional delta-methods are desirable. While conditional functional delta-methods for some special cases already exist, there is a lack of more general conditional functional delta-methods for resampling procedures as the permutation and pooled bootstrap method. This gap is addressed in the present paper. Thereby, a general multiple sample problem is considered. The flexible application of the developed conditional delta-method is shown in various relevant examples.</div></div>","PeriodicalId":51160,"journal":{"name":"Stochastic Processes and their Applications","volume":"195 ","pages":"Article 104885"},"PeriodicalIF":1.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146023276","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Tube formula for spherically contoured random fields with subexponential marginals","authors":"Satoshi Kuriki , Evgeny Spodarev","doi":"10.1016/j.spa.2025.104858","DOIUrl":"10.1016/j.spa.2025.104858","url":null,"abstract":"<div><div>It is widely known that the tube method, or equivalently the Euler characteristic heuristic, provides a very accurate approximation for the tail probability that the supremum of a smooth Gaussian random field exceeds a threshold value <em>c</em>. The relative approximation error Δ(<em>c</em>) is exponentially small as a function of <em>c</em> when <em>c</em> tends to infinity. On the other hand, little is known about non-Gaussian random fields.</div><div>In this paper, we obtain the approximation error of the tube method applied to the canonical isotropic random fields on a unit sphere defined by <em>u</em>↦⟨<em>u, ξ</em>⟩, <span><math><mrow><mi>u</mi><mo>∈</mo><mi>M</mi><mo>⊂</mo><msup><mi>S</mi><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow></msup></mrow></math></span>, where <em>ξ</em> is a spherically contoured random vector. These random fields have statistical applications in multiple testing and simultaneous regression inference when the unknown variance is estimated. The decay rate of the relative error Δ(<em>c</em>) depends on the tail of the distribution of ‖<em>ξ</em>‖<sup>2</sup> and the critical radius of the index set <em>M</em>. If this distribution is subexponential but not regularly varying, Δ(<em>c</em>) → 0 as <em>c</em> → ∞. However, in the regularly varying case, Δ(<em>c</em>) does not vanish and hence is not negligible. To address this limitation, we provide simple upper and lower bounds for Δ(<em>c</em>) and for the tube formula itself. Numerical studies are conducted to assess the accuracy of the asymptotic approximation.</div></div>","PeriodicalId":51160,"journal":{"name":"Stochastic Processes and their Applications","volume":"195 ","pages":"Article 104858"},"PeriodicalIF":1.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145886100","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A law of large numbers concerning the distribution of critical points of random Fourier series","authors":"Qiangang “Brandon” Fu, Liviu I. Nicolaescu","doi":"10.1016/j.spa.2026.104899","DOIUrl":"10.1016/j.spa.2026.104899","url":null,"abstract":"<div><div>On the flat torus <span><math><mrow><msup><mi>T</mi><mi>m</mi></msup><mo>=</mo><msup><mi>R</mi><mi>m</mi></msup><mo>/</mo><msup><mi>Z</mi><mi>m</mi></msup></mrow></math></span> we consider the Gaussian random function <span><math><msubsup><mi>F</mi><mi>a</mi><mi>R</mi></msubsup></math></span> defined as a random Fourier series (1.1). The Fourier coefficients are mean zero independent normal variables whose variances depend on the frequencies via an even Schwartz function <span><math><mi>a</mi></math></span> on <span><math><mi>R</mi></math></span> and large rescaling parameter <em>R</em>. For any open subset <em>U</em> of the torus denote by <em>Z<sub>R</sub></em>(<em>U</em>) the number of critical points of <span><math><msubsup><mi>F</mi><mi>a</mi><mi>R</mi></msubsup></math></span> in <em>U</em>. We prove that if <em>U</em> is contained in a geodesic ball, then the variance of <em>Z<sub>R</sub></em>(<em>U</em>) is asymptotic to <em>const</em> × <em>R<sup>m</sup>vol</em>[<em>U</em>] as <em>R</em> → ∞. We use this to prove that if <em>m</em> ≥ 2, then as <em>N</em> → ∞, the random measures <span><math><mrow><msup><mi>N</mi><mrow><mo>−</mo><mi>m</mi></mrow></msup><msub><mi>Z</mi><mi>N</mi></msub><mrow><mo>(</mo><mo>−</mo><mo>)</mo></mrow></mrow></math></span> converge a.s. to an explicit multiple of the volume measure on the flat torus.</div></div>","PeriodicalId":51160,"journal":{"name":"Stochastic Processes and their Applications","volume":"195 ","pages":"Article 104899"},"PeriodicalIF":1.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146173740","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Large-sample analysis of cost functionals for inference under the coalescent","authors":"Martina Favero , Jere Koskela","doi":"10.1016/j.spa.2026.104894","DOIUrl":"10.1016/j.spa.2026.104894","url":null,"abstract":"<div><div>The coalescent is a foundational model of latent genealogical trees under neutral evolution, but suffers from intractable sampling probabilities. Methods for approximating these sampling probabilities either introduce bias or fail to scale to large sample sizes. We show that a class of cost functionals of the coalescent with recurrent mutation and a finite number of alleles converge to tractable processes in the infinite-sample limit. A particular choice of costs yields insight about importance sampling methods, which are a classical tool for coalescent sampling probability approximation. These insights reveal that the behaviour of coalescent importance sampling algorithms differs markedly from standard sequential importance samplers, with or without resampling. We conduct a simulation study to verify that our asymptotics are accurate for algorithms with finite (and moderate) sample sizes. Our results constitute the first theoretical description of large-sample importance sampling algorithms for the coalescent, provide heuristics for the a priori optimisation of computational effort, and identify settings where resampling is harmful for algorithm performance. We observe strikingly different behaviour for importance sampling methods under the infinite sites model of mutation, which is regarded as a good and more tractable approximation of finite alleles mutation in most respects.</div></div>","PeriodicalId":51160,"journal":{"name":"Stochastic Processes and their Applications","volume":"195 ","pages":"Article 104894"},"PeriodicalIF":1.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146078392","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"On the 1/H-variation of the divergence integral with respect to a Hermite process","authors":"Petr Čoupek, Pavel Kříž, Matěj Svoboda","doi":"10.1016/j.spa.2026.104891","DOIUrl":"10.1016/j.spa.2026.104891","url":null,"abstract":"<div><div>In this paper, a divergence-type integral of a random integrand with respect to the Hermite process of order <span><math><mrow><mi>k</mi><mo>∈</mo><mi>N</mi></mrow></math></span> with Hurst parameter <em>H</em> ∈ (1/2, 1) is defined and it is shown that the integral is of finite 1/<em>H</em>-variation.</div></div>","PeriodicalId":51160,"journal":{"name":"Stochastic Processes and their Applications","volume":"195 ","pages":"Article 104891"},"PeriodicalIF":1.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146078432","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Numerical approximation of ergodic BSDEs using non linear Feynman-Kac formulas","authors":"Emmanuel Gobet , Adrien Richou , Lukasz Szpruch","doi":"10.1016/j.spa.2026.104871","DOIUrl":"10.1016/j.spa.2026.104871","url":null,"abstract":"<div><div>In this work we study the numerical approximation of a class of ergodic Backward Stochastic Differential Equations. These equations are formulated in an infinite horizon framework and provide a probabilistic representation for elliptic Partial Differential Equations of ergodic type. In order to build our numerical scheme, we put forward a new representation of the PDE solution by using a classical probabilistic representation of the gradient. Then, based on this representation, we propose a fully implementable numerical scheme using a Picard iteration procedure, a grid space discretization and a Monte-Carlo approximation. Up to a limiting technical condition that guarantees the contraction of the Picard procedure, we obtain an upper bound for the numerical error. We also provide some numerical experiments that show the efficiency of this approach for small dimensions.</div></div>","PeriodicalId":51160,"journal":{"name":"Stochastic Processes and their Applications","volume":"195 ","pages":"Article 104871"},"PeriodicalIF":1.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145978772","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Sharp asymptotics for N-point correlation functions of coalescing heavy-tailed random walks","authors":"Jinjiong Yu","doi":"10.1016/j.spa.2026.104897","DOIUrl":"10.1016/j.spa.2026.104897","url":null,"abstract":"<div><div>We study a system of coalescing continuous-time random walks starting from every site on <span><math><mi>Z</mi></math></span>, where the jump increments lie in the domain of attraction of an <em>α</em>-stable distribution with <em>α</em> ∈ (0, 1]. We establish sharp asymptotics for the <em>N</em>-point correlation function of the system. Our analysis relies on two precise tail estimates for the system density, as well as the non-collision probability of <em>N</em> independent random walks with arbitrary fixed initial configurations. In addition, we derive refined estimates for heavy-tailed random walks, which may be of independent interest.</div></div>","PeriodicalId":51160,"journal":{"name":"Stochastic Processes and their Applications","volume":"195 ","pages":"Article 104897"},"PeriodicalIF":1.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146173741","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Moments of generalized fractional polynomial processes","authors":"Johannes Assefa, Martin Keller-Ressel","doi":"10.1016/j.spa.2026.104901","DOIUrl":"10.1016/j.spa.2026.104901","url":null,"abstract":"<div><div>We derive a moment formula for generalized fractional polynomial processes, i.e., for polynomial-preserving Markov processes time-changed by an inverse Lévy-subordinator. If the time change is inverse <em>α</em>-stable, the time-derivative of the Kolmogorov backward equation is replaced by a Caputo fractional derivative of order <em>α</em>, and we demonstrate that moments of such processes are computable, in a closed form, using matrix Mittag-Leffler functions. The same holds true for cross-moments in equilibrium, generalizing results of Leonenko, Meerschaert and Sikorskii from the one-dimensional diffusive case of second-order moments to the multivariate, jump-diffusive case of moments of arbitrary order. We show that also in this more general setting, fractional polynomial processes exhibit long-range dependence, with correlations decaying as a power law with exponent <em>α</em>.</div></div>","PeriodicalId":51160,"journal":{"name":"Stochastic Processes and their Applications","volume":"195 ","pages":"Article 104901"},"PeriodicalIF":1.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146173739","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A new stochastic SIS-type modelling framework for analysing epidemic dynamics in continuous space","authors":"Apolline Louvet , Bastian Wiederhold","doi":"10.1016/j.spa.2026.104896","DOIUrl":"10.1016/j.spa.2026.104896","url":null,"abstract":"<div><div>We propose a new stochastic epidemiological model defined in a continuous space of arbitrary dimension, based on SIS dynamics implemented in a spatial Λ-Fleming-Viot (SLFV) process. The model can be described by as little as three parameters, and is dual to a spatial branching process with competition linked to genealogies of infected individuals. Therefore, it is a possible modelling framework to develop computationally tractable inference tools for epidemics in a continuous space using demographic and genetic data.</div><div>We provide mathematical constructions of the process based on well-posed martingale problems as well as driving space-time Poisson point processes. With these devices and the duality relation in hand, we unveil some of the drivers of the transition between extinction and survival of the epidemic. In particular, we show that extinction is in large parts independent of the initial condition, and identify a strong candidate for the reproduction number R<sub>0</sub> of the epidemic in such a model.</div></div>","PeriodicalId":51160,"journal":{"name":"Stochastic Processes and their Applications","volume":"195 ","pages":"Article 104896"},"PeriodicalIF":1.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146173743","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}