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

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

व्युत्क्रम परिवर्तन प्रतिचयन एक संख्या $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}$ $X$ की संयोजी प्रसरण फलन $F_{X}$ का सामान्यीकृत व्युत्क्रम है और $U$ $[0,1]$ पर समरूप है।

वास्तविकता विवरण के लिए, निर्दिष्ट प्रायिकता अंकित्र का विपरीत अनुरूप सांख्यिकीय ज्ञानकोश यह स्थापित करता है कि प्रायिकता वितरण के अनुरूप एक निर्दिष्ट प्रायिकता अंकित्र $X$ के लिए, जिसका संयोजी वितरण फलन $F_{X}$ है, यादृच्छिक चर $U=F_{X}(X)$ समरूप $[0,1]$ पर होता है।

$x$ से $F(x)$ तक के लिए व्युत्क्रम की तकनीक का आरेख। नीचे दाएं कोने में हम सतत फलन देखते हैं और ऊपर बाएं कोने में इसका यूतक्रम।

अंतर्बोध

यदि $U\sim \mathrm {Unif} [0,1]$ है, तो हम $CDF$ $F_{X}(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} ,$ यहां अंतिम चरण में उपयोग किया गया कि जब $U$ $[0,1]$ पर समरूप होता है, तो $\Pr(U\leq y)=y$ होता है।

हमने $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 \ toMathrm{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)$.

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

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

यह भी देखें

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