भिन्नात्मक SDE में ड्रिफ्ट फलन के एक गणनीय स्कोरोखोड समाकल आधारित प्रक्षेपण अनुमानक पर

मूल जानकारी

• पेपर ID: 2307.04949
• शीर्षक: भिन्नात्मक SDE में ड्रिफ्ट फलन के एक गणनीय स्कोरोखोड समाकल आधारित प्रक्षेपण अनुमानक पर
• लेखक: निकोलस मैरी (लेबोरेटोयर मोडल'एक्स, विश्वविद्यालय पेरिस नांटेरे)
• वर्गीकरण: math.ST (सांख्यिकी सिद्धांत), math.PR (संभाव्यता), stat.TH (सांख्यिकी सिद्धांत)
• प्रकाशन समय: 15 जनवरी 2024 (arXiv v2)
• पेपर लिंक: https://arxiv.org/abs/2307.04949

सारांश

यह पेपर स्कोरोखोड समाकल आधारित प्रक्षेपण अनुमानक $\hat{b}_m$ का अध्ययन करता है जो भिन्नात्मक ब्राउनियन गति द्वारा संचालित यादृच्छिक अवकल समीकरण में ड्रिफ्ट फलन $b_0$ का अनुमान लगाता है। यह अनुमानक समाधान $X$ की $N$ स्वतंत्र प्रतियों $X^1,\ldots,X^N$ के आधार पर गणना की जाती है, जहाँ $dX_t = b_0(X_t)dt + \sigma dB_t$, $B$ हर्स्ट सूचकांक $H \in (1/2,1)$ के साथ भिन्नात्मक ब्राउनियन गति है। यद्यपि स्कोरोखोड समाकल आधारित अनुमानक को प्रेक्षित डेटा से सीधे गणना नहीं की जा सकती, यह पेपर इसके गणनीय सन्निकटन के लिए $L^2$ त्रुटि सीमाएँ स्थापित करता है।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या परिभाषा

यह पेपर भिन्नात्मक यादृच्छिक अवकल समीकरण में ड्रिफ्ट फलन के अ-पैरामीट्रिक अनुमान की मूल समस्या को हल करता है। निम्नलिखित भिन्नात्मक SDE पर विचार करें: $X_t = X_0 + \int_0^t b_0(X_s)ds + \sigma B_t, \quad t \in [0,T]$

जहाँ $B$ द्विपक्षीय भिन्नात्मक ब्राउनियन गति है, $H \in (1/2,1)$ हर्स्ट सूचकांक है, $b_0 \in C^1(\mathbb{R})$ अनुमानित किया जाने वाला ड्रिफ्ट फलन है।

अनुसंधान प्रेरणा

1. व्यावहारिक अनुप्रयोग की आवश्यकता: औषध गतिविज्ञान जैसे क्षेत्रों में, अक्सर एकल नमूने के दीर्घकालीन व्यवहार के बजाय कई स्वतंत्र नमूनों (जैसे नैदानिक परीक्षणों में कई रोगियों) के आधार पर मापदंडों का अनुमान लगाने की आवश्यकता होती है।
2. गणनीय चुनौती: पारंपरिक स्कोरोखोड समाकल आधारित अनुमानक सैद्धांतिक रूप से व्यावहारिक हैं, लेकिन प्रेक्षित डेटा से सीधे गणना नहीं की जा सकती, जो इसकी व्यावहारिकता को गंभीर रूप से सीमित करता है।
3. सैद्धांतिक अंतराल: मौजूदा अनुसंधान मुख्य रूप से पैरामीट्रिक अनुमान या दीर्घकालीन व्यवहार आधारित अनुमानकों पर केंद्रित है, स्वतंत्र प्रतियों के आधार पर अ-पैरामीट्रिक अनुमान पर अनुसंधान अपेक्षाकृत कम है।

मूल योगदान

1. गणनीय अनुमानक का प्रस्ताव: निश्चित बिंदु रणनीति के माध्यम से, गैर-गणनीय स्कोरोखोड समाकल अनुमानक को गणनीय निश्चित बिंदु अनुमानक $\tilde{b}_m$ में परिवर्तित करता है।
2. सैद्धांतिक त्रुटि सीमाएँ स्थापित करना: सहायक अनुमानक $\hat{b}_m$ और निश्चित बिंदु अनुमानक $\tilde{b}_m$ के लिए $L^2$ त्रुटि सीमाएँ स्थापित करता है।
3. अभिसरण दर विश्लेषण: स्थिर स्थिति में, अभिसरण दर $N^{-(2H-1)/(10H)}$ को सिद्ध करता है, और नमूना आकार $N$, समय विंडो $T$ और प्रक्षेपण आयाम $m$ के बीच इष्टतम संतुलन संबंध प्रदान करता है।
4. अ-पैरामीट्रिक ढाँचे का विस्तार: मैरी 12 की पैरामीट्रिक अनुमान निश्चित बिंदु रणनीति को अ-पैरामीट्रिक अनुमान ढाँचे तक सफलतापूर्वक विस्तारित करता है।

विधि विवरण

कार्य परिभाषा

भिन्नात्मक SDE की $N$ स्वतंत्र समाधान पथों $X^1,\ldots,X^N$ को देखते हुए, समय अंतराल $[0,T]$ पर प्रेक्षित, लक्ष्य ड्रिफ्ट फलन $b_0$ का अनुमान लगाना है।

सैद्धांतिक आधार: स्कोरोखोड समाकल

मान लीजिए $H = \{h: \|h\|_H < \infty\}$ $B_T$ का पुनरुत्पादक कर्नल हिल्बर्ट स्पेस है, जिसका आंतरिक गुणनफल है: $\langle h,\eta \rangle_H = \alpha_H \int_0^T \int_0^T h(s)\eta(t)|t-s|^{2H-2}dsdt$ जहाँ $\alpha_H = H(2H-1)$।

स्कोरोखोड समाकल को मल्लियाविन व्युत्पन्न के संलग्न ऑपरेटर के रूप में परिभाषित किया गया है। प्रक्रिया $Y$ के लिए, इसका स्कोरोखोड समाकल है: $\int_0^t Y_s \delta X_s = \int_0^t Y_s b_0(X_s)ds + \sigma \int_0^t Y_s \delta B_s$

प्रक्षेपण अनुमानक निर्माण

अंतराल $I$ पर एक ऑर्थोनॉर्मल परिवार $(\phi_1,\ldots,\phi_m)$ पर विचार करें, आदर्श प्रक्षेपण अनुमानक है: $\hat{b}_m(x) = \frac{1}{f(x)} \sum_{j=1}^m \left(\frac{1}{NT} \sum_{i=1}^N \int_0^T \phi_j(X^i_s)\delta X^i_s\right) \phi_j(x)$

जहाँ $f(x) = \frac{1}{T}\int_0^T f_s(x)ds$ औसत घनत्व फलन है।

निश्चित बिंदु विधि

चूंकि स्कोरोखोड समाकल को सीधे गणना नहीं किया जा सकता, नुआलर्ट 15 के परिणाम का उपयोग करते हुए: $\int_0^T \phi_j(X^i_s)\delta X^i_s = I_{i,j} - a\int_0^T \int_0^t \phi'_j(X^i_t)\exp\left(\int_s^t b'_0(X^i_u)du\right)|t-s|^{2H-2}dsdt$

जहाँ $I_{i,j} = \int_0^T \phi_j(X^i_s)dX^i_s$ गणनीय यंग समाकल है, $a = \sigma^2 H(2H-1)$।

यह संबंध $\hat{b}_m = \Phi_m(b_0) \approx \Phi_m(\hat{b}_m)$ की ओर ले जाता है, जहाँ $\Phi_m$ मानचित्र है: $\Phi_m(\phi)(x) = \frac{1}{f(x)} \sum_{j=1}^m \left[\frac{1}{NT}\sum_{i=1}^N \left[I_{i,j} - a\int_0^T \int_0^t \phi'_j(X^i_t)\exp\left(\int_s^t \phi'(X^i_u)du\right)|t-s|^{2H-2}dsdt\right]\right]\phi_j(x)$

निश्चित बिंदु अनुमानक $\tilde{b}_m$ को समुच्चय $S_{m,c} = \{\phi \in S_m: \|\phi'\|_\infty \leq c\}$ में $\Phi_m$ के निश्चित बिंदु के रूप में परिभाषित किया गया है, जहाँ $c > \|b'_0\|_\infty$।

प्रायोगिक सेटअप

सैद्धांतिक विश्लेषण ढाँचा

पेपर मुख्य रूप से सैद्धांतिक विश्लेषण करता है, दो मुख्य स्थितियों पर विचार करते हुए:

1. निर्धारणीय प्रारंभिक मान स्थिति: $X_0(\cdot) = x_0 \in \mathbb{R}$
2. स्थिर स्थिति: $b_0$ विघटनकारी शर्त को संतुष्ट करता है $\exists m > 0: \forall x \in \mathbb{R}, b'_0(x) \leq -m$

आधार फलन चयन

मुख्य रूप से त्रिकोणमितीय आधार फलन स्थिति पर विचार करें: कॉम्पैक्ट अंतराल $I = [\ell,r]$ पर, $\phi_1(x) = \sqrt{\frac{1}{r-\ell}}1_I(x)$$\phi_{2j+1}(x) = \sqrt{\frac{2}{r-\ell}}\sin\left(\frac{2\pi j(x-\ell)}{r-\ell}\right)1_I(x)$$\phi_{2j}(x) = \sqrt{\frac{2}{r-\ell}}\cos\left(\frac{2\pi j(x-\ell)}{r-\ell}\right)1_I(x)$

मुख्य मान्यताएँ

• मान्यता 3.1: $X_t$ के वितरण में घनत्व $f_t$ है, और $t \mapsto f_t(x) \in L^1([0,T])$
• मान्यता 4.1: $\exists m_f > 0$ ऐसा कि $f(x) \geq m_f, \forall x \in I$
• मान्यता 4.4: $\exists$ स्थिरांक $c_\phi > 0$ ऐसा कि $I(m) \vee L(m) \leq c_\phi R(m)$

जहाँ $L(m) = \sum_{j=1}^m \|\phi_j\|_\infty^2$, $R(m) = \sum_{j=1}^m \|\phi'_j\|_\infty^2$।

प्रायोगिक परिणाम

मुख्य सैद्धांतिक परिणाम

प्रमेय 3.3 (सहायक अनुमानक त्रुटि सीमा): $E(\|\hat{b}_m - b_0\|_{f^2}^2) \leq \|b_m - b_0\|_{f^2}^2 + \frac{2}{N}\left(\|b_0\|_f^2 L(m) + c_{2.5}\sigma^2 \frac{m_T}{T^{2-2H}}(L(m) + R(m))\right)$

प्रमेय 4.7 (निश्चित बिंदु अनुमानक त्रुटि सीमा): स्थिर स्थिति में, $E(\|\tilde{b}_m^{c,l} - b_0\|_{f^2}^2) \leq 2\|b_m - b_0\|_{f^2}^2 + c_{4.7}R(m)V(N,T)$ जहाँ $V(N,T) = N^{-1/2}T^{-1} + T^{2H-1}$।

अभिसरण दर विश्लेषण

त्रिकोणमितीय आधार फलनों के लिए, जब $(b_0f)|_I \in W_2^\beta(I)$ ($\beta$ क्रम सोबोलेव स्पेस) हो: $E(\|\tilde{b}_m^{c,l} - b_0\|_{f^2}^2) \leq c_{4.7}(m^{-2\beta} + m^3 N^{-(2H-1)/(4H)})$

इष्टतम चयन:

• समय विंडो: $T(N) = N^{-1/(4H)}$
• प्रक्षेपण आयाम: $m \asymp N^{(2H-1)/(4H(3+2\beta))}$
• अभिसरण दर: $N^{-(2H-1)/(10H)}$

संभाव्यता सीमा अनुमान

लेम्मा 4.6: उपयुक्त शर्तों के तहत, $P(\Delta_m^c) \leq c_{4.6,2}R(m)(N^{-1/2}T^{-1} + T^{2H-1})$ जहाँ $\Delta_m$ वह घटना है जहाँ $\Phi_m$ एक संकुचन मानचित्र है।

संबंधित कार्य

ऐतिहासिक विकास

1. दीर्घकालीन व्यवहार अनुमानक: क्लेप्त्सिना और ले ब्रेटन 9, ट्यूडर और विएन्स 17, हू और नुआलर्ट 7 आदि ने एकल नमूने के दीर्घकालीन अवलोकन के आधार पर पैरामीट्रिक अनुमान का अध्ययन किया।
2. अ-पैरामीट्रिक अनुमान: सॉसेरॉ 16, कॉम्टे और मैरी 2 ने भिन्नात्मक SDE के अ-पैरामीट्रिक अनुमान का अध्ययन किया।
3. स्वतंत्र प्रतियों की विधि: कॉम्टे और जेनॉन-कैटलॉट 1 ने पहली बार शास्त्रीय विसरण प्रक्रियाओं में स्वतंत्र प्रतियों के आधार पर अनुमान का अध्ययन किया, मैरी 12 ने इसे भिन्नात्मक स्थिति में पैरामीट्रिक अनुमान तक विस्तारित किया।

इस पेपर का योगदान

यह पेपर पहली बार स्वतंत्र प्रतियों के आधार पर विधि को भिन्नात्मक SDE के अ-पैरामीट्रिक अनुमान तक विस्तारित करता है, और स्कोरोखोड समाकल की गणनीय समस्या को हल करता है।

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

1. भिन्नात्मक SDE ड्रिफ्ट फलन के गणनीय अ-पैरामीट्रिक अनुमानक का सफलतापूर्वक निर्माण
2. अस्तित्व, अद्वितीयता और अभिसरण दर सहित पूर्ण सैद्धांतिक विश्लेषण ढाँचा स्थापित करना
3. नमूना आकार, अवलोकन समय और प्रक्षेपण आयाम के बीच जटिल संतुलन संबंध को उजागर करना

सीमाएँ

1. गणनीय जटिलता: निश्चित बिंदु विधि को पुनरावृत्ति गणना की आवश्यकता है, व्यावहारिक कार्यान्वयन काफी जटिल है
2. घनत्व अनुमान निर्भरता: व्यावहारिक अनुप्रयोग में अज्ञात घनत्व फलन $f$ का एक साथ अनुमान लगाने की आवश्यकता है
3. मान्यता शर्तें: मजबूत नियमितता मान्यताओं और सीमित व्युत्पन्न शर्तों की आवश्यकता है

भविष्य की दिशाएँ

1. अज्ञात घनत्व फलन स्थिति में संयुक्त अनुमान समस्या का अध्ययन करना
2. अधिक कुशल संख्यात्मक एल्गोरिदम विकसित करना
3. बहु-आयामी स्थिति और अधिक सामान्य भिन्नात्मक प्रक्रियाओं तक विस्तार करना

गहन मूल्यांकन

शक्तियाँ

1. सैद्धांतिक नवीनता: पहली बार भिन्नात्मक SDE में स्कोरोखोड समाकल अनुमानक की गणनीय समस्या को हल करता है
2. विधि कठोरता: पूर्ण गणितीय प्रमाण और त्रुटि विश्लेषण
3. व्यावहारिक मूल्य: भिन्नात्मक प्रक्रिया के सांख्यिकीय अनुमान के लिए नई गणनीय ढाँचा प्रदान करता है
4. तकनीकी गहराई: मल्लियाविन विश्लेषण और निश्चित बिंदु सिद्धांत को चतुराई से जोड़ता है

कमियाँ

1. व्यावहारिक सत्यापन की कमी: पेपर विशुद्ध सैद्धांतिक है, संख्यात्मक प्रयोग सत्यापन की कमी है
2. गणनीय दक्षता: निश्चित बिंदु पुनरावृत्ति के अभिसरण गति और गणनीय जटिलता का गहन विश्लेषण नहीं है
3. मान्यता सीमाएँ: कुछ तकनीकी मान्यताएँ व्यावहारिक अनुप्रयोग में सत्यापित करना कठिन हो सकती हैं

प्रभाव

1. शैक्षणिक योगदान: भिन्नात्मक यादृच्छिक प्रक्रिया सांख्यिकी के लिए नई दिशा खोलता है
2. पद्धति मूल्य: निश्चित बिंदु रणनीति अन्य समान समस्याओं पर लागू हो सकती है
3. अनुप्रयोग संभावनाएँ: वित्त, जीव विज्ञान आदि क्षेत्रों में संभावित अनुप्रयोग मूल्य

लागू परिदृश्य

• वित्तीय बाजार में भिन्नात्मक ब्राउनियन गति मॉडल पैरामीटर अनुमान
• जैव चिकित्सा में कई व्यक्तिगत डेटा के आधार पर पैरामीटर अनुमान
• भौतिक प्रणाली में दीर्घ स्मृति प्रक्रिया मॉडलिंग

संदर्भ

पेपर 17 संबंधित संदर्भों का हवाला देता है, मुख्य रूप से:

• मल्लियाविन विश्लेषण मूल सिद्धांत (नुआलर्ट 15, डिक्रीसेफॉन्ड 4)
• भिन्नात्मक SDE सिद्धांत (हेयरर और ओहाशी 6)
• सांख्यिकीय अनुमान विधि (कॉम्टे और मैरी 2,3, मैरी 12)
• सन्निकटन सिद्धांत (डेवोर और लोरेंत्ज़ 5)

समग्र मूल्यांकन: यह भिन्नात्मक यादृच्छिक अवकल समीकरण के अ-पैरामीट्रिक अनुमान क्षेत्र में एक उच्च गुणवत्ता का सैद्धांतिक सांख्यिकी पेपर है जो महत्वपूर्ण योगदान देता है। यद्यपि संख्यात्मक प्रयोगों की कमी है, लेकिन सैद्धांतिक विश्लेषण गहन और कठोर है, जो इस क्षेत्र के आगे विकास के लिए एक मजबूत आधार स्थापित करता है।