We study the local asymptotic normality (LAN) property for the likelihood function associated with discretely observed $d$-dimensional McKean-Vlasov stochastic differential equations over a fixed time interval. The model involves a joint parameter in both the drift and diffusion coefficients, introducing challenges due to its dependence on the process distribution. We derive a stochastic expansion of the log-likelihood ratio using Malliavin calculus techniques and establish the LAN property under appropriate conditions. The main technical challenge arises from the implicit nature of the transition densities, which we address through integration by parts and Gaussian-type bounds. This work extends existing LAN results for interacting particle systems to the mean-field regime, contributing to statistical inference in non-linear stochastic models
معرّف الورقة : 2511.13366العنوان : Local asymptotic normality for discretely observed McKean-Vlasov diffusionsالمؤلفون : Akram Heidari, Mark Podolskij (جامعة لوكسمبرغ)التصنيف : math.ST, stat.TH (النظرية الإحصائية)تاريخ الإرسال : 17 نوفمبر 2025رابط الورقة : https://arxiv.org/abs/2511.13366 تدرس هذه الورقة خاصية الحالة الطبيعية المقاربة المحلية (LAN) لدالة الاحتمالية لمعادلات McKean-Vlasov العشوائية التفاضلية ذات البعد d المرصودة بشكل منفصل على فترة زمنية ثابتة. يتضمن النموذج معاملات مشتركة في كل من معاملات الانجراف والانتشار، مما يقدم تحديات بسبب الاعتماد على توزيع العملية. يستخدم المؤلفون تقنيات حساب Malliavin لاشتقاق توسع عشوائي لنسبة اللوغاريتم الاحتمالي وإنشاء خاصية LAN في ظل شروط مناسبة. تنشأ التحديات التقنية الرئيسية من الطبيعة الضمنية لكثافة الانتقال، والتي يتم حلها من خلال التكامل بالأجزاء والحدود من النوع الغاوسي. يوسع هذا العمل نتائج LAN الموجودة للأنظمة الجزيئية التفاعلية إلى آليات المجال المتوسط، مما يساهم في الاستدلال الإحصائي للنماذج العشوائية غير الخطية.
تدرس هذه الورقة مشكلة تقدير المعاملات لمعادلات McKean-Vlasov العشوائية التفاضلية (SDEs)، وتركز بشكل محدد على إنشاء الحالة الطبيعية المقاربة المحلية (LAN) لدالة الاحتمالية بيانات الملاحظة في الوقت المنفصل. تأخذ معادلة McKean-Vlasov الشكل:
d X t i , θ = b θ 1 ( X t i , θ , μ t θ ) d t + a θ 2 ( X t i , θ ) d W t i dX^{i,\theta}_t = b_{\theta_1}(X^{i,\theta}_t, \mu^\theta_t)dt + a_{\theta_2}(X^{i,\theta}_t)dW^i_t d X t i , θ = b θ 1 ( X t i , θ , μ t θ ) d t + a θ 2 ( X t i , θ ) d W t i
حيث μ t θ \mu^\theta_t μ t θ هو توزيع X t i , θ X^{i,\theta}_t X t i , θ ، مما يجعل المعادلة ذات خاصية لاخطية متأصلة.
الانتشار الواسع للتطبيقات : معادلات McKean-Vlasov لها تطبيقات واسعة في الفيزياء الإحصائية والتمويل وألعاب المجال المتوسط وغيرهاالأساس النظري : خاصية LAN هي أداة أساسية للاستدلال الإحصائي المقارب، قدمها Le Cam، ويمكن استخدامها لاشتقاق الحدود الدنيا لتباين المقدرات المقاربةنظرية المجال المتوسط : ربط الاستدلال الإحصائي للأنظمة الجزيئية الدقيقة بحدود المجال المتوسط الكليةالملاحظة المستمرة مقابل المنفصلة : تركز نتائج LAN الموجودة بشكل أساسي على حالة الملاحظة المستمرة 13 ، حيث يمكن الحصول على تعبير مغلق لدالة الاحتمالية من خلال نظرية Girsanovالطبيعة الضمنية لكثافة الانتقال : في حالة الملاحظة المنفصلة، لا تحتوي كثافة الانتقال على تعبير صريح، مما يتطلب تقنيات جديدةتحديات الأنظمة الجزيئية التفاعلية : بالنسبة للأنظمة الجزيئية التفاعلية (3.18)، يتطلب التعامل مع كثافة انتقال مشتركة ذات بعد dN، وهناك نقص في الأدب الحالي لحدود مناسبةملء الفراغ في نظرية LAN لمعادلات McKean-Vlasov المرصودة بشكل منفصل تطوير تقنيات حساب Malliavin للتعامل مع كثافات الانتقال الضمنية توفير أساس نظري للاستدلال الإحصائي لنماذج المجال المتوسط، وإنشاء اتصال نظري مع الطرق الحديثة للتقدير 1 إنشاء خاصية LAN : للمرة الأولى، يتم إنشاء خاصية LAN لمعادلات McKean-Vlasov المرصودة بشكل منفصل، تحت آلية التقارب Δ n → 0 , N → ∞ \Delta_n \to 0, N \to \infty Δ n → 0 , N → ∞ مع الفترة الزمنية T T T ثابتةتقنيات حساب Malliavin : استخدام صيغة التكامل بالأجزاء في حساب Malliavin لاشتقاق تمثيل صريح لمشتق لوغاريتم كثافة الانتقال (الاقتراح 3.1)التوسع العشوائي : إنشاء توسع عشوائي دقيق لنسبة اللوغاريتم الاحتمالي (الاقتراح 3.2)، مع تحديد الحد الرئيسي والحدود المتبقيةمصفوفة التباين المقاربة الصريحة : توفير مصفوفة تباين مقاربة صريحة Σ θ 0 \Sigma^{\theta_0} Σ θ 0 ، والتي تتضمن حد المشتق الدالي ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 ، وهو فريد لنماذج McKean-Vlasovمعدلات التقدير المختلفة : إثبات أن معاملات الانجراف يتم تقديرها بمعدل N \sqrt{N} N ، ومعاملات الانتشار بمعدل N / Δ n \sqrt{N/\Delta_n} N / Δ n ، وهو متسق مع طرق التقدير المقارنة في 1 الابتكار التقني : التعامل مع كثافات الانتقال الضمنية من خلال الحدود من النوع الغاوسي (الاقتراح 4.2) وتقديرات اللحظات (اللمة 4.1)، مما يتغلب على العقبات التقنية الرئيسيةبيانات الملاحظة :
{ X t k i , θ } i = 1 , … , N k = 1 , … , n \{X^{i,\theta}_{t_k}\}_{i=1,\ldots,N}^{k=1,\ldots,n} { X t k i , θ } i = 1 , … , N k = 1 , … , n
حيث t k = T k / n t_k = Tk/n t k = T k / n ، Δ n = T / n \Delta_n = T/n Δ n = T / n هي خطوة التقسيم
اضطراب المعاملات :
θ + = ( θ 1 + , θ 2 + ) = ( θ 1 0 + u N , θ 2 0 + v N / Δ n ) \theta^+ = (\theta_1^+, \theta_2^+) = \left(\theta_1^0 + \frac{u}{\sqrt{N}}, \theta_2^0 + \frac{v}{\sqrt{N/\Delta_n}}\right) θ + = ( θ 1 + , θ 2 + ) = ( θ 1 0 + N u , θ 2 0 + N / Δ n v )
الهدف : إثبات أن نسبة اللوغاريتم الاحتمالي
z ( θ 0 , θ + ) : = log d P θ + d P θ 0 z(\theta_0, \theta^+) := \log \frac{dP_{\theta^+}}{dP_{\theta_0}} z ( θ 0 , θ + ) := log d P θ 0 d P θ +
تحقق خاصية LAN، أي
z ( θ 0 , θ + ) → l a w ( u v ) ⊤ N θ 0 − 1 2 ( u v ) ⊤ Σ θ 0 ( u v ) z(\theta_0, \theta^+) \xrightarrow{law} \begin{pmatrix} u \\ v \end{pmatrix}^\top N_{\theta_0} - \frac{1}{2}\begin{pmatrix} u \\ v \end{pmatrix}^\top \Sigma_{\theta_0} \begin{pmatrix} u \\ v \end{pmatrix} z ( θ 0 , θ + ) l a w ( u v ) ⊤ N θ 0 − 2 1 ( u v ) ⊤ Σ θ 0 ( u v )
حيث N θ 0 ∼ N ( 0 , Σ θ 0 ) N_{\theta_0} \sim N(0, \Sigma_{\theta_0}) N θ 0 ∼ N ( 0 , Σ θ 0 ) .
يفترض النموذج N N N جزيء مستقل وموزع بشكل متطابق، كل منها يحقق:
d X t i , θ = b θ 1 ( X t i , θ , μ t θ ) d t + a θ 2 ( X t i , θ ) d W t i dX^{i,\theta}_t = b_{\theta_1}(X^{i,\theta}_t, \mu^\theta_t)dt + a_{\theta_2}(X^{i,\theta}_t)dW^i_t d X t i , θ = b θ 1 ( X t i , θ , μ t θ ) d t + a θ 2 ( X t i , θ ) d W t i
الخصائص الرئيسية:
الاعتماد على التوزيع : يعتمد حد الانجراف على التوزيع الهامشي μ t θ = Law ( X t i , θ ) \mu^\theta_t = \text{Law}(X^{i,\theta}_t) μ t θ = Law ( X t i , θ ) فصل المعاملات : يظهر معامل الانجراف θ 1 \theta_1 θ 1 ومعامل الانتشار θ 2 \theta_2 θ 2 في معاملات مختلفةالاستقلالية : حركات براونية مختلفة الجزيئات ( W i ) 1 ≤ i ≤ N (W^i)_{1\leq i \leq N} ( W i ) 1 ≤ i ≤ N مستقلة عن بعضها البعضباستخدام خاصية Markov:
z ( θ 0 , θ + ) = ∑ k = 1 n ∑ i = 1 N log p θ + p θ 0 ( t k , t k + 1 , X t k i , X t k + 1 i ) z(\theta_0, \theta^+) = \sum_{k=1}^n \sum_{i=1}^N \log \frac{p^{\theta^+}}{p^{\theta_0}}(t_k, t_{k+1}, X^i_{t_k}, X^i_{t_{k+1}}) z ( θ 0 , θ + ) = ∑ k = 1 n ∑ i = 1 N log p θ 0 p θ + ( t k , t k + 1 , X t k i , X t k + 1 i )
يتم تحليله بشكل أكبر إلى أجزاء الانجراف والانتشار:
z ( θ 0 , θ + ) = ∑ k = 1 n ∑ i = 1 N ( ζ k i , θ 1 + ζ k i , θ 2 ) z(\theta_0, \theta^+) = \sum_{k=1}^n \sum_{i=1}^N (\zeta^{i,\theta_1}_k + \zeta^{i,\theta_2}_k) z ( θ 0 , θ + ) = ∑ k = 1 n ∑ i = 1 N ( ζ k i , θ 1 + ζ k i , θ 2 )
حيث
ζ k i , θ 1 = u N ∫ 0 1 ∂ θ 1 p θ 1 ( l ) , θ 2 + p θ 1 ( l ) , θ 2 + ( t k , t k + 1 , X t k i , X t k + 1 i ) d l \zeta^{i,\theta_1}_k = \frac{u}{\sqrt{N}} \int_0^1 \frac{\partial_{\theta_1} p^{\theta_1(l), \theta_2^+}}{p^{\theta_1(l), \theta_2^+}}(t_k, t_{k+1}, X^i_{t_k}, X^i_{t_{k+1}})dl ζ k i , θ 1 = N u ∫ 0 1 p θ 1 ( l ) , θ 2 + ∂ θ 1 p θ 1 ( l ) , θ 2 + ( t k , t k + 1 , X t k i , X t k + 1 i ) d l
ζ k i , θ 2 = v N / Δ n ∫ 0 1 ∂ θ 2 p θ 1 0 , θ 2 ( l ) p θ 1 0 , θ 2 ( l ) ( t k , t k + 1 , X t k i , X t k + 1 i ) d l \zeta^{i,\theta_2}_k = \frac{v}{\sqrt{N/\Delta_n}} \int_0^1 \frac{\partial_{\theta_2} p^{\theta_1^0, \theta_2(l)}}{p^{\theta_1^0, \theta_2(l)}}(t_k, t_{k+1}, X^i_{t_k}, X^i_{t_{k+1}})dl ζ k i , θ 2 = N / Δ n v ∫ 0 1 p θ 1 0 , θ 2 ( l ) ∂ θ 2 p θ 1 0 , θ 2 ( l ) ( t k , t k + 1 , X t k i , X t k + 1 i ) d l
التقنية الرئيسية : بالنسبة لمشتقات معاملات كثافة الانتقال، استخدام حساب Malliavin للحصول على:
∂ θ 1 p θ p θ ( t k , t k + 1 , x , y ) = 1 Δ n E t k , x θ [ ∑ r = 1 d δ ( ∂ θ 1 X r , Δ n i , θ U r i ) ∣ X t k + 1 i , θ = y ] \frac{\partial_{\theta_1} p^\theta}{p^\theta}(t_k, t_{k+1}, x, y) = \frac{1}{\Delta_n} E^\theta_{t_k,x}\left[\sum_{r=1}^d \delta(\partial_{\theta_1} X^{i,\theta}_{r,\Delta_n} U^i_r) \Big| X^{i,\theta}_{t_{k+1}} = y\right] p θ ∂ θ 1 p θ ( t k , t k + 1 , x , y ) = Δ n 1 E t k , x θ [ ∑ r = 1 d δ ( ∂ θ 1 X r , Δ n i , θ U r i ) X t k + 1 i , θ = y ]
حيث:
δ \delta δ هو التكامل Skorohod (المؤثر الثنائي لمشتق Malliavin)U s i = a θ 2 − 1 ( X t k + s i , θ ) Y s i , θ ( Y Δ n i , θ ) − 1 U^i_s = a^{-1}_{\theta_2}(X^{i,\theta}_{t_k+s}) Y^{i,\theta}_s (Y^{i,\theta}_{\Delta_n})^{-1} U s i = a θ 2 − 1 ( X t k + s i , θ ) Y s i , θ ( Y Δ n i , θ ) − 1 Y t i , θ Y^{i,\theta}_t Y t i , θ هي مصفوفة عملية تحقق معادلة SDE خطية (3.14)عملية مشتق المعامل : ∂ θ 1 X t i , θ \partial_{\theta_1} X^{i,\theta}_t ∂ θ 1 X t i , θ تحقق SDE:
∂ θ 1 X t i , θ = ∫ 0 t ( ∂ θ 1 b θ 1 + ∇ x b θ 1 ∂ θ 1 X s i , θ + ∫ R d ∂ μ b θ 1 ( X t k + s i , θ , y , μ t k + s θ ) ∂ θ 1 μ t k + s θ ( d y ) ) d s + … \partial_{\theta_1} X^{i,\theta}_t = \int_0^t \left(\partial_{\theta_1} b_{\theta_1} + \nabla_x b_{\theta_1} \partial_{\theta_1} X^{i,\theta}_s + \int_{\mathbb{R}^d} \partial_\mu b_{\theta_1}(X^{i,\theta}_{t_k+s}, y, \mu^\theta_{t_k+s}) \partial_{\theta_1}\mu^\theta_{t_k+s}(dy)\right)ds + \ldots ∂ θ 1 X t i , θ = ∫ 0 t ( ∂ θ 1 b θ 1 + ∇ x b θ 1 ∂ θ 1 X s i , θ + ∫ R d ∂ μ b θ 1 ( X t k + s i , θ , y , μ t k + s θ ) ∂ θ 1 μ t k + s θ ( d y ) ) d s + …
لاحظ أن الحد الثالث يتضمن المشتق الدالي ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 ، وهو فريد لنماذج McKean-Vlasov.
جزء الانجراف : إثبات
δ ( ∂ θ 1 X r , Δ n i , θ U r i ) = Δ n z r , θ 1 θ ( X t k i , θ ) [ a θ 2 − 2 ( X t k i , θ ) ( X t k + 1 i , θ − m t k , t k + 1 θ ( X t k i , θ ) ) ] r + H t k + 1 i \delta(\partial_{\theta_1} X^{i,\theta}_{r,\Delta_n} U^i_r) = \Delta_n z^{\theta}_{r,\theta_1}(X^{i,\theta}_{t_k}) [a^{-2}_{\theta_2}(X^{i,\theta}_{t_k})(X^{i,\theta}_{t_{k+1}} - m^\theta_{t_k,t_{k+1}}(X^{i,\theta}_{t_k}))]_r + H^i_{t_{k+1}} δ ( ∂ θ 1 X r , Δ n i , θ U r i ) = Δ n z r , θ 1 θ ( X t k i , θ ) [ a θ 2 − 2 ( X t k i , θ ) ( X t k + 1 i , θ − m t k , t k + 1 θ ( X t k i , θ )) ] r + H t k + 1 i
حيث H t k + 1 i H^i_{t_{k+1}} H t k + 1 i هي الحد المتبقي، والذي يحقق ( E t k , x θ ∣ H t k + 1 i ∣ τ ) 1 / τ = R t k i ( Δ n 2 ) (E^\theta_{t_k,x}|H^i_{t_{k+1}}|^\tau)^{1/\tau} = R^i_{t_k}(\Delta_n^2) ( E t k , x θ ∣ H t k + 1 i ∣ τ ) 1/ τ = R t k i ( Δ n 2 ) .
الكمية الرئيسية :
z t θ ( x ) : = ∂ θ 1 b θ 1 ( x , μ t θ ) + ∫ R d ∂ μ b θ 1 ( x , y , μ t θ ) ∂ θ 1 μ t θ ( d y ) z^\theta_t(x) := \partial_{\theta_1} b_{\theta_1}(x, \mu^\theta_t) + \int_{\mathbb{R}^d} \partial_\mu b_{\theta_1}(x, y, \mu^\theta_t) \partial_{\theta_1}\mu^\theta_t(dy) z t θ ( x ) := ∂ θ 1 b θ 1 ( x , μ t θ ) + ∫ R d ∂ μ b θ 1 ( x , y , μ t θ ) ∂ θ 1 μ t θ ( d y )
هذه الكمية تلعب دوراً أساسياً في مصفوفة التباين المقاربة.
المسار التقني :
استخدام صيغة التكامل بالأجزاء (2.7): δ ( F u ) = F δ ( u ) − ⟨ D F , u ⟩ H \delta(Fu) = F\delta(u) - \langle DF, u\rangle_H δ ( F u ) = F δ ( u ) − ⟨ D F , u ⟩ H تقريب U r i U^i_r U r i بـ U ^ r i = a θ 2 − 1 ( X t k + r i , θ ) \hat{U}^i_r = a^{-1}_{\theta_2}(X^{i,\theta}_{t_k+r}) U ^ r i = a θ 2 − 1 ( X t k + r i , θ ) إثبات أن كل الحدود المتبقية H n i , j , j = 1 , 2 , 3 H^{i,j}_n, j=1,2,3 H n i , j , j = 1 , 2 , 3 هي من رتبة Δ n 2 \Delta_n^2 Δ n 2 جزء الانتشار : إثبات بشكل مشابه
δ ( ∂ θ 2 X r , Δ n i , θ U r i ) = [ ∂ θ 2 a θ 2 ( X t k i ) a θ 2 − 1 ( X t k i ) ( X t k + 1 i − m t k , t k + 1 θ ( X t k i ) ) ] r × [ … ] + حد متبقي \delta(\partial_{\theta_2} X^{i,\theta}_{r,\Delta_n} U^i_r) = [\partial_{\theta_2} a_{\theta_2}(X^i_{t_k}) a^{-1}_{\theta_2}(X^i_{t_k})(X^i_{t_{k+1}} - m^\theta_{t_k,t_{k+1}}(X^i_{t_k}))]_r \times [\ldots] + \text{حد متبقي} δ ( ∂ θ 2 X r , Δ n i , θ U r i ) = [ ∂ θ 2 a θ 2 ( X t k i ) a θ 2 − 1 ( X t k i ) ( X t k + 1 i − m t k , t k + 1 θ ( X t k i )) ] r × [ … ] + حد متبقي
الحد المتبقي يحقق رتبة Δ n 3 / 2 \Delta_n^{3/2} Δ n 3/2 .
طريقة Gobet 19,20 : الطريقة الأصلية موجهة نحو عمليات الانتشار المتوازنة، وتعتمد على الخصائص المتوازنة للملاحظة طويلة الأجلطريقة هذه الورقة : لا تتطلب افتراضات التوازن، والتقارب المقارب يقوده عدد الجزيئات N → ∞ N \to \infty N → ∞ المشتق الدالي : تظهر حدود ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 في نماذج McKean-Vlasov، وهي غير موجودة في SDE الكلاسيكيةحدود من نوع Aronson :
1 L Δ n d / 2 exp ( − c ∥ x − y ∥ 2 Δ n ) exp ( − c Δ n ∥ x ∥ 2 ) ≤ p θ ( t k , t k + 1 , x , y ) \frac{1}{L\Delta_n^{d/2}} \exp\left(-c\frac{\|x-y\|^2}{\Delta_n}\right) \exp(-c\Delta_n\|x\|^2) \leq p^\theta(t_k, t_{k+1}, x, y) L Δ n d /2 1 exp ( − c Δ n ∥ x − y ∥ 2 ) exp ( − c Δ n ∥ x ∥ 2 ) ≤ p θ ( t k , t k + 1 , x , y ) ≤ L Δ n d / 2 exp ( − ∥ x − y ∥ 2 c Δ n ) exp ( c Δ n ∥ x ∥ 2 ) \leq \frac{L}{\Delta_n^{d/2}} \exp\left(-\frac{\|x-y\|^2}{c\Delta_n}\right) \exp(c\Delta_n\|x\|^2) ≤ Δ n d /2 L exp ( − c Δ n ∥ x − y ∥ 2 ) exp ( c Δ n ∥ x ∥ 2 )
حدود مشتقات المعاملات :
E t k , x θ ˉ [ ∣ ∂ θ 1 p θ p θ ( t k , t k + 1 , x , X t k + 1 i ) ∣ m ] ≤ L Δ n m / 2 exp ( c Δ n ∥ x ∥ 2 ) ( 1 + ∥ x ∥ ) q E^{\bar{\theta}}_{t_k,x}\left[\left|\frac{\partial_{\theta_1} p^\theta}{p^\theta}(t_k, t_{k+1}, x, X^i_{t_{k+1}})\right|^m\right] \leq \frac{L}{\Delta_n^{m/2}} \exp(c\Delta_n\|x\|^2)(1+\|x\|)^q E t k , x θ ˉ [ p θ ∂ θ 1 p θ ( t k , t k + 1 , x , X t k + 1 i ) m ] ≤ Δ n m /2 L exp ( c Δ n ∥ x ∥ 2 ) ( 1 + ∥ x ∥ ) q
هذه الحدود حاسمة لإثبات قابلية تجاهل الحدود المتبقية (الاقتراح 4.4).
ملاحظة : هذه ورقة نظرية بحتة ولا تتضمن تجارب عددية. النتائج الرئيسية هي إثبات النظريات النظرية.
على الرغم من عدم وجود تجارب عددية، تتحقق الورقة من معقولية النظرية بالطرق التالية:
الاتساق مع النتائج الموجودة : مصفوفة التباين المقاربة Σ θ 0 \Sigma^{\theta_0} Σ θ 0 تحت الشرط N Δ n → 0 N\Delta_n \to 0 N Δ n → 0 متسقة مع طريقة التقدير المقارنة في 1 معدلات التقدير :معامل الانجراف: معدل N \sqrt{N} N معامل الانتشار: معدل N / Δ n \sqrt{N/\Delta_n} N / Δ n هذا متسق مع نظرية SDE الكلاسيكية والأدب الحديث 1 الحالات الخاصة : عندما ∂ μ b θ 1 = 0 \partial_\mu b_{\theta_1} = 0 ∂ μ b θ 1 = 0 (بدون اعتماد على التوزيع)، تتحلل النتائج إلى نتائج LAN لعمليات الانتشار الكلاسيكيةخاصية LAN : تحت الافتراضات A1-A5،
z ( θ 0 , θ + ) → P θ 0 − l a w ( u v ) ⊤ N θ 0 − 1 2 ( u v ) ⊤ Σ θ 0 ( u v ) z(\theta_0, \theta^+) \xrightarrow{P^{\theta_0}-law} \begin{pmatrix} u \\ v \end{pmatrix}^\top N_{\theta_0} - \frac{1}{2}\begin{pmatrix} u \\ v \end{pmatrix}^\top \Sigma_{\theta_0} \begin{pmatrix} u \\ v \end{pmatrix} z ( θ 0 , θ + ) P θ 0 − l a w ( u v ) ⊤ N θ 0 − 2 1 ( u v ) ⊤ Σ θ 0 ( u v )
مصفوفة التباين المقاربة :
Σ θ 0 = ( Σ b θ 0 0 0 Σ a θ 0 ) \Sigma_{\theta_0} = \begin{pmatrix} \Sigma^{\theta_0}_b & 0 \\ 0 & \Sigma^{\theta_0}_a \end{pmatrix} Σ θ 0 = ( Σ b θ 0 0 0 Σ a θ 0 )
حيث
Σ b θ 0 = ∫ 0 T ∫ R d z s θ 0 ( x ) ⊤ a θ 2 0 − 2 ( x ) z s θ 0 ( x ) μ s θ 0 ( d x ) d s \Sigma^{\theta_0}_b = \int_0^T \int_{\mathbb{R}^d} z^{\theta_0}_s(x)^\top a^{-2}_{\theta_2^0}(x) z^{\theta_0}_s(x) \mu^{\theta_0}_s(dx)ds Σ b θ 0 = ∫ 0 T ∫ R d z s θ 0 ( x ) ⊤ a θ 2 0 − 2 ( x ) z s θ 0 ( x ) μ s θ 0 ( d x ) d s
Σ a θ 0 = 2 ∫ 0 T ∫ R d tr ( ∂ θ 2 a θ 2 0 ( x ) a θ 2 0 − 1 ( x ) ∂ θ 2 a θ 2 0 ( x ) a θ 2 0 − 1 ( x ) ) μ s θ 0 ( d x ) d s \Sigma^{\theta_0}_a = 2\int_0^T \int_{\mathbb{R}^d} \text{tr}(\partial_{\theta_2} a_{\theta_2^0}(x) a^{-1}_{\theta_2^0}(x) \partial_{\theta_2} a_{\theta_2^0}(x) a^{-1}_{\theta_2^0}(x)) \mu^{\theta_0}_s(dx)ds Σ a θ 0 = 2 ∫ 0 T ∫ R d tr ( ∂ θ 2 a θ 2 0 ( x ) a θ 2 0 − 1 ( x ) ∂ θ 2 a θ 2 0 ( x ) a θ 2 0 − 1 ( x )) μ s θ 0 ( d x ) d s
البنية القطرية : Σ θ 0 \Sigma_{\theta_0} Σ θ 0 هي مصفوفة قطرية، مما يشير إلى الاستقلالية المقاربة لمعاملات الانجراف والانتشاردور المشتق الدالي : z s θ 0 ( x ) z^{\theta_0}_s(x) z s θ 0 ( x ) يتضمن حد ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 ، وهو فريد لنماذج McKean-Vlasov، مما يعكس تأثير الاعتماد على التوزيعالفرق عن الأنظمة الجزيئية التفاعلية :نموذج McKean-Vlasov (1.1): التباين يتضمن ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 نظام الجزيئات التفاعلية (3.18): التباين في z t θ ( x ) z^\theta_t(x) z t θ ( x ) يبسط إلى ∂ θ 1 b θ 1 ( x , μ t θ ) \partial_{\theta_1} b_{\theta_1}(x, \mu^\theta_t) ∂ θ 1 b θ 1 ( x , μ t θ ) يتم إثبات النظرية 3.4 من خلال التحقق من ستة شروط تقارب (4.30)-(4.36):
الشرط (4.30) : اللحظة الأولى لجزء الانجراف
∑ k = 1 n ∑ i = 1 N E t k θ 0 [ ζ ^ k i , θ 1 ] → P θ 0 − 1 2 u 2 Σ b θ 0 \sum_{k=1}^n \sum_{i=1}^N E^{\theta_0}_{t_k}[\hat{\zeta}^{i,\theta_1}_k] \xrightarrow{P^{\theta_0}} -\frac{1}{2}u^2 \Sigma^{\theta_0}_b ∑ k = 1 n ∑ i = 1 N E t k θ 0 [ ζ ^ k i , θ 1 ] P θ 0 − 2 1 u 2 Σ b θ 0
الخطوة الرئيسية: استخدام توسع Taylor
m t k , t k + 1 θ 0 ( X t k i ) − m t k , t k + 1 θ 1 ( l ) , θ 2 + ( X t k i ) = − l u Δ n N z t k θ 0 ( X t k i ) + R t k i ( ε n , N Δ n / N ) m^{\theta_0}_{t_k,t_{k+1}}(X^i_{t_k}) - m^{\theta_1(l),\theta_2^+}_{t_k,t_{k+1}}(X^i_{t_k}) = -\frac{lu\Delta_n}{\sqrt{N}} z^{\theta_0}_{t_k}(X^i_{t_k}) + R^i_{t_k}(\varepsilon_{n,N}\Delta_n/\sqrt{N}) m t k , t k + 1 θ 0 ( X t k i ) − m t k , t k + 1 θ 1 ( l ) , θ 2 + ( X t k i ) = − N l u Δ n z t k θ 0 ( X t k i ) + R t k i ( ε n , N Δ n / N )
الشرط (4.31) : تقارب اللحظة الثانية لجزء الانجراف إلى u 2 Σ b θ 0 u^2\Sigma^{\theta_0}_b u 2 Σ b θ 0
الحيلة:
الحد المتقاطع ∑ i 1 ≠ i 2 ζ ^ k i 1 , θ 1 ζ ^ k i 2 , θ 1 \sum_{i_1 \neq i_2} \hat{\zeta}^{i_1,\theta_1}_k \hat{\zeta}^{i_2,\theta_1}_k ∑ i 1 = i 2 ζ ^ k i 1 , θ 1 ζ ^ k i 2 , θ 1 قابل للتجاهل بشكل مقارب الحد الرئيسي يأتي من ∑ i = 1 N ( ζ ^ k i , θ 1 ) 2 \sum_{i=1}^N (\hat{\zeta}^{i,\theta_1}_k)^2 ∑ i = 1 N ( ζ ^ k i , θ 1 ) 2 الشرط (4.32) : شرط اللحظة الرابعة
∑ k = 1 n E t k θ 0 [ ∣ ∑ i = 1 N ζ ^ k i , θ 1 ∣ 4 ] → P θ 0 0 \sum_{k=1}^n E^{\theta_0}_{t_k}\left[\left|\sum_{i=1}^N \hat{\zeta}^{i,\theta_1}_k\right|^4\right] \xrightarrow{P^{\theta_0}} 0 ∑ k = 1 n E t k θ 0 [ ∑ i = 1 N ζ ^ k i , θ 1 4 ] P θ 0 0
من خلال تحليل الحد الرابع إلى مجموعات فهرس مختلفة (جميع مختلفة، زوجان، جميع متطابقة) وإثبات أن كل جزء يميل إلى 0.
الشروط (4.33)-(4.35) : شروط مشابهة لجزء الانتشار، باستخدام توسع التباين الشرطي
V t k , t k + 1 θ 0 ( x ) − V t k , t k + 1 θ 1 0 , θ 2 ( l ) ( x ) = − 2 l v Δ n 3 / 2 N ∂ θ 2 a θ 2 0 ( X t k i ) a θ 2 0 ( X t k i ) + … V^{\theta_0}_{t_k,t_{k+1}}(x) - V^{\theta_1^0,\theta_2(l)}_{t_k,t_{k+1}}(x) = -\frac{2lv\Delta_n^{3/2}}{\sqrt{N}} \partial_{\theta_2} a_{\theta_2^0}(X^i_{t_k}) a_{\theta_2^0}(X^i_{t_k}) + \ldots V t k , t k + 1 θ 0 ( x ) − V t k , t k + 1 θ 1 0 , θ 2 ( l ) ( x ) = − N 2 l v Δ n 3/2 ∂ θ 2 a θ 2 0 ( X t k i ) a θ 2 0 ( X t k i ) + …
الشرط (4.36) : الاستقلالية المقاربة للانجراف والانتشار، من خلال إثبات أن الحد المتقاطع يميل بشكل مقارب إلى 0
الملاحظة المنفصلة :
1 Amorino وآخرون (2023): تقدير مقارن للأنظمة الجزيئية التفاعلية، كائن المقارنة المباشر لهذه الورقة6 Bishwal (2011): تقدير الانتشار التفاعلي9 Chen (2021): تقدير الاحتمالية القصوى للجهد من بيانات المسار الفردي16,17 Genon-Catalot & Larédo (2021): نماذج McKean-Vlasov ذات التباين الصغير والوقت الطويل27 Liu & Qiao (2022): معادلات SDE لـ McKean-Vlasov التي تعتمد على المسار31 Sharrock وآخرون (2021): تقدير المعاملات عبر الإنترنتالملاحظة المستمرة :
13 Della Maestra & Hoffmann (2023): خاصية LAN لنماذج McKean-Vlasov تحت آلية المجال المتوسط (ذات صلة مباشرة بهذه الورقة)
الفرق: الملاحظة المستمرة يمكن استخدام نظرية Girsanov، ودالة الاحتمالية لها تعبير مغلق 2 Amorino وآخرون (2024): معدل متعدد الحدود من خلال إلغاء الالتفاف4 Belomestny وآخرون (2022): تقدير شبه بارامتري لمعادلات SDE لـ McKean-Vlasov11 Comte وآخرون (2024): طريقة لحظات غير بارامترية12 Della Maestra & Hoffmann (2022): تقدير غير بارامتري للأنظمة الجزيئية التفاعلية29 Nickl وآخرون (2025): استدلال بايزي غير بارامتري لنماذج McKean-Vlasov19,20 Gobet (2001,2002):
الحالة الطبيعية المقاربة المختلطة المحلية للانتشار الإهليلجي خاصية LAN للانتشار المتوازن من الملاحظة المنفصلة أساس طريقة هذه الورقة : استخدام حساب Malliavin لاشتقاق تمثيل مشتقات كثافة الانتقالمقابل 13 (الملاحظة المستمرة) :التعامل مع التحديات التقنية للملاحظة المنفصلة عدم الاعتماد على نظرية Girsanov مقابل 1 (التقدير المقارن) :توفير أساس نظري لطريقة الاحتمالية إنشاء خاصية LAN، والتي يمكن استخدامها لاشتقاق الأمثلية المقاربة للمقدرات مقابل 20 (SDE الكلاسيكية) :التوسع إلى إعداد McKean-Vlasov عدم الحاجة إلى افتراضات التوازن التعامل مع المشتق الدالي ∂ μ b \partial_\mu b ∂ μ b مقابل الأنظمة الجزيئية التفاعلية :تجنب حدود كثافة الانتقال المشتركة عالية الأبعاد (تشير الملاحظة 3.3 إلى أن هذا هو العقبة الرئيسية لـ LAN للأنظمة الجزيئية التفاعلية) استخدام البنية i.i.d. لتبسيط التحليل إنشاء خاصية LAN : للمرة الأولى، يتم إنشاء خاصية LAN لمعادلات McKean-Vlasov المرصودة بشكل منفصل، مما يملأ الفراغ النظري في هذا المجالالشكل الصريح لمصفوفة التباين المقاربة :
Σ θ 0 = diag ( Σ b θ 0 , Σ a θ 0 ) \Sigma_{\theta_0} = \text{diag}(\Sigma^{\theta_0}_b, \Sigma^{\theta_0}_a) Σ θ 0 = diag ( Σ b θ 0 , Σ a θ 0 )
حيث يتضمن جزء الانجراف المشتق الدالي ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 ، مما يعكس الاعتماد على التوزيعتأكيد معدلات التقدير :الانجراف: N \sqrt{N} N الانتشار: N / Δ n \sqrt{N/\Delta_n} N / Δ n متسق مع طرق التقدير المقارنة الحديثة 1 المساهمة التقنية : تطوير تقنيات حساب Malliavin للتعامل مع كثافات الانتقال الضمنية، مع دمج الحدود من النوع الغاوسي والتكامل بالأجزاءشروط الافتراضات قوية :A2 : المعاملات محدودة و Lipschitz مستمرةA3 : السلاسة من الدرجة العالية للمعاملات (C 2 C^2 C 2 والمشتقات نمو متعدد الحدود)A5 : الإهليلجية المنتظمة لمصفوفة الانتشار قد لا تكون هذه الشروط مرضية في التطبيقات العمليةقيد المعامل أحادي البعد : على الرغم من أن الورقة تشير إلى إمكانية التوسع إلى حالة متعددة المعاملات، إلا أنها تتعامل بالتفصيل فقط مع حالة θ 1 , θ 2 ∈ R \theta_1, \theta_2 \in \mathbb{R} θ 1 , θ 2 ∈ R الفراغ في الأنظمة الجزيئية التفاعلية :تشير الملاحظة 3.3 إلى أنه بالنسبة للنظام الجزيئي التفاعلي (3.18)، بسبب نقص الحدود لكثافة الانتقال المشتركة ذات البعد dN، لا يمكن إنشاء خاصية LAN هذه مشكلة مفتوحة مهمة آلية التقارب : يتطلب Δ n → 0 , N → ∞ \Delta_n \to 0, N \to \infty Δ n → 0 , N → ∞ بشكل متزامن، مع متطلبات معينة على السرعة النسبية لـ N Δ n N\Delta_n N Δ n التوزيع الأولي : يتطلب الافتراض A1 أن يكون التوزيع الأولي μ 0 \mu_0 μ 0 من النوع sub-Gaussian، مما يحد من نطاق التطبيقLAN للأنظمة الجزيئية التفاعلية : تطوير حدود كثافة الانتقال عالية الأبعاد، وإنشاء خاصية LAN للنموذج (3.18)تخفيف شروط الافتراضات :دراسة حالة الانتشار غير الإهليلجي السماح بمعاملات غير محدودة أو Lipschitz محلية فقط توسع متعدد المعاملات : معالجة كاملة لحالة θ 1 ∈ R p , θ 2 ∈ R q \theta_1 \in \mathbb{R}^p, \theta_2 \in \mathbb{R}^q θ 1 ∈ R p , θ 2 ∈ R q بناء المقدرات المثلى : استخدام خاصية LAN لبناء مقدرات مقاربة فعالةاختبار الفرضيات : تطوير نظرية اختبار الفرضيات لنماذج McKean-Vlasov بناءً على خاصية LANالحالات غير المتوازنة : التوسع إلى عمليات McKean-Vlasov غير المتوازنةالبيانات عالية التردد : دراسة النظرية المقاربة عندما يكون Δ n → 0 \Delta_n \to 0 Δ n → 0 أسرعالصرامة النظرية :إثبات كامل وتفصيلي (القسم 4 يشغل نصف الورقة) كل خطوة تقنية لها حجة واضحة استخدام مناسب لأدوات التحليل العشوائي الحديثة (حساب Malliavin) الابتكار في الطريقة :تطبيق ذكي لحساب Malliavin : تقنية توسيع تكامل Skorohod إلى حد رئيسي زائد حد متبقي (الاقتراح 3.2) هي الابتكار الأساسيالتعامل مع المشتق الدالي : التعرف الصحيح والتعامل مع حد ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 ، وهو فريد لنماذج McKean-Vlasovالتحكم في الحدود المتبقية : معالجة موحدة لقابلية تجاهل أنواع مختلفة من الحدود المتبقية من خلال الاقتراح 4.4المساهمة النظرية :ملء الفراغ في نظرية LAN لمعادلات McKean-Vlasov المرصودة بشكل منفصل ربط طريقة الاحتمالية وطريقة التقدير المقارنة (الاتصال مع 1 ) توفير أساس نظري للاستدلال الإحصائي المقارب لنماذج المجال المتوسط وضوح الكتابة :بنية واضحة: الافتراضات → النتائج الرئيسية → الإثبات نظام تدوين شامل (القسم 2.1) شرح كافٍ للصعوبات الرئيسية وطرق الحل (المقدمة والملاحظة 3.3) مراجعة الأدب الشاملة : تحديد دقيق لموقع هذه الورقة في أدب الاستدلال الإحصائي لنماذج McKean-Vlasovالفائدة العملية محدودة :شروط الافتراضات قوية، قد لا تكون مرضية في البيانات الفعلية لا توجد محاكاة عددية للتحقق من النتائج النظرية لا توجد مناقشة حول كيفية التحقق من الافتراضات في الممارسة العملية قابلية قراءة التفاصيل التقنية :إثبات القسم 4 معقد جداً من الناحية التقنية، غير ودود للمتخصصين غير المتخصصين بعض عدم المساواة الرئيسية (مثل إثبات الاقتراح 4.2) تشير إلى 20 دون شرح كامل لكيفية التكيف مع إعداد McKean-Vlasov حدود النتائج :خاصية LAN للنظام الجزيئي التفاعلي (3.18) لا تزال مفتوحة (الملاحظة 3.3) يتم التعامل فقط مع الفترة الزمنية الثابتة T T T ، بدون مناقشة حالة T → ∞ T \to \infty T → ∞ معالجة الحالة متعددة المعاملات :على الرغم من الادعاء بإمكانية التوسع إلى حالة متعددة المعاملات، إلا أن الإطار فقط يتم تقديمه التفاصيل التقنية للحالة متعددة المعاملات (خاصة العناصر غير القطرية لـ Σ θ 0 \Sigma_{\theta_0} Σ θ 0 ) لم يتم توسيعها بالكامل الفصل عن التطبيقات :لا توجد أمثلة تطبيقية محددة لا توجد مناقشة حول كيفية استخدام النتائج في التمويل والعلوم العصبية وغيرها المساهمة في المجال :الأساس النظري : توفير أساس نظري صارم للاستدلال الإحصائي لنماذج McKean-Vlasovالمنهجية : تطبيق منهجي لحساب Malliavin في الإحصاء لنماذج McKean-Vlasovالمشاكل المفتوحة : تحديد واضح للعقبات التقنية في LAN للأنظمة الجزيئية التفاعلية (كثافة الانتقال عالية الأبعاد)، مما يوجه البحث المستقبليالقيمة العملية :تقييم المقدرات : يمكن استخدامها لتقييم الكفاءة المقاربة للمقدرات الموجودة (مثل التقدير المقارن في 1 )الحد الأدنى : توفير خاصية LAN حد Cramér-Rao النوع الأدنى لتباين المقدراتبناء المقدرات المثلى : يمكن أن توجه بناء مقدرات مقاربة فعالةقابلية إعادة الإنتاج :✅ النتائج النظرية قابلة للتحقق بالكامل (الإثبات كامل) ❌ لا توجد أكواد أو تجارب عددية ✅ الافتراضات واضحة وصريحة ⚠️ بعض التفاصيل التقنية تتطلب الرجوع إلى المراجع 19,20,30 وغيرها الحالة المتوقعة للاستشهاد :قصير الأجل : سيتم الاستشهاد به من قبل الخبراء في مجال الاستدلال الإحصائي للمجال المتوسطمتوسط الأجل : قد يصبح مرجعاً معيارياً للاستدلال الإحصائي لنماذج McKean-Vlasovطويل الأجل : إذا تم حل مشكلة الأنظمة الجزيئية التفاعلية، سيزداد التأثير بشكل أكبرالبحث النظري :النظرية الإحصائية لنماذج McKean-Vlasov تقدير المعاملات في ألعاب المجال المتوسط الإحصاء المقارب للمعادلات التفاضلية العشوائية غير الخطية مجالات التطبيق المحتملة :التمويل : نماذج المخاطر النظامية 18 ، تسعير الخيارات 21 العلوم العصبية : نماذج شبكات الخلايا العصبية 3 الفيزياء الإحصائية : حد المجال المتوسط لأنظمة الجزيئاتالديناميكا الاجتماعية : نماذج ديناميكا الآراء 8 قابلية تطبيق الطريقة :✅ عينة كبيرة (N N N كبيرة) ✅ ملاحظة عالية التردد (Δ n \Delta_n Δ n صغيرة) ✅ فترة زمنية ثابتة ✅ معاملات سلسة ❌ عينة صغيرة أو ملاحظة منخفضة التردد ❌ انتشار غير إهليلجي المقارنة مع الطرق الأخرى :مقابل التقدير المقارن 1 : LAN يوفر الأمثلية النظرية، التقدير المقارن أسهل في الحسابمقابل الطريقة البايزية 29 : LAN هي طريقة تكرارية، الطريقة البايزية أكثر مرونة لكن حسابياً أكثر كثافةمقابل الطرق غير البارامترية 12 : LAN موجهة نحو النماذج البارامترية، الطرق غير البارامترية مناسبة عندما يكون هناك عدم يقين في النموذج1 Amorino وآخرون (2023) : التقدير المقارن للأنظمة الجزيئية التفاعلية، كائن المقارنة المباشر لهذه الورقة13 Della Maestra & Hoffmann (2023) : LAN لـ McKean-Vlasov تحت الملاحظة المستمرة، السلف المباشر لهذه الورقة19,20 Gobet (2001,2002) : المصدر الأصلي لطريقة حساب Malliavin30 Nualart (1995) : المرجع المعياري لحساب Malliavin22,25,26 سلسلة Le Cam : الأدب الأساسي لنظرية LANهذه الورقة هي مساهمة نظرية مهمة في مجال الاستدلال الإحصائي لمعادلات McKean-Vlasov العشوائية التفاضلية. من خلال التطبيق الذكي لحساب Malliavin، نجح المؤلفون في إنشاء الحالة الطبيعية المقاربة المحلية في حالة الملاحظة المنفصلة، مما يملأ الفراغ النظري في هذا المجال. مستوى التقنية في الورقة عالٍ، والإثبات صارم، مما يوفر أساساً نظرياً صارماً للاستدلال الإحصائي المقارب لنماذج المجال المتوسط.
القيمة الرئيسية تكمن في: (1) الاكتمال النظري: إنشاء منهجي لنظرية LAN لنماذج McKean-Vlasov؛ (2) الابتكار في الطريقة: تطوير تقنيات للتعامل مع كثافات الانتقال الضمنية؛ (3) التوجيه النظري: توفير معيار للأمثلية المقاربة للمقدرات.
القيود الرئيسية: (1) شروط الافتراضات قوية؛ (2) نقص التحقق العددي؛ (3) LAN للأنظمة الجزيئية التفاعلية لا تزال مفتوحة.
بالنسبة للباحثين الذين يعملون على الاستدلال الإحصائي لنماذج McKean-Vlasov، هذه ورقة يجب قراءتها. بالنسبة للباحثين التطبيقيين، يجب تقييم ما إذا كانت شروط الافتراضات مرضية للمشكلة المحددة.