आबादी से नमूना करते समय UMV का


10

Let घनत्व f _ {\ theta} (x) = \ the थीटा x ^ {\ theta-1} \ mathbf1_ {0 <<1} \ _ \ _ से एक यादृच्छिक नमूना हो। \ थीटा> 0(X1,X2,,Xn)

fθ(x)=θxθ110<x<1,θ>0

मैं θ1+θ के UMVUE को खोजने का प्रयास कर रहा हूं ।

(X1,,Xn) का संयुक्त घनत्व है

fθ(x1,,xn)=θn(i=1nxi)θ110<x1,,xn<1=exp[(θ1)i=1nlnxi+nlnθ+ln(10<x1,,xn<1)],θ>0

जनसंख्या पीडीएफ के रूप में fθ एक पैरामीटर घातीय परिवार से है, यह दिखाता है कि के लिए एक पूर्ण पर्याप्त आंकड़ा θ है

T(X1,,Xn)=i=1nlnXi

चूंकि E(X1)=θ1+θ , पहले इस विचार से, E(X1T) की UMVUE मुझे देना होगा θ1+θ द्वारा लेहमैन-शेफ़ी प्रमेय। निश्चित नहीं है कि यह सशर्त अपेक्षा सीधे मिल सकती है या किसी को सशर्त वितरण X1i=1nlnXi

दूसरी ओर, मैंने निम्नलिखित दृष्टिकोण पर विचार किया:

हमारे पास Xii.i.dBeta(θ,1)2θlnXii.i.dχ22 , ताकि 2θTχ2n2

इसलिए th ऑर्डर कच्चे पल के बारे में शून्य है, जैसा कि ची-स्क्वायर पीडीएफ का उपयोग करके गणना की गई हैr2θT

E(2θT)r=2rΓ(n+r)Γ(n),n+r>0

तो ऐसा लगता है कि विभिन्न पूर्णांक विकल्पों के लिए , मुझे के विभिन्न पूर्णांक शक्तियों के निष्पक्ष अनुमानक (और UMVUE) मिलेंगे । उदाहरण के लिए, और सीधे मुझे UMVUE के और क्रमशः दें।rθE(Tn)=1θE(1nT)=θ1θθ

अब, जब हमारे पास ।θ>1θ1+θ=(1+1θ)1=11θ+1θ21θ3+

मुझे निश्चित रूप से UMVUE की _ फ़्रेक इत्यादि मिल सकते हैं। इसलिए इन UMVUE के I को संयोजित करने पर आपको आवश्यक _ _ UMVUE मिल सकता है । क्या यह विधि मान्य है या मुझे पहली विधि के साथ आगे बढ़ना चाहिए? जब यह मौजूद है तो UMVUE अद्वितीय है, दोनों को मुझे एक ही उत्तर देना चाहिए।1θ,1θ2,1θ3θ1+θ

स्पष्ट होने के लिए, मुझे

E(1+Tn+T2n(n+1)+T3n(n+1)(n+2)+)=11θ+1θ21θ3+

अर्थात,

E(r=0Trn(n+1)...(n+r1))=θ1+θ

क्या यह संभव है कि मेरी अपेक्षित UMVUE जब ?r=0Trn(n+1)...(n+r1)θ>1

के लिए , मैं मिलेगा , और UMVUE अलग होगा तो।0<θ<1g(θ)=θ(1+θ+θ2+)


यह मानते हुए कि पहले दृष्टिकोण में सशर्त अपेक्षा सीधे नहीं पाई जा सकती थी, और चूंकि , मैं आगे बढ़ चुका था। सशर्त वितरण को खोजने के लिए । उसके लिए, मुझे के संयुक्त घनत्व की आवश्यकता थी ।E(X1lnXi=t)=E(X1Xi=et)X1Xi(X1,Xi)

मैंने चर जैसे कि for all । के संयुक्त समर्थन करने के लिए यह नेतृत्व किया जा रहा है ।(X1,,Xn)(Y1,,Yn)Yi=j=1iXji=1,2,,n(Y1,,Yn)S={(y1,,yn):0<y1<1,0<yj<yj1 for j=2,3,,n}

जैकोबियन निर्धारक ।J=(i=1n1yi)1

इसलिए मुझे रूप में का संयुक्त घनत्व मिला(Y1,,Yn)

fY(y1,y2,,yn)=θnynθ1i=1n1yi1S

का संयुक्त घनत्व इसलिए(Y1,Yn)

fY1,Yn(y1,yn)=θnynθ1y10yn20yn30y11y3y4...yn1dy2y2dyn2dyn1

क्या मैं यहां एक अलग परिवर्तन कर सकता हूं जो संयुक्त घनत्व की व्युत्पत्ति को कम बोझिल बना देगा? मुझे यकीन नहीं है कि मैंने यहां सही परिवर्तन लिया है।


टिप्पणी अनुभाग में कुछ उत्कृष्ट सुझावों के आधार पर, मैंने संयुक्त घनत्व के स्थान पर का संयुक्त घनत्व जहां और ।(U,U+V)(X1,Xi)U=lnX1V=i=2nlnXi

यह तुरंत देखा जाता है कि और स्वतंत्र हैं।UExp(θ)VGamma(n1,θ)

और वास्तव में, ।U+VGamma(n,θ)

के लिए , के संयुक्त घनत्व हैn>1(U,V)

fU,V(u,v)=θeθu1u>0θn1Γ(n1)eθvvn21v>0

चर बदलने, मैं के संयुक्त घनत्व मिल गया के रूप में(U,U+V)

fU,U+V(u,z)=θnΓ(n1)eθz(zu)n210<u<z

तो, का सशर्त घनत्वUU+V=z

fUU+V(uz)=(n1)(zu)n2zn110<u<z

अब, मेरा UMVUE बिल्कुल , जैसा कि मैंने सही उल्लेख किया था इस पोस्ट की शुरुआत में।E(eUU+V=z)=E(X1i=1nlnXi=z)

इसलिए सभी को करने के लिए

E(eUU+V=z)=n1zn10zeu(zu)n2du

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


आपको के सशर्त वितरण को खोजने के लिए पहली विधि के साथ आगे बढ़ना चाहिए इस एप्लिकेशन के साथ काम करने के लिए पर्याप्त आँकड़ा का कौन सा रूप , को आसान बनाता है। X[1]|Xi
जूलमैन

1
उस बिंदु पर (शुरुआत में) जहाँ आप परिचय देते हैं, आपको चर संदर्भ में काम करने का मन होना चाहिएयह लगभग तत्काल है कि वे वितरणों के समानुपाती हैं , जो कि और के संयुक्त वितरण पर विचार करने में आपकी समस्या को कम करता है। यह गणित के शेष दो पृष्ठों को शार्टकट करेगा और आपको समाधान का त्वरित मार्ग देगा। TYi=logXi.Γ(1)(U,U+V)UΓ(1)VΓ(n1).
whuber

स्पष्ट होने के लिए, क्या आप सुझाव दे रहे हैं कि मैं का घनत्व पहले और घनत्व का पता ? मैंने देखा था कि की घातीय चर दर के साथ (जो कि आपके कहे अनुसार एक गामा चर भी है), लेकिन उसके साथ काम करने के बारे में नहीं सोचा था। ( एक्स 1 , Π एक्स मैं ) - ln एक्स मैं θ(lnX1,lnX1i=2nlnXi)(X1,Xi)lnXiθ
स्टबबोर्नटॉम

1
@whuber लेकिन सीधे से मुझे कैसे मिलेगा ? ( ln एक्स 1 | )E(X1...)E(lnX1...)
StubbornAtom

1
@whuber कृपया मेरे संपादन पर एक नज़र है। मैंने इसे लगभग पूरा कर लिया है, लेकिन निश्चित नहीं है कि उस अभिन्न के साथ क्या किया जाए। मुझे पूरा विश्वास है कि मेरी गणना सही है।
स्टबबोर्नटॉम

जवाबों:


5

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

यहाँ, अर्थ है गामा घनत्व जहां , और मतलब , ( ) के साथ एक घातांक वितरण को दर्शाता है । स्पष्ट रूप से, ।f ( y ) = θ nGamma(n,θ)θ,n>0ऍक्स्प(θ)1/θθ>0ऍक्स्प(θ)गामा(1,θ)f(y)=θnΓ(n)eθyyn11y>0θ,n>0Exp(θ)1/θθ>0Exp(θ)Gamma(1,θ)

के बाद से के लिए पूरा पर्याप्त है और , द्वारा लेहमैन-Scheffe प्रमेय the का UMVUE है । इसलिए हमें इस सशर्त अपेक्षा को खोजना होगा। θ ( एक्स 1 ) = θT=i=1nlnXiθ(एक्स1|टी)θE(X1)=θ1+θE(X1T)θ1+θ

हम ध्यान दें कि ।Xii.i.dBeta(θ,1)lnXii.i.dExp(θ)TGamma(n,θ)

विधि I:

चलो और , ताकि और स्वतंत्र हैं। दरअसल, और , का अर्थ है । वी = - Σ n मैं = 2 ln एक्स मैं यू वी यू ~ ऍक्स्प ( θ ) वी ~ गामा ( n - 1 , θ ) यू + वी ~ गामा ( n , θ )U=lnX1V=i=2nlnXiUVUExp(θ)VGamma(n1,θ)U+VGamma(n,θ)

तो, ।E(X1i=1nlnXi=t)=E(eUU+V=t)

अब हम का सशर्त वितरण पाते हैं ।UU+V

के लिए और , के संयुक्त घनत्व हैθ > 0 ( यू , वी )n>1θ>0(U,V)

fU,V(u,v)=θeθu1u>0θn1Γ(n1)eθvvn21v>0=θnΓ(n1)eθ(u+v)vn21u,v>0

परिवर्तनशील चर, यह तत्काल है कि का संयुक्त घनत्वएफ यू , यू + वी(U,U+V)

fU,U+V(u,z)=θnΓ(n1)eθz(zu)n210<u<z

चलो के घनत्व हो । इस प्रकार का सशर्त घनत्वयू + वी यू | यू + वी = z यू | यू + वी ( यू | z )fU+V()U+VUU+V=z

fUU+V(uz)=fU,U+V(u,z)fU+V(z)=(n1)(zu)n2zn110<u<z

इसलिए, ।E(eUU+V=z)=n1zn10zeu(zu)n2du

अर्थात्, frak का ( एक्स 1 |θ1+θ(1)E(X1T)=n1(T)n10Teu(Tu)n2du

विधि II:

के रूप में के लिए एक पूर्ण पर्याप्त आंकड़ा है , के किसी भी निष्पक्ष आकलनकर्ता जो है के एक समारोह की UMVUE हो जाएगा लेहमैन-शेफ़े प्रमेय द्वारा। तो हम के क्षणों को खोजने के लिए आगे बढ़ते हैं , जिसका वितरण हमें ज्ञात है। हमारे पास है,θ θTθ टीθθ1+θT -टीθ1+θT

E(T)r=0yrθneθyyn1Γ(n)dy=Γ(n+r)θrΓ(n),n+r>0

इस समीकरण का उपयोग करके हम प्रत्येक पूर्णांक लिए निष्पक्ष अनुमानक (और UMVUE के) प्राप्त । आर 11/θrr1

अब , हमारे पासθθ>1θ1+θ=(1+1θ)1=11θ+1θ21θ3+

के निष्पक्ष अनुमानकों को मिलाकर हम( 1 + टी)1/θr

E(1+Tn+T2n(n+1)+T3n(n+1)(n+2)+)=11θ+1θ21θ3+

अर्थात,

E(r=0Trn(n+1)...(n+r1))=θ1+θ

इसलिए _ , UMVUE ऑफ _ फ़्रेक कोθθ>1 जी ( टी ) = Σ आर = 0θ1+θ(2)g(T)=r=0Trn(n+1)...(n+r1)


मैं दूसरी विधि में बारे में निश्चित नहीं हूं ।0<θ<1

मैथेमैटिका के अनुसार , अधूरा गामा फ़ंक्शन के संदर्भ में समीकरण एक बंद रूप है। और समीकरण , हम उत्पाद को सामान्य गामा फ़ंक्शन के रूप में रूप में व्यक्त कर सकते हैं । यह शायद और बीच स्पष्ट संबंध प्रदान करता है ।( 2 ) n ( n + 1 ) ( n + 2 ) ( n + r - 1 ) n ( n + 1 ) ( n + 2 ) ( n + r - 1 ) = Γ ( n + r )(1)(2)n(n+1)(n+2)...(n+r1) (1)(2)n(n+1)(n+2)...(n+r1)=Γ(n+r)Γ(n)(1)(2)

Mathematica का उपयोग करके मैं यह सत्यापित कर सकता था कि और वास्तव में एक ही चीज हैं।( 2 )(1)(2)


वास्तव में, और बीच समानता in की शक्ति श्रृंखला विस्तार लिखकर और फिर अभिन्न और योग को इंटरचेंज करके लिखती है । ( 2 ) - यू ( 2 )(1)(2)eu(2)
स्टबबोर्नटॉम

1

मुझे लगता है कि ऊपरी अधूरा गामा फ़ंक्शन के संबंध में आपको और अधिक कॉम्पैक्ट उत्तर मिल सकता है। पहली विधि का उपयोग करते हुए, मुझे अभिव्यक्ति मिली

E[X1|X1X2Xn=eT]=(n1)01zr(1r)nzdr,
जेड = टी जहाँz=eT.

वोल्फ्राम अल्फा को प्राप्त करने के लिए इसे एकीकृत करता है

E[X1|X1X2Xn=eT]=eT(n1)Tn1[(n2)!Γ(n1,T)]

जब पूर्णांक होता है तो अपूर्ण गामा फ़ंक्शन शब्द का एक बंद रूप होता है। यह हैn

Γ(n1,T)=Γ(n1)eTj=0n2Tjj!

अपेक्षा को फिर से देखना और सरल करना, हम पाते हैं

E[X1|X1X2Xn=eT]=Γ(n)Tn1[eTj=0n2Tjj!]

मेरे पास उस सॉफ़्टवेयर तक पहुंच नहीं है जो आपके परिणामों और साथ तुल्यता को सत्यापित करेगा , लेकिन और लिए हाथ की गणना आपके साथ मेल खाते हैं ।(1)(2)n=2n=3(1)


जहाँ आपने लिखा है आपको[ एक्स 1 | एक्स 1 एक्स 2एक्स एन = टी ]E[X1x1x2xn=eT],E[X1X1X2Xn=eT].
माइकल हार्डी
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.