क्षणों की सामान्यीकृत विधि: Difference between revisions
No edit summary |
No edit summary |
||
(18 intermediate revisions by 3 users not shown) | |||
Line 1: | Line 1: | ||
{{Short description|Parameter estimation technique in statistics, particularly econometrics}}[[अर्थमिति]] और सांख्यिकी में, क्षणों की सामान्यीकृत विधि (जीएमएम) [[सांख्यिकीय मॉडल|सांख्यिकीय प्रारूपों]] में [[सांख्यिकीय पैरामीटर|मानदंडों]] का अनुमान लगाने के लिए एक सामान्य विधि है। आमतौर पर इसे [[अर्धपैरामीट्रिक मॉडल|अर्धपैरामीट्रिक प्रारूप]] के संदर्भ में लागू किया जाता है, जहां पैरामीटर मुख्यतः परिमित-आयामी होता है, जबकि डेटा के वितरण फलन का पूर्ण आकार ज्ञात नहीं हो सकता है, और इसलिए [[अधिकतम संभावना अनुमान|अधिकतम प्रायिकता अनुमान]] लागू नहीं होता है। | {{Short description|Parameter estimation technique in statistics, particularly econometrics}}[[अर्थमिति]] और सांख्यिकी में, '''क्षणों की सामान्यीकृत विधि (जीएमएम),''' [[सांख्यिकीय मॉडल|सांख्यिकीय प्रारूपों]] में [[सांख्यिकीय पैरामीटर|मानदंडों]] का अनुमान लगाने के लिए एक सामान्य विधि है। आमतौर पर इसे [[अर्धपैरामीट्रिक मॉडल|अर्धपैरामीट्रिक प्रारूप]] के संदर्भ में लागू किया जाता है, जहां पैरामीटर मुख्यतः परिमित-आयामी होता है, जबकि डेटा के वितरण फलन का पूर्ण आकार ज्ञात नहीं हो सकता है, और इसलिए [[अधिकतम संभावना अनुमान|अधिकतम प्रायिकता अनुमान]] लागू नहीं होता है। | ||
इस विधि में यह आवश्यक है कि प्रारूप के लिए निर्दिष्ट संख्या के "क्षण स्थितियों" को निर्दिष्ट किया जाए। ये क्षण स्थिति प्रारूप, पैरामीटर और डेटा के फलन होते हैं, जिनका [[प्रत्याशित मान|अपेक्षित मान]] पैरामीटर के वास्तविक मान पर शून्य होता है। जीएमएम विधि तब क्षण स्थितियों के प्रारूप औसत के किसी परिमित [[मान]] को कम करती है, और इसलिए इसे न्यूनतम-दूरी अनुमान के एक विशेष परिप्रेक्ष्य के रूप में सोचा जा सकता है।<ref>{{cite book |first=Fumio |last=Hayashi |title=अर्थमिति|publisher=Princeton University Press |year=2000 |isbn=0-691-01018-8 |page=206 |url=https://books.google.com/books?id=QyIW8WUIyzcC&pg=PA206 }}</ref> | इस विधि में यह आवश्यक है कि प्रारूप के लिए निर्दिष्ट संख्या के "क्षण स्थितियों" को निर्दिष्ट किया जाए। ये क्षण स्थिति प्रारूप, पैरामीटर और डेटा के फलन होते हैं, जिनका [[प्रत्याशित मान|अपेक्षित मान]] पैरामीटर के वास्तविक मान पर शून्य होता है। जीएमएम विधि तब क्षण स्थितियों के प्रारूप औसत के किसी परिमित [[मान]] को कम करती है, और इसलिए इसे न्यूनतम-दूरी अनुमान के एक विशेष परिप्रेक्ष्य के रूप में सोचा जा सकता है।<ref>{{cite book |first=Fumio |last=Hayashi |title=अर्थमिति|publisher=Princeton University Press |year=2000 |isbn=0-691-01018-8 |page=206 |url=https://books.google.com/books?id=QyIW8WUIyzcC&pg=PA206 }}</ref> | ||
Line 5: | Line 5: | ||
== विवरण == | == विवरण == | ||
मान लीजिए कि उपलब्ध डेटा में टी अवलोकन | मान लीजिए कि उपलब्ध डेटा में टी अवलोकन {{nowrap|{''Y<sub>t</sub>'' }<sub> ''t'' {{=}} 1,...,''T''</sub>}}, सम्मिलित है जहां प्रत्येक अवलोकन Y<sub>t</sub>एक एन-आयामी [[बहुभिन्नरूपी यादृच्छिक चर]] है। हम मान लेते हैं कि डेटा एक अपरिमित सांख्यिकीय प्रारूप से प्राप्त होता है, जिसे एक अज्ञात [[पैरामीटर]] {{nowrap|''θ'' ∈ Θ}} तक परिभाषित किया गया है। अनुमान समस्या का लक्ष्य इस पैरामीटर का "सही" मान, θ<sub>0</sub>, या कम से कम एक यथोचित निकटतम अनुमान खोजना है। | ||
जीएमएम की एक सामान्य धारणा यह है कि डेटा ''Y<sub>t</sub>'' एक [[स्थिर प्रक्रिया|कमजोर स्थिर]] [[अभ्यतिप्राय प्रक्रिया|अभ्यतिप्राय]] [[प्रसंभाव्य प्रक्रिया]] द्वारा उत्पन्न किया जाता है। [[स्वतंत्र और समान रूप से वितरित यादृच्छिक चर|स्वतंत्र और समान रूप से वितरित]] चर ''Y<sub>t</sub>'' का परिप्रेक्ष्य इस स्थिति का एक विशेष उदाहरण है। | |||
जीएमएम लागू करने के लिए, हमें क्षण | जीएमएम लागू करने के लिए, हमें क्षण स्थितियों की आवश्यकता है, अर्थात, हमें एक [[सदिश-मान फलन]] g(Y,θ) को जानना होगा जैसे कि | ||
: <math> | : <math> | ||
m(\theta_0) \equiv \operatorname{E}[\,g(Y_t,\theta_0)\,]=0, | m(\theta_0) \equiv \operatorname{E}[\,g(Y_t,\theta_0)\,]=0, | ||
</math> | </math> | ||
जहां E अपेक्षित | जहां E अपेक्षित मान दर्शाता है, और Y<sub>t</sub>एक सामान्य अवलोकन है. इसके अतिरिक्त, फ़ंक्शन m(θ) शून्य से भिन्न होना चाहिए {{nowrap|''θ'' ≠ ''θ''<sub>0</sub>}}, अन्यथा पैरामीटर θ बिंदु-[[पहचान]] नहीं होगा। | ||
जीएमएम के पीछे मूल विचार सैद्धांतिक अपेक्षित | जीएमएम के पीछे मूल विचार सैद्धांतिक अपेक्षित मान ई[⋅] को उसके अनुभवजन्य एनालॉग-प्रारूप औसत से परिवर्तित करना है: | ||
: <math> | : <math> | ||
\hat{m}(\theta) \equiv \frac{1}{T}\sum_{t=1}^T g(Y_t,\theta) | \hat{m}(\theta) \equiv \frac{1}{T}\sum_{t=1}^T g(Y_t,\theta) | ||
</math> | </math> | ||
और फिर θ के संबंध में इस | और फिर θ के संबंध में इस व्यंजक के मानदंड को कम करने के लिए। θ का न्यूनतम मान θ के लिए हमारा अनुमान <sub>0</sub> है. | ||
बड़ी संख्या के नियम | |||
''T'' के अत्यधिक बड़े मानों के लिए [[बड़ी संख्या के नियम]] के अनुसार, <math style="vertical-align:-.37em">\scriptstyle\hat{m}(\theta),\approx;\operatorname{E}[g(Y_t,\theta)],=,m(\theta)</math> होता है, और इसलिए हम अपेक्षा करते हैं कि <math style="vertical-align:-.3em">\scriptstyle\hat{m}(\theta_0);\approx;m(\theta_0);=;0</math> होगा। क्षणों की सामान्यीकृत विधि एक ऐसे संख्या <math style="vertical-align:0">\scriptstyle\hat\theta</math> की खोज करती है जो <math style="vertical-align:-.3em">\scriptstyle\hat{m}(\;\!\hat\theta\;\!)</math> को शून्य के निकट स्थापित करेगा। गणितीय रूप से, यह एक निश्चित मानक को न्यूनतम करने के बराबर है <math style="vertical-align:-.3em">\scriptstyle\hat{m}(\theta)</math> (m का मान, जिसे ||m|| के रूप में दर्शाया गया है, m और शून्य के बीच की दूरी को मापता है)। ये परिणामी अनुमानक के गुण मानक फलन की विशेष पसंद पर निर्भर होंगे, और इसलिए जीएमएम का सिद्धांत मानदंडों के एक पूरे परिवार पर विचार करता है, जिसे इस प्रकार परिभाषित किया गया है | |||
: <math> | : <math> | ||
\| \hat{m}(\theta) \|^2_{W} = \hat{m}(\theta)^{\mathsf{T}}\,W\hat{m}(\theta), | \| \hat{m}(\theta) \|^2_{W} = \hat{m}(\theta)^{\mathsf{T}}\,W\hat{m}(\theta), | ||
</math> | </math> | ||
जहां W एक [[ | जहां W एक [[धनात्मक-परिमित आव्यूह]] है | धनात्मक-परिमित भार आव्यूह, और <math>m^{\mathsf{T}}</math> स्थानान्तरण को दर्शाता है। व्यवहार में, भार आव्यूह डब्ल्यू की गणना उपलब्ध डेटा समुच्चयों के आधार पर की जाती है, जिसे इस प्रकार दर्शाया जाता है <math style="vertical-align:0em">\scriptstyle\hat{W}</math>. इस प्रकार, जीएमएम अनुमानकों को इस प्रकार लिखा जा सकता है | ||
: <math> | : <math> | ||
\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 = \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) | ||
</math> | </math> | ||
उपयुक्त परिस्थितियों में यह अनुमानक सुसंगत अनुमानक, [[स्पर्शोन्मुख सामान्यता]] और भार | उपयुक्त परिस्थितियों में यह अनुमानक सुसंगत अनुमानक, [[स्पर्शोन्मुख सामान्यता]] और भार आव्यूह के सही विकल्प के साथ <math style="vertical-align:0em">\scriptstyle\hat{W}</math> कुशल अनुमानक भी है। | ||
== गुण == | == गुण == | ||
=== | === संगतता === | ||
सुसंगत अनुमानक | सुसंगत अनुमानक किसी अनुमानक का सांख्यिकीय गुण है जिसमें कहा गया है कि, पर्याप्त संख्या में अवलोकन होने पर, अनुमानक पैरामीटर के वास्तविक मूल्य की प्रायिकता में अभिसरण करेगा: | ||
: <math>\hat\theta \xrightarrow{p} \theta_0\ \text{as}\ T\to\infty.</math> | : <math>\hat\theta \xrightarrow{p} \theta_0\ \text{as}\ T\to\infty.</math> | ||
जीएमएम अनुमानक के संगतता के लिए पर्याप्त स्थितियाँ इस प्रकार हैं: | |||
# <math>\hat{W}_T \xrightarrow{p} W,</math> जहां W एक सकारात्मक | # <math>\hat{W}_T \xrightarrow{p} W,</math> जहां ''W'' एक [[सकारात्मक प्रायिकता योग्य मैट्रिक्स|सकारात्मक समीभूत आव्यूह]] है, | ||
# <math>\,W\operatorname{E}[\,g(Y_t,\theta)\,]=0</math>केवल <math>\,\theta=\theta_0,</math> | # <math>\,W\operatorname{E}[\,g(Y_t,\theta)\,]=0</math> केवल <math>\,\theta=\theta_0,</math> | ||
# संभावित मापदंडों का स्थान | # संभावित मापदंडों का स्थान <math>\Theta \subset \mathbb{R}^{k}</math> [[ सघन स्थान ]] है। | ||
# <math>\,g(Y,\theta)</math>प्रायिकता एक के साथ प्रत्येक θ पर | # <math>\,g(Y,\theta)</math> प्रायिकता एक के साथ प्रत्येक θ पर सतत है। | ||
# <math>\operatorname{E}[\,\textstyle\sup_{\theta\in\Theta} \lVert g(Y,\theta)\rVert\,]<\infty.</math> | # <math>\operatorname{E}[\,\textstyle\sup_{\theta\in\Theta} \lVert g(Y,\theta)\rVert\,]<\infty.</math> | ||
यहां दूसरी स्थिति (तथाकथित वैश्विक पहचान स्थिति) को सत्यापित करना | यहां दूसरी स्थिति (तथाकथित वैश्विक पहचान स्थिति) को सत्यापित करना प्रायः विशेष रूप से कठिन होता है। वहाँ सरल आवश्यक परंतु पर्याप्त स्थितियाँ उपलब्ध नहीं हैं, जिनका उपयोग गैर-पहचान समस्या का पता लगाने के लिए किया जा सकता है: | ||
* | *क्रम स्थिति. क्षण फलन ''m(θ)'' का आयाम कम से कम पैरामीटर सदिश ''θ'' के आयाम जितना बड़ा होना चाहिए। | ||
* | *'''स्थानिक पहचान'''। यदि <math>\theta_0</math> के पास में ''g(Y,θ)'' सतत रूप से अविभाज्य है, तो आरेख <math>W\operatorname{E}[\nabla_\theta g(Y_t,\theta_0)]</math> का पूर्ण [[क्रमशः पंक्ति-क्रमशः|पंक्ति क्रमशः]] होना चाहिए। | ||
व्यावहारिक रूप से लागू अर्थशास्त्री | व्यावहारिक रूप से लागू अर्थशास्त्री प्रायः वास्तव में इसे सिद्ध किए बिना यह मान लेते हैं कि वैश्विक पहचान मान्य है।<ref>{{cite book |last1=Newey |first1=W. |last2=McFadden |first2=D. |year=1994 |chapter=Large sample estimation and hypothesis testing |title=अर्थमिति की पुस्तिका|volume=4 |publisher=Elsevier Science |pages=2111–2245 |doi=10.1016/S1573-4412(05)80005-4 |citeseerx=10.1.1.724.4480 |isbn=9780444887665 }}</ref>{{rp|2127}} | ||
=== | |||
=== उपगामी सामान्यता === | |||
उपगामी सामान्यता एक उपयोगी गुण है, क्योंकि यह हमें अनुमानक के लिए [[परिमित अंतराल]] बनाने और विभिन्न परीक्षण करने की अनुमति देता है। इससे पहले कि हम जीएमएम अनुमानक के उपगामी वितरण के बारे में एक बयान दे सकें, हमें दो सहायक आव्यूह को परिभाषित करने की आवश्यकता है: | |||
: <math>G = \operatorname{E}[\,\nabla_{\!\theta}\,g(Y_t,\theta_0)\,], \qquad | : <math>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}}\,]</math> | \Omega = \operatorname{E}[\,g(Y_t,\theta_0)g(Y_t,\theta_0)^{\mathsf{T}}\,]</math> | ||
फिर नीचे सूचीबद्ध | फिर नीचे सूचीबद्ध स्थितियों 1-6 के अंतर्गत, जीएमएम अनुमानक [[वितरण में अभिसरण]] के साथ असम्बद्ध रूप से सामान्य होगा: | ||
<math>\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].</math> | <math>\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].</math> | ||
स्थितियाँ: | स्थितियाँ: | ||
# <math>\hat\theta</math> सुसंगत है (पिछला अनुभाग देखें), | # <math>\hat\theta</math> सुसंगत है (पिछला अनुभाग देखें), | ||
# संभावित मापदंडों का | # संभावित मापदंडों का समुच्चय <math>\Theta \subset \mathbb{R}^{k}</math> संयोजी समुच्चय है, | ||
# <math>\,g(Y,\theta)</math> | # <math>\,g(Y,\theta)</math> अपने निकट N में <math>\theta_0</math> प्रायिकता के साथ सतत भिन्न होता है , | ||
# <math>\operatorname{E}[\,\lVert g(Y_t,\theta) \rVert^2\,]<\infty,</math> | # <math>\operatorname{E}[\,\lVert g(Y_t,\theta) \rVert^2\,]<\infty,</math> | ||
# <math>\operatorname{E}[\,\textstyle\sup_{\theta\in N}\lVert \nabla_\theta g(Y_t,\theta) \rVert\,]<\infty,</math> | # <math>\operatorname{E}[\,\textstyle\sup_{\theta\in N}\lVert \nabla_\theta g(Y_t,\theta) \rVert\,]<\infty,</math> | ||
# गणित का | # गणित का प्रश्न <math>G'WG</math> निरर्थक है। | ||
=== सापेक्ष दक्षता === | === सापेक्ष दक्षता === | ||
अब तक हमने | अब तक हमने आव्यूह <nowiki>''W''</nowiki> के चयन के बारे में कुछ नहीं कहा है, केवल इतना कहा है कि यह सकारात्मक समीभूत होना चाहिए। वास्तव में, ऐसी कोई भी आव्यूह एक सुसंगत और अनंतिक रूप से सामान्य जीएमएम अनुमापक का निर्माण करेगी, एकमात्र अंतर होगा कि उस अनुमापक की उपगामी चारधिकता में होगा। इसे निम्नलिखित रूप से सिद्ध जा सकता है | ||
: <math> W \propto\ \Omega^{-1} </math> | : <math> W \propto\ \Omega^{-1} </math> | ||
क्षण अनुमानकों की सभी | जो क्षण अनुमानकों की सभी सामान्यीकृत विधियों की श्रेणी में सबसे कुशल अनुमानक का परिणाम होगा। केवल अनंत संख्या में ऑर्थोगोनल स्थितियों से सबसे छोटा विचरण, क्रैमर-राव विचरण प्राप्त होता है। | ||
इस | इस परिप्रेक्ष में जीएमएम अनुमानक के उपगामी वितरण का सूत्र सरल हो जाता है | ||
: <math>\sqrt{T}\big(\hat\theta - \theta_0\big)\ \xrightarrow{d}\ \mathcal{N}\big[0, (G^{\mathsf{T}}\,\Omega^{-1}G)^{-1}\big]</math> | : <math>\sqrt{T}\big(\hat\theta - \theta_0\big)\ \xrightarrow{d}\ \mathcal{N}\big[0, (G^{\mathsf{T}}\,\Omega^{-1}G)^{-1}\big]</math> | ||
यह प्रमाण कि | यह प्रमाण कि भार आव्यूह का ऐसा विकल्प वास्तव में स्थानीय रूप से इष्टतम है, अन्य अनुमानकों की दक्षता स्थापित करते समय प्रायः मामूली संशोधनों के साथ अपनाया जाता है। एक नियम के रूप में, जब एक वेटिंग आव्यूह के क्रम के निकट रूप में परिवर्तित होती है, तब वह अनुकरणता की ओर आगे बढ़ती है। | ||
{| style="margin-left:30pt" | {| style="margin-left:30pt" | ||
| colspan="2" | ''''' | | colspan="2" | '''''प्रमाण'''''. हम अपरिमित चरणीय विसरण के बीच का अंतर यादृच्छिक ''W'' के साथ और <math>W=\Omega^{-1}</math> के साथ की उपगामी विसरण के बीच की तुलना करेंगे। यदि हम इस अंतर को ''CC''' के रूप में किसी आव्यूह ''C'' के साथ गुणन कर सकते हैं, तो यह यह दावा सिद्ध करेगा कि यह अंतर अविनाशीय असकारात्मक होगा, और इस प्रकार <math>W=\Omega^{-1}</math> परिभाषा के अनुसार सर्वोत्तम होगा। | ||
|- | |- | ||
| <math>\,V(W)-V(\Omega^{-1})</math> | | <math>\,V(W)-V(\Omega^{-1})</math> | ||
Line 84: | Line 118: | ||
| || <math>\,=A(I-B)A^{\mathsf{T}},</math> | | || <math>\,=A(I-B)A^{\mathsf{T}},</math> | ||
|- | |- | ||
| colspan="2" | | | colspan="2" | जहां हमने अक्षरण को थोड़ा सरल बनाने के लिए आव्यूह ''A'' और ''B'' को प्रस्तुत किया है; ''I'' एक [[identity matrix|इकाई आव्यूह]] है। हम देख सकते हैं कि यहां आव्यूह ''B'' सममित और [[Projection (linear algebra)|आइडेम्पोटेंट]] है: <math>B^2=B</math>। इसका अर्थ है कि ''I−B'' भी सममित और आइडेम्पोटेंट होगा: <math>I-B=(I-B)(I-B)^{\mathsf{T}}</math>। इस प्रकार हम पिछले व्यंजकों को गुणा करना जारी रख सकते हैं | ||
|- | |- | ||
| || <math>\,=A(I-B)(I-B)^{\mathsf{T}}A^{\mathsf{T}} = \Big(A(I-B)\Big)\Big(A(I-B)\Big)^{\mathsf{T}} \geq 0</math> | | || <math>\,=A(I-B)(I-B)^{\mathsf{T}}A^{\mathsf{T}} = \Big(A(I-B)\Big)\Big(A(I-B)\Big)^{\mathsf{T}} \geq 0</math> | ||
|} | |} | ||
== कार्यान्वयन == | == कार्यान्वयन == | ||
आरेखित विधि को लागू करने में एक कठिनाई यह है कि हम {{nowrap|''W'' {{=}} Ω<sup>−1</sup>}} नहीं ले सकते हैं क्योंकि, आव्यूह Ω की परिभाषा के अनुसार, हमें इस आव्यूह को गणना करने के लिए ''θ''<sub>0</sub> की मान को जानने की आवश्यकता होती है, और ''θ''<sub>0</sub> वही मात्रा है जिसे हम नहीं जानते हैं और पहले ही अपेक्षा कर रहे हैं कि हम उसे अनुमानित करेंगे। वाई के परिप्रेक्ष्य में <sub>t</sub> आईआईडी होने के कारण हम डब्ल्यू का अनुमान लगा सकते हैं | |||
: <math>\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}.</math> | : <math>\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}.</math> | ||
इस समस्या | इस समस्या के समाधान के लिए कई दृष्टिकोण उपलब्ध हैं, जिनमें से पहला सबसे लोकप्रिय है: | ||
{{unordered list | {{unordered list | ||
|1= ''' | |1= '''दो-चरणीय संभव जीएमएम''': | ||
{{unordered list | {{unordered list | ||
|1= '' | |1= ''चरण 1'': ''W = I'' (इकाई आव्यूह) या कोई अन्य सकारात्मक परिभाषित आव्यूह लें और प्राथमिक जीएमएम अनुमापन <math style="vertical-align:-.5em">\scriptstyle\hat\theta_{(1)}</math> की गणना करें। यह अनुमापक ''θ''<sub>0</sub> के लिए सुसंगत है, यद्यपि प्रभावी नहीं है। | ||
|2= '' | |2= ''चरण 2'': <math style="vertical-align:-.5em">\hat{W}T(\hat\theta{(1)})</math> प्रायिकता में Ω<sup>−1</sup> के प्रति संघटित होता है, और इसलिए यदि हम इस वेटिंग मैट्रिक्स के साथ <math style="vertical-align:0">\scriptstyle\hat\theta</math> की गणना करें, तो अनुमापक [[Efficiency (statistics)#Asymptotic efficiency|उपगामी रूप से प्रभावी]] होगा। | ||
}} | }} | ||
|2= ''' | |2= '''अनुक्रमित जीएमएम'''. मूलतः दो-चरणीय जीएमएम की तरहीं प्रक्रिया होती है, बस इसके बाद आव्यूह <math>\hat{W}_T</math> कई बार पुनर्गणना की जाती है। अर्थात, चरण 2 में प्राप्त अनुमापित मान का उपयोग चरण 3 के लिए भार आव्यूह की गणना के लिए किया जाता है, और इसी प्रकार चरण 4 तक जारी रखा जाता है, जब तक कि कुछ संघटन मापदंड पूरे नहीं हो जाते हैं। | ||
: <math> | : <math> | ||
\hat\theta_{(i+1)} = \operatorname{arg}\min_{\theta\in\Theta}\bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg)^{\mathsf{T}} \hat{W}_T(\hat\theta_{(i)}) \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg) | \hat\theta_{(i+1)} = \operatorname{arg}\min_{\theta\in\Theta}\bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg)^{\mathsf{T}} \hat{W}_T(\hat\theta_{(i)}) \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg) | ||
</math> | </math> | ||
ऐसे अनुक्रमितता के माध्यम से असिम्प्टोटिक रूप से कोई सुधार नहीं हो सकता है, यद्यपि कुछ मॉन्टे-कार्लो प्रयोग इस अनुमापक की सीमित-संख्याक गुणों को थोड़ा बेहतर दिखाते हैं।{{Citation needed|date=June 2009}} | |||
|3= ''' | |3= '''सतत अपडेटिंग जीएमएम''' (सीयूजीएमएम, या सीयूई). अनुमापन <math style="vertical-align:0">\scriptstyle\hat\theta</math> और भार आव्यूह ''W'' का एक साथ अनुमापन करता है: | ||
: <math> | : <math> | ||
\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}_T(\theta) \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg) | \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}_T(\theta) \bigg(\frac{1}{T}\sum_{t=1}^T g(Y_t,\theta)\bigg) | ||
</math> | </math> | ||
मॉन्टे-कार्लो प्रयोगों में, यह विधि पारंपरिक दो-चरणीय जीएमएम से बेहतर प्रदर्शन करती है: अनुमापक में छोटी माध्यिक भूल और अधिकाधिक स्थितियों में अधिक विश्वसनीय परिमिति बाधाओं के लिए जे-परीक्षण कार्यक्षमता प्रदान करती है।<ref>{{cite journal | last1 = Hansen | first1 = Lars Peter | last2 = Heaton | first2 = John | last3 = Yaron | first3 = Amir | year = 1996 | title = Finite-sample properties of some alternative GMM estimators | journal = Journal of Business & Economic Statistics | volume = 14 | issue = 3 | pages = 262–280 | jstor = 1392442 | doi=10.1080/07350015.1996.10524656 | url = http://dspace.mit.edu/bitstream/1721.1/47970/1/finitesampleprop00hans.pdf| hdl = 1721.1/47970 | hdl-access = free }}</ref> | |||
}} | }} | ||
न्यूनतमकरण प्रक्रिया के कार्यान्वयन में एक और महत्वपूर्ण | न्यूनतमकरण प्रक्रिया के कार्यान्वयन में एक और महत्वपूर्ण विषय यह है कि फलन को (संभवतः उच्च-आयामी) पैरामीटर स्पेस Θ के माध्यम से खोजना होता है और θ का मान ढूंढना होता है जो उद्देश्य फलन को न्यूनतम करता है। ऐसी प्रक्रिया के लिए कोई सामान्य अनुशंसा उपलब्ध नहीं है, यह अपने स्वयं के क्षेत्र, [[संख्यात्मक अनुकूलन]] का विषय है। | ||
== सर्गन-हैनसेन जे परीक्षण == | == सर्गन-हैनसेन जे परीक्षण == | ||
{{main article| | {{main article|सर्गन-हैनसेन जे परीक्षण}} | ||
जब क्षण स्थितियों की संख्या पैरामीटर वेक्टर θ के आयाम से अधिक होती है, तो | जब क्षण स्थितियों की संख्या पैरामीटर वेक्टर θ के आयाम से अधिक होती है, तो प्रारूप को अति-पहचानित कहा जाता है। सरगन (1958) ने वाद्य चर अनुमानकों के आधार पर अति-पहचान प्रतिबंधों के लिए परीक्षणों का प्रस्ताव रखा, जिन्हें स्वतंत्रता के क्रम के साथ ची-वर्ग चर के रूप में बड़े प्रारूपों में वितरित किया जाता है जो अति-पहचान प्रतिबंधों की संख्या पर निर्भर करते हैं। इसके बाद, हैनसेन (1982) ने इस परीक्षण को जीएमएम अनुमानकों के गणितीय समकक्ष सूत्रण पर लागू किया। यद्यपि, ध्यान दें कि ऐसे आँकड़े अनुभवजन्य अनुप्रयोगों में जहाँ प्रारूप गलत निर्दिष्ट हैं, नकारात्मक हो सकते हैं, और प्रायिकता अनुपात परीक्षण अंतर्दृष्टि प्राप्त कर सकते हैं क्योंकि प्रारूप का अनुमान शून्य और वैकल्पिक दोनों परिकल्पनाओं (भार्गव और सरगन, 1983) के अंतर्गत लगाया गया है। | ||
संकल्पनात्मक रूप से हम जाँच सकते हैं कि क्या <math>\hat{m}(\hat\theta)</math> यह सुझाव देने के लिए शून्य के पर्याप्त | संकल्पनात्मक रूप से हम जाँच सकते हैं कि क्या <math>\hat{m}(\hat\theta)</math> यह सुझाव देने के लिए शून्य के पर्याप्त निकट है कि प्रारूप डेटा को अच्छी तरह से संयोजित करता है। फिर जीएमएम विधि ने समीकरण <math>\hat{m}(\theta)=0</math> को हल करने की समस्या को, जिसमें <math>\theta</math> को पूरी तरह से प्रतिबंधों के साथ मेल खाता है, एक कम से कमीकरण के लिए गणना के द्वारा परिवर्तित कर दिया है। न्यूनतमकरण <math>m(\theta_0)=0</math> के लिए कोई ऐसा <math>\theta_0</math> उपलब्ध न होने के बाद भी, कम से कमीकरण सदैव की जा सकती है। जे-परीक्षण यही करता है। जे-परीक्षण को प्रायः ''प्रतिबंधों के लिए परीक्षण'' कहा जाता है। | ||
औपचारिक रूप से हम दो [[सांख्यिकीय परिकल्पना परीक्षण]] पर विचार करते हैं: | औपचारिक रूप से हम दो [[सांख्यिकीय परिकल्पना परीक्षण]] पर विचार करते हैं: | ||
* <math>H_0:\ m(\theta_0)=0</math>([[शून्य परिकल्पना]] | * <math>H_0:\ m(\theta_0)=0</math>([[शून्य परिकल्पना]] का प्रारूप "वैध" है), और | ||
* <math>H_1:\ m(\theta)\neq 0,\ \forall \theta\in\Theta</math>([[वैकल्पिक परिकल्पना]] | * <math>H_1:\ m(\theta)\neq 0,\ \forall \theta\in\Theta</math>([[वैकल्पिक परिकल्पना]] का प्रारूप "अमान्य" है; जो डेटा प्रतिबंधों को पूरा करने के निकट नहीं आता है) | ||
परिकल्पना के अंतर्गत <math>H_0</math>, निम्नलिखित तथाकथित जे- | परिकल्पना के अंतर्गत <math>H_0</math>, निम्नलिखित तथाकथित जे-डाटा उपगामी रूप से [[ची-वर्ग वितरण]] है | ची-वर्ग स्वतंत्रता के k-l क्रम के साथ वितरित किया जाता है। J को निम्नलिखित रूप से परिभाषित किया जा सकता है: | ||
:<math>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}</math> अंतर्गत <math>H_0,</math> | :<math>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}</math> अंतर्गत <math>H_0,</math> | ||
जहाँ <math>\hat\theta</math> पैरामीटर का जीएमएम अनुमानक <math>\theta_0</math> है, k क्षण स्थितियों की संख्या (सदिश g का आयाम) है, और l अनुमानित मापदंडों की संख्या (सदिश θ का आयाम) है। मानचित्र <math>\hat{W}_T</math> की प्रायिकता में <math>\Omega^{-1}</math> की ओर संघटित होनी चाहिए, यह प्रभावी भार आव्यूह है। ध्यान दें कि पहले हमने सिर्फ यह चयन किया था कि ''W'' <math>\Omega^{-1}</math> के आनुपातिक होना चाहिए जिससे अनुमापक प्रभावी हो; यद्यपि जे-परीक्षण को आयोजित करने के लिए ''W'' बिल्कुल <math>\Omega^{-1}</math> के बराबर होना चाहिए तथा सिर्फ़ आनुपातिक नहीं होना चाहिए। | |||
वैकल्पिक परिकल्पना के अंतर्गत <math>H_1</math>, जे- | वैकल्पिक परिकल्पना के अंतर्गत <math>H_1</math>, जे-डाटा उपगामी रूप से अपरिमित है: | ||
: <math>J\ \xrightarrow{p}\ \infty</math>अंतर्गत <math>H_1</math> | : <math>J\ \xrightarrow{p}\ \infty</math>अंतर्गत <math>H_1</math> | ||
परीक्षण करने के लिए हम डेटा से J के मान की गणना करते हैं। यह एक अऋणात्मक संख्या है. हम इसकी तुलना (उदाहरण के लिए) 0.95 [[ मात्रात्मक ]] से करते | परीक्षण करने के लिए हम डेटा से J के मान की गणना करते हैं। यह एक अऋणात्मक संख्या है. हम इसकी तुलना (उदाहरण के लिए) 0.95 [[ मात्रात्मक ]] रूप से करते हैं। | ||
<math>\chi^2_{k-\ell}</math> वितरण: | <math>\chi^2_{k-\ell}</math> वितरण: | ||
* <math>H_0</math> | * <math>H_0</math> को 95% विश्वसनीयता स्तर पर अस्वीकार कर दिया जाता है यदि <math> J > q_{0.95}^{\chi^2_{k-\ell}} </math> | ||
* <math>H_0</math> | * <math>H_0</math> को 95% विश्वसनीयता स्तर पर अस्वीकार नहीं किया जा सकता है यदि <math> J < q_{0.95}^{\chi^2_{k-\ell}} </math> | ||
== | == विस्तार == | ||
जीएमएम अनुकूलन के संदर्भ में कई अन्य लोकप्रिय अनुमान | जीएमएम अनुकूलन के संदर्भ में कई अन्य लोकप्रिय अनुमान विधियों को अपनाया जा सकता है: | ||
{{unordered list | {{unordered list | ||
|1= [[Ordinary least squares#Generalized method of moments| | |1= [[Ordinary least squares#Generalized method of moments|साधारण अत्यंत न्यून वर्ग]] (ओएलएस) क्षण स्थितियों के साथ जीएमएम के समान होती है: | ||
: <math>\operatorname{E}[ | : <math>\operatorname{E}[,x_t(y_t - x_t^{\mathsf{T}}\beta),]=0</math> | ||
|2= [[ | |2= [[भारित अत्यंत न्यून वर्ग]] (डब्ल्यूएलएस) | ||
: <math>\operatorname{E}[ | : <math>\operatorname{E}[,x_t(y_t - x_t^{\mathsf{T}}\beta)/\sigma^2(x_t),]=0</math> | ||
|3= [[ | |3= [[साधारण समकर्मी चर]] रेखांकन (आइवी) | ||
: <math>\operatorname{E}[ | : <math>\operatorname{E}[,z_t(y_t - x_t^{\mathsf{T}}\beta),]=0</math> | ||
|4= [[ | |4= [[गैर-रैखिक अत्यंत न्यून वर्ग]] (एनएलएलएस): | ||
: <math>\operatorname{E}[ | : <math>\operatorname{E}[,\nabla_{!\beta}, g(x_t,\beta)\cdot(y_t - g(x_t,\beta)),]=0</math> | ||
|5= [[ | |5= [[अधिकतम प्रायिकता]] अनुमापन (एमएलई): | ||
: <math>\operatorname{E}[ | : <math>\operatorname{E}[,\nabla_{!\theta} \ln f(x_t,\theta) ,]=0</math> | ||
}} | }} | ||
== कार्यान्वयन == | == कार्यान्वयन == | ||
*विकीबुक:आर प्रोग्रामिंग/ | *विकीबुक:आर प्रोग्रामिंग/क्षण की विधियाँ|आर प्रोग्रामिंग विकिबुक, क्षण की विधियाँ | ||
*[https://cran.r-project.org/web/packages/gmm/gmm.pdf आर] | *[https://cran.r-project.org/web/packages/gmm/gmm.pdf आर] | ||
*[https://www.stata.com/manuals15/rgmm.pdf स्टेटा] | *[https://www.stata.com/manuals15/rgmm.pdf स्टेटा] | ||
*[https://web.archive.org/web/20100327141233/http://www.eviews.com/EViews7/ev7features.html | *[https://web.archive.org/web/20100327141233/http://www.eviews.com/EViews7/ev7features.html ई व्यू] | ||
*[http://support.sas.com/documentation/cdl/en/etsug/60372/HTML/default/viewer.htm#etsug_model_sect035.htm | *[http://support.sas.com/documentation/cdl/en/etsug/60372/HTML/default/viewer.htm#etsug_model_sect035.htm एसएएस ] | ||
*[http://gretl.sourceforge.net/gretl-help/cmdref.html#gmm ग्रेटल] | *[http://gretl.sourceforge.net/gretl-help/cmdref.html#gmm ग्रेटल] | ||
== यह भी देखें == | == यह भी देखें == | ||
Line 200: | Line 289: | ||
[[Category:Articles with unsourced statements from June 2009]] | [[Category:Articles with unsourced statements from June 2009]] | ||
[[Category:Created On 06/07/2023]] | [[Category:Created On 06/07/2023]] | ||
[[Category:Lua-based templates]] | |||
[[Category:Machine Translated Page]] | [[Category:Machine Translated Page]] | ||
[[Category:Pages with script errors]] | [[Category:Pages with script errors]] | ||
Line 205: | Line 295: | ||
[[Category:Template documentation pages|Short description/doc]] | [[Category:Template documentation pages|Short description/doc]] | ||
[[Category:Templates Vigyan Ready]] | [[Category:Templates Vigyan Ready]] | ||
[[Category:Templates that add a tracking category]] | |||
[[Category:Templates that generate short descriptions]] | |||
[[Category:Templates using TemplateData]] |
Latest revision as of 21:32, 15 July 2023
अर्थमिति और सांख्यिकी में, क्षणों की सामान्यीकृत विधि (जीएमएम), सांख्यिकीय प्रारूपों में मानदंडों का अनुमान लगाने के लिए एक सामान्य विधि है। आमतौर पर इसे अर्धपैरामीट्रिक प्रारूप के संदर्भ में लागू किया जाता है, जहां पैरामीटर मुख्यतः परिमित-आयामी होता है, जबकि डेटा के वितरण फलन का पूर्ण आकार ज्ञात नहीं हो सकता है, और इसलिए अधिकतम प्रायिकता अनुमान लागू नहीं होता है।
इस विधि में यह आवश्यक है कि प्रारूप के लिए निर्दिष्ट संख्या के "क्षण स्थितियों" को निर्दिष्ट किया जाए। ये क्षण स्थिति प्रारूप, पैरामीटर और डेटा के फलन होते हैं, जिनका अपेक्षित मान पैरामीटर के वास्तविक मान पर शून्य होता है। जीएमएम विधि तब क्षण स्थितियों के प्रारूप औसत के किसी परिमित मान को कम करती है, और इसलिए इसे न्यूनतम-दूरी अनुमान के एक विशेष परिप्रेक्ष्य के रूप में सोचा जा सकता है।[1] जीएमएम अनुमानकों को सभी अनुमानकों की श्रेणी में सुसंगत अनुमानक, स्पर्शोन्मुख वितरण और सबसे कुशल अनुमानक के रूप में जाना जाता है जो क्षण स्थितियों में निहित जानकारी के अतिरिक्त किसी भी अन्य जानकारी का उपयोग नहीं करते हैं। जीएमएम को 1982 में लार्स पीटर हैंसन द्वारा प्रस्तावित किया गया था, जो क्षण विधि का विस्तार है,[2] जिसे 1894 में कार्ल पियरसन ने प्रस्तुत किया था। यद्यपि, ये अनुमापक गणनाओं के साथ "लंबकोणीयता अवस्था" (सारगान, 1958, 1959) या "अमुख्य अनुमापक समीकरण" (ह्यूबर, 1967; वांग आदि, 1997) पर आधारित उन्नत अनुमापकों के लिए गणितीय रूप से समान होते हैं।
विवरण
मान लीजिए कि उपलब्ध डेटा में टी अवलोकन {Yt } t = 1,...,T, सम्मिलित है जहां प्रत्येक अवलोकन Ytएक एन-आयामी बहुभिन्नरूपी यादृच्छिक चर है। हम मान लेते हैं कि डेटा एक अपरिमित सांख्यिकीय प्रारूप से प्राप्त होता है, जिसे एक अज्ञात पैरामीटर θ ∈ Θ तक परिभाषित किया गया है। अनुमान समस्या का लक्ष्य इस पैरामीटर का "सही" मान, θ0, या कम से कम एक यथोचित निकटतम अनुमान खोजना है।
जीएमएम की एक सामान्य धारणा यह है कि डेटा Yt एक कमजोर स्थिर अभ्यतिप्राय प्रसंभाव्य प्रक्रिया द्वारा उत्पन्न किया जाता है। स्वतंत्र और समान रूप से वितरित चर Yt का परिप्रेक्ष्य इस स्थिति का एक विशेष उदाहरण है।
जीएमएम लागू करने के लिए, हमें क्षण स्थितियों की आवश्यकता है, अर्थात, हमें एक सदिश-मान फलन g(Y,θ) को जानना होगा जैसे कि
जहां E अपेक्षित मान दर्शाता है, और Ytएक सामान्य अवलोकन है. इसके अतिरिक्त, फ़ंक्शन m(θ) शून्य से भिन्न होना चाहिए θ ≠ θ0, अन्यथा पैरामीटर θ बिंदु-पहचान नहीं होगा।
जीएमएम के पीछे मूल विचार सैद्धांतिक अपेक्षित मान ई[⋅] को उसके अनुभवजन्य एनालॉग-प्रारूप औसत से परिवर्तित करना है:
और फिर θ के संबंध में इस व्यंजक के मानदंड को कम करने के लिए। θ का न्यूनतम मान θ के लिए हमारा अनुमान 0 है.
T के अत्यधिक बड़े मानों के लिए बड़ी संख्या के नियम के अनुसार, होता है, और इसलिए हम अपेक्षा करते हैं कि होगा। क्षणों की सामान्यीकृत विधि एक ऐसे संख्या की खोज करती है जो को शून्य के निकट स्थापित करेगा। गणितीय रूप से, यह एक निश्चित मानक को न्यूनतम करने के बराबर है (m का मान, जिसे ||m|| के रूप में दर्शाया गया है, m और शून्य के बीच की दूरी को मापता है)। ये परिणामी अनुमानक के गुण मानक फलन की विशेष पसंद पर निर्भर होंगे, और इसलिए जीएमएम का सिद्धांत मानदंडों के एक पूरे परिवार पर विचार करता है, जिसे इस प्रकार परिभाषित किया गया है
जहां W एक धनात्मक-परिमित आव्यूह है | धनात्मक-परिमित भार आव्यूह, और स्थानान्तरण को दर्शाता है। व्यवहार में, भार आव्यूह डब्ल्यू की गणना उपलब्ध डेटा समुच्चयों के आधार पर की जाती है, जिसे इस प्रकार दर्शाया जाता है . इस प्रकार, जीएमएम अनुमानकों को इस प्रकार लिखा जा सकता है
उपयुक्त परिस्थितियों में यह अनुमानक सुसंगत अनुमानक, स्पर्शोन्मुख सामान्यता और भार आव्यूह के सही विकल्प के साथ कुशल अनुमानक भी है।
गुण
संगतता
सुसंगत अनुमानक किसी अनुमानक का सांख्यिकीय गुण है जिसमें कहा गया है कि, पर्याप्त संख्या में अवलोकन होने पर, अनुमानक पैरामीटर के वास्तविक मूल्य की प्रायिकता में अभिसरण करेगा:
जीएमएम अनुमानक के संगतता के लिए पर्याप्त स्थितियाँ इस प्रकार हैं:
- जहां W एक सकारात्मक समीभूत आव्यूह है,
- केवल
- संभावित मापदंडों का स्थान सघन स्थान है।
- प्रायिकता एक के साथ प्रत्येक θ पर सतत है।
यहां दूसरी स्थिति (तथाकथित वैश्विक पहचान स्थिति) को सत्यापित करना प्रायः विशेष रूप से कठिन होता है। वहाँ सरल आवश्यक परंतु पर्याप्त स्थितियाँ उपलब्ध नहीं हैं, जिनका उपयोग गैर-पहचान समस्या का पता लगाने के लिए किया जा सकता है:
- क्रम स्थिति. क्षण फलन m(θ) का आयाम कम से कम पैरामीटर सदिश θ के आयाम जितना बड़ा होना चाहिए।
- स्थानिक पहचान। यदि के पास में g(Y,θ) सतत रूप से अविभाज्य है, तो आरेख का पूर्ण पंक्ति क्रमशः होना चाहिए।
व्यावहारिक रूप से लागू अर्थशास्त्री प्रायः वास्तव में इसे सिद्ध किए बिना यह मान लेते हैं कि वैश्विक पहचान मान्य है।[3]: 2127
उपगामी सामान्यता
उपगामी सामान्यता एक उपयोगी गुण है, क्योंकि यह हमें अनुमानक के लिए परिमित अंतराल बनाने और विभिन्न परीक्षण करने की अनुमति देता है। इससे पहले कि हम जीएमएम अनुमानक के उपगामी वितरण के बारे में एक बयान दे सकें, हमें दो सहायक आव्यूह को परिभाषित करने की आवश्यकता है:
फिर नीचे सूचीबद्ध स्थितियों 1-6 के अंतर्गत, जीएमएम अनुमानक वितरण में अभिसरण के साथ असम्बद्ध रूप से सामान्य होगा:
स्थितियाँ:
- सुसंगत है (पिछला अनुभाग देखें),
- संभावित मापदंडों का समुच्चय संयोजी समुच्चय है,
- अपने निकट N में प्रायिकता के साथ सतत भिन्न होता है ,
- गणित का प्रश्न निरर्थक है।
सापेक्ष दक्षता
अब तक हमने आव्यूह ''W'' के चयन के बारे में कुछ नहीं कहा है, केवल इतना कहा है कि यह सकारात्मक समीभूत होना चाहिए। वास्तव में, ऐसी कोई भी आव्यूह एक सुसंगत और अनंतिक रूप से सामान्य जीएमएम अनुमापक का निर्माण करेगी, एकमात्र अंतर होगा कि उस अनुमापक की उपगामी चारधिकता में होगा। इसे निम्नलिखित रूप से सिद्ध जा सकता है
जो क्षण अनुमानकों की सभी सामान्यीकृत विधियों की श्रेणी में सबसे कुशल अनुमानक का परिणाम होगा। केवल अनंत संख्या में ऑर्थोगोनल स्थितियों से सबसे छोटा विचरण, क्रैमर-राव विचरण प्राप्त होता है।
इस परिप्रेक्ष में जीएमएम अनुमानक के उपगामी वितरण का सूत्र सरल हो जाता है
यह प्रमाण कि भार आव्यूह का ऐसा विकल्प वास्तव में स्थानीय रूप से इष्टतम है, अन्य अनुमानकों की दक्षता स्थापित करते समय प्रायः मामूली संशोधनों के साथ अपनाया जाता है। एक नियम के रूप में, जब एक वेटिंग आव्यूह के क्रम के निकट रूप में परिवर्तित होती है, तब वह अनुकरणता की ओर आगे बढ़ती है।
प्रमाण. हम अपरिमित चरणीय विसरण के बीच का अंतर यादृच्छिक W के साथ और के साथ की उपगामी विसरण के बीच की तुलना करेंगे। यदि हम इस अंतर को CC' के रूप में किसी आव्यूह C के साथ गुणन कर सकते हैं, तो यह यह दावा सिद्ध करेगा कि यह अंतर अविनाशीय असकारात्मक होगा, और इस प्रकार परिभाषा के अनुसार सर्वोत्तम होगा। | |
जहां हमने अक्षरण को थोड़ा सरल बनाने के लिए आव्यूह A और B को प्रस्तुत किया है; I एक इकाई आव्यूह है। हम देख सकते हैं कि यहां आव्यूह B सममित और आइडेम्पोटेंट है: । इसका अर्थ है कि I−B भी सममित और आइडेम्पोटेंट होगा: । इस प्रकार हम पिछले व्यंजकों को गुणा करना जारी रख सकते हैं | |
कार्यान्वयन
आरेखित विधि को लागू करने में एक कठिनाई यह है कि हम W = Ω−1 नहीं ले सकते हैं क्योंकि, आव्यूह Ω की परिभाषा के अनुसार, हमें इस आव्यूह को गणना करने के लिए θ0 की मान को जानने की आवश्यकता होती है, और θ0 वही मात्रा है जिसे हम नहीं जानते हैं और पहले ही अपेक्षा कर रहे हैं कि हम उसे अनुमानित करेंगे। वाई के परिप्रेक्ष्य में t आईआईडी होने के कारण हम डब्ल्यू का अनुमान लगा सकते हैं
इस समस्या के समाधान के लिए कई दृष्टिकोण उपलब्ध हैं, जिनमें से पहला सबसे लोकप्रिय है:
- दो-चरणीय संभव जीएमएम:
- चरण 1: W = I (इकाई आव्यूह) या कोई अन्य सकारात्मक परिभाषित आव्यूह लें और प्राथमिक जीएमएम अनुमापन की गणना करें। यह अनुमापक θ0 के लिए सुसंगत है, यद्यपि प्रभावी नहीं है।
- चरण 2: प्रायिकता में Ω−1 के प्रति संघटित होता है, और इसलिए यदि हम इस वेटिंग मैट्रिक्स के साथ की गणना करें, तो अनुमापक उपगामी रूप से प्रभावी होगा।
- अनुक्रमित जीएमएम. मूलतः दो-चरणीय जीएमएम की तरहीं प्रक्रिया होती है, बस इसके बाद आव्यूह कई बार पुनर्गणना की जाती है। अर्थात, चरण 2 में प्राप्त अनुमापित मान का उपयोग चरण 3 के लिए भार आव्यूह की गणना के लिए किया जाता है, और इसी प्रकार चरण 4 तक जारी रखा जाता है, जब तक कि कुछ संघटन मापदंड पूरे नहीं हो जाते हैं।
- सतत अपडेटिंग जीएमएम (सीयूजीएमएम, या सीयूई). अनुमापन और भार आव्यूह W का एक साथ अनुमापन करता है:
न्यूनतमकरण प्रक्रिया के कार्यान्वयन में एक और महत्वपूर्ण विषय यह है कि फलन को (संभवतः उच्च-आयामी) पैरामीटर स्पेस Θ के माध्यम से खोजना होता है और θ का मान ढूंढना होता है जो उद्देश्य फलन को न्यूनतम करता है। ऐसी प्रक्रिया के लिए कोई सामान्य अनुशंसा उपलब्ध नहीं है, यह अपने स्वयं के क्षेत्र, संख्यात्मक अनुकूलन का विषय है।
सर्गन-हैनसेन जे परीक्षण
जब क्षण स्थितियों की संख्या पैरामीटर वेक्टर θ के आयाम से अधिक होती है, तो प्रारूप को अति-पहचानित कहा जाता है। सरगन (1958) ने वाद्य चर अनुमानकों के आधार पर अति-पहचान प्रतिबंधों के लिए परीक्षणों का प्रस्ताव रखा, जिन्हें स्वतंत्रता के क्रम के साथ ची-वर्ग चर के रूप में बड़े प्रारूपों में वितरित किया जाता है जो अति-पहचान प्रतिबंधों की संख्या पर निर्भर करते हैं। इसके बाद, हैनसेन (1982) ने इस परीक्षण को जीएमएम अनुमानकों के गणितीय समकक्ष सूत्रण पर लागू किया। यद्यपि, ध्यान दें कि ऐसे आँकड़े अनुभवजन्य अनुप्रयोगों में जहाँ प्रारूप गलत निर्दिष्ट हैं, नकारात्मक हो सकते हैं, और प्रायिकता अनुपात परीक्षण अंतर्दृष्टि प्राप्त कर सकते हैं क्योंकि प्रारूप का अनुमान शून्य और वैकल्पिक दोनों परिकल्पनाओं (भार्गव और सरगन, 1983) के अंतर्गत लगाया गया है।
संकल्पनात्मक रूप से हम जाँच सकते हैं कि क्या यह सुझाव देने के लिए शून्य के पर्याप्त निकट है कि प्रारूप डेटा को अच्छी तरह से संयोजित करता है। फिर जीएमएम विधि ने समीकरण को हल करने की समस्या को, जिसमें को पूरी तरह से प्रतिबंधों के साथ मेल खाता है, एक कम से कमीकरण के लिए गणना के द्वारा परिवर्तित कर दिया है। न्यूनतमकरण के लिए कोई ऐसा उपलब्ध न होने के बाद भी, कम से कमीकरण सदैव की जा सकती है। जे-परीक्षण यही करता है। जे-परीक्षण को प्रायः प्रतिबंधों के लिए परीक्षण कहा जाता है।
औपचारिक रूप से हम दो सांख्यिकीय परिकल्पना परीक्षण पर विचार करते हैं:
- (शून्य परिकल्पना का प्रारूप "वैध" है), और
- (वैकल्पिक परिकल्पना का प्रारूप "अमान्य" है; जो डेटा प्रतिबंधों को पूरा करने के निकट नहीं आता है)
परिकल्पना के अंतर्गत , निम्नलिखित तथाकथित जे-डाटा उपगामी रूप से ची-वर्ग वितरण है | ची-वर्ग स्वतंत्रता के k-l क्रम के साथ वितरित किया जाता है। J को निम्नलिखित रूप से परिभाषित किया जा सकता है:
- अंतर्गत
जहाँ पैरामीटर का जीएमएम अनुमानक है, k क्षण स्थितियों की संख्या (सदिश g का आयाम) है, और l अनुमानित मापदंडों की संख्या (सदिश θ का आयाम) है। मानचित्र की प्रायिकता में की ओर संघटित होनी चाहिए, यह प्रभावी भार आव्यूह है। ध्यान दें कि पहले हमने सिर्फ यह चयन किया था कि W के आनुपातिक होना चाहिए जिससे अनुमापक प्रभावी हो; यद्यपि जे-परीक्षण को आयोजित करने के लिए W बिल्कुल के बराबर होना चाहिए तथा सिर्फ़ आनुपातिक नहीं होना चाहिए।
वैकल्पिक परिकल्पना के अंतर्गत , जे-डाटा उपगामी रूप से अपरिमित है:
- अंतर्गत
परीक्षण करने के लिए हम डेटा से J के मान की गणना करते हैं। यह एक अऋणात्मक संख्या है. हम इसकी तुलना (उदाहरण के लिए) 0.95 मात्रात्मक रूप से करते हैं।
वितरण:
- को 95% विश्वसनीयता स्तर पर अस्वीकार कर दिया जाता है यदि
- को 95% विश्वसनीयता स्तर पर अस्वीकार नहीं किया जा सकता है यदि
विस्तार
जीएमएम अनुकूलन के संदर्भ में कई अन्य लोकप्रिय अनुमान विधियों को अपनाया जा सकता है:
- साधारण अत्यंत न्यून वर्ग (ओएलएस) क्षण स्थितियों के साथ जीएमएम के समान होती है:
- भारित अत्यंत न्यून वर्ग (डब्ल्यूएलएस)
- साधारण समकर्मी चर रेखांकन (आइवी)
- गैर-रैखिक अत्यंत न्यून वर्ग (एनएलएलएस):
- अधिकतम प्रायिकता अनुमापन (एमएलई):
कार्यान्वयन
- विकीबुक:आर प्रोग्रामिंग/क्षण की विधियाँ|आर प्रोग्रामिंग विकिबुक, क्षण की विधियाँ
- आर
- स्टेटा
- ई व्यू
- एसएएस
- ग्रेटल
यह भी देखें
- अधिकतम प्रायिकता की विधि
- सामान्यीकृत अनुभवजन्य प्रायिकता
- अरेलानो-बॉन्ड अनुमानक
- अनुमानित बायेसियन गणना
संदर्भ
- ↑ Hayashi, Fumio (2000). अर्थमिति. Princeton University Press. p. 206. ISBN 0-691-01018-8.
- ↑ Hansen, Lars Peter (1982). "Large Sample Properties of Generalized Method of Moments Estimators". Econometrica. 50 (4): 1029–1054. doi:10.2307/1912775. JSTOR 1912775.
- ↑ Newey, W.; McFadden, D. (1994). "Large sample estimation and hypothesis testing". अर्थमिति की पुस्तिका. Vol. 4. Elsevier Science. pp. 2111–2245. CiteSeerX 10.1.1.724.4480. doi:10.1016/S1573-4412(05)80005-4. ISBN 9780444887665.
- ↑ Hansen, Lars Peter; Heaton, John; Yaron, Amir (1996). "Finite-sample properties of some alternative GMM estimators" (PDF). Journal of Business & Economic Statistics. 14 (3): 262–280. doi:10.1080/07350015.1996.10524656. hdl:1721.1/47970. JSTOR 1392442.
अग्रिम पठन
- 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.
- Imbens, Guido W.; Spady, Richard H.; Johnson, Phillip (1998). "Information theoretic approaches to inference in moment condition models" (PDF). Econometrica. 66 (2): 333–357. doi:10.2307/2998561. JSTOR 2998561.
- 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.
- Hayashi, Fumio (2000). Econometrics. Princeton: Princeton University Press. ISBN 0-691-01018-8.
- Hansen, Lars Peter (2002). "Method of Moments". In Smelser, N. J.; Bates, P. B. (eds.). International Encyclopedia of the Social and Behavior Sciences. Oxford: Pergamon.
- Hall, Alastair R. (2005). Generalized Method of Moments. Advanced Texts in Econometrics. Oxford University Press. ISBN 0-19-877520-2.
- Faciane, Kirby Adam Jr. (2006). Statistics for Empirical and Quantitative Finance. Statistics for Empirical and Quantitative Finance. H.C. Baird. ISBN 0-9788208-9-4.
- Special issues of Journal of Business and Economic Statistics: vol. 14, no. 3 and vol. 20, no. 4.