चरघातांकी आनमन

चरघातांकी आनमन (ET), चरघातांकी व्यावर्तन, या चरघातांकी माप का परिवर्तन (ECM) एक वितरण स्थानांतरण तकनीक है जिसका उपयोग गणित के कई हिस्सों में किया जाता है। एक यादृच्छिक चर $$X$$ के विभिन्न चरघातांकी आनमन को $$X$$ के प्राकृतिक घातीय समूह के रूप में जाना जाता है।

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

सिंहावलोकन
एक यादृच्छिक चर दिया गया $$X$$ संभाव्यता वितरण के साथ $$\mathbb{P}$$, घनत्व $$f$$, और मोमेंट-जनरेटिंग_फंक्शन (एमजीएफ) $$M_X(\theta) = \mathbb{E}[e^{\theta X}] < \infty$$, चरचरघातांकी रूप से झुका हुआ माप $$\mathbb{P}_\theta$$ को इस प्रकार परिभाषित किया गया है:


 * $$\mathbb{P}_\theta(X \in dx) = \frac{\mathbb{E}[e^{\theta X}\mathbb{I}[X \in dx]]}{M_X(\theta)}=e^{\theta x - \kappa(\theta)}\mathbb{P}(X\in dx),$$

कहाँ $$\kappa(\theta)$$ संचयी  (सीजीएफ) के रूप में परिभाषित किया गया है


 * $$\kappa(\theta) = \log\mathbb{E}[e^{\theta X}] = \log M_X(\theta).$$ हम बुलाते है


 * $$\mathbb{P}_\theta(X \in dx)=f_\theta(x)$$ $$\theta$$-झुका हुआ संभाव्यता घनत्व फ़ंक्शन $$X$$. यह संतुष्ट करता है $$f_\theta(x) \propto e^{\theta x}f(x)$$.

एक यादृच्छिक वेक्टर का घातीय आनमन $$X$$ एक समान परिभाषा है:


 * $$\mathbb{P}_{\theta}(X\in dx) = e^{\theta^T x - \kappa(\theta)}\mathbb{P}(X\in dx),$$

कहाँ $$\kappa(\theta) = \log \mathbb{E}[\exp\{\theta^T X\}]$$.

उदाहरण
कई मामलों में चरचरघातांकी रूप से झुके हुए माप का पैरामीट्रिक रूप भी वैसा ही होता है $$X$$. एक-आयामी उदाहरणों में सामान्य वितरण, घातीय वितरण, द्विपद वितरण और पॉइसन वितरण शामिल हैं।

उदाहरण के लिए, सामान्य वितरण के मामले में, $$N( \mu, \sigma ^2)$$ झुका हुआ घनत्व $$f_\theta(x)$$ है $$N( \mu + \theta \sigma ^2, \sigma ^2)$$ घनत्व। नीचे दी गई तालिका झुके हुए घनत्व के अधिक उदाहरण प्रदान करती है।

हालाँकि, कुछ वितरणों के लिए, घातीय रूप से झुका हुआ वितरण उसी पैरामीट्रिक समूह से संबंधित नहीं है $$f$$. इसका एक उदाहरण पेरेटो वितरण है $$f(x) = \alpha /(1 + x) ^\alpha, x > 0$$, कहाँ $$f_\theta(x)$$ के लिए अच्छी तरह से परिभाषित है $$ \theta < 0 $$ लेकिन यह मानक वितरण नहीं है। ऐसे उदाहरणों में, यादृच्छिक परिवर्तनीय पीढ़ी हमेशा सीधी नहीं हो सकती है।

फायदे
कई मामलों में, झुका हुआ वितरण मूल के समान पैरामीट्रिक समूह से संबंधित होता है। यह विशेष रूप से सच है जब मूल घनत्व वितरण के घातीय समूह से संबंधित है। यह मोंटे-कार्लो सिमुलेशन के दौरान यादृच्छिक चर पीढ़ी को सरल बनाता है। यदि यह मामला नहीं है तो घातीय आनमन अभी भी उपयोगी हो सकता है, हालांकि सामान्यीकरण संभव होना चाहिए और अतिरिक्त नमूना एल्गोरिदम की आवश्यकता हो सकती है।

इसके अलावा, मूल और झुके हुए सीएफजी के बीच एक सरल संबंध मौजूद है,


 * $$\kappa_\theta(\eta) = \log(\mathbb{E}_\theta[e^{\eta X}]) = \kappa(\theta + \eta) - \kappa(\theta).$$ इसका अवलोकन हम कर सकते हैं


 * $$F_\theta(x) = \int\limits_{\infty}^x\exp\{\theta y - \kappa(\theta)\}f(y)dy.$$ इस प्रकार,



\begin{align} \kappa_{\theta}(\eta) &= \log \int e^{\eta x}dF_{\theta}(x) \\ &= \log \int e^{\eta x}e^{\theta x - \kappa(\theta)}dF(x) \\ &= \log\mathbb{E}[e^{(\eta+\theta)X-\kappa(\theta)}] \\ &= \log(e^{\kappa(\eta+\theta)-\kappa(\theta)}) \\ &= \kappa(\eta+\theta)-\kappa(\theta) \end{align} $$.

स्पष्ट रूप से, यह संबंध झुके हुए वितरण के सीजीएफ और इस प्रकार वितरण क्षणों की आसान गणना की अनुमति देता है। इसके अलावा, इसका परिणाम संभावना अनुपात का एक सरल रूप है। विशेष रूप से,


 * $$\ell = \frac{d\mathbb{P}}{d\mathbb{P}_\theta} = \frac{f(x)}{f_{\theta}(x)} = e^{- \theta x + \kappa(\theta)}$$.

गुण

 * अगर $$\kappa(\eta)=\log \mathrm{E}[\exp(\eta X)]$$ का सीजीएफ है $$X$$, फिर का सी.जी.एफ $$\theta$$-झुका हुआ $$X$$ है


 * $$\kappa_\theta(\eta) = \kappa(\theta + \eta) - \kappa(\theta).$$ :इसका मतलब यह है कि $$i$$झुका हुआ का -वाँ संचयी $$X$$ है $$\kappa^{(i)}(\theta)$$. खास तौर पर झुके हुए वितरण की अपेक्षा है


 * $$\mathrm{E}_\theta[X]=\tfrac{d}{d\eta}\kappa_\theta(\eta)|_{\eta=0} = \kappa'(\theta)$$.


 * झुके हुए वितरण का विचरण है


 * $$\mathrm{Var}_\theta[X]=\tfrac{d^2}{d\eta^2}\kappa_\theta(\eta)|_{\eta=0} = \kappa''(\theta)$$.


 * बार-बार झुकना योगात्मक है। यानी सबसे पहले झुकना $$\theta_1$$ और तब $$\theta_2$$ एक बार झुकने के समान है $$\theta_1+\theta_2$$.


 * अगर $$X$$ स्वतंत्र, लेकिन जरूरी नहीं कि समान यादृच्छिक चर का योग है $$X_1, X_2, \dots$$, फिर $$\theta$$- का झुका हुआ वितरण $$X$$ का योग है $$X_1, X_2, \dots$$ प्रत्येक $$\theta$$-व्यक्तिगत रूप से झुका हुआ.


 * अगर $$\mu=\mathrm{E}[X]$$, तब $$\kappa(\theta)-\theta \mu$$ कुल्बैक-लीब्लर विचलन है


 * $$D_\text{KL}(P \parallel P_\theta)=\mathrm{E} \left[\log\tfrac{P}{P_\theta}\right]$$ :झुके हुए वितरण के बीच $$P_\theta$$ और मूल वितरण $$P$$ का $$X$$.


 * इसी प्रकार, चूँकि $$\mathrm{E}_{\theta}[X]=\kappa'(\theta)$$, हमारे पास कुल्बैक-लीब्लर विचलन है


 * $$D_\text{KL}(P_\theta \parallel P) = \mathrm{E}_\theta \left[\log\tfrac{P_\theta}{P} \right] = \theta \kappa'(\theta) - \kappa(\theta)$$.

दुर्लभ-घटना अनुकरण
का घातीय आनमन $$X$$यह मानते हुए कि यह अस्तित्व में है, वितरण के एक समूह की आपूर्ति करता है जिसका उपयोग अस्वीकृति नमूने के लिए प्रस्ताव वितरण के रूप में किया जा सकता है। स्वीकृति-अस्वीकृति नमूनाकरण या महत्व नमूने के लिए महत्व वितरण। एक सामान्य अनुप्रयोग डोमेन के उप-क्षेत्र पर सशर्त वितरण से नमूना लेना है, अर्थात। $$X|X\in A$$. के उचित विकल्प के साथ $$\theta$$, से नमूनाकरण $$\mathbb{P}_\theta$$ नमूने की आवश्यक मात्रा या अनुमानक के विचरण को सार्थक रूप से कम कर सकता है।

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


 * $$f_{\theta}(\bar{x}) = f(\bar{x})\exp\{n(\theta \bar{x} - \kappa(\theta))\}$$.

के लिए एजवर्थ श्रृंखला को लागू करना $$f_{\theta}(\bar{x})$$, अपने पास


 * $$f_{\theta}(\bar{x}) = \psi(z)(\mathrm{Var}[\bar{X}])^{-1/2}\left\{1 + \frac{\rho_3(\theta)h_3(z)}{6} + \frac{\rho_4(\theta)h_4(z)}{24}\dots\right\},$$

कहाँ $$\psi(z)$$ का मानक सामान्य घनत्व है


 * $$z = \frac{\bar{x} - \kappa_{\bar{x}} ' (\theta)}{\kappa_{\bar{x}} ''(\theta)}$$,
 * $$\rho_n(\theta) = \kappa^{(n)}(\theta)\{\kappa ''(\theta)^{n/2}\}$$,

और $$h_n$$ हर्मिट बहुपद हैं.

के मूल्यों पर विचार करते समय $$\bar{x}$$ वितरण के केंद्र से उत्तरोत्तर दूर, $$|z|\rightarrow \infty $$ और यह $$h_n(z)$$ शर्तें असीमित हो जाती हैं। हालाँकि, प्रत्येक मान के लिए $$\bar{x}$$, हम चुन सकते हैं $$\theta$$ ऐसा है कि


 * $$\kappa '(\theta) = \bar{x}.$$

का यह मान $$\theta$$ इसे सैडल-पॉइंट के रूप में जाना जाता है, और उपरोक्त विस्तार का मूल्यांकन हमेशा झुके हुए वितरण की अपेक्षा पर किया जाता है। इस विकल्प का $$\theta$$ द्वारा दिए गए सन्निकटन के अंतिम प्रतिनिधित्व की ओर ले जाता है


 * $$f(\bar{x}) \approx \left(\frac{n}{2\pi\kappa ''(\theta)}\right)^{1/2}\exp\{n(\kappa(\theta) - \theta\bar{x})\}.$$

अस्वीकृति नमूनाकरण
झुके हुए वितरण का उपयोग करना $$\mathbb{P}_{\theta}$$ प्रस्ताव के रूप में, अस्वीकृति नमूनाकरण एल्गोरिदम से नमूनाकरण निर्धारित करता है $$f_\theta(x)$$ और संभाव्यता के साथ स्वीकार करना


 * $$\frac{1}{c} \exp(-\theta x + \kappa(\theta)),$$ कहाँ


 * $$c = \sup\limits_{x\in X}\frac{d\mathbb{P}}{d\mathbb{P}_{\theta}}(x).$$

अर्थात् एक समान रूप से वितरित यादृच्छिक चर $$p \sim \mbox{Unif}(0,1)$$ उत्पन्न होता है, और से नमूना $$f_\theta(x)$$ स्वीकार किया जाता है यदि


 * $$p \leq \frac{1}{c} \exp(-\theta x + \kappa(\theta)).$$

महत्वपूर्ण नमूनाकरण
घातीय रूप से झुके हुए वितरण को महत्व वितरण के रूप में लागू करने से समीकरण प्राप्त होता है


 * $$\mathbb{E}(h(X)) = \mathbb{E}_{\theta}[\ell(X)h(X)]$$,

कहाँ


 * $$\ell(X) = \frac{d\mathbb{P}}{d\mathbb{P}_{\theta}}$$

संभाव्यता फलन है. तो, से एक नमूना $$f_{\theta}$$ महत्व वितरण के अंतर्गत संभाव्यता का अनुमान लगाना $$\mathbb{P}(dX)$$ और फिर इसे संभावना अनुपात से गुणा कर देता है। इसके अलावा, हमारे पास इसके द्वारा दिया गया विचरण है


 * $$\mbox{Var}(X) = \mathbb{E}[(\ell(X)h(X)^2]$$.

उदाहरण
स्वतंत्र और समान रूप से वितरित मान लें $$\{X_i\}$$ ऐसा है कि $$\kappa(\theta) < \infty$$. अनुमान लगाने के लिए $$\mathbb{P}(X_1 + \cdots + X_n > c)$$, हम महत्व का नमूना लेकर उसे नियोजित कर सकते हैं


 * $$h(X) = \mathbb{I}(\sum_{i = 1}^n X_i > c)$$.

अटल $$c$$ के रूप में पुनः लिखा जा सकता है $$na$$ किसी अन्य स्थिरांक के लिए $$a$$. तब,


 * $$\mathbb{P}(\sum_{i = 1}^n X_i > na) = \mathbb{E}_{\theta_a} \left[\exp\{-\theta_a\sum_{i = 1}^n X_i + n\kappa(\theta_a)\}\mathbb{I}(\sum_{i = 1}^n X_i > na) \right]$$,

कहाँ $$\theta_a$$ को दर्शाता है $$\theta$$ सैडल-पॉइंट समीकरण द्वारा परिभाषित


 * $$\kappa '(\theta_a) = a$$.

स्टोकेस्टिक प्रक्रियाएं
एक सामान्य आर.वी. के आनमन को देखते हुए, यह सहज है कि घातीय आनमन $$X_t$$, बहाव के साथ एक एक प्रकार कि गति $$\mu$$ और विचरण $$\sigma^2$$, बहाव के साथ एक ब्राउनियन गति है $$\mu + \theta\sigma^2$$ और विचरण $$\sigma^2$$. इस प्रकार, बहाव के साथ कोई भी ब्राउनियन गति $$\mathbb{P}$$ बिना किसी बहाव के ब्राउनियन गति के रूप में सोचा जा सकता है $$\mathbb{P}_{\theta^*}$$. इसे देखने के लिए प्रक्रिया पर विचार करें $$X_t = B_t + \mu_t$$. $$f(X_t) = f_{\theta^*}(X_t)\frac{d\mathbb{P}}{d\mathbb{P}_{\theta^*}} = f(B_t)\exp\{\mu B_T - \frac{1}{2}\mu^2T\}$$. संभाव्यता अनुपात पद, $$\exp\{\mu B_{T} - \frac{1}{2}\mu^2 T\}$$, एक मार्टिंगेल (संभावना सिद्धांत) है और आमतौर पर निरूपित किया जाता है $$M_T$$. इस प्रकार, बहाव प्रक्रिया के साथ एक ब्राउनियन गति (साथ ही ब्राउनियन निस्पंदन के लिए अनुकूलित कई अन्य निरंतर प्रक्रियाएं) एक है $$\mathbb{P}_{\theta^*}$$-मार्टिंगेल.

स्टोकेस्टिक विभेदक समीकरण
उपरोक्त स्टोकेस्टिक विभेदक समीकरण के वैकल्पिक प्रतिनिधित्व की ओर ले जाता है $$dX(t) = \mu(t)dt + \sigma(t) dB(t)$$: $$dX_{\theta}(t) = \mu_{\theta}(t) dt + \sigma(t) dB(t)$$, कहाँ $$\mu_{\theta}(t)$$ = $$\mu(t) + \theta\sigma(t)$$. गिरसानोव का फॉर्मूला संभावना अनुपात बताता है $$\frac{d\mathbb{P}}{d\mathbb{P}_{\theta}} = \exp\{-\int\limits_0^T\frac{\mu_{\theta}(t) - \mu(t)}{\sigma^2(t)}dB(t) + \int\limits_0^T(\frac{\sigma^2(t)}{2})dt\}$$. इसलिए, गिरसानोव के फॉर्मूला का उपयोग कुछ एसडीई के लिए महत्व के नमूने को लागू करने के लिए किया जा सकता है।

किसी प्रक्रिया का अनुकरण करने के लिए आनमन भी उपयोगी हो सकता है $$X(t)$$ एसडीई के अस्वीकृति नमूने के माध्यम से $$dX(t) = \mu(X(t))dt+ dB(t)$$. हम एसडीई पर ध्यान केंद्रित कर सकते हैं क्योंकि हम यह जानते हैं $$X(t)$$ लिखा जा सकता है $$\int\limits_0^t dX(t) + X(0)$$. जैसा कि पहले कहा गया है, बहाव के साथ ब्राउनियन गति को बहाव के बिना ब्राउनियन गति में झुकाया जा सकता है। इसलिए, हम चुनते हैं $$\mathbb{P}_{proposal}=\mathbb{P}_{\theta^*}$$. संभाव्यता अनुपात $$\frac{d\mathbb{P}_{\theta^*}}{d\mathbb{P}}(dX(s): 0 \leq s \leq t) =$$ $$\prod\limits_{\tau\geq t}\exp\{\mu(X(\tau))dX(\tau) - \frac{\mu(X(\tau))^2}{2}\}dt = \exp\{\int\limits_0^t\mu(X(\tau))dX(\tau) - \int\limits_0^t\frac{\mu(X(s))^2}{2}\}dt$$. इस संभावना अनुपात को दर्शाया जाएगा $$M(t)$$. यह सुनिश्चित करने के लिए कि यह एक वास्तविक संभावना अनुपात है, इसे दिखाया जाना चाहिए $$\mathbb{E}[M(t)] = 1$$. यह स्थिति मानते हुए, यह दिखाया जा सकता है $$f_{X(t)}(y) = f_{X(t)}^{\theta^*}(y)\mathbb{E}_{\theta ^*}[M(t)|X(t) = y]$$. इसलिए, अस्वीकृति नमूनाकरण निर्धारित करता है कि एक मानक ब्राउनियन गति से नमूना लें और संभाव्यता के साथ स्वीकार करें $$\frac{f_{X(t)}(y)}{f_{X(t)}^{\theta ^*}(y)}\frac{1}{c} = \frac{1}{c}\mathbb{E}_{\theta ^*}[M(t)|X(t) = y]$$.

सिगमंड का एल्गोरिदम
मान लीजिए आई.आई.डी. एक्स लाइट टेल्ड डिस्ट्रीब्यूशन के साथ और $$\mathbb{E}[X] > 0$$. अनुमान लगाने के लिए $$\psi(c) = \mathbb{P}(\tau(c) < \infty)$$ कहाँ $$\tau(c) = \inf\{t:\sum\limits_{i=1}^t X_i> c\}$$, कब  $$c$$ बड़ा है और इसलिए  $$\psi(c)$$ छोटा, एल्गोरिथ्म महत्व वितरण प्राप्त करने के लिए घातीय आनमन का उपयोग करता है। एल्गोरिदम का उपयोग कई पहलुओं में किया जाता है, जैसे अनुक्रमिक परीक्षण, जी/जी/1 कतार प्रतीक्षा समय, और $$\psi$$ बर्बाद सिद्धांत में अंतिम बर्बादी की संभावना के रूप में उपयोग किया जाता है। इस संदर्भ में, यह सुनिश्चित करना तर्कसंगत है $$\mathbb{P}_\theta(\tau(c) < \infty) = 1$$. कसौटी $$\theta > \theta_0$$, कहाँ $$\theta_0$$ एस.टी. है $$\kappa'(\theta_0) = 0$$ इसे हासिल करता है. सिगमंड के एल्गोरिदम का उपयोग करता है $$\theta = \theta^*$$, यदि यह मौजूद है, तो कहां $$\theta^*$$ निम्नलिखित प्रकार से परिभाषित किया गया है: $$\kappa(\theta^*) = 0$$. ऐसा दिखाया गया है $$\theta^*$$ सीमित सापेक्ष त्रुटि उत्पन्न करने वाला एकमात्र आनमन पैरामीटर है ($$\underset{x \rightarrow \infty}{\lim\sup}\frac{Var\mathbb{I}_{A(x)}}{\mathbb{P}A(x)^2} < \infty$$).

ब्लैक-बॉक्स एल्गोरिदम
हम ब्लैक बॉक्स की संरचना को जाने बिना केवल उसके इनपुट और आउटपुट को देख सकते हैं। एल्गोरिदम को इसकी संरचना पर केवल न्यूनतम जानकारी का उपयोग करना है। जब हम यादृच्छिक संख्याएँ उत्पन्न करते हैं, तो आउटपुट नहीं हो सकता है समान सामान्य पैरामीट्रिक वर्ग के भीतर, जैसे सामान्य या चरघातांकी वितरण। ईसीएम करने के लिए स्वचालित तरीके का उपयोग किया जा सकता है। होने देना $$X_1, X_2,...$$आई.आई.डी. हो वितरण के साथ आर.वी $$G$$; सरलता के लिए हम मान लेते हैं $$X\geq 0$$. परिभाषित करना $$ \mathfrak{F}_n = \sigma(X_1,...,X_n,U_1,..., U_n) $$, कहाँ $$U_1, U_2$$,. . . स्वतंत्र (0, 1) वर्दी हैं। के लिए एक यादृच्छिक रुकने का समय $$X_1, X_2$$,. . . तब रुकने का समय w.r.t. है निस्पंदन $$ \{\mathfrak{F}_n\}$$,. . . आगे चलो $$ \mathfrak{G}$$ वितरण का एक वर्ग बनें $$G$$ पर $$ [0, \infty)$$ साथ $$ k_G = \int_0^\infty e^{\theta x}G(dx) < \infty$$ और परिभाषित करें $$G_\theta$$ द्वारा $$\frac{dG_\theta}{dG(x)} = e^{\theta x - k_G}$$. हम दिए गए के लिए ईसीएम के लिए एक ब्लैक-बॉक्स एल्गोरिदम परिभाषित करते हैं $$\theta$$ और दी गई कक्षा $$\mathfrak{G}$$यादृच्छिक रोक समय की एक जोड़ी के रूप में वितरण का $$\tau$$ और एक $$ \mathfrak{F}_\tau- $$ मापने योग्य आर.वी. $$Z $$ ऐसा है कि $$Z $$ के अनुसार वितरित किया जाता है $$G_\theta $$ किसी के लिए $$ G \in \mathfrak{G}$$. औपचारिक रूप से, हम इसे इस प्रकार लिखते हैं $$ \mathbb{P}_G (Z<x) = G_\theta (x) $$ सभी के लिए  $$x $$. दूसरे शब्दों में, गेम के नियम यह हैं कि एल्गोरिदम का उपयोग किया जा सकता है से सिम्युलेटेड मान  $$G $$ और आर.वी. तैयार करने के लिए अतिरिक्त वर्दी। से  $$G_\theta $$.

यह भी देखें

 * महत्व नमूनाकरण
 * अस्वीकृति नमूनाकरण
 * मोंटे कार्लो विधि
 * घातीय समूह
 * एस्चेर परिवर्तन