व्युत्क्रम परिवर्तन नमूनाकरण

व्युत्क्रम परिवर्तन प्रतिचयन जिसे व्युत्क्रम प्रतिचयन, व्युत्क्रम संभाव्यता अभिन्न परिवर्तन, व्युत्क्रम परिवर्तन विधि, निकोलाई स्मिरनोव परिवर्तन, या स्वर्ण नियम के नाम से भी जाना जाता है एक मूलभूत विधि है जो छद्म-यादृच्छिक संख्या प्रतिचयन के लिए काम में लाई जाती है, अर्थात किसी भी प्रायिकता वितरण से उसकी करगणना संचार फलन देते हुए किसी भी यादृच्छिक प्रतिचयन संख्याओं को उत्पन्न करने के लिए इसका प्रयोग किया जाता है।

व्युत्क्रम परिवर्तन प्रतिचयन एक संख्या $u$ के 0 और 1 के बीच के संख्या प्रतिचयन (जिनकी व्याख्या प्रायिकता के रूप में की जाती हैं) का उपयोग करता है, और पुनः एक ऐसी सबसे छोटी संख्या $x\in \mathbb {R}$ देता है जिसके लिए $F(x)\geq u$ होता है, यहां $F$ एक यादृच्छिक चर के लिए संयोजी वितरण फलन है। उदाहरण के लिए, मान लीजिए कि $F$ माध्य शून्य और मानक प्रसरण एक मानक सामान्य वितरण है। नीचे दी गई तालिका समान वितरण से लिए गए प्रारूप तथा मानक सामान्य वितरण पर उनका प्रतिनिधित्व दर्शाती है।

समरूप प्रतिचयन से सामान्य प्रतिचयन में परिवर्तन

$u$	$F^{-1}(u)$
.5	0
.975	1.95996
.995	2.5758
.999999	4.75342
1-2−52	8.12589

सामान्य वितरण के लिए व्युत्क्रम परिवर्तन प्रतिचयन

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

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

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

औपचारिक कथन

किसी भी यादृच्छिक चर के लिए $X\in \mathbb {R}$, यादृच्छिक चर $F_{X}^{-1}(U)$ एक ही कानून है जैसा $X$, कहाँ $F_{X}^{-1}$ संचयी वितरण फलन का संचयी वितरण फलन # व्युत्क्रम_वितरण_फलन_(विभाजक_फलन) है $F_{X}$ का $X$ और $U$ पर एक समान है $[0,1]$. रैंडम_वेरिएबल#कंटीन्युअस_रैंडम_वेरिएबल के लिए, व्युत्क्रम संभाव्यता इंटीग्रल ट्रांसफॉर्म वास्तव में संभाव्यता अभिन्न परिवर्तन का व्युत्क्रम है, जो बताता है कि एक सतत यादृच्छिक चर के लिए $X$ संचयी वितरण फलन के साथ $F_{X}$, यादृच्छिक चर $U=F_{X}(X)$ पर एकसमान वितरण (निरंतर) है $[0,1]$.

उलटा तकनीक का ग्राफ़ से $x$ को $F(x)$. नीचे दाईं ओर हम नियमित फलन देखते हैं और ऊपर बाईं ओर इसका उलटा देखते हैं।

अंतर्ज्ञान

से $U\sim \mathrm {Unif} [0,1]$, हम उत्पन्न करना चाहते हैं $X$ संचयी वितरण फलन के साथ $F_{X}(x).$ हम यह मानते है कि $F_{X}(x)$ एक सतत, सख्ती से बढ़ता कार्य होना, जो अच्छा अंतर्ज्ञान प्रदान करता है।

हम यह देखना चाहते हैं कि क्या हम कुछ सख्ती से नीरस परिवर्तन पा सकते हैं $T:[0,1]\mapsto \mathbb {R}$, ऐसा है कि $T(U){\overset {d}{=}}X$. हमारे पास होगा

$F_{X}(x)=\Pr(X\leq x)=\Pr(T(U)\leq x)=\Pr(U\leq T^{-1}(x))=T^{-1}(x),{\text{ for }}x\in \mathbb {R} ,$ जहां अंतिम चरण में उसका उपयोग किया गया $\Pr(U\leq y)=y$ कब $U$ पर एक समान है $[0,1]$.

तो हमें मिल गया $F_{X}$ का व्युत्क्रम कार्य होना $T$, या, समकक्ष $T(u)=F_{X}^{-1}(u),u\in [0,1].$ इसलिए, हम उत्पन्न कर सकते हैं $X$ से $F_{X}^{-1}(U).$

विधि

व्युत्क्रम परिवर्तन नमूने का योजनाबद्ध। का उलटा कार्य $y=F_{X}(x)$ द्वारा परिभाषित किया जा सकता है $F_{X}^{-1}(y)=\mathrm {inf} \{x|F_{X}(x)\geq y\}$.

व्युत्क्रम परिवर्तन प्रतिचयन कैसे समान रूप से वितरित यादृच्छिक मानों से सामान्य रूप से वितरित यादृच्छिक मान उत्पन्न करता है, इसका एक एनीमेशन

व्युत्क्रम परिवर्तन प्रतिचयन विधि द्वारा हल की जाने वाली समस्या इस प्रकार है:

• होने देना $X$ एक यादृच्छिक चर हो जिसका वितरण संचयी वितरण फलन द्वारा वर्णित किया जा सकता है $F_{X}$.
• हम मूल्य उत्पन्न करना चाहते हैं $X$ जो इस वितरण के अनुसार वितरित किये जाते हैं।

व्युत्क्रम परिवर्तन प्रतिचयन विधि निम्नानुसार काम करती है:

1. छद्म यादृच्छिक संख्या जनरेटर $u$ अंतराल में मानक समान वितरण से $[0,1]$, यानी से $U\sim \mathrm {Unif} [0,1].$
2. वांछित सीडीएफ का संचयी वितरण फलन # व्युत्क्रम_वितरण_फलन_(विभाजक_फलन) ढूंढें, यानी। $F_{X}^{-1}(u)$.
3. गणना करें $X'(u)=F_{X}^{-1}(u)$. परिकलित यादृच्छिक चर $X'(U)$ वितरण है $F_{X}$ और इस प्रकार वही कानून है $X$.

संचयी वितरण फलन को देखते हुए, अलग-अलग तरीके से व्यक्त किया गया $F_{X}$ और एक समान चर $U\in [0,1]$, यादृच्छिक चर $X=F_{X}^{-1}(U)$ वितरण है $F_{X}$.

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

उदाहरण

• उदाहरण के तौर पर, मान लीजिए कि हमारे पास एक यादृच्छिक चर है $U\sim \mathrm {Unif} (0,1)$ और एक संचयी वितरण फलन

${\begin{aligned}F(x)=1-\exp(-{\sqrt {x}})\end{aligned}}$

एक व्युत्क्रम निष्पादित करने के लिए जिसे हम हल करना चाहते हैं $F(F^{-1}(u))=u$

${\begin{aligned}F(F^{-1}(u))&=u\\1-\exp \left(-{\sqrt {F^{-1}(u)}}\right)&=u\\F^{-1}(u)&=(-\log(1-u))^{2}\\&=(\log(1-u))^{2}\end{aligned}}$

यहां से हम चरण एक, दो और तीन करेंगे।

• एक अन्य उदाहरण के रूप में, हम घातीय वितरण का उपयोग करते हैं $F_{X}(x)=1-e^{-\lambda x}$ x ≥ 0 के लिए (और अन्यथा 0)। y=F(x) को हल करके हम व्युत्क्रम फलन प्राप्त करते हैं

$x=F^{-1}(y)=-{\frac {1}{\lambda }}\ln(1-y).$

इसका मतलब यह है कि यदि हम कुछ चित्र बनाते हैं $y_{0}$एक से $U\sim \mathrm {Unif} (0,1)$ और गणना करें $x_{0}=F_{X}^{-1}(y_{0})=-{\frac {1}{\lambda }}\ln(1-y_{0}),$ यह $x_{0}$ घातीय वितरण है.

यह विचार निम्नलिखित ग्राफ़ में दर्शाया गया है:

यादृच्छिक संख्या y_i 0 और 1, यानी Y ~ U(0, 1) के बीच एक समान वितरण से उत्पन्न होते हैं। उन्हें y-अक्ष पर रंगीन बिंदुओं के रूप में चित्रित किया गया है। प्रत्येक बिंदु को x=F के अनुसार मैप किया गया है⁻¹(y), जिसे दो उदाहरण बिंदुओं के लिए ग्रे तीरों के साथ दिखाया गया है। इस उदाहरण में, हमने एक घातीय वितरण का उपयोग किया है। इसलिए, x ≥ 0 के लिए, संभाव्यता घनत्व है $\varrho _{X}(x)=\lambda e^{-\lambda \,x}$ और संचयी वितरण फलन है $F(x)=1-e^{-\lambda \,x}$. इसलिए, $x=F^{-1}(y)=-{\frac {\ln(1-y)}{\lambda }}$. हम देख सकते हैं कि इस पद्धति का उपयोग करते हुए, कई बिंदु 0 के करीब पहुंच जाते हैं और केवल कुछ बिंदुओं पर उच्च x-मान होते हैं - जैसा कि एक घातांकीय वितरण के लिए अपेक्षित है।

: ध्यान दें कि यदि हम y के बजाय 1-y से शुरू करते हैं तो वितरण नहीं बदलता है। संगणनीय उद्देश्यों के लिए, इसलिए [0, 1] में यादृच्छिक संख्या y उत्पन्न करना और फिर बस गणना करना पर्याप्त है

$x=F^{-1}(y)=-{\frac {1}{\lambda }}\ln(y).$

सटीकता का प्रमाण

होने देना $F$ एक संचयी वितरण फलन बनें, और चलो $F^{-1}$ इसका संचयी वितरण फलन हो#इनवर्स_डिस्ट्रिब्यूशन_फंक्शन_(विभाजक_फंक्शन) (न्यूनतम का उपयोग करते हुए क्योंकि सीडीएफ कमजोर रूप से मोनोटोनिक और कैडलैग|राइट-कंटीन्यूअस हैं):

$F^{-1}(u)=\inf \;\{x\mid F(x)\geq u\}\qquad (0<u<1).$

दावा: यदि $U$ एक समान वितरण (निरंतर) यादृच्छिक चर है $[0,1]$ तब $F^{-1}(U)$ है $F$ इसके सीडीएफ के रूप में।

सबूत:

${\begin{aligned}&\Pr(F^{-1}(U)\leq x)\\&{}=\Pr(U\leq F(x))\quad &(F{\text{ is right-continuous, so }}\{u:F^{-1}(u)\leq x\}=\{u:u\leq F(x)\})\\&{}=F(x)\quad &({\text{because }}\Pr(U\leq u)=u,{\text{ when U is uniform on}}[0,1])\\\end{aligned}}$

कटा हुआ वितरण

व्युत्क्रम परिवर्तन प्रतिचयन को अंतराल पर काटे गए वितरण के मामलों तक आसानी से बढ़ाया जा सकता है $(a,b]$ अस्वीकृति नमूने की लागत के बिना: समान विधिकलन का पालन किया जा सकता है, लेकिन एक यादृच्छिक संख्या उत्पन्न करने के बजाय $u$ 0 और 1 के बीच समान रूप से वितरित, उत्पन्न करें $u$ के बीच समान रूप से वितरित $F(a)$ और $F(b)$, और फिर दोबारा ले लो $F^{-1}(u)$.

व्युत्क्रमों की संख्या में कमी

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

यह भी देखें

• संभाव्यता अभिन्न परिवर्तन
• कोपुला (सांख्यिकी), संभाव्यता अभिन्न परिवर्तन के माध्यम से परिभाषित।
• व्युत्क्रम सीडीएफ के स्पष्ट निर्माण के लिए विभाजक फलन।
• असतत घटकों के साथ वितरण के लिए सटीक गणितीय परिभाषा के लिए संचयी वितरण फलन # व्युत्क्रम।