আমি অন্যান্য উত্তরগুলি পছন্দ করি তবে এখনও কেউ নিচেরটি উল্লেখ করেনি। ইভেন্ট ঘটে এবং যদি কেবলমাত্র { m a x ( U , V ) ≤ t } হয় , তাই যদি U এবং V স্বতন্ত্র এবং ডাব্লু = মি একটি x ( U , V ) হয় তবে F ডাব্লু ( টি ) = এফ ইউ ( টি ) ∗। ইউ≤ টি , ভি ≤ টি }{ m a x ( ইউ, ভি) ≤ টি }ইউভীওয়াট= মি a এক্স ( ইউ, ভি) তাই জন্য α একটি ধনাত্মক পূর্ণসংখ্যা (বলুন, α = ঢ ) নিতে এক্স = মি একটি এক্স ( জেড 1 , । । । জেড এন ) যেখানে জেড 'র দ্বারা IIDএফওয়াট( টি ) = এফইউ( টি ) ∗ চভী( টি )α। = এনএক্স= মি একটি এক্স ( টু Z1,...Zn)জেড
জন্য আমরা পেতে switcheroo করতে এফ জেড = এফ এন এক্স , তাই এক্স যে দৈব চলক যেমন যে এর সর্বোচ্চ হবে এন স্বাধীন কপি হিসাবে একই ডিস্ট্রিবিউশন আছে জেড (এবং এই আমাদের পরিচিত এক বন্ধু হবে না , সাধারণভাবে)। α = 1 / এনএফজেড= এফএনএক্সএক্সএনজেড
ক্ষেত্রে একটি ইতিবাচক মূলদ সংখ্যা (বলুন, α = মি / এন ) যেহেতু পূর্ববর্তী থেকে অনুসরণ করে
( এফ জেড ) মি / এন = ( এফ 1 / এন জেড ) মি ।αα = মি / এন
( চজেড)মি / এন= ( চ1 / এনজেড)মি।
জন্য একটি অযৌক্তিক ইতিবাচক rationals একটা ক্রম চয়ন একটি ট সমকেন্দ্রি করার α ; তারপরে ক্রম এক্স কে (যেখানে আমরা প্রতিটি কে জন্য আমাদের উপরের কৌশলগুলি ব্যবহার করতে পারি ) এক্স পছন্দসই বিতরণে রূপান্তরিত হবে ।αএকটিটαএক্সটটএক্স
এটি আপনি যে বৈশিষ্ট্যটির সন্ধান করছেন তা নাও হতে পারে তবে এটি কীভাবে ভাবেন সে সম্পর্কে কিছুটা ধারণা দেয় জন্য α উপযুক্ত সুন্দর। অন্যদিকে, আমি সত্যিই নিশ্চিত হতে পারি না যে এটি আসলেই কতটা ভাল পেতে পারে: আপনার কাছে ইতিমধ্যে সিডিএফ রয়েছে, তাই শৃঙ্খলা নিয়ম আপনাকে পিডিএফ দেয় এবং আপনি সূর্য অস্তমিত হওয়া অবধি মুহুর্ত গণনা করতে পারবেন ...? এটা সত্য সবচেয়ে যে জেড 'র একটি থাকবে না এক্স যে জন্য পরিচিত α = √এফαজেডαজেডএক্স , তবে আমি আকর্ষণীয় কিছুসন্ধানেরজন্য উদাহরণের সাথে খেলতে চাইলে আমিএফ(জেড)=জেড,0<জেড<1 এরসাথে ইউনিট ব্যবধানেজেডটিসমানভাবে বিতরণ করারচেষ্টা করতে পারি।। = 2-√জেডএফ( জেড)) = জেড0 < জেড< 1
সম্পাদনা: আমি @ জেএমএস উত্তরে কিছু মন্তব্য লিখেছিলাম, এবং আমার পাটিগণিত সম্পর্কে একটি প্রশ্ন ছিল, তাই আমি আরও আশাবাদী হয়ে কী বোঝাতে চেয়েছি তা লিখব।
@ জেডএমএস উত্তরের মন্তব্যে @ কার্ডিনাল সঠিকভাবে লিখেছেন যে সমস্যাটি
বা আরও সাধারণভাবে যখন জেড অগত্যা এন ( 0 , 1 ) না হয় , তখন সহজ হয় আছে
এক্স = ছ - 1 ( Y ) = এফ - 1 ( এফ α ( Y ) ) ।
ছ- 1( y)) = Φ- 1( Φ)α( y)) ) ,
জেডএন( 0 , 1 )x = g- 1( y)) = চ- 1( চα( y)) ) ।
আমার বক্তব্যটি ছিল যে যখন
একটি সুন্দর বিপরীত ফাংশন থাকে আমরা কেবলমাত্র বুনিয়াদি বীজগণিত সহ
y = g ( x ) ফাংশনটি সমাধান করতে পারি । আমি মন্তব্যে লিখেছেন
ছ হওয়া উচিত
Y = ছ ( এক্স ) = এফ - 1 ( এফ 1 / α ( এক্স ) ) ।এফY= জি( এক্স )ছY= জি( এক্স ) = এফ- 1( চ1 / α( এক্স ) ) ।
আসুন একটি বিশেষ কেস নেওয়া যাক, জিনিসগুলি প্লাগ ইন করুন এবং দেখুন এটি কীভাবে কাজ করে। সিডিএফ এফ ( এক্স ) = ( 1 - ই - এক্স ) , এক্স > 0 ,
এবং বিপরীতমুখী সিডিএফ
এফ - 1 ( y ) = - এলএন ( 1 - y ) সহ এক্স (1) বিতরণ করা
যাক । জি
সন্ধানের জন্য সমস্ত কিছু প্লাগ করা সহজ ; কাজ শেষ হয়ে গেলে আমরা y = g ( x ) = - পাই
এক্স
এফ( এক্স ) = ( 1 - ই)- এক্স) , x > 0 ,
এফ- 1( y)) = - ln( 1 - y)).
g
সুতরাং, সংক্ষেপে, আমার দাবিটি যদি
X ∼ E x p ( 1 ) এবং আমরা
Y = - ln ( 1 - ( 1 - e - X ) সংজ্ঞায়িত করি তবে
) 1 / α ) ,
তারপরে
Y এর সিডিএফ থাকবে যা দেখতে
এফ ওয়াই ( y ) = (y=g(x)=−ln(1−(1−e−x)1/α)
X∼Exp(1)Y=−ln(1−(1−e−X)1/α),
Y
আমরা এটি সরাসরি প্রমাণ করতে পারি (
পি(Y≤y দেখুন)এবং এক্সপ্রেশনটি পেতে বীজগণিত ব্যবহার করতে পারি, শেষ ধাপের পরবর্তীটিতে আমাদের সম্ভাব্য ইন্টিগ্রাল ট্রান্সফর্ম দরকার)। আমি যে পাগল হয়ে গেছি (প্রায়শই পুনরাবৃত্তি হয়) সে ক্ষেত্রে ডাবল-চেক করার জন্য কিছু সিমুলেশন চালিয়েছি যা এটি কাজ করে, ... এবং তা ঘটে। নিচে দেখ. কোডটি সহজ করার জন্য আমি দুটি তথ্য ব্যবহার করেছি:
যদি X ∼ F হয় তবে ইউ = এফ ( এক্স ) ∼ ইউ এন আই এফ ( 0 , 1 )FY(y)=(1−e−y)α.
P(Y≤y)If X∼F then U=F(X)∼Unif(0,1).
যদি ইউ~ U n i f( 0 , 1 ) তারপর ইউ1/ α∼ বি ই টি এ ( α , ১ ) ।
সিমুলেশন ফলাফলের প্লটটি অনুসরণ করে।
প্লট (বিয়োগের লেবেল) তৈরি করতে ব্যবহৃত আর কোডটি
n <- 10000; alpha <- 0.7
z <- rbeta(n, shape1 = alpha, shape2 = 1)
y <- -log(1 - z)
plot(ecdf(y))
f <- function(x) (pexp(x, rate = 1))^alpha
curve(f, add = TRUE, lty = 2, lwd = 2)
ফিট বেশ ভাল দেখাচ্ছে, আমার মনে হয়? হয়তো আমি পাগল নই (এবার)?