Θ এর UMVUE


10

যাক (X1,X2,,Xn) ঘনত্ব থেকে একটি র্যান্ডম নমুনা হতে

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

আমি M এর UMVUE সন্ধান করার চেষ্টা করছি θ1+θ

যুগ্ম ঘনত্ব (X1,,Xn) হয়

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

জনসংখ্যা PDF হিসেবে fθ এক প্যারামিটার সূচকীয় পরিবারের জন্যে, এই শো জন্য একটি সম্পূর্ণ যথেষ্ট পরিসংখ্যাত যে θ হল

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

যেহেতু E(X1)=θ1+θ , প্রথমভাবাতে,E(X1T)আমাকেθএর UMVUE দেবেθ1+θ লেহম্যান-Scheffe উপপাদ্য দ্বারা। না নিশ্চিত এই শর্তসাপেক্ষ প্রত্যাশা সরাসরি পাওয়া যাবে যদি বা এক শর্তাধীন বিতরণ খুঁজে পেতে হয়েছে X1i=1nlnXi

অন্যদিকে, আমি নিম্নলিখিত পদ্ধতির বিষয়টি বিবেচনা করেছি:

আমাদের Xii.i.dBeta(θ,1)2θlnXii.i.dχ22 , যাতে2θTχ2n2

তাই r ম অর্ডার কাঁচা মুহুর্ত 2θT প্রায় শূন্য, চি-স্কোয়ার পিডিএফ ব্যবহার করে গণনা করা হয়

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

তাই মনে হচ্ছে যে বিভিন্ন পূর্ণসংখ্যা পছন্দ জন্য r , আমি বিভিন্ন পূর্ণসংখ্যা ক্ষমতা পক্ষপাতিত্বহীন estimators (এবং UMVUEs) পেতে হবে θ । উদাহরণস্বরূপ, E(Tn)=1θ এবংE(1nT)=θআমাকে সরাসরি1এর ইউএমভিউ দেয়1θ এবংθযথাক্রমে।

এখন, যখন θ>1 আমাদের θ থাকে θ1+θ=(1+1θ)1=11θ+1θ21θ3+

আমি অবশ্যই 1 এর UMVUE পেতে পারি1θ,1θ2,1θ3 এবং আরও। মিশ্রন তাই এই UMVUE আমি প্রয়োজনীয় UMVUE পেতে পারেনθ1+θ । এই পদ্ধতিটি বৈধ নাকি আমার প্রথম পদ্ধতিটি নিয়ে এগিয়ে যাওয়া উচিত? যখন এটি উপস্থিত থাকে তখন UMVUE অনন্য, সুতরাং উভয়েরই আমাকে একই উত্তর দেওয়া উচিত।

স্পষ্টভাবে বলতে গেলে, আমি E ( 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?

জন্য 0<θ<1 , আমি পেতে হবে g(θ)=θ(1+θ+θ2+) , এবং UMVUE ভিন্ন হবে তাই।


প্রথম দৃষ্টিভঙ্গিতে শর্তাধীন প্রত্যাশা সরাসরি খুঁজে পাওয়া যায়নি এবং E(X1lnXi=t)=E(X1Xi=et) , আমি এটি অগ্রসর ছিল শর্তাধীন বিতরণX1Xi । যে জন্য, আমি যুগ্ম ঘনত্ব প্রয়োজন(X1,Xi)

আমি ভেরিয়েবলের পরিবর্তন ব্যবহার করেছি (X1,,Xn)(Y1,,Yn) যেমন Yi=j=1iXjসবার জন্যআমি=1,2,,এন। যুগ্ম সমর্থন এই নেতৃত্ব(ওয়াই1,,ওয়াইএন)হচ্ছেএস={ ( y 1 ,i=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) এর যৌথ ঘনত্বটি f y ( y 1 , y 2 , , y n ) = θ n হিসাবে পেয়েছি

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

(Y1,Yn) এর যৌথ ঘনত্ব তাই

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

আমি এখানে আলাদা আলাদা রূপান্তর ব্যবহার করতে পারি যা যৌথ ঘনত্বের উত্সকে কম জটিল করে তুলবে? আমি এখানে সঠিক রূপান্তর গ্রহণ করেছি কিনা তা নিশ্চিত নই।


মন্তব্য বিভাগে কয়েকটি দুর্দান্ত পরামর্শের ভিত্তিতে, আমি যৌথ ঘনত্বের পরিবর্তে (U,U+V) এর যৌথ ঘনত্ব (X1,Xi) পেয়েছি যেখানে U=lnX1 এবং V=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
jboman

1
পয়েন্ট (ভোরে) যেখানে আপনাকে পরিচয় করিয়ে এ আপনি ভেরিয়েবল পরিপ্রেক্ষিতে কাজ সংস্কারমুক্ত থাকা উচিত ওয়াই আমি = - লগ এক্স আমিএটা প্রায় অবিলম্বে তারা সমানুপাতিক হয় এর Γ ( 1 ) ডিস্ট্রিবিউশন, যা দ্রুত যুগ্ম বিতরণ বিবেচনা করার জন্য আপনার সমস্যা হ্রাস করা ( ইউ , ইউ + + ভি ) যেখানে ইউ ~ Γ ( 1 ) এবং ভী ~ Γ ( এন - 1 ) TYi=logXi.Γ(1)(U,U+V)UΓ(1)VΓ(n1). এটি গণিতের বাকী দুটি পৃষ্ঠাকে শর্টকাট দেবে এবং আপনাকে সমাধানের দ্রুত পথ দেবে।
whuber

@ স্পষ্ট করে বলতে গেলে, আপনি কি পরামর্শ দিচ্ছেন যে আমি এর ঘনত্বটি খুঁজে পাই এবং সেখান থেকে ( এক্স 1 , এক্স i ) এর ঘনত্বটি খুঁজে পাই ) ? আমি লক্ষ করেছি যে - ln এক্স i(lnX1,lnX1i=2nlnXi)(X1,Xi)lnXi হ'ল রেট সহ এক্সপেনশনাল ভেরিয়েবল (যা আপনি বলছেন তেমন গামা ভেরিয়েবলও) তবে এর সাথে কাজ করার কথা ভাবেনি। θ
জেদীআটম

1
@whuber কিন্তু কিভাবে আমি পেতে হবে থেকে ( Ln এক্স 1 | ) সরাসরি? E(X1...)E(lnX1...)
জেদীআটম

1
@ হুবার দয়া করে আমার সম্পাদনাটি একবার দেখুন। আমি প্রায় এটি সম্পন্ন করেছি তবে সেই অবিচ্ছেদ্য সাথে কী করব তা নিশ্চিত নই। আমি আমার হিসাব সঠিক বলে মোটামুটি আত্মবিশ্বাসী।
জেদীআটম

উত্তর:


5

দেখা যাচ্ছে যে আমার মূল পোস্টে উভয় পন্থা (আমার প্রাথমিক প্রচেষ্টা এবং মন্তব্য বিভাগে পরামর্শের ভিত্তিতে অন্য) একই উত্তর দেয়। প্রশ্নের সম্পূর্ণ উত্তরের জন্য আমি উভয় পদ্ধতির রূপরেখা করব।

এখানে, অর্থ গামা ঘনত্ব f ( y ) = θ nGamma(n,θ)যেখানেθ,এন>0, এবংমেপুঃ(θ)গড় সঙ্গে একটি সূচকীয় বণ্টনের উল্লেখ করে1/θ(θ>0)। স্পষ্টতই, Exp(θ)amগামা(1)f(y)=θnΓ(n)eθyyn11y>0θ,n>0Exp(θ)1/θθ>0Exp(θ)Gamma(1,θ)

যেহেতু θ এবং E ( X 1 ) = এর জন্য যথেষ্ট যথেষ্টT=i=1nlnXiθ , লেহম্যান-শেফি উপপাদক ই দ্বারা(এক্স1টি)এর ইউএমভিউE(X1)=θ1+θE(X1T) । সুতরাং আমাদের এই শর্তসাপেক্ষ প্রত্যাশা খুঁজে বের করতে হবে।θ1+θ

আমরা লক্ষ করি যে Xii.i.dBeta(θ,1)lnXii.i.dExp(θ)TGamma(n,θ)

পদ্ধতি আমি:

এবং V = - n i = 2 ln X i , যাতে U এবং V স্বতন্ত্র হয় Let প্রকৃতপক্ষে, ইউ এক্সপ ( θ ) এবং ভি গামা ( এন - 1 , θ ) , ইউ + ভি গামা ( এন , θ ) বোঝায়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 , z- র ) = θ এন(U,U+V)

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

যাক ঘনত্ব হতে ইউ + + ভি । সুতরাং এর শর্তাধীন ঘনত্ব ইউ | ইউ + + ভী = z- র হয় ইউ | ইউ + + ভি ( U | 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

অর্থাৎ θ এর UMVUE হল ( এক্স 1 | টি ) = - 1θ1+θ(1)E(X1T)=n1(T)n10Teu(Tu)n2du

পদ্ধতি দ্বিতীয়:

যেহেতু θ এর জন্য সম্পূর্ণ পর্যাপ্ত পরিসংখ্যান , তাই কোনও any এর নিরপেক্ষ অনুমানক θTθ যাটি এরএকটি ক্রিয়াকলাপহবেθএর UMVUEθθ1+θT লেহম্যান-Scheffe উপপাদ্য দ্বারা। সুতরাং আমরা-টিএর মুহুর্তগুলি সন্ধান করতে এগিয়ে যাই, যার বিতরণটি আমাদের জানা। আমাদের আছে,θ1+θT

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

এই সমীকরণটি ব্যবহার করে আমরা প্রতিটি পূর্ণসংখ্যার r 1 এর জন্য এর নিরপেক্ষ অনুমানক (এবং UMVUE এর) পাই ।1/θrr1

এখন , আমাদের কাছে have রয়েছে θθ>1θ1+θ=(1+1θ)1=11θ+1θ21θ3+

নিরপেক্ষ অনুমানের সংমিশ্রণে আমরা E ( 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+θ

θ>1θ1+θ(2)g(T)=r=0Trn(n+1)...(n+r1)


0<θ<1

(1)(2)n(n+1)(n+2)...(n+r1)n(n+1)(n+2)...(n+r1)=Γ(n+r)Γ(n)(1)(2)

(1)(2)


(1)(2)eu(2)

1

আমি মনে করি আপনি আরও কমপ্যাক্ট উত্তর পেতে পারেন যা আপনি উপরের অসম্পূর্ণ গামা ফাংশনের ক্ষেত্রে সম্মতি দিয়েছেন। প্রথম পদ্ধতিটি ব্যবহার করে আমি প্রকাশটি পেয়েছি

[এক্স1|এক্স1এক্স2এক্সএন=টি]=(এন-1)01z- রR(1-R)এন-z- রR,
কোথায় z- র=টি

ওল্ফ্রাম আলফা এটি পেতে সংহত করে

[এক্স1|এক্স1এক্স2এক্সএন=টি]=টি(এন-1)টিএন-1[(এন-2)!-Γ(এন-1,টি)]

এখন অসম্পূর্ণ গামা ফাংশন টার্মটি যখন বন্ধ থাকে তখন এনএকটি পূর্ণসংখ্যা এটাই

Γ(এন-1,টি)=Γ(এন-1)টিΣ=0এন-2টি!

Rewriting the expectation and simplifying, we find

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

I don't have access to the software that will verify equivalence with your results (1) and (2), but hand calculations for n=2 and n=3 do match with your (1).


Where you wrote E[X1x1x2xn=eT], you ought to have E[X1X1X2Xn=eT].
Michael Hardy
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.