क्षणों की सामान्यीकृत विधि

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

इस विधि में यह आवश्यक है कि प्रारूप के लिए निर्दिष्ट संख्या के "क्षण स्थितियों" को निर्दिष्ट किया जाए। ये क्षण स्थिति प्रारूप, पैरामीटर और डेटा के फलन होते हैं, जिनका अपेक्षित मान पैरामीटर के वास्तविक मान पर शून्य होता है। जीएमएम विधि तब क्षण स्थितियों के प्रारूप औसत के किसी परिमित मान को कम करती है, और इसलिए इसे न्यूनतम-दूरी अनुमान के एक विशेष परिप्रेक्ष्य के रूप में सोचा जा सकता है। जीएमएम अनुमानकों को सभी अनुमानकों की श्रेणी में सुसंगत अनुमानक, स्पर्शोन्मुख वितरण और सबसे कुशल अनुमानक के रूप में जाना जाता है जो क्षण स्थितियों में निहित जानकारी के अतिरिक्त किसी भी अन्य जानकारी का उपयोग नहीं करते हैं। जीएमएम को 1982 में लार्स पीटर हैंसन द्वारा प्रस्तावित किया गया था, जो क्षण विधि का विस्तार है, जिसे 1894 में कार्ल पियरसन ने प्रस्तुत किया था। यद्यपि, ये अनुमापक गणनाओं के साथ "लंबकोणीयता अवस्था" (सारगान, 1958, 1959) या "अमुख्य अनुमापक समीकरण" (ह्यूबर, 1967; वांग आदि, 1997) पर आधारित उन्नत अनुमापकों के लिए गणितीय रूप से समान होते हैं।

विवरण
मान लीजिए कि उपलब्ध डेटा में टी अवलोकन {Yt }t = 1,...,T, सम्मिलित है जहां प्रत्येक अवलोकन Ytएक एन-आयामी बहुभिन्नरूपी यादृच्छिक चर है। हम मान लेते हैं कि डेटा एक अपरिमित सांख्यिकीय प्रारूप से प्राप्त होता है, जिसे एक अज्ञात पैरामीटर θ ∈ Θ तक परिभाषित किया गया है। अनुमान समस्या का लक्ष्य इस पैरामीटर का "सही" मान, θ0, या कम से कम एक यथोचित निकटतम अनुमान खोजना है।

जीएमएम की एक सामान्य धारणा यह है कि डेटा Yt एक कमजोर स्थिर अभ्यतिप्राय प्रसंभाव्य प्रक्रिया द्वारा उत्पन्न किया जाता है। स्वतंत्र और समान रूप से वितरित चर Yt का परिप्रेक्ष्य इस स्थिति का एक विशेष उदाहरण है।

जीएमएम लागू करने के लिए, हमें क्षण स्थितियों की आवश्यकता है, अर्थात, हमें एक सदिश-मान फलन g(Y,θ) को जानना होगा जैसे कि

m(\theta_0) \equiv \operatorname{E}[\,g(Y_t,\theta_0)\,]=0, $$ जहां E अपेक्षित मान दर्शाता है, और Ytएक सामान्य अवलोकन है. इसके अतिरिक्त, फ़ंक्शन m(θ) शून्य से भिन्न होना चाहिए θ ≠ θ0, अन्यथा पैरामीटर θ बिंदु-पहचान नहीं होगा।

जीएमएम के पीछे मूल विचार सैद्धांतिक अपेक्षित मान ई[⋅] को उसके अनुभवजन्य एनालॉग-प्रारूप औसत से परिवर्तित करना है:

\hat{m}(\theta) \equiv \frac{1}{T}\sum_{t=1}^T g(Y_t,\theta) $$ और फिर θ के संबंध में इस व्यंजक के मानदंड को कम करने के लिए। θ का न्यूनतम मान θ के लिए हमारा अनुमान 0 है.

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

\| \hat{m}(\theta) \|^2_{W} = \hat{m}(\theta)^{\mathsf{T}}\,W\hat{m}(\theta), $$ जहां W एक धनात्मक-परिमित आव्यूह है | धनात्मक-परिमित भार आव्यूह, और $$m^{\mathsf{T}}$$ स्थानान्तरण को दर्शाता है। व्यवहार में, भार आव्यूह डब्ल्यू की गणना उपलब्ध डेटा समुच्चयों के आधार पर की जाती है, जिसे इस प्रकार दर्शाया जाता है. इस प्रकार, जीएमएम अनुमानकों को इस प्रकार लिखा जा सकता है

\hat\theta = \operatorname{arg}\min_{\theta\in\Theta} \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg)^{\mathsf{T}} \hat{W} \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg) $$ उपयुक्त परिस्थितियों में यह अनुमानक सुसंगत अनुमानक, स्पर्शोन्मुख सामान्यता और भार आव्यूह के सही विकल्प के साथ कुशल अनुमानक भी है।

संगतता
सुसंगत अनुमानक किसी अनुमानक का सांख्यिकीय गुण है जिसमें कहा गया है कि, पर्याप्त संख्या में अवलोकन होने पर, अनुमानक पैरामीटर के वास्तविक मूल्य की प्रायिकता में अभिसरण करेगा:
 * $$\hat\theta \xrightarrow{p} \theta_0\ \text{as}\ T\to\infty.$$

जीएमएम अनुमानक के संगतता के लिए पर्याप्त स्थितियाँ इस प्रकार हैं: यहां दूसरी स्थिति (तथाकथित वैश्विक पहचान स्थिति) को सत्यापित करना प्रायः विशेष रूप से कठिन होता है। वहाँ सरल आवश्यक परंतु पर्याप्त स्थितियाँ उपलब्ध नहीं हैं, जिनका उपयोग गैर-पहचान समस्या का पता लगाने के लिए किया जा सकता है:
 * 1) $$\hat{W}_T \xrightarrow{p} W,$$ जहां W एक सकारात्मक समीभूत आव्यूह है,
 * 2) $$\,W\operatorname{E}[\,g(Y_t,\theta)\,]=0$$ केवल $$\,\theta=\theta_0,$$
 * 3) संभावित मापदंडों का स्थान $$\Theta \subset \mathbb{R}^{k}$$  सघन स्थान  है।
 * 4) $$\,g(Y,\theta)$$ प्रायिकता एक के साथ प्रत्येक θ पर सतत है।
 * 5) $$\operatorname{E}[\,\textstyle\sup_{\theta\in\Theta} \lVert g(Y,\theta)\rVert\,]<\infty.$$
 * क्रम स्थिति. क्षण फलन m(θ) का आयाम कम से कम पैरामीटर सदिश θ के आयाम जितना बड़ा होना चाहिए।
 * स्थानिक पहचान। यदि $$\theta_0$$ के पास में g(Y,θ) सतत रूप से अविभाज्य है, तो आरेख $$W\operatorname{E}[\nabla_\theta g(Y_t,\theta_0)]$$ का पूर्ण पंक्ति क्रमशः होना चाहिए।

व्यावहारिक रूप से लागू अर्थशास्त्री प्रायः वास्तव में इसे सिद्ध किए बिना यह मान लेते हैं कि वैश्विक पहचान मान्य है।

उपगामी सामान्यता
उपगामी सामान्यता एक उपयोगी गुण है, क्योंकि यह हमें अनुमानक के लिए परिमित अंतराल बनाने और विभिन्न परीक्षण करने की अनुमति देता है। इससे पहले कि हम जीएमएम अनुमानक के उपगामी वितरण के बारे में एक बयान दे सकें, हमें दो सहायक आव्यूह को परिभाषित करने की आवश्यकता है:
 * $$G = \operatorname{E}[\,\nabla_{\!\theta}\,g(Y_t,\theta_0)\,], \qquad

\Omega = \operatorname{E}[\,g(Y_t,\theta_0)g(Y_t,\theta_0)^{\mathsf{T}}\,]$$ फिर नीचे सूचीबद्ध स्थितियों 1-6 के अंतर्गत, जीएमएम अनुमानक वितरण में अभिसरण के साथ असम्बद्ध रूप से सामान्य होगा:

$$\sqrt{T}\big(\hat\theta - \theta_0\big)\ \xrightarrow{d}\ \mathcal{N}\big[0, (G^{\mathsf{T}}WG)^{-1}G^{\mathsf{T}}W\Omega W^{\mathsf{T}}G(G^{\mathsf{T}}W^{\mathsf{T}}G)^{-1}\big].$$ स्थितियाँ:
 * 1) $$\hat\theta$$ सुसंगत है (पिछला अनुभाग देखें),
 * 2) संभावित मापदंडों का समुच्चय $$\Theta \subset \mathbb{R}^{k}$$ संयोजी समुच्चय है,
 * 3) $$\,g(Y,\theta)$$ अपने निकट N में $$\theta_0$$ प्रायिकता के साथ सतत भिन्न होता है ,
 * 4) $$\operatorname{E}[\,\lVert g(Y_t,\theta) \rVert^2\,]<\infty,$$
 * 5) $$\operatorname{E}[\,\textstyle\sup_{\theta\in N}\lVert \nabla_\theta g(Y_t,\theta) \rVert\,]<\infty,$$
 * 6) गणित का प्रश्न $$G'WG$$ निरर्थक है।

सापेक्ष दक्षता
अब तक हमने आव्यूह W के चयन के बारे में कुछ नहीं कहा है, केवल इतना कहा है कि यह सकारात्मक समीभूत होना चाहिए। वास्तव में, ऐसी कोई भी आव्यूह एक सुसंगत और अनंतिक रूप से सामान्य जीएमएम अनुमापक का निर्माण करेगी, एकमात्र अंतर होगा कि उस अनुमापक की उपगामी चारधिकता में होगा। इसे निम्नलिखित रूप से सिद्ध जा सकता है
 * $$ W \propto\ \Omega^{-1} $$

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

इस परिप्रेक्ष में जीएमएम अनुमानक के उपगामी वितरण का सूत्र सरल हो जाता है
 * $$\sqrt{T}\big(\hat\theta - \theta_0\big)\ \xrightarrow{d}\ \mathcal{N}\big[0, (G^{\mathsf{T}}\,\Omega^{-1}G)^{-1}\big]$$

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

कार्यान्वयन
आरेखित विधि को लागू करने में एक कठिनाई यह है कि हम नहीं ले सकते हैं क्योंकि, आव्यूह Ω की परिभाषा के अनुसार, हमें इस आव्यूह को गणना करने के लिए θ0 की मान को जानने की आवश्यकता होती है, और θ0 वही मात्रा है जिसे हम नहीं जानते हैं और पहले ही अपेक्षा कर रहे हैं कि हम उसे अनुमानित करेंगे। वाई के परिप्रेक्ष्य में t आईआईडी होने के कारण हम डब्ल्यू का अनुमान लगा सकते हैं
 * $$\hat{W}_T(\hat\theta) = \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\hat\theta)g(Y_t,\hat\theta)^{\mathsf{T}}\bigg)^{-1}.$$

इस समस्या के समाधान के लिए कई दृष्टिकोण उपलब्ध हैं, जिनमें से पहला सबसे लोकप्रिय है:

न्यूनतमकरण प्रक्रिया के कार्यान्वयन में एक और महत्वपूर्ण विषय यह है कि फलन को (संभवतः उच्च-आयामी) पैरामीटर स्पेस Θ के माध्यम से खोजना होता है और θ का मान ढूंढना होता है जो उद्देश्य फलन को न्यूनतम करता है। ऐसी प्रक्रिया के लिए कोई सामान्य अनुशंसा उपलब्ध नहीं है, यह अपने स्वयं के क्षेत्र, संख्यात्मक अनुकूलन का विषय है।

सर्गन-हैनसेन जे परीक्षण
जब क्षण स्थितियों की संख्या पैरामीटर वेक्टर θ के आयाम से अधिक होती है, तो मॉडल को अति-पहचानित कहा जाता है। सरगन (1958) ने वाद्य चर अनुमानकों के आधार पर अति-पहचान प्रतिबंधों के लिए परीक्षणों का प्रस्ताव रखा, जिन्हें स्वतंत्रता की डिग्री के साथ ची-स्क्वायर चर के रूप में बड़े नमूनों में वितरित किया जाता है जो अति-पहचान प्रतिबंधों की संख्या पर निर्भर करते हैं। इसके बाद, हैनसेन (1982) ने इस परीक्षण को जीएमएम अनुमानकों के गणितीय समकक्ष फॉर्मूलेशन पर लागू किया। हालाँकि, ध्यान दें कि ऐसे आँकड़े अनुभवजन्य अनुप्रयोगों में नकारात्मक हो सकते हैं जहाँ मॉडल गलत निर्दिष्ट हैं, और प्रायिकता अनुपात परीक्षण अंतर्दृष्टि प्राप्त कर सकते हैं क्योंकि मॉडल का अनुमान शून्य और वैकल्पिक दोनों परिकल्पनाओं (भार्गव और सरगन, 1983) के तहत लगाया गया है।

संकल्पनात्मक रूप से हम जाँच सकते हैं कि क्या $$\hat{m}(\hat\theta)$$ यह सुझाव देने के लिए शून्य के पर्याप्त करीब है कि मॉडल डेटा को अच्छी तरह से फिट बैठता है। जीएमएम विधि ने समीकरण को हल करने की समस्या का स्थान ले लिया है $$\hat{m}(\theta)=0$$, जो चुनता है $$\theta$$ न्यूनतमकरण गणना द्वारा, प्रतिबंधों का सटीक मिलान करना। न्यूनतमकरण हमेशा नहीं होने पर भी किया जा सकता है $$\theta_0$$ ऐसा मौजूद है $$m(\theta_0)=0$$. जे-टेस्ट यही करता है। जे-परीक्षण को प्रतिबंधों की अधिक पहचान के लिए परीक्षण भी कहा जाता है।

औपचारिक रूप से हम दो सांख्यिकीय परिकल्पना परीक्षण पर विचार करते हैं:
 * $$H_0:\ m(\theta_0)=0$$(शून्य परिकल्पना कि मॉडल "वैध" है), और
 * $$H_1:\ m(\theta)\neq 0,\ \forall \theta\in\Theta$$(वैकल्पिक परिकल्पना कि मॉडल "अमान्य" है; डेटा प्रतिबंधों को पूरा करने के करीब नहीं आता है)

परिकल्पना के अंतर्गत $$H_0$$, निम्नलिखित तथाकथित जे-आँकड़ा स्पर्शोन्मुख रूप से ची-वर्ग वितरण है | ची-वर्ग स्वतंत्रता की k-l डिग्री के साथ वितरित किया जाता है। J को परिभाषित करें:
 * $$J \equiv T \cdot \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\hat\theta)\bigg)^{\mathsf{T}} \hat{W}_T \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\hat\theta)\bigg)\ \xrightarrow{d}\ \chi^2_{k-\ell}$$ अंतर्गत $$H_0,$$

कहाँ $$\hat\theta$$ पैरामीटर का GMM अनुमानक है $$\theta_0$$, k क्षण स्थितियों की संख्या (वेक्टर g का आयाम) है, और l अनुमानित मापदंडों की संख्या (वेक्टर θ का आयाम) है। आव्यूह $$\hat{W}_T$$ की संभाव्यता में अभिसरण होना चाहिए $$\Omega^{-1}$$, कुशल भार आव्यूह (ध्यान दें कि पहले हमें केवल यह आवश्यक था कि W के समानुपाती हो $$\Omega^{-1}$$ अनुमानक के कुशल होने के लिए; हालाँकि, J-परीक्षण करने के लिए W के बिल्कुल बराबर होना चाहिए $$\Omega^{-1}$$, केवल आनुपातिक नहीं)।

वैकल्पिक परिकल्पना के अंतर्गत $$H_1$$, जे-आँकड़ा स्पर्शोन्मुख रूप से असीमित है:
 * $$J\ \xrightarrow{p}\ \infty$$अंतर्गत $$H_1$$

परीक्षण करने के लिए हम डेटा से J के मान की गणना करते हैं। यह एक अऋणात्मक संख्या है. हम इसकी तुलना (उदाहरण के लिए) 0.95 मात्रात्मक  से करते हैं $$\chi^2_{k-\ell}$$ वितरण:
 * $$H_0$$ यदि 95% विश्वास स्तर पर खारिज कर दिया जाता है $$ J > q_{0.95}^{\chi^2_{k-\ell}} $$
 * $$H_0$$ यदि 95% विश्वास स्तर पर अस्वीकार नहीं किया जा सकता है $$ J < q_{0.95}^{\chi^2_{k-\ell}} $$

दायरा
जीएमएम अनुकूलन के संदर्भ में कई अन्य लोकप्रिय अनुमान तकनीकों को अपनाया जा सकता है:

कार्यान्वयन

 * विकीबुक:आर प्रोग्रामिंग/मेथड ऑफ मोमेंट्स|आर प्रोग्रामिंग विकिबुक, मेथड ऑफ मोमेंट्स
 * आर
 * स्टेटा
 * EViews
 * SAS
 * ग्रेटल

यह भी देखें

 * अधिकतम प्रायिकता की विधि
 * सामान्यीकृत अनुभवजन्य प्रायिकता
 * अरेलानो-बॉन्ड अनुमानक
 * अनुमानित बायेसियन गणना

अग्रिम पठन

 * Huber, P. (1967). The behavior of maximum likelihood estimates under nonstandard conditions. Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics and Probability 1, 221-233.


 * Newey W., McFadden D. (1994). Large sample estimation and hypothesis testing, in Handbook of Econometrics, Ch.36. Elsevier Science.




 * Sargan, J.D. (1958). The estimation of economic relationships using instrumental variables. Econometrica, 26, 393-415.


 * Sargan, J.D. (1959). The estimation of relationships with autocorrelated residuals by the use on instrumental variables. Journal of the Royal Statistical Society B, 21, 91-105.


 * Wang, C.Y., Wang, S., and Carroll, R. (1997). Estimation in choice-based sampling with measurement error and bootstrap analysis. Journal of Econometrics, 77, 65-86.


 * Bhargava, A., and Sargan, J.D. (1983). Estimating dynamic random effects from panel data covering short time periods. Econometrica, 51, 6, 1635-1659.


 * Special issues of Journal of Business and Economic Statistics: vol. 14, no. 3 and vol. 20, no. 4.
 * Special issues of Journal of Business and Economic Statistics: vol. 14, no. 3 and vol. 20, no. 4.
 * Special issues of Journal of Business and Economic Statistics: vol. 14, no. 3 and vol. 20, no. 4.
 * Special issues of Journal of Business and Economic Statistics: vol. 14, no. 3 and vol. 20, no. 4.
 * Special issues of Journal of Business and Economic Statistics: vol. 14, no. 3 and vol. 20, no. 4.


 * Short Introduction to the Generalized Method of Moments