স্বাভাবিক এলোমেলো ভেরিয়েবলের জন্য আনুমানিক অর্ডার পরিসংখ্যান


38

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

সম্পাদনা: স্পষ্ট করার জন্য, আমি আনুমানিক সূত্রগুলির সন্ধান করছি যা আরও-বা-কম স্পষ্টভাবে মূল্যায়ন করা যায়, সঠিক অখণ্ডতা প্রকাশ নয়।

উদাহরণস্বরূপ, আমি একটি সাধারণ আরভি'র প্রথম অর্ডার পরিসংখ্যানের জন্য (যেমন সর্বনিম্ন) নীচের দুটি অনুমানটি দেখেছি:

e1:nμn12n1σ

এবং

e1:nμ+Φ1(1n+1)σ

এগুলোর মধ্যে প্রথম, জন্য n=200 , দেয় প্রায় e1:200μ10σ যা দুর্দান্তভাবে আলগা আবদ্ধ মত মনে হয়।

দ্বিতীয় দেয় e1:200μ2.58σ যেহেতু একটি দ্রুত মন্টে কার্লো দেয় e1:200μ2.75σ , তাই এটি একটি খারাপ পড়তা না কিন্তু না মহান পারেন, এবং আরো গুরুত্বপূর্ণ আমি কোন অনুভূতি নেই ওটা কোথা থেকে এসেছে.

কোন সাহায্য?


4
আপনি যদি আর ব্যবহার করেন তবে পয়েন্টস ফাংশনটি দেখুন।
কার্ডিনাল

1
আপনার সম্ভাব্য তালিকাগুলির জন্য @ প্রোব্যাবিলিটিস্লোগিক কিছু ভাল অন্তর্দৃষ্টি দিয়েছে। আমি যদি বিকল্প দৃষ্টিকোণ থেকে আরও কিছু দিয়ে থাকি, বা আপনি কি এই বিষয়ে আপনার কৌতূহলকে সন্তুষ্ট করেছেন তবে কি আদৌ সহায়ক হবে?
কার্ডিনাল

উত্তর:


31

ক্লাসিক রেফারেন্স হ'ল রায়স্টন (1982) [1] যার সুস্পষ্ট সূত্রের বাইরে অ্যালগরিদম রয়েছে। এটি ব্লম (1958) দ্বারা একটি সুপরিচিত সূত্রটি উদ্ধৃত করেছে: সঙ্গেα=0,375। এই সূত্রটিn=200,r=1 এরজন্য -2.73 এর গুণক দেয়।E(r:n)μ+Φ1(rαn2α+1)σα=0.375n=200,r=1

[1]: অ্যালগরিদম এএস 177: প্রত্যাশিত সাধারণ আদেশের পরিসংখ্যান (সঠিক এবং আনুমানিক) জেপি রয়স্টন। রয়্যাল স্ট্যাটিস্টিকাল সোসাইটির জার্নাল। সিরিজ সি (প্রয়োগ পরিসংখ্যান) খণ্ড। 31, নং 2 (1982), পৃষ্ঠা 161-165


21

continuous কোনও ক্রমাগত র‌্যান্ডম এর আইথ অর্ডার পরিসংখ্যান বিতরণ পিডিএফ সহ ভেরিয়েবল "বিটা-এফ" যৌগিক বিতরণ দ্বারা দেওয়া হয়। এই বিতরণ সম্পর্কে ভাবার স্বজ্ঞাত উপায় হ'ল একটি নমুনায় বিথ অর্ডার পরিসংখ্যান বিবেচনা করা । এখন একটি র্যান্ডম ভেরিয়েবল এর ইথ অর্ডার পরিসংখ্যানের সমান হতে আমাদের 3 টি শর্ত প্রয়োজন:NXx
  1. i1 মান নীচে , প্রতিটি পর্যবেক্ষণের জন্য এটির সম্ভাব্যতা , যেখানে এলোমেলো পরিবর্তনশীল এক্স এর সিডিএফ isxFX(x)FX(x)=Pr(X<x)
  2. Ni উপরে মান , এটির সম্ভাব্যতাx1FX(x)
  3. একটি অন্তর অন্তর অন্তর 1 মান , এর সম্ভাব্যতা যেখানে হয় র্যান্ডম ভেরিয়েবল এর পিডিএফxfX(x)dxfX(x)dx=dFX(x)=Pr(x<X<x+dx)X

এই পছন্দটি করার জন্য এখানে , সুতরাং আমাদের কাছে রয়েছে:(N1)(N1i1)

fi(xi)=N!(i1)!(Ni)!fX(xi)[1FX(xi)]Ni[FX(xi)]i1dx

সম্পাদনা আমার মূল পোস্টে, আমি এই বিন্দু থেকে যাচ্ছে আরো একটি অত্যন্ত দরিদ্র চেষ্টা করেন, এবং নীচে মন্তব্য এই প্রতিফলিত করে। আমি নীচে এটি সংশোধন করার চেষ্টা করেছি

আমরা যদি এই পিডিএফটির গড় মূল্য গ্রহণ করি তবে আমরা পাই:

E(Xi)=xifi(xi)dxi

এবং এই অবিচ্ছেদে, আমরা নিম্নলিখিত পরিবর্তনশীল (@ ইঙ্গিত গ্রহণ) এর পরিবর্তন করব এবং অবিচ্ছেদ্য হয়ে যায়:pi=FX(xi)

E(Xi)=01FX1(pi)Beta(pi|i,Ni+1)dpi=EBeta(pi|i,Ni+1)[FX1(pi)]

সুতরাং এটি বিপরীতমুখী সিডিএফের প্রত্যাশিত মান, যা ডেল্টা পদ্ধতিটি দেওয়ার জন্য ভালভাবে অনুমান করা যায়:

EBeta(pi|i,Ni+1)[FX1(pi)]FX1[EBeta(pi|i,Ni+1)]=FX1[iN+1]

আরও ভাল অনুমান করার জন্য, আমরা ২ য় অর্ডারে প্রাইম প্রসারিত করতে পারি (প্রাইম ডোনটিং পৃথককরণ) এবং লক্ষ্য করে যে একটি বিপরীতের দ্বিতীয় ব্যয়টি হ'ল:

2a2FX1(a)=FX(FX1(a))[FX(FX1(a))]3=fX(FX1(a))[fX(FX1(a))]3

Let । তারপর আমাদের আছে:νi=FX1[iN+1]

EBeta(pi|i,Ni+1)[FX1(pi)]FX1[νi]VarBeta(pi|i,Ni+1)[pi]2fX(νi)[fX(νi)]3
=νi(iN+1)(1iN+1)2(N+2)fX(νi)[fX(νi)]3

এখন, সাধারণ ক্ষেত্রে বিশেষীকরণ করে আমাদের কাছে

fX(x)=1σϕ(xμσ)fX(x)=xμσ3ϕ(xμσ)=xμσ2fX(x)
FX(x)=Φ(xμσ)FX1(x)=μ+σΦ1(x)

দ্রষ্টব্য যে এবং প্রত্যাশা প্রায় হয়ে যায়:fX(νi)=1σϕ[Φ1(iN+1)]

E[xi]μ+σΦ1(iN+1)+(iN+1)(1iN+1)2(N+2)σΦ1(iN+1)[ϕ[Φ1(iN+1)]]2

এবং পরিশেষে:

E[xi]μ+σΦ1(iN+1)[1+(iN+1)(1iN+1)2(N+2)[ϕ[Φ1(iN+1)]]2]

যদিও @ হুবার উল্লেখ করেছেন, এটি লেজগুলিতে সঠিক হবে না। আসলে আমি মনে করি এটি আরও খারাপ হতে পারে, কারণ বিভিন্ন পরামিতিগুলির সাথে একটি বিটার স্কিউনেস


1
"একটি এলোমেলো পরিবর্তনশীল সর্বাধিক সম্ভাবনা অনুমানক "? এটি কী তা নিশ্চিত নয় তবে আমি মনে করি আপনি (প্রায়) মোডটি গণনা করেছেন ।
কার্ডিনাল

1
হঠাৎ সতর্কতা বা সংজ্ঞা ছাড়াই এবং হাজির হওয়ার পথে প্রায় দুই-তৃতীয়াংশ রহস্যজনক কিছু ঘটে । μσ
whuber

2
আমি "গাদা" করা বলতে চাইছি না, তবে বন্ধনীর পরিমাণ কীভাবে aণাত্মক সংখ্যার দ্বারা প্রায় অনুমান করা যায় তা আমার পক্ষেও শক্ত।
কার্ডিনাল

1
@ প্রব্যাবিলিসিলোগিক, ক্যালকুলাসের স্তরে থাকাকালীন, আপনি বলতে পারেন যে এই ক্ষেত্রে আমরা একটি দ্বিচারিত ফাংশন বিবেচনা করছি এবং কেবল অন্যটির পরিবর্তে একটি পরিবর্তনশীলকে সর্বোচ্চ করে তুলছি, আমি মনে করি গাণিতিক, পরিসংখ্যানগত এবং পাঠ্যক্রমিক কারণগুলির জন্য আপনাকে কল না করার কারণ রয়েছে 'সর্বাধিক সম্ভাবনার প্রাক্কলন' সম্পন্ন হয়েছে। তারা এই স্থানটিতে গণনা করার জন্য অনেক বেশি, তবে একটি সহজ যেটি আমি মনে করি যে যথেষ্ট পরিমাণে বাধ্যযোগ্য তা হ'ল আমরা কোনও কারণে কোনও কারণে পরিসংখ্যানগুলিতে একটি নির্দিষ্ট, আরকেন শব্দভাণ্ডার ব্যবহার করি। একক সমস্যার জন্য কৌতূহলে পরিবর্তন করা ভুল বোঝাবুঝির কারণ হতে পারে ... / ...
কার্ডিনাল

2
সংশোধিত উত্তরের জন্য @probabilityislogic (+1)। এক পরামর্শ, হয়তো বেশী ভালো মানে "বোঝা"। আপনি কিছু কনভার্জেন্স দাবি করছেন না তা বুঝতে কয়েক সেকেন্ডের জন্য কয়েকটি লাইনে ঘুরে দেখেছিল।
কার্ডিনাল

13

অনিকোর উত্তর ব্লমের সুপরিচিত সূত্রে নির্ভর করে যা একটি পছন্দ জড়িত । দেখা যাচ্ছে যে এই সূত্রটি নিজেই জি। এলফভিং (১৯৪,) এর কারণে একটি সঠিক উত্তরের নিছক অনুমান, সাধারণ জনসংখ্যার বায়োমেট্রিকা, খণ্ডের নমুনায় পরিসীমাটির বিস্ময়কর বন্টন । 34, পৃষ্ঠা 111-119। এল্ফিংয়ের সূত্রটি ন্যূনতম এবং সর্বাধিক স্যাম্পলকে লক্ষ্য করে, যার জন্য আলফার সঠিক পছন্দটি । Blom এর সূত্র ফলাফল যখন আমরা আনুমানিক দ্বারা ।α=3/8π/8π3

ব্লমের সান্নিধ্যের চেয়ে এলিভিং সূত্রটি ব্যবহার করে আমরা -2.744165 এর গুণক পাই। এই সংখ্যাটি এরিক পি এর সঠিক উত্তরের (-২.7466) এবং মন্টে কার্লো অনুমানের (-২.75)) এর নিকটতর, ব্লুমের কাছাকাছি (-২.73)) তুলনায়, সঠিক সূত্রের তুলনায় কার্যকর করা সহজ।


এল্ফিংয়ের মাধ্যমে কীভাবে আগমন হয় সে সম্পর্কে আপনি আরও কিছু বিশদ সরবরাহ করতে পারেন ? এটি নিবন্ধে স্পষ্ট নয়। α=π/8
অ্যান্টনি

1
অ্যান্টনি - আমি গণিত সংক্রান্ত পরিসংখ্যানের পাঠ্যপুস্তকের উপর নির্ভর করছি, স্যামুয়েল উইলকস, পাব দ্বারা। উইলে (1962)। পি তে 8.21 অনুশীলন করুন। 249 বলেছে: "যদি x_ (1), x_ (n) হল ক্রমাগত সিডিএফ এফ (এক্স) থেকে আকার n এর নমুনার বৃহত্তম এবং বৃহত্তম অর্ডার পরিসংখ্যান ... র্যান্ডম ভেরিয়েবল 2n * স্কয়ার্ট {[এফ (এক্স_ ( 1))] [1-এফ (x_ (n))]} এর পি / 2 এবং গড় 4- (পাই ^ 2) / 4 এর সাথে এন -> অনন্ত হিসাবে সীমাবদ্ধতা রয়েছে "" (দুঃখিত, আমি মার্কআপ কোড জানি না!) প্রতিসাম্য বিতরণের জন্য, এফ (এক্স_ (1)) = 1-এফ (এক্স_ (এন))। সুতরাং এফ (x_ (n)) পাই / (4 এন), বা x_ (n) প্রায় F ^ (- 1) (পাই / (4 এন)) এর মধ্যে। ব্লম সূত্রটি 3 / (4 এন) এর অনুমান ব্যবহার করে।
হাল এম এম সুইটকে

এটি আমাকে ইন্ডিয়ানা রাজ্য আইনসভায় দায়ী কুখ্যাত " " বিলের কথা মনে করিয়ে দেয় । (যদিও উইকিপিডিয়া নিবন্ধটি সুপারিশ করেছে যে গল্পটির জনপ্রিয় সংস্করণটি সঠিক নয়))π=3
স্টিভো'আমেরিকা

7

আপনি যা করতে চান তার উপর নির্ভর করে এই উত্তরটি সাহায্য করতে পারে বা নাও পারে - আমি ম্যাপেলের পরিসংখ্যান প্যাকেজ থেকে নিম্নলিখিত সঠিক সূত্রটি পেয়েছি ।

with(Statistics):
X := OrderStatistic(Normal(0, 1), 1, n):
m := Mean(X):
m;

1/2_t0n!2e1/2_t02(1/21/2erf(1/2_t02))1+n(1+n)!πd_t0

নিজে থেকে এটি খুব কার্যকর নয় (এবং এটি সম্ভবত হাতে থেকে সহজেই উত্পন্ন করা যেতে পারে, যেহেতু এটি সর্বনিম্ন এলোমেলো ভেরিয়েবল), তবে এটি প্রদত্ত মানগুলির জন্য দ্রুত এবং খুব নির্ভুল সান্নিধ্যের অনুমতি দেয় - এর চেয়ে অনেক বেশি নির্ভুল মন্টে কার্লো:nn

evalf(eval(m, n = 200));
evalf[25](eval(m, n = 200));

যথাক্রমে -2.746042447 এবং -2.746042447451154492412344 দেয়।

(সম্পূর্ণ প্রকাশ - আমি এই প্যাকেজটি বজায় রাখি))


1
@ প্রব্যাবিলিটিআইসলোগিক তার জবাবের প্রথমার্ধে সমস্ত অর্ডার পরিসংখ্যানের জন্য এই অবিচ্ছেদ্য উদ্ভূত হয়েছে।
হুবুহু
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.