Multivariate change estimation for a stochastic heat equation from local measurements

Anton Tiepner,Lukas Trottner
2024-09-23
Abstract:We study a stochastic heat equation with piecewise constant diffusivity $\theta$ having a jump at a hypersurface $\Gamma$ that splits the underlying space $[0,1]^d$, $d\geq2,$ into two disjoint sets $\Lambda_-\cup\Lambda_+.$ Based on multiple spatially localized measurement observations on a regular $\delta$-grid of $[0,1]^d$, we propose a joint M-estimator for the diffusivity values and the set $\Lambda_+$ that is inspired by statistical image reconstruction methods. We study convergence of the domain estimator $\hat{\Lambda}_+$ in the vanishing resolution level regime $\delta \to 0$ and with respect to the expected symmetric difference pseudometric. Our main finding is a characterization of the convergence rate for $\hat{\Lambda}_+$ in terms of the complexity of $\Gamma$ measured by the number of intersecting hypercubes from the regular $\delta$-grid. Implications of our general result are discussed under two specific structural assumptions on $\Lambda_+$. For a $\beta$-Hölder smooth boundary fragment $\Gamma$, the set $\Lambda_+$ is estimated with rate $\delta^\beta$. If we assume $\Lambda_+$ to be convex, we obtain a $\delta$-rate. While our approach only aims at optimal domain estimation rates, we also demonstrate consistency of our diffusivity estimators.
Statistics Theory
What problem does this paper attempt to address?