चरघातांकी आनमन: Difference between revisions

From Vigyanwiki
(Created page with "{{short description|Monte Carlo distribution shifting technique}} एक्सपोनेंशियल टिल्टिंग (ईटी), एक्सपोनेंश...")
 
No edit summary
 
(23 intermediate revisions by 4 users not shown)
Line 1: Line 1:
{{short description|Monte Carlo distribution shifting technique}}
{{short description|Monte Carlo distribution shifting technique}}
एक्सपोनेंशियल टिल्टिंग (ईटी), एक्सपोनेंशियल ट्विस्टिंग, या एक्सपोनेंशियल चेंज ऑफ मेजरमेंट (ईसीएम) एक वितरण स्थानांतरण तकनीक है जिसका उपयोग गणित के कई हिस्सों में किया जाता है।
'''चरघातांकी आनमन''' (ET), '''चरघातांकी व्यावर्तन''', या '''चरघातांकी माप का परिवर्तन''' (ECM) एक वितरण स्थानांतरण तकनीक है जिसका उपयोग गणित के कई हिस्सों में किया जाता है। एक यादृच्छिक चर <math>X</math> के विभिन्न चरघातांकी आनमन को <math>X</math> के [[प्राकृतिक घातीय परिवार|प्राकृतिक घातीय समूह]] के रूप में जाना जाता है।
एक यादृच्छिक चर के विभिन्न घातांकीय झुकाव <math>X</math> के [[प्राकृतिक घातीय परिवार]] के रूप में जाना जाता है <math>X</math>.


एक्सपोनेंशियल टिल्टिंग का उपयोग [[मोंटे कार्लो विधि]] में दुर्लभ-घटना सिमुलेशन और विशेष रूप से [[अस्वीकृति नमूनाकरण]] और [[महत्व नमूनाकरण]] के लिए किया जाता है।
चरघातांकी आनमन का उपयोग [[मोंटे कार्लो विधि|मोंटे कार्लो]] [[अनुमान]] में दुर्लभ-घटना अनुकरण और विशेष रूप से [[अस्वीकृति नमूनाकरण|अस्वीकृति]] और [[महत्व नमूनाकरण|महत्व प्रतिदर्श]] के लिए किया जाता है। गणितीय वित्त में <ref>{{Cite journal|last=H.U. Gerber & E.S.W. Shiu|date=1994|title=Esscher द्वारा विकल्प मूल्य निर्धारण परिवर्तन|journal=Transactions of the Society of Actuaries|volume=46|pages=99–191}}</ref> चरघातांकी आनमन को एस्चेर आनमन (या [[एस्चर परिवर्तन]]) के रूप में भी जाना जाता है, और इसे प्रायः अप्रत्यक्ष [[एजवर्थ श्रृंखला]] के साथ जोड़ा जाता है और इसका उपयोग बीमा वायदा मूल्य निर्धारण जैसे संदर्भों में किया जाता है।<ref>{{Cite book
गणितीय वित्त में <ref>{{Cite journal|last=H.U. Gerber & E.S.W. Shiu|date=1994|title=Esscher द्वारा विकल्प मूल्य निर्धारण परिवर्तन|journal=Transactions of the Society of Actuaries|volume=46|pages=99–191}}</ref> एक्सपोनेंशियल टिल्टिंग को एस्चेर टिल्टिंग (या [[एस्चर परिवर्तन]]) के रूप में भी जाना जाता है, और इसे अक्सर अप्रत्यक्ष [[एजवर्थ श्रृंखला]] के साथ जोड़ा जाता है और इसका उपयोग बीमा वायदा मूल्य निर्धारण जैसे संदर्भों में किया जाता है।<ref>{{Cite book
  |title = परिचालन जोखिम और बीमा विश्लेषण के मौलिक पहलू|last=Cruz|first=Marcelo|publisher=Wiley|year=2015|isbn=978-1-118-11839-9
  |title = परिचालन जोखिम और बीमा विश्लेषण के मौलिक पहलू|last=Cruz|first=Marcelo|publisher=Wiley|year=2015|isbn=978-1-118-11839-9
  |pages= 784–796
  |pages= 784–796
}}</ref>
}}</ref>
एक्सपोनेंशियल टिल्टिंग की प्रारंभिक औपचारिकता का श्रेय अक्सर [[Esscher]] को दिया जाता है<ref>{{Cite book|title=अनुप्रयोगों के साथ सैडल प्वाइंट सन्निकटन|url=https://archive.org/details/saddlepointappro00butl|url-access=limited|last=Butler|first=Ronald|publisher=Cambridge University Press|year=2007|isbn=9780521872508|pages=[https://archive.org/details/saddlepointappro00butl/page/n169 156]}}</ref> महत्व के नमूने में इसके उपयोग का श्रेय [[डेविड सिगमंड]] को दिया जाता है।<ref>{{cite journal
 
चरघातांकी आनमन की प्रारंभिक औपचारिकता का श्रेय प्रायः [[Esscher|एस्चेर]] को दिया जाता है<ref>{{Cite book|title=अनुप्रयोगों के साथ सैडल प्वाइंट सन्निकटन|url=https://archive.org/details/saddlepointappro00butl|url-access=limited|last=Butler|first=Ronald|publisher=Cambridge University Press|year=2007|isbn=9780521872508|pages=[https://archive.org/details/saddlepointappro00butl/page/n169 156]}}</ref> जबकि महत्व प्रतिदर्श में इसके उपयोग का श्रेय [[डेविड सिगमंड]] को दिया जाता है।<ref>{{cite journal
  |last1=Siegmund|first1=D.  
  |last1=Siegmund|first1=D.  
  |year = 1976
  |year = 1976
Line 17: Line 16:
|doi=10.1214/aos/1176343541  
|doi=10.1214/aos/1176343541  
  |doi-access= free}}</ref>
  |doi-access= free}}</ref>
 
== अवलोकन ==
 
प्रायिकता वितरण <math>\mathbb{P}</math>, घनत्व <math>f</math>, और [[आघुर्णजनक फलन]] (एमजीएफ) <math>M_X(\theta) = \mathbb{E}[e^{\theta X}] < \infty</math> के साथ एक यादृच्छिक चर <math>X</math> को देखते हुए, चरघातांकी रूप से आनत माप <math>\mathbb{P}_\theta</math> को इस प्रकार परिभाषित किया गया है,
== सिंहावलोकन ==
एक यादृच्छिक चर दिया गया <math>X</math> संभाव्यता वितरण के साथ <math>\mathbb{P}</math>, घनत्व <math>f</math>, और मोमेंट-जनरेटिंग_फंक्शन (एमजीएफ) <math>M_X(\theta) = \mathbb{E}[e^{\theta X}] < \infty</math>, चरघातांकीय रूप से झुका हुआ माप <math>\mathbb{P}_\theta</math> को इस प्रकार परिभाषित किया गया है:


:<math>\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),</math>
:<math>\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),</math>
कहाँ <math>\kappa(\theta)</math> [[ संचयी ]] (सीजीएफ) के रूप में परिभाषित किया गया है
जहां <math>\kappa(\theta)</math>[[ संचयी | संचयी जनक फलन]] (सीजीएफ) है जिसे
 
:<math>\kappa(\theta) = \log\mathbb{E}[e^{\theta X}] = \log M_X(\theta).</math> हम बुलाते है


:<math>\mathbb{P}_\theta(X \in dx)=f_\theta(x)</math> <math>\theta</math>-झुका हुआ संभाव्यता घनत्व फ़ंक्शन <math>X</math>. यह संतुष्ट करता है <math>f_\theta(x) \propto e^{\theta x}f(x)</math>.
:<math>\kappa(\theta) = \log\mathbb{E}[e^{\theta X}] = \log M_X(\theta).</math> के रूप में परिभाषित किया गया है।
हम <math>\mathbb{P}_\theta(X \in dx)=f_\theta(x)</math> को <math>\theta</math>-का आनत [[घनत्व]] <math>X</math> कहते हैं। यह <math>f_\theta(x) \propto e^{\theta x}f(x)</math>. को संतुष्ट करता है।


एक यादृच्छिक वेक्टर का घातीय झुकाव <math>X</math> एक समान परिभाषा है:
एक यादृच्छिक सदिश <math>X</math> के घातीय आनमन की एक समान परिभाषा है,


:<math>\mathbb{P}_{\theta}(X\in dx) = e^{\theta^T x - \kappa(\theta)}\mathbb{P}(X\in dx),</math>
:<math>\mathbb{P}_{\theta}(X\in dx) = e^{\theta^T x - \kappa(\theta)}\mathbb{P}(X\in dx),</math>
कहाँ <math>\kappa(\theta) = \log \mathbb{E}[\exp\{\theta^T X\}]</math>.
जहां <math>\kappa(\theta) = \log \mathbb{E}[\exp\{\theta^T X\}]</math> दिया गया है।


===उदाहरण===
===उदाहरण===


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


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


{| class="wikitable"
{| class="wikitable"
|-
|-
! Original distribution<ref>{{Cite book|last=Asmussen Soren & Glynn Peter|title=Stochastic Simulation|publisher=Springer|year=2007|isbn=978-0-387-30679-7|pages=130}}</ref><ref>{{cite journal|last1=Fuh|first1=Cheng-Der|last2=Teng|first2=Huei-Wen|last3=Wang|first3=Ren-Her|title=Efficient Importance Sampling for Rare Event Simulation with Applications|url=https://archive.org/details/arxiv-1302.0583|date=2013}}</ref>
! मूल वितरण<ref>{{Cite book|last=Asmussen Soren & Glynn Peter|title=Stochastic Simulation|publisher=Springer|year=2007|isbn=978-0-387-30679-7|pages=130}}</ref><ref>{{cite journal|last1=Fuh|first1=Cheng-Der|last2=Teng|first2=Huei-Wen|last3=Wang|first3=Ren-Her|title=Efficient Importance Sampling for Rare Event Simulation with Applications|url=https://archive.org/details/arxiv-1302.0583|date=2013}}</ref>
! θ-Tilted distribution
! θ-आनत वितरण
|-
|-
|  <math>\mathrm{Gamma} (\alpha, \beta)</math>
|  <math>\mathrm{Gamma} (\alpha, \beta)</math>
Line 66: Line 62:
|<math>\mathrm{Gamma}\left(\frac{\kappa}{2}, \frac{2}{1-2\theta}\right)</math>
|<math>\mathrm{Gamma}\left(\frac{\kappa}{2}, \frac{2}{1-2\theta}\right)</math>
|}
|}
हालाँकि, कुछ वितरणों के लिए, घातीय रूप से झुका हुआ वितरण उसी पैरामीट्रिक परिवार से संबंधित नहीं है <math>f</math>. इसका एक उदाहरण [[पेरेटो वितरण]] है <math>f(x) = \alpha /(1 + x) ^\alpha, x > 0</math>, कहाँ <math>f_\theta(x)</math> के लिए अच्छी तरह से परिभाषित है <math> \theta < 0 </math> लेकिन यह मानक वितरण नहीं है। ऐसे उदाहरणों में, यादृच्छिक परिवर्तनीय पीढ़ी हमेशा सीधी नहीं हो सकती है।<ref>Asmussen, Soren & Glynn, Peter (2007). Stochastic Simulation. Springer. pp. 164–167. {{ISBN|978-0-387-30679-7}}</ref>
हालाँकि, कुछ वितरणों के लिए, घातीय रूप से आनत वितरण <math>f</math> के समान प्राचलिक समूह से संबंधित नहीं है। इसका एक उदाहरण <math>f(x) = \alpha /(1 + x) ^\alpha, x > 0</math> [[पेरेटो वितरण]] है, जहां <math>f_\theta(x)</math> को <math> \theta < 0 </math> के लिए अच्छी तरह से परिभाषित किया गया है लेकिन यह एक मानक वितरण नहीं है। ऐसे उदाहरणों में, यादृच्छिक परिवर्तनीय पीढ़ी हमेशा स्पष्ट नहीं हो सकती है।<ref>Asmussen, Soren & Glynn, Peter (2007). Stochastic Simulation. Springer. pp. 164–167. {{ISBN|978-0-387-30679-7}}</ref>
 
===लाभ===
 
कई स्थितियों में, आनत वितरण मूल के समान प्राचलिक समूह से संबंधित होता है। यह विशेष रूप से सच है कि एक मूल घनत्व वितरण [[घातीय परिवार|घातीय समूह]] से संबंधित होता है। यह मोंटे-कार्लो अनुकरण के दौरान यादृच्छिक चर पीढ़ी को सरल बनाता है। यदि यह स्थिति नहीं है तो घातीय आनमन अभी भी उपयोगी हो सकता है, हालांकि सामान्यीकरण संभव होना चाहिए क्योकि अतिरिक्त प्रतिदर्श कलन विधि की आवश्यकता हो सकती है।
===फायदे===
कई मामलों में, झुका हुआ वितरण मूल के समान पैरामीट्रिक परिवार से संबंधित होता है। यह विशेष रूप से सच है जब मूल घनत्व वितरण के [[घातीय परिवार]] से संबंधित है। यह मोंटे-कार्लो सिमुलेशन के दौरान यादृच्छिक चर पीढ़ी को सरल बनाता है। यदि यह मामला नहीं है तो घातीय झुकाव अभी भी उपयोगी हो सकता है, हालांकि सामान्यीकरण संभव होना चाहिए और अतिरिक्त नमूना एल्गोरिदम की आवश्यकता हो सकती है।
 
इसके अलावा, मूल और झुके हुए सीएफजी के बीच एक सरल संबंध मौजूद है,


:<math>\kappa_\theta(\eta) = \log(\mathbb{E}_\theta[e^{\eta X}]) = \kappa(\theta + \eta) - \kappa(\theta).</math> इसका अवलोकन हम कर सकते हैं
इसके अलावा, मूल और आनत सीएफजी,
 
:<math>F_\theta(x) = \int\limits_{\infty}^x\exp\{\theta y - \kappa(\theta)\}f(y)dy.</math> इस प्रकार,


:<math>\kappa_\theta(\eta) = \log(\mathbb{E}_\theta[e^{\eta X}]) = \kappa(\theta + \eta) - \kappa(\theta).</math> के बीच एक सरल संबंध उपस्थित है। इसका अवलोकन हम इस प्रकार कर सकते हैं,<math>F_\theta(x) = \int\limits_{\infty}^x\exp\{\theta y - \kappa(\theta)\}f(y)dy.</math>।
:इस प्रकार से,
:<math>
:<math>
\begin{align}
\begin{align}
Line 88: Line 80:
</math>.
</math>.


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


:<math>\ell = \frac{d\mathbb{P}}{d\mathbb{P}_\theta} = \frac{f(x)}{f_{\theta}(x)} = e^{- \theta x + \kappa(\theta)}</math>.
:<math>\ell = \frac{d\mathbb{P}}{d\mathbb{P}_\theta} = \frac{f(x)}{f_{\theta}(x)} = e^{- \theta x + \kappa(\theta)}</math>. सरल रूप है।


== गुण ==
== गुण ==


* अगर <math>\kappa(\eta)=\log \mathrm{E}[\exp(\eta X)]</math> का सीजीएफ है <math>X</math>, फिर का सी.जी.एफ <math>\theta</math>-झुका हुआ <math>X</math> है
* यदि <math>\kappa(\eta)=\log \mathrm{E}[\exp(\eta X)]</math>, <math>X</math> का सीजीएफ है, तो <math>X</math> आनत <math>\theta</math>- का सीजीएफ
 
::<math>\kappa_\theta(\eta) = \kappa(\theta + \eta) - \kappa(\theta).</math> :इसका मतलब यह है कि <math>i</math>झुका हुआ का -वाँ संचयी <math>X</math> है <math>\kappa^{(i)}(\theta)</math>. खास तौर पर झुके हुए वितरण की अपेक्षा है
 
::<math>\mathrm{E}_\theta[X]=\tfrac{d}{d\eta}\kappa_\theta(\eta)|_{\eta=0} = \kappa'(\theta)</math>.


:झुके हुए वितरण का विचरण है
::<math>\kappa_\theta(\eta) = \kappa(\theta + \eta) - \kappa(\theta).</math>है। इसका मतलब यह है कि आनत <math>X</math> का <math>i</math> -वाँ [[संचयी]] <math>\kappa^{(i)}(\theta)</math> है। विशेष रूप से, आनत वितरण की अपेक्षा<math>\mathrm{E}_\theta[X]=\tfrac{d}{d\eta}\kappa_\theta(\eta)|_{\eta=0} = \kappa'(\theta)</math> है।
::आनत वितरण का विचरण  
::<math>\mathrm{Var}_\theta[X]=\tfrac{d^2}{d\eta^2}\kappa_\theta(\eta)|_{\eta=0} = \kappa''(\theta)</math>. है।


::<math>\mathrm{Var}_\theta[X]=\tfrac{d^2}{d\eta^2}\kappa_\theta(\eta)|_{\eta=0} = \kappa''(\theta)</math>.
* पुनरावर्ती आनत योगात्मक है। अर्थात् पहले <math>\theta_1</math> और फिर <math>\theta_2</math> से आनत एक बार <math>\theta_1+\theta_2</math> के आनत के समान है।


* बार-बार झुकना योगात्मक है। यानी सबसे पहले झुकना <math>\theta_1</math> और तब <math>\theta_2</math> एक बार झुकने के समान है <math>\theta_1+\theta_2</math>.
* यदि <math>X</math> स्वतंत्र, लेकिन आवश्यक रूप से समान यादृच्छिक चर <math>X_1, X_2, \dots</math> का योग नहीं है, तो <math>X</math> का <math>\theta</math>- आनत वितरण प्रत्येक -<math>\theta</math> को व्यक्तिगत रूप से आनत <math>X_1, X_2, \dots</math> का योग है।


* अगर <math>X</math> स्वतंत्र, लेकिन जरूरी नहीं कि समान यादृच्छिक चर का योग है <math>X_1, X_2, \dots</math>, फिर <math>\theta</math>- का झुका हुआ वितरण <math>X</math> का योग है <math>X_1, X_2, \dots</math> प्रत्येक <math>\theta</math>-व्यक्तिगत रूप से झुका हुआ.
* यदि <math>\mu=\mathrm{E}[X]</math>, तो <math>\kappa(\theta)-\theta \mu</math> आनत वितरण <math>P_\theta</math> और <math>X</math> के मूल वितरण <math>P</math> के बीच [[कुल्बैक-लीबलर विचलन]] <math>D_\text{KL}(P \parallel P_\theta)=\mathrm{E} \left[\log\tfrac{P}{P_\theta}\right]</math>है।


* अगर <math>\mu=\mathrm{E}[X]</math>, तब <math>\kappa(\theta)-\theta \mu</math> कुल्बैक-लीब्लर विचलन है
* इसी प्रकार, <math>\mathrm{E}_{\theta}[X]=\kappa'(\theta)</math> के बाद से, हमारे पास,


::<math>D_\text{KL}(P \parallel P_\theta)=\mathrm{E} \left[\log\tfrac{P}{P_\theta}\right]</math> :झुके हुए वितरण के बीच <math>P_\theta</math> और मूल वितरण <math>P</math> का <math>X</math>.
::<math>D_\text{KL}(P_\theta \parallel P) = \mathrm{E}_\theta \left[\log\tfrac{P_\theta}{P} \right] = \theta \kappa'(\theta) - \kappa(\theta)</math> के रूप में कुल्बैक-लीबलर विचलन है।
 
* इसी प्रकार, चूँकि <math>\mathrm{E}_{\theta}[X]=\kappa'(\theta)</math>, हमारे पास कुल्बैक-लीब्लर विचलन है
 
::<math>D_\text{KL}(P_\theta \parallel P) = \mathrm{E}_\theta \left[\log\tfrac{P_\theta}{P} \right] = \theta \kappa'(\theta) - \kappa(\theta)</math>.


==अनुप्रयोग==
==अनुप्रयोग==


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


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


:<math>f_{\theta}(\bar{x}) = f(\bar{x})\exp\{n(\theta \bar{x} - \kappa(\theta))\}</math>.
:<math>f_{\theta}(\bar{x}) = f(\bar{x})\exp\{n(\theta \bar{x} - \kappa(\theta))\}</math> है।


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


:<math>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\},</math>
:<math>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\},</math>प्राप्त कर सकते है,
कहाँ <math>\psi(z)</math> का मानक सामान्य घनत्व है
जहां <math>\psi(z)</math> ,


:<math>z = \frac{\bar{x} - \kappa_{\bar{x}} ' (\theta)}{\kappa_{\bar{x}} ''(\theta)}</math>,  
:<math>z = \frac{\bar{x} - \kappa_{\bar{x}} ' (\theta)}{\kappa_{\bar{x}} ''(\theta)}</math> का मानक सामान्य घनत्व है,
:<math>\rho_n(\theta) = \kappa^{(n)}(\theta)\{\kappa ''(\theta)^{n/2}\}</math>,
:<math>\rho_n(\theta) = \kappa^{(n)}(\theta)\{\kappa ''(\theta)^{n/2}\}</math>,


और <math>h_n</math> हर्मिट बहुपद हैं.
और <math>h_n</math> [[हर्मिट बहुपद]] हैं।


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


::<math>\kappa '(\theta) = \bar{x}.</math>
::<math>\kappa '(\theta) = \bar{x}.</math>
का यह मान <math>\theta</math> इसे सैडल-पॉइंट के रूप में जाना जाता है, और उपरोक्त विस्तार का मूल्यांकन हमेशा झुके हुए वितरण की अपेक्षा पर किया जाता है। इस विकल्प का <math>\theta</math> द्वारा दिए गए सन्निकटन के अंतिम प्रतिनिधित्व की ओर ले जाता है
<math>\theta</math> के इस मान को सैडल-बिंदु के रूप में जाना जाता है, और उपरोक्त विस्तार का मूल्यांकन हमेशा आनत वितरण की अपेक्षा पर किया जाता है। <math>\theta</math> का यह विकल्प


:<math>f(\bar{x}) \approx \left(\frac{n}{2\pi\kappa ''(\theta)}\right)^{1/2}\exp\{n(\kappa(\theta) - \theta\bar{x})\}.</math><ref>{{Cite book|title=अनुप्रयोगों के साथ सैडल प्वाइंट सन्निकटन|url=https://archive.org/details/saddlepointappro00butl|url-access=limited|last=Butler|first=Ronald|publisher=Cambridge University Press|year=2007|isbn=9780521872508|pages=[https://archive.org/details/saddlepointappro00butl/page/n169 156]–157}}</ref><ref>{{Cite book|title=जीएलआईएम और सांख्यिकीय मॉडलिंग में प्रगति|last=Seeber|first=G.U.H.|publisher=Springer|year=1992|isbn=978-0-387-97873-4 |pages=195–200}}</ref>
:<math>f(\bar{x}) \approx \left(\frac{n}{2\pi\kappa ''(\theta)}\right)^{1/2}\exp\{n(\kappa(\theta) - \theta\bar{x})\}.</math> द्वारा दिए गए सन्निकटन के अंतिम प्रतिनिधित्व की ओर ले जाता है।<ref>{{Cite book|title=अनुप्रयोगों के साथ सैडल प्वाइंट सन्निकटन|url=https://archive.org/details/saddlepointappro00butl|url-access=limited|last=Butler|first=Ronald|publisher=Cambridge University Press|year=2007|isbn=9780521872508|pages=[https://archive.org/details/saddlepointappro00butl/page/n169 156]–157}}</ref><ref>{{Cite book|title=जीएलआईएम और सांख्यिकीय मॉडलिंग में प्रगति|last=Seeber|first=G.U.H.|publisher=Springer|year=1992|isbn=978-0-387-97873-4 |pages=195–200}}</ref>


===अस्वीकृति प्रतिदर्श===
प्रस्ताव के रूप में आनत वितरण <math>\mathbb{P}_{\theta}</math> का उपयोग करते हुए, अस्वीकृति प्रतिदर्श कलन विधि <math>f_\theta(x)</math> से प्रतिदर्श लेने और प्रायिकता


===अस्वीकृति नमूनाकरण===
:<math>\frac{1}{c} \exp(-\theta x + \kappa(\theta)),</math> के साथ स्वीकार करने को निर्धारित करता है जहां
झुके हुए वितरण का उपयोग करना <math>\mathbb{P}_{\theta}</math> प्रस्ताव के रूप में, अस्वीकृति नमूनाकरण एल्गोरिदम से नमूनाकरण निर्धारित करता है <math>f_\theta(x)</math> और संभाव्यता के साथ स्वीकार करना


:<math>\frac{1}{c} \exp(-\theta x + \kappa(\theta)),</math> कहाँ
:<math>c = \sup\limits_{x\in X}\frac{d\mathbb{P}}{d\mathbb{P}_{\theta}}(x).</math> दिया गया है।
अर्थात् एक समान रूप से वितरित यादृच्छिक चर <math>p \sim \mbox{Unif}(0,1)</math> उत्पन्न होता है, और  <math>f_\theta(x)</math> से प्रतिदर्श स्वीकार किया जाता है यदि


:<math>c = \sup\limits_{x\in X}\frac{d\mathbb{P}}{d\mathbb{P}_{\theta}}(x).</math>
:<math>p \leq \frac{1}{c} \exp(-\theta x + \kappa(\theta)).</math>दिया गया हो।
अर्थात् एक समान रूप से वितरित यादृच्छिक चर <math>p \sim \mbox{Unif}(0,1)</math> उत्पन्न होता है, और से नमूना <math>f_\theta(x)</math> स्वीकार किया जाता है यदि
===महत्व प्रतिदर्श===
घातीय रूप से आनत वितरण को महत्व वितरण के रूप में लागू करने से समीकरण


:<math>p \leq \frac{1}{c} \exp(-\theta x + \kappa(\theta)).</math>
:<math>\mathbb{E}(h(X)) = \mathbb{E}_{\theta}[\ell(X)h(X)]</math> प्राप्त होता है,


 
जहां
===महत्वपूर्ण नमूनाकरण===
घातीय रूप से झुके हुए वितरण को महत्व वितरण के रूप में लागू करने से समीकरण प्राप्त होता है
 
:<math>\mathbb{E}(h(X)) = \mathbb{E}_{\theta}[\ell(X)h(X)]</math>,
 
कहाँ


:<math>\ell(X) = \frac{d\mathbb{P}}{d\mathbb{P}_{\theta}}</math>
:<math>\ell(X) = \frac{d\mathbb{P}}{d\mathbb{P}_{\theta}}</math>
संभाव्यता फलन है. तो, से एक नमूना <math>f_{\theta}</math> महत्व वितरण के अंतर्गत संभाव्यता का अनुमान लगाना <math>\mathbb{P}(dX)</math> और फिर इसे संभावना अनुपात से गुणा कर देता है। इसके अलावा, हमारे पास इसके द्वारा दिया गया विचरण है
संभाव्यता फलन है। तो, महत्व वितरण <math>\mathbb{P}(dX)</math> के तहत प्रायिकता का अनुमान लगाने के लिए <math>f_{\theta}</math> से एक प्रतिदर्श लें और फिर इसे संभावना अनुपात से गुणा करें। इसके अलावा, हमारे पास  


:<math>\mbox{Var}(X) = \mathbb{E}[(\ell(X)h(X)^2]</math>.
:<math>\mbox{Var}(X) = \mathbb{E}[(\ell(X)h(X)^2]</math> द्वारा दिया गया विचरण है।


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


:<math>h(X) = \mathbb{I}(\sum_{i = 1}^n X_i > c)</math>.
:<math>h(X) = \mathbb{I}(\sum_{i = 1}^n X_i > c)</math> लेकर महत्व प्रतिदर्श का उपयोग कर सकते हैं।


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


:<math>\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]</math>,
:<math>\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]</math>,


कहाँ <math>\theta_a</math> को दर्शाता है <math>\theta</math> सैडल-पॉइंट समीकरण द्वारा परिभाषित
जहां <math>\theta_a</math> सैडल-बिंदु समीकरण <math>\kappa '(\theta_a) = a</math> द्वारा परिभाषित  <math>\theta</math> को दर्शाता है।
 
:<math>\kappa '(\theta_a) = a</math>.


===स्टोकेस्टिक प्रक्रियाएं===
===प्रसंभाव्य प्रक्रम===
एक सामान्य आर.वी. के झुकाव को देखते हुए, यह सहज है कि घातीय झुकाव <math>X_t</math>, बहाव के साथ एक [[एक प्रकार कि गति]] <math>\mu</math> और विचरण <math>\sigma^2</math>, बहाव के साथ एक ब्राउनियन गति है <math>\mu + \theta\sigma^2</math> और विचरण <math>\sigma^2</math>. इस प्रकार, बहाव के साथ कोई भी ब्राउनियन गति <math>\mathbb{P}</math> बिना किसी बहाव के ब्राउनियन गति के रूप में सोचा जा सकता है <math>\mathbb{P}_{\theta^*}</math>. इसे देखने के लिए प्रक्रिया पर विचार करें <math>X_t = B_t + \mu_t</math>. <math>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\}</math>. संभाव्यता अनुपात पद, <math>\exp\{\mu B_{T} - \frac{1}{2}\mu^2 T\}</math>, एक [[मार्टिंगेल (संभावना सिद्धांत)]] है और आमतौर पर निरूपित किया जाता है <math>M_T</math>. इस प्रकार, बहाव प्रक्रिया के साथ एक ब्राउनियन गति (साथ ही ब्राउनियन निस्पंदन के लिए अनुकूलित कई अन्य निरंतर प्रक्रियाएं) एक है <math>\mathbb{P}_{\theta^*}</math>-मार्टिंगेल.<ref>{{Cite book|last=Asmussen Soren & Glynn Peter|title=स्टोकेस्टिक सिमुलेशन|publisher=Springer|year=2007|isbn=978-0-387-30679-7|pages=407}}</ref><ref>{{Cite book|title=स्टोकेस्टिक कैलकुलस और वित्तीय अनुप्रयोग|url=https://archive.org/details/stochasticcalcul00stee|url-access=limited|last=Steele|first=J. Michael|publisher=Springer|year=2001|isbn=978-1-4419-2862-7 |pages=[https://archive.org/details/stochasticcalcul00stee/page/n224 213]–229}}</ref>
एक सामान्य आर.वी. के आनमन को देखते हुए, यह सहज है कि <math>X_t</math> का घातांकीय आनमन, प्रक्षेप <math>\mu</math> विचरण <math>\sigma^2</math>, प्रक्षेप के साथ एक [[ब्राउनियन गति]] है, और प्रक्षेप <math>\mu + \theta\sigma^2</math> और विचरण <math>\sigma^2</math> के साथ एक ब्राउनियन गति है। इस प्रकार, <math>\mathbb{P}</math> के नीचे प्रक्षेप के साथ किसी भी ब्राउनियन गति को <math>\mathbb{P}_{\theta^*}</math>के नीचे प्रक्षेप के बिना ब्राउनियन गति के रूप में माना जा सकता है। इसे देखने के लिए प्रक्रिया <math>X_t = B_t + \mu_t</math> पर विचार करें। <math>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\}</math>. संभाव्यता अनुपात शब्द, <math>\exp\{\mu B_{T} - \frac{1}{2}\mu^2 T\}</math>, एक [[मार्टिंगेल (संभावना सिद्धांत)|मार्टिंगेल]] है और आमतौर पर <math>M_T</math> द्वारा निरूपित किया जाता है। इस प्रकार, प्रक्षेप प्रक्रिया (साथ ही ब्राउनियन निस्पंदन के लिए अनुकूलित कई अन्य निरंतर प्रक्रियाएं) के साथ एक ब्राउनियन गति एक <math>\mathbb{P}_{\theta^*}</math>-मार्टिंगेल है ।<ref>{{Cite book|last=Asmussen Soren & Glynn Peter|title=स्टोकेस्टिक सिमुलेशन|publisher=Springer|year=2007|isbn=978-0-387-30679-7|pages=407}}</ref><ref>{{Cite book|title=स्टोकेस्टिक कैलकुलस और वित्तीय अनुप्रयोग|url=https://archive.org/details/stochasticcalcul00stee|url-access=limited|last=Steele|first=J. Michael|publisher=Springer|year=2001|isbn=978-1-4419-2862-7 |pages=[https://archive.org/details/stochasticcalcul00stee/page/n224 213]–229}}</ref>
====प्रसंभाव्य विभेदक समीकरण====
उपरोक्त [[प्रसंभाव्य विभेदक समीकरण]] <math>dX(t) = \mu(t)dt + \sigma(t) dB(t)</math> के वैकल्पिक प्रतिनिधित्व की ओर ले जाता है, जहाँ <math>dX_{\theta}(t) = \mu_{\theta}(t) dt + \sigma(t) dB(t)</math> और <math>\mu_{\theta}(t)</math> = <math>\mu(t) + \theta\sigma(t)</math> दिया गया है।गिरसानोव का सूत्र संभावना अनुपात <math>\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\}</math> बताता है। इसलिए, गिरसानोव के सूत्र का उपयोग कुछ एसडीई के महत्व प्रतिदर्श को लागू करने के लिए किया जा सकता है।


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


====[[स्टोकेस्टिक विभेदक समीकरण]]====
==आनमन प्राचल का विकल्प==
उपरोक्त स्टोकेस्टिक विभेदक समीकरण के वैकल्पिक प्रतिनिधित्व की ओर ले जाता है <math>dX(t) = \mu(t)dt + \sigma(t) dB(t)</math>: <math>dX_{\theta}(t) = \mu_{\theta}(t) dt + \sigma(t) dB(t)</math>, कहाँ <math>\mu_{\theta}(t)</math> = <math>\mu(t) + \theta\sigma(t)</math>. गिरसानोव का फॉर्मूला संभावना अनुपात बताता है <math>\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\}</math>. इसलिए, गिरसानोव के फॉर्मूला का उपयोग कुछ एसडीई के लिए महत्व के नमूने को लागू करने के लिए किया जा सकता है।
 
किसी प्रक्रिया का अनुकरण करने के लिए झुकाव भी उपयोगी हो सकता है <math>X(t)</math> एसडीई के अस्वीकृति नमूने के माध्यम से <math>dX(t) = \mu(X(t))dt+ dB(t)</math>. हम एसडीई पर ध्यान केंद्रित कर सकते हैं क्योंकि हम यह जानते हैं <math>X(t)</math> लिखा जा सकता है <math>\int\limits_0^t dX(t) + X(0)</math>. जैसा कि पहले कहा गया है, बहाव के साथ ब्राउनियन गति को बहाव के बिना ब्राउनियन गति में झुकाया जा सकता है। इसलिए, हम चुनते हैं <math>\mathbb{P}_{proposal}=\mathbb{P}_{\theta^*}</math>. संभाव्यता अनुपात <math>\frac{d\mathbb{P}_{\theta^*}}{d\mathbb{P}}(dX(s): 0 \leq s \leq t) =</math>
<math>\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</math>. इस संभावना अनुपात को दर्शाया जाएगा <math>M(t)</math>. यह सुनिश्चित करने के लिए कि यह एक वास्तविक संभावना अनुपात है, इसे दिखाया जाना चाहिए <math>\mathbb{E}[M(t)] = 1</math>. यह स्थिति मानते हुए, यह दिखाया जा सकता है <math>f_{X(t)}(y) = f_{X(t)}^{\theta^*}(y)\mathbb{E}_{\theta ^*}[M(t)|X(t) = y]</math>. इसलिए, अस्वीकृति नमूनाकरण निर्धारित करता है कि एक मानक ब्राउनियन गति से नमूना लें और संभाव्यता के साथ स्वीकार करें <math>\frac{f_{X(t)}(y)}{f_{X(t)}^{\theta ^*}(y)}\frac{1}{c} = \frac{1}{c}\mathbb{E}_{\theta ^*}[M(t)|X(t) = y]</math>.
 
==झुकाव पैरामीटर का विकल्प==
 
===सिगमंड का एल्गोरिदम===
मान लीजिए आई.आई.डी. एक्स लाइट टेल्ड डिस्ट्रीब्यूशन के साथ और <math>\mathbb{E}[X] > 0</math>. अनुमान लगाने के लिए  <math>\psi(c) = \mathbb{P}(\tau(c) < \infty)</math> कहाँ <math>\tau(c) = \inf\{t:\sum\limits_{i=1}^t X_i> c\}</math>, कब  <math>c</math> बड़ा है और इसलिए  <math>\psi(c)</math> छोटा, एल्गोरिथ्म महत्व वितरण प्राप्त करने के लिए घातीय झुकाव का उपयोग करता है। एल्गोरिदम का उपयोग कई पहलुओं में किया जाता है, जैसे अनुक्रमिक परीक्षण,<ref>D. Siegmund (1985) Sequential Analysis. Springer-Verlag</ref> जी/जी/1 कतार प्रतीक्षा समय, और <math>\psi</math> [[बर्बाद सिद्धांत]] में अंतिम बर्बादी की संभावना के रूप में उपयोग किया जाता है। इस संदर्भ में, यह सुनिश्चित करना तर्कसंगत है <math>\mathbb{P}_\theta(\tau(c) < \infty) = 1</math>. कसौटी <math>\theta > \theta_0</math>, कहाँ <math>\theta_0</math> एस.टी. है <math>\kappa'(\theta_0) = 0</math> इसे हासिल करता है. सिगमंड के एल्गोरिदम का उपयोग करता है <math>\theta = \theta^*</math>, यदि यह मौजूद है, तो कहां <math>\theta^*</math> निम्नलिखित प्रकार से परिभाषित किया गया है:
<math>\kappa(\theta^*) = 0</math>.
ऐसा दिखाया गया है <math>\theta^*</math> सीमित सापेक्ष त्रुटि उत्पन्न करने वाला एकमात्र झुकाव पैरामीटर है (<math>\underset{x \rightarrow \infty}{\lim\sup}\frac{Var\mathbb{I}_{A(x)}}{\mathbb{P}A(x)^2} < \infty</math>).<ref>{{Cite book|last=Asmussen Soren & Glynn Peter|first=Peter|title=स्टोकेस्टिक सिमुलेशन|publisher=Springer|year=2007|isbn=978-0-387-30679-7|pages=164–167}}</ref>
 
 
===ब्लैक-बॉक्स एल्गोरिदम===
हम ब्लैक बॉक्स की संरचना को जाने बिना केवल उसके इनपुट और आउटपुट को देख सकते हैं। एल्गोरिदम को इसकी संरचना पर केवल न्यूनतम जानकारी का उपयोग करना है। जब हम यादृच्छिक संख्याएँ उत्पन्न करते हैं, तो आउटपुट नहीं हो सकता है
समान सामान्य पैरामीट्रिक वर्ग के भीतर, जैसे सामान्य या घातांकीय वितरण। ईसीएम करने के लिए स्वचालित तरीके का उपयोग किया जा सकता है। होने देना <math>X_1, X_2,...</math>आई.आई.डी. हो वितरण के साथ आर.वी <math>G</math>; सरलता के लिए हम मान लेते हैं <math>X\geq 0</math>. परिभाषित करना <math> \mathfrak{F}_n = \sigma(X_1,...,X_n,U_1,..., U_n) </math>, कहाँ <math>U_1, U_2</math>, . . . स्वतंत्र (0, 1) वर्दी हैं। के लिए एक यादृच्छिक रुकने का समय <math>X_1, X_2</math>, . . . तब रुकने का समय w.r.t. है निस्पंदन <math> \{\mathfrak{F}_n\}</math>, . . . आगे चलो <math> \mathfrak{G}</math> वितरण का एक वर्ग बनें <math>G</math> पर <math> [0, \infty)</math> साथ <math> k_G = \int_0^\infty e^{\theta x}G(dx) < \infty</math> और परिभाषित करें <math>G_\theta</math> द्वारा <math>\frac{dG_\theta}{dG(x)} = e^{\theta x - k_G}</math>. हम दिए गए के लिए ईसीएम के लिए एक ब्लैक-बॉक्स एल्गोरिदम परिभाषित करते हैं <math>\theta</math> और दी गई कक्षा  <math>\mathfrak{G}</math>यादृच्छिक रोक समय की एक जोड़ी के रूप में वितरण का <math>\tau</math> और एक <math> \mathfrak{F}_\tau- </math> मापने योग्य आर.वी. <math>Z </math> ऐसा है कि <math>Z </math> के अनुसार वितरित किया जाता है <math>G_\theta </math> किसी के लिए <math> G \in \mathfrak{G}</math>. औपचारिक रूप से, हम इसे इस प्रकार लिखते हैं <math> \mathbb{P}_G (Z<x) = G_\theta (x) </math> सभी के लिए  <math>x </math>. दूसरे शब्दों में, गेम के नियम यह हैं कि एल्गोरिदम का उपयोग किया जा सकता है
से सिम्युलेटेड मान  <math>G </math> और आर.वी. तैयार करने के लिए अतिरिक्त वर्दी। से  <math>G_\theta </math>.<ref>Asmussen, Soren & Glynn, Peter (2007). Stochastic Simulation. Springer. pp. 416–420. {{ISBN|978-0-387-30679-7}}</ref>


===सिगमंड का कलन विधि===
मान लीजिए आई.आई.डी. एक्स हल्के पुच्छल वाले वितरण के साथ <math>\mathbb{E}[X] > 0</math> है।  <math>\psi(c) = \mathbb{P}(\tau(c) < \infty)</math> का अनुमान लगाने के लिए जहां<math>\tau(c) = \inf\{t:\sum\limits_{i=1}^t X_i> c\}</math>, जब <math>c</math> बड़ा है और इसलिए <math>\psi(c)</math> छोटा है, कलन विधि महत्व वितरण प्राप्त करने के लिए घातीय आनमन का उपयोग करता है। कलन विधि का उपयोग कई पहलुओं में किया जाता है, जैसे अनुक्रमिक परीक्षण,<ref>D. Siegmund (1985) Sequential Analysis. Springer-Verlag</ref> [[जी/जी/1 कतार]] प्रतीक्षा समय, और <math>\psi</math> का उपयोग [[बर्बाद सिद्धांत|रुइन सिद्धांत]] में अंतिम रुइन की प्रायिकता के रूप में उपयोग किया जाता है। इस संदर्भ में, यह सुनिश्चित करना तर्कसंगत है कि <math>\mathbb{P}_\theta(\tau(c) < \infty) = 1</math> होगा।  मानदंड <math>\theta > \theta_0</math>, जहां <math>\theta_0</math> एस.टी. है इसीलिए <math>\kappa'(\theta_0) = 0</math> इसे प्राप्त करता है। सिगमंड का कलन विधि <math>\theta = \theta^*</math> का उपयोग करता है, यदि यह उपस्थित है, तो <math>\theta^*</math> को निम्नलिखित तरीके से परिभाषित किया गया है,
<math>\kappa(\theta^*) = 0</math>।
यह दिखाया गया है कि <math>\theta^*</math> एकमात्र आनमन वाला प्राचल है जो सीमित सापेक्ष त्रुटि (<math>\underset{x \rightarrow \infty}{\lim\sup}\frac{Var\mathbb{I}_{A(x)}}{\mathbb{P}A(x)^2} < \infty</math>) उत्पन्न करता है।<ref>{{Cite book|last=Asmussen Soren & Glynn Peter|first=Peter|title=स्टोकेस्टिक सिमुलेशन|publisher=Springer|year=2007|isbn=978-0-387-30679-7|pages=164–167}}</ref>
===ब्लैक-बॉक्स कलन विधि===
हम ब्लैक बॉक्स की संरचना को जाने बिना केवल उसके निविष्ट और निर्गत को देख सकते हैं। कलन विधि को इसकी संरचना पर केवल न्यूनतम जानकारी का उपयोग करना है। जब हम यादृच्छिक संख्याएँ उत्पन्न करते हैं, तो निर्गत सामान्य या चरघातांकी वितरण  जैसे समान सामान्य प्राचलिक वर्ग के भीतर नहीं हो सकता है। ईसीएम करने के लिए स्वचालित तरीके का उपयोग किया जा सकता है। मान लीजिये <math>X_1, X_2,...</math>आई.आई.डी. वितरण <math>G</math> के साथआर.वी हो, इसीलिए सरलता के लिए हम <math>X\geq 0</math> मानते हैं। <math> \mathfrak{F}_n = \sigma(X_1,...,X_n,U_1,..., U_n) </math> को परिभाषित करने पर, हमे<math>U_1, U_2</math>, . . . स्वतंत्र (0, 1) एकसमान प्राप्त होते हैं। <math>X_1, X_2</math>, . . . के लिए एक यादृच्छिक रुकने का समय फिर निस्पंदन <math> \{\mathfrak{F}_n\}</math>, . . .  के संबंध में एक रुकने का समय है। मान लीजिए  <math> \mathfrak{G}</math>, <math> [0, \infty)</math> के साथ वितरण <math>G</math> का एक वर्ग है और<math> k_G = \int_0^\infty e^{\theta x}G(dx) < \infty</math> को  <math>G_\theta</math> द्वारा <math>\frac{dG_\theta}{dG(x)} = e^{\theta x - k_G}</math>से परिभाषित करता है। हम दिए गए <math>\theta</math> और वितरण के दिए गए वर्ग <math>\mathfrak{G}</math> के लिए ईसीएम के लिए एक ब्लैक-बॉक्स कलन विधि को यादृच्छिक रोक समय <math>\tau</math> और <math> \mathfrak{F}_\tau- </math> मापने योग्य आर.वी. की एक जोड़ी के रूप में परिभाषित करते हैं, <math>Z </math> इस प्रकार है कि <math>Z </math> को किसी भी <math> G \in \mathfrak{G}</math>  के लिए <math>G_\theta </math> के अनुसार वितरित किया जाता है। औपचारिक रूप से, हम इसे सभी <math>x </math> के लिए <math> \mathbb{P}_G (Z<x) = G_\theta (x) </math> के रूप में लिखते हैं।


दूसरे शब्दों में, गेम के नियम यह हैं कि कलन विधि <math>G </math> से अनुकारित मानों और <math>G_\theta </math> से आर.वी. उत्पन्न करने के लिए अतिरिक्त एकसमान का उपयोग कर सकता है।<ref>Asmussen, Soren & Glynn, Peter (2007). Stochastic Simulation. Springer. pp. 416–420. {{ISBN|978-0-387-30679-7}}</ref>
== यह भी देखें ==
== यह भी देखें ==
* महत्व नमूनाकरण
* [[महत्व प्रतिदर्श]]
* अस्वीकृति नमूनाकरण
* [[अस्वीकृति प्रतिदर्श]]
*मोंटे कार्लो विधि
*[[मोंटे कार्लो विधि]]
* घातीय परिवार
* [[घातीय समूह]]
* एस्चेर परिवर्तन
* [[एस्चेर परिवर्तन]]


==संदर्भ==
==संदर्भ==
{{Reflist|30em}}
{{Reflist|30em}}
[[Category: नमूनाकरण तकनीक]]


[[Category: Machine Translated Page]]
[[Category:CS1 errors]]
[[Category:Created On 08/07/2023]]
[[Category:Created On 08/07/2023]]
[[Category:Lua-based templates]]
[[Category:Machine Translated Page]]
[[Category:Pages with script errors]]
[[Category:Templates Vigyan Ready]]
[[Category:Templates that add a tracking category]]
[[Category:Templates that generate short descriptions]]
[[Category:Templates using TemplateData]]
[[Category:नमूनाकरण तकनीक]]

Latest revision as of 19:25, 21 July 2023

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

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

चरघातांकी आनमन की प्रारंभिक औपचारिकता का श्रेय प्रायः एस्चेर को दिया जाता है[3] जबकि महत्व प्रतिदर्श में इसके उपयोग का श्रेय डेविड सिगमंड को दिया जाता है।[4]

अवलोकन

प्रायिकता वितरण , घनत्व , और आघुर्णजनक फलन (एमजीएफ) के साथ एक यादृच्छिक चर को देखते हुए, चरघातांकी रूप से आनत माप को इस प्रकार परिभाषित किया गया है,

जहां संचयी जनक फलन (सीजीएफ) है जिसे

के रूप में परिभाषित किया गया है।

हम को -का आनत घनत्व कहते हैं। यह . को संतुष्ट करता है।

एक यादृच्छिक सदिश के घातीय आनमन की एक समान परिभाषा है,

जहां दिया गया है।

उदाहरण

कई स्थितियों में चरचरघातांकी रूप से आनत माप का प्राचलिक रूप के समान होता है। एक-आयामी उदाहरणों में सामान्य वितरण, घातीय वितरण, द्विपद वितरण और पॉइसन वितरण सम्मिलित हैं।

उदाहरण के लिए, सामान्य वितरण की स्थिति में, आनत घनत्व , घनत्व है। नीचे दी गई तालिका आनत घनत्व के अधिक उदाहरण प्रदान करती है।

मूल वितरण[5][6] θ-आनत वितरण

हालाँकि, कुछ वितरणों के लिए, घातीय रूप से आनत वितरण के समान प्राचलिक समूह से संबंधित नहीं है। इसका एक उदाहरण पेरेटो वितरण है, जहां को के लिए अच्छी तरह से परिभाषित किया गया है लेकिन यह एक मानक वितरण नहीं है। ऐसे उदाहरणों में, यादृच्छिक परिवर्तनीय पीढ़ी हमेशा स्पष्ट नहीं हो सकती है।[7]

लाभ

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

इसके अलावा, मूल और आनत सीएफजी,

के बीच एक सरल संबंध उपस्थित है। इसका अवलोकन हम इस प्रकार कर सकते हैं,
इस प्रकार से,
.

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

. सरल रूप है।

गुण

  • यदि , का सीजीएफ है, तो आनत - का सीजीएफ
है। इसका मतलब यह है कि आनत का -वाँ संचयी है। विशेष रूप से, आनत वितरण की अपेक्षा है।
आनत वितरण का विचरण
. है।
  • पुनरावर्ती आनत योगात्मक है। अर्थात् पहले और फिर से आनत एक बार के आनत के समान है।
  • यदि स्वतंत्र, लेकिन आवश्यक रूप से समान यादृच्छिक चर का योग नहीं है, तो का - आनत वितरण प्रत्येक - को व्यक्तिगत रूप से आनत का योग है।
  • यदि , तो आनत वितरण और के मूल वितरण के बीच कुल्बैक-लीबलर विचलन है।
  • इसी प्रकार, के बाद से, हमारे पास,
के रूप में कुल्बैक-लीबलर विचलन है।

अनुप्रयोग

दुर्लभ-घटना अनुकरण

का घातीय आनमन, यह मानते हुए कि यह उपस्थित है, यह वितरण के एक समूह की आपूर्ति करता है जिसका उपयोग स्वीकृति-अस्वीकृति प्रतिदर्श के लिए प्रस्ताव वितरण या महत्व प्रतिदर्श के लिए महत्व वितरण के रूप में किया जा सकता है। एक सामान्य अनुप्रयोग प्रक्षेत्र के उप-क्षेत्र पर सशर्त वितरण से प्रतिदर्श, अर्थात लेना है। के उचित विकल्प के साथ, के प्रतिदर्श सार्थक रूप से प्रतिदर्श की आवश्यक मात्रा या अनुमानक के विचरण को कम कर सकता है।

सैडलबिंदु सन्निकटन

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

है।

के लिए एजवर्थ विस्तार लागू करने पर, हमा

प्राप्त कर सकते है,

जहां ,

का मानक सामान्य घनत्व है,
,

और हर्मिट बहुपद हैं।

वितरण के केंद्र से उत्तरोत्तर दूर के मूल्यों पर विचार करते समय, और पद अपरिबद्ध हो जाते हैं। हालाँकि, के प्रत्येक मान के लिए , हम को इस प्रकार चुन सकते हैं जैसे कि

के इस मान को सैडल-बिंदु के रूप में जाना जाता है, और उपरोक्त विस्तार का मूल्यांकन हमेशा आनत वितरण की अपेक्षा पर किया जाता है। का यह विकल्प

द्वारा दिए गए सन्निकटन के अंतिम प्रतिनिधित्व की ओर ले जाता है।[8][9]

अस्वीकृति प्रतिदर्श

प्रस्ताव के रूप में आनत वितरण का उपयोग करते हुए, अस्वीकृति प्रतिदर्श कलन विधि से प्रतिदर्श लेने और प्रायिकता

के साथ स्वीकार करने को निर्धारित करता है जहां
दिया गया है।

अर्थात् एक समान रूप से वितरित यादृच्छिक चर उत्पन्न होता है, और से प्रतिदर्श स्वीकार किया जाता है यदि

दिया गया हो।

महत्व प्रतिदर्श

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

प्राप्त होता है,

जहां

संभाव्यता फलन है। तो, महत्व वितरण के तहत प्रायिकता का अनुमान लगाने के लिए से एक प्रतिदर्श लें और फिर इसे संभावना अनुपात से गुणा करें। इसके अलावा, हमारे पास

द्वारा दिया गया विचरण है।

उदाहरण

स्वतंत्र और समान रूप से वितरित मान लें कि होगा। का अनुमान लगाने के लिए, हम

लेकर महत्व प्रतिदर्श का उपयोग कर सकते हैं।

किसी अन्य स्थिरांक के लिए स्थिरांक को के रूप में पुनः लिखा जा सकता है।तब,

,

जहां सैडल-बिंदु समीकरण द्वारा परिभाषित को दर्शाता है।

प्रसंभाव्य प्रक्रम

एक सामान्य आर.वी. के आनमन को देखते हुए, यह सहज है कि का घातांकीय आनमन, प्रक्षेप विचरण , प्रक्षेप के साथ एक ब्राउनियन गति है, और प्रक्षेप और विचरण के साथ एक ब्राउनियन गति है। इस प्रकार, के नीचे प्रक्षेप के साथ किसी भी ब्राउनियन गति को के नीचे प्रक्षेप के बिना ब्राउनियन गति के रूप में माना जा सकता है। इसे देखने के लिए प्रक्रिया पर विचार करें। . संभाव्यता अनुपात शब्द, , एक मार्टिंगेल है और आमतौर पर द्वारा निरूपित किया जाता है। इस प्रकार, प्रक्षेप प्रक्रिया (साथ ही ब्राउनियन निस्पंदन के लिए अनुकूलित कई अन्य निरंतर प्रक्रियाएं) के साथ एक ब्राउनियन गति एक -मार्टिंगेल है ।[10][11]

प्रसंभाव्य विभेदक समीकरण

उपरोक्त प्रसंभाव्य विभेदक समीकरण के वैकल्पिक प्रतिनिधित्व की ओर ले जाता है, जहाँ और = दिया गया है।गिरसानोव का सूत्र संभावना अनुपात बताता है। इसलिए, गिरसानोव के सूत्र का उपयोग कुछ एसडीई के महत्व प्रतिदर्श को लागू करने के लिए किया जा सकता है।

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

आनमन प्राचल का विकल्प

सिगमंड का कलन विधि

मान लीजिए आई.आई.डी. एक्स हल्के पुच्छल वाले वितरण के साथ है। का अनुमान लगाने के लिए जहां, जब बड़ा है और इसलिए छोटा है, कलन विधि महत्व वितरण प्राप्त करने के लिए घातीय आनमन का उपयोग करता है। कलन विधि का उपयोग कई पहलुओं में किया जाता है, जैसे अनुक्रमिक परीक्षण,[12] जी/जी/1 कतार प्रतीक्षा समय, और का उपयोग रुइन सिद्धांत में अंतिम रुइन की प्रायिकता के रूप में उपयोग किया जाता है। इस संदर्भ में, यह सुनिश्चित करना तर्कसंगत है कि होगा। मानदंड , जहां एस.टी. है इसीलिए इसे प्राप्त करता है। सिगमंड का कलन विधि का उपयोग करता है, यदि यह उपस्थित है, तो को निम्नलिखित तरीके से परिभाषित किया गया है,

यह दिखाया गया है कि एकमात्र आनमन वाला प्राचल है जो सीमित सापेक्ष त्रुटि () उत्पन्न करता है।[13]

ब्लैक-बॉक्स कलन विधि

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

दूसरे शब्दों में, गेम के नियम यह हैं कि कलन विधि से अनुकारित मानों और से आर.वी. उत्पन्न करने के लिए अतिरिक्त एकसमान का उपयोग कर सकता है।[14]

यह भी देखें

संदर्भ

  1. H.U. Gerber & E.S.W. Shiu (1994). "Esscher द्वारा विकल्प मूल्य निर्धारण परिवर्तन". Transactions of the Society of Actuaries. 46: 99–191.
  2. Cruz, Marcelo (2015). परिचालन जोखिम और बीमा विश्लेषण के मौलिक पहलू. Wiley. pp. 784–796. ISBN 978-1-118-11839-9.
  3. Butler, Ronald (2007). अनुप्रयोगों के साथ सैडल प्वाइंट सन्निकटन. Cambridge University Press. pp. 156. ISBN 9780521872508.
  4. Siegmund, D. (1976). "Importance Sampling in the Monte Carlo Study of Sequential Tests". The Annals of Statistics. 4 (4): 673–684. doi:10.1214/aos/1176343541.
  5. Asmussen Soren & Glynn Peter (2007). Stochastic Simulation. Springer. p. 130. ISBN 978-0-387-30679-7.
  6. Fuh, Cheng-Der; Teng, Huei-Wen; Wang, Ren-Her (2013). "Efficient Importance Sampling for Rare Event Simulation with Applications". {{cite journal}}: Cite journal requires |journal= (help)
  7. Asmussen, Soren & Glynn, Peter (2007). Stochastic Simulation. Springer. pp. 164–167. ISBN 978-0-387-30679-7
  8. Butler, Ronald (2007). अनुप्रयोगों के साथ सैडल प्वाइंट सन्निकटन. Cambridge University Press. pp. 156–157. ISBN 9780521872508.
  9. Seeber, G.U.H. (1992). जीएलआईएम और सांख्यिकीय मॉडलिंग में प्रगति. Springer. pp. 195–200. ISBN 978-0-387-97873-4.
  10. Asmussen Soren & Glynn Peter (2007). स्टोकेस्टिक सिमुलेशन. Springer. p. 407. ISBN 978-0-387-30679-7.
  11. Steele, J. Michael (2001). स्टोकेस्टिक कैलकुलस और वित्तीय अनुप्रयोग. Springer. pp. 213–229. ISBN 978-1-4419-2862-7.
  12. D. Siegmund (1985) Sequential Analysis. Springer-Verlag
  13. Asmussen Soren & Glynn Peter, Peter (2007). स्टोकेस्टिक सिमुलेशन. Springer. pp. 164–167. ISBN 978-0-387-30679-7.
  14. Asmussen, Soren & Glynn, Peter (2007). Stochastic Simulation. Springer. pp. 416–420. ISBN 978-0-387-30679-7