<p>From the observation of a diffusion path <InlineEquation ID="IEq1"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq1.gif" Format="GIF" Height="21" Rendition="HTML" Resolution="72" Type="Linedraw" Width="69" /> </InlineMediaObject> <EquationSource Format="TEX">\((X_t)_{t\in [0,T]}\)</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mrow> <mo stretchy="false">(</mo> <msub> <mi>X</mi> <mi>t</mi> </msub> <mo stretchy="false">)</mo> </mrow> <mrow> <mi>t</mi> <mo>∈</mo> <mo stretchy="false">[</mo> <mn>0</mn> <mo>,</mo> <mi>T</mi> <mo stretchy="false">]</mo> </mrow> </msub> </math></EquationSource> </InlineEquation> on a compact connected <i>d</i>-dimensional manifold <InlineEquation ID="IEq2"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq2.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="25" /> </InlineMediaObject> <EquationSource Format="TEX">\(\mathcal {M}\)</EquationSource> <EquationSource Format="MATHML"><math> <mi mathvariant="script">M</mi> </math></EquationSource> </InlineEquation> without boundary, we consider the problem of estimating the stationary measure <InlineEquation ID="IEq3"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq3.gif" Format="GIF" Height="12" Rendition="HTML" Resolution="72" Type="Linedraw" Width="15" /> </InlineMediaObject> <EquationSource Format="TEX">\(\mu \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>μ</mi> </math></EquationSource> </InlineEquation> of the process. Wang and Zhu (2023) showed that for the Wasserstein metric <InlineEquation ID="IEq4"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq4.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="25" /> </InlineMediaObject> <EquationSource Format="TEX">\(\mathcal {W}_2\)</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi mathvariant="script">W</mi> <mn>2</mn> </msub> </math></EquationSource> </InlineEquation> and for <InlineEquation ID="IEq5"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq5.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="41" /> </InlineMediaObject> <EquationSource Format="TEX">\(d\ge 5\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>d</mi> <mo>≥</mo> <mn>5</mn> </mrow> </math></EquationSource> </InlineEquation>, the convergence rate of <InlineEquation ID="IEq6"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq6.gif" Format="GIF" Height="18" Rendition="HTML" Resolution="72" Type="Linedraw" Width="65" /> </InlineMediaObject> <EquationSource Format="TEX">\(T^{-1/(d-2)}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>T</mi> <mrow> <mo>-</mo> <mn>1</mn> <mo stretchy="false">/</mo> <mo stretchy="false">(</mo> <mi>d</mi> <mo>-</mo> <mn>2</mn> <mo stretchy="false">)</mo> </mrow> </msup> </math></EquationSource> </InlineEquation> is attained by the occupation measure of the path <InlineEquation ID="IEq7"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq1.gif" Format="GIF" Height="21" Rendition="HTML" Resolution="72" Type="Linedraw" Width="69" /> </InlineMediaObject> <EquationSource Format="TEX">\((X_t)_{t\in [0,T]}\)</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mrow> <mo stretchy="false">(</mo> <msub> <mi>X</mi> <mi>t</mi> </msub> <mo stretchy="false">)</mo> </mrow> <mrow> <mi>t</mi> <mo>∈</mo> <mo stretchy="false">[</mo> <mn>0</mn> <mo>,</mo> <mi>T</mi> <mo stretchy="false">]</mo> </mrow> </msub> </math></EquationSource> </InlineEquation> when <InlineEquation ID="IEq8"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq1.gif" Format="GIF" Height="21" Rendition="HTML" Resolution="72" Type="Linedraw" Width="69" /> </InlineMediaObject> <EquationSource Format="TEX">\((X_t)_{t\in [0,T]}\)</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mrow> <mo stretchy="false">(</mo> <msub> <mi>X</mi> <mi>t</mi> </msub> <mo stretchy="false">)</mo> </mrow> <mrow> <mi>t</mi> <mo>∈</mo> <mo stretchy="false">[</mo> <mn>0</mn> <mo>,</mo> <mi>T</mi> <mo stretchy="false">]</mo> </mrow> </msub> </math></EquationSource> </InlineEquation> is a Langevin diffusion. We extend their result in several directions. First, we show that the rate of convergence holds for a large class of diffusion paths, whose generators are uniformly elliptic. Second, the regularity of the density <i>p</i> of the stationary measure <InlineEquation ID="IEq9"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq3.gif" Format="GIF" Height="12" Rendition="HTML" Resolution="72" Type="Linedraw" Width="15" /> </InlineMediaObject> <EquationSource Format="TEX">\(\mu \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>μ</mi> </math></EquationSource> </InlineEquation> with respect to the volume measure of <InlineEquation ID="IEq10"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq2.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="25" /> </InlineMediaObject> <EquationSource Format="TEX">\(\mathcal {M}\)</EquationSource> <EquationSource Format="MATHML"><math> <mi mathvariant="script">M</mi> </math></EquationSource> </InlineEquation> can be leveraged to obtain faster estimators: when <i>p</i> belongs to a Sobolev space of order <InlineEquation ID="IEq11"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq11.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="40" /> </InlineMediaObject> <EquationSource Format="TEX">\(\ell \ge 2\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>ℓ</mi> <mo>≥</mo> <mn>2</mn> </mrow> </math></EquationSource> </InlineEquation>, smoothing the occupation measure by convolution with a kernel yields an estimator whose rate of convergence is of order <InlineEquation ID="IEq12"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="440_2025_1437_Article_IEq12.gif" Format="GIF" Height="18" Rendition="HTML" Resolution="72" Type="Linedraw" Width="108" /> </InlineMediaObject> <EquationSource Format="TEX">\(T^{-(\ell +1)/(2\ell +d-2)}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>T</mi> <mrow> <mo>-</mo> <mo stretchy="false">(</mo> <mi>ℓ</mi> <mo>+</mo> <mn>1</mn> <mo stretchy="false">)</mo> <mo stretchy="false">/</mo> <mo stretchy="false">(</mo> <mn>2</mn> <mi>ℓ</mi> <mo>+</mo> <mi>d</mi> <mo>-</mo> <mn>2</mn> <mo stretchy="false">)</mo> </mrow> </msup> </math></EquationSource> </InlineEquation>. We further show that this rate is the minimax rate of estimation for this problem.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Measure estimation on a manifold explored by a diffusion process

  • Vincent Divol,
  • Hélène Guérin,
  • Dinh–Toan Nguyen,
  • Viet Chi Tran

摘要

From the observation of a diffusion path \((X_t)_{t\in [0,T]}\) ( X t ) t [ 0 , T ] on a compact connected d-dimensional manifold \(\mathcal {M}\) M without boundary, we consider the problem of estimating the stationary measure \(\mu \) μ of the process. Wang and Zhu (2023) showed that for the Wasserstein metric \(\mathcal {W}_2\) W 2 and for \(d\ge 5\) d 5 , the convergence rate of \(T^{-1/(d-2)}\) T - 1 / ( d - 2 ) is attained by the occupation measure of the path \((X_t)_{t\in [0,T]}\) ( X t ) t [ 0 , T ] when \((X_t)_{t\in [0,T]}\) ( X t ) t [ 0 , T ] is a Langevin diffusion. We extend their result in several directions. First, we show that the rate of convergence holds for a large class of diffusion paths, whose generators are uniformly elliptic. Second, the regularity of the density p of the stationary measure \(\mu \) μ with respect to the volume measure of \(\mathcal {M}\) M can be leveraged to obtain faster estimators: when p belongs to a Sobolev space of order \(\ell \ge 2\) 2 , smoothing the occupation measure by convolution with a kernel yields an estimator whose rate of convergence is of order \(T^{-(\ell +1)/(2\ell +d-2)}\) T - ( + 1 ) / ( 2 + d - 2 ) . We further show that this rate is the minimax rate of estimation for this problem.