সিডিএফ একটি শক্তি উত্থাপিত?


15

যদি FZ কোনও সিডিএফ হয় তবে মনে হয় FZ(z)α ( α>0 ) এছাড়াও একটি সিডিএফ।

প্রশ্ন: এটি কি স্ট্যান্ডার্ড ফলাফল?

প্রশ্ন: একটি ভালো উপায় একটি ফাংশন খুঁজে পেতে g সঙ্গে Xg(Z) St FX(x)=FZ(z)α , যেখানে xg(z)

মূলত, আমি হাতের কাছে পাইনে সিডিএফ, have FZ(z)α । কিছু হ্রাস আকারে আমি যে সিডিএফ উত্পাদন করে র্যান্ডম ভেরিয়েবল বৈশিষ্ট্যযুক্ত করতে চাই।

সম্পাদনা: বিশেষত জন্য বিশ্লেষণাত্মক ফলাফল পেতে পারলে আমি খুশি হব জেড~এন(0,1)। বা কমপক্ষে জেনে থাকুন যে এই জাতীয় ফলাফলটি অক্ষম।


2
হ্যাঁ, এটি একটি বেশ সুপরিচিত ফলাফল এবং সাধারণীকরণ করা সহজ। (কীভাবে?) আপনি g , কমপক্ষে অন্তর্নিহিতভাবেও খুঁজে পেতে পারেন । এটি মূলত একটি স্বেচ্ছাসেবী বিতরণের এলোমেলো রূপগুলি উত্পন্ন করতে সাধারণত ব্যবহৃত বিপরীতমুখী কৌশলটির একটি প্রয়োগ।
কার্ডিনাল

2
দয়া করে উত্তর দিন। দলটি পরে অভিযোগ করছে যে আমরা কম জবাব অনুপাতের সাথে লড়াই করছি না।

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

2
@ কার্ডিনাল আমাদের কারও কারও জন্য একই কারণ, একই কারণে দোষী ...
হুঙ্কার

2
@ ব্রায়ানজিডি হ্যাঁ, এটি একটি সুপরিচিত ফলাফল যা শিল্পজাতভাবে "সাধারণীকরণ" বিতরণ উত্পাদন করতে ব্যবহৃত হয়েছে, দেখুন । এটির মতো অনেকগুলি রূপান্তর রয়েছে এবং লোকেরা তাদের এ উদ্দেশ্যে ব্যবহার করে: তারা একটি প্যারাম্যাট্রিক রূপান্তর খুঁজে পায়, এটি একটি বিতরণ এবং voilá এ প্রয়োগ করে, আপনার কেবল একটি কাগজ রয়েছে কেবলমাত্র এর বৈশিষ্ট্য গণনা করে। এবং অবশ্যই, স্বাভাবিকটিই প্রথম 'শিকার'।

উত্তর:


11

আমি অন্যান্য উত্তরগুলি পছন্দ করি তবে এখনও কেউ নিচেরটি উল্লেখ করেনি। ইভেন্ট ঘটে এবং যদি কেবলমাত্র { m a x ( U , V ) t } হয় , তাই যদি U এবং V স্বতন্ত্র এবং ডাব্লু = মি একটি x ( U , V ) হয় তবে F ডাব্লু ( টি ) = এফ ইউ ( টি ) {ইউটি, ভীটি}{মিএকটিএক্স(ইউ,ভী)টি}ইউভীওয়াট=মিএকটিএক্স(ইউ,ভী) তাই জন্য α একটি ধনাত্মক পূর্ণসংখ্যা (বলুন, α = ) নিতে এক্স = মি একটি এক্স ( জেড 1 , জেড এন ) যেখানে জেড 'র দ্বারা IIDএফওয়াট(টি)=এফইউ(টি)*এফভী(টি)αα=এনX=max(Z1,...Zn)জেড

জন্য আমরা পেতে switcheroo করতে এফ জেড = এফ এন এক্স , তাই এক্স যে দৈব চলক যেমন যে এর সর্বোচ্চ হবে এন স্বাধীন কপি হিসাবে একই ডিস্ট্রিবিউশন আছে জেড (এবং এই আমাদের পরিচিত এক বন্ধু হবে না , সাধারণভাবে)। α=1/এনএফজেড=এফএক্সএনএক্সএনজেড

ক্ষেত্রে একটি ইতিবাচক মূলদ সংখ্যা (বলুন, α = মি / এন ) যেহেতু পূর্ববর্তী থেকে অনুসরণ করে ( এফ জেড ) মি / এন = ( এফ 1 / এন জেড ) মিαα=মি/এন

(এফজেড)মি/এন=(এফজেড1/এন)মি

জন্য একটি অযৌক্তিক ইতিবাচক rationals একটা ক্রম চয়ন একটি সমকেন্দ্রি করার α ; তারপরে ক্রম এক্স কে (যেখানে আমরা প্রতিটি কে জন্য আমাদের উপরের কৌশলগুলি ব্যবহার করতে পারি ) এক্স পছন্দসই বিতরণে রূপান্তরিত হবে ।αএকটিαএক্সএক্স

এটি আপনি যে বৈশিষ্ট্যটির সন্ধান করছেন তা নাও হতে পারে তবে এটি কীভাবে ভাবেন সে সম্পর্কে কিছুটা ধারণা দেয় জন্য α উপযুক্ত সুন্দর। অন্যদিকে, আমি সত্যিই নিশ্চিত হতে পারি না যে এটি আসলেই কতটা ভাল পেতে পারে: আপনার কাছে ইতিমধ্যে সিডিএফ রয়েছে, তাই শৃঙ্খলা নিয়ম আপনাকে পিডিএফ দেয় এবং আপনি সূর্য অস্তমিত হওয়া অবধি মুহুর্ত গণনা করতে পারবেন ...? এটা সত্য সবচেয়ে যে জেড 'র একটি থাকবে না এক্স যে জন্য পরিচিত α = এফজেডααজেডএক্স , তবে আমি আকর্ষণীয় কিছুসন্ধানেরজন্য উদাহরণের সাথে খেলতে চাইলে আমিএফ(জেড)=জেড,0<জেড<1 এরসাথে ইউনিট ব্যবধানেজেডটিসমানভাবে বিতরণ করারচেষ্টা করতে পারি।α=2জেডএফ(z- র)=z- র0<z- র<1


সম্পাদনা: আমি @ জেএমএস উত্তরে কিছু মন্তব্য লিখেছিলাম, এবং আমার পাটিগণিত সম্পর্কে একটি প্রশ্ন ছিল, তাই আমি আরও আশাবাদী হয়ে কী বোঝাতে চেয়েছি তা লিখব।

@ জেডএমএস উত্তরের মন্তব্যে @ কার্ডিনাল সঠিকভাবে লিখেছেন যে সমস্যাটি বা আরও সাধারণভাবে যখন জেড অগত্যা এন ( 0 , 1 ) না হয় , তখন সহজ হয় আছে এক্স = - 1 ( Y ) = এফ - 1 ( এফ α ( Y ) )

-1(Y)=Φ-1(Φα(Y)),
জেডএন(0,1)
এক্স=-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--এক্স), এক্স>0,
F1(y)=ln(1y).
g সুতরাং, সংক্ষেপে, আমার দাবিটি যদি X E x p ( 1 ) এবং আমরা Y = - ln ( 1 - ( 1 - e - X ) সংজ্ঞায়িত করি তবে ) 1 / α ) , তারপরে Y এর সিডিএফ থাকবে যা দেখতে এফ ওয়াই ( y ) = (
y=g(x)=ln(1(1ex)1/α)
XExp(1)
Y=ln(1(1eX)1/α),
Y আমরা এটি সরাসরি প্রমাণ করতে পারি (পি(Yy দেখুন)এবং এক্সপ্রেশনটি পেতে বীজগণিত ব্যবহার করতে পারি, শেষ ধাপের পরবর্তীটিতে আমাদের সম্ভাব্য ইন্টিগ্রাল ট্রান্সফর্ম দরকার)। আমি যে পাগল হয়ে গেছি (প্রায়শই পুনরাবৃত্তি হয়) সে ক্ষেত্রে ডাবল-চেক করার জন্য কিছু সিমুলেশন চালিয়েছি যা এটি কাজ করে, ... এবং তা ঘটে। নিচে দেখ. কোডটি সহজ করার জন্য আমি দুটি তথ্য ব্যবহার করেছি: যদি  X F হয়  তবে  ইউ = এফ ( এক্স ) ইউ এন আই এফ ( 0 , 1 )
FY(y)=(1ey)α.
P(Yy)
If XF then U=F(X)Unif(0,1).
যদি ইউ~ইউএনআমি(0,1) তারপর ইউ1/α~বিটিএকটি(α,1)

সিমুলেশন ফলাফলের প্লটটি অনুসরণ করে।

ECDF এবং F আলফাতে

প্লট (বিয়োগের লেবেল) তৈরি করতে ব্যবহৃত আর কোডটি

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)

ফিট বেশ ভাল দেখাচ্ছে, আমার মনে হয়? হয়তো আমি পাগল নই (এবার)?


@ জেএমএস উত্তরে আমার মন্তব্য দেখুন। জন্যজেড~এন(0,1) দেখে মনে হচ্ছে উত্তরটি আছে (z- র)=Φ-1(Φ1/α(z- র))যা বন্ধ রূপ নয় তবে সহজেই যথেষ্ট পরিমাণে গণনা করা যায়। এবং বিপরীত সিডিএফ-এ ইনপুটটি উপযুক্তভাবে নির্বাচিত বিটা বিতরণ তা স্বীকৃতি দিয়ে আপনি এটি আরও সহজ করে তুলতে পারেন। বিপরীতমুখী সিডিএফ সুন্দর ক্ষেত্রে উত্তরগুলি দুর্দান্ত হবে এবং আশেপাশে চলছে এমন কয়েকটি রয়েছে।

আপনার পাটিগণিত ডাবল চেক করা ভাল হবে।
কার্ডিনাল

@ কার্ডিনাল এরর ... ঠিক আছে, আমি করেছি, ... এবং ঠিক আছে? আপনি দয়া করে ত্রুটিটি নির্দেশ করুন?

(+1) দুঃখিত। আমি যখন এটি প্রথম দেখলাম তখন আমার মাথা কোথায় ছিল তা আমি নিশ্চিত নই। এটি স্পষ্টতই (ভাল, হওয়া উচিত ছিল!) সঠিক।
কার্ডিনাল

পছন্দ করেছেন আমি স্বীকার করি, যদিও, আপনি সত্যিই আমাকে এক মিনিটের জন্য ঘামছিলেন! :-)

14

কথা ছাড়া প্রমাণ

এখানে চিত্র বর্ণনা লিখুন

নীচের নীল বক্ররেখা হয় এফ, উপরের লাল বক্ররেখা হয় এফα (কেস টাইপ করা হচ্ছে) α<1), এবং তীরগুলি কীভাবে যেতে হবে তা দেখায় z- র প্রতি এক্স=(z- র)


সুন্দর ছবি! প্রশ্ন: এটি কি আঁকা ছিল? TikZ?
নিম্নমানের

1
@ ব্রায়ানজিডি: যদি আমি মনে করি, @ হুবুথ ম্যাথমেটিকা ​​ব্যবহার করে তার অনেক প্লট করেন।
কার্ডিনাল

3
@ কার্ডিনালাল আপনি ঠিক বলেছেন আসলে, আমি যা কিছু সুবিধাজনক তা ব্যবহার করি এবং দেখে মনে হয় এটি দ্রুত কোনও ভাল কাজ করবে। এফডব্লিউআইডাব্লু, এখানে কোডটি রয়েছে:Module[ {y, w, a = 0.1, z = 3.24, f = ChiDistribution[7.6], xmin=0, xmax=5}, y = CDF[f,z]; w = InverseCDF[f, y^(1/a)]; Show[ Plot[{CDF[f, x],CDF[f,x]^a} , {x, xmin, xmax}, Filling->{1->{2}}], Graphics[{ Dashed, Arrow[{{z,0}, {z,y}}], Arrow[{{z,y}, {w,y}}], Arrow[{{w,y}, {w,0}}] }] ] ]
হুবুহু

6

প্র 1) হ্যাঁ। এটি ভেরিয়েবলগুলি উত্পন্ন করার জন্য দরকারী যা স্টোকাস্টিকিক অর্ডার দেওয়া হয়; আপনি এটি @ ভুবারের সুন্দর ছবি :) থেকে দেখতে পারেন।α>1 স্টোকাস্টিক অর্ডার অদলবদল করে।

এটি একটি বৈধ সিডিএফ প্রয়োজনীয় শর্তাদি যাচাই করার বিষয় মাত্র: এফz- র(z- র)αহতে হয়েছে cadlag , nondecreasing এবং সীমা1 অনন্ত এবং 0 নেতিবাচক অনন্ত সময়ে। এফz- র এই বৈশিষ্ট্যগুলি তাই এটি দেখানো সহজ।

প্রশ্ন 2) বিশ্লেষণাত্মকভাবে এটির মতো বেশ শক্ত হবে বলে মনে হচ্ছে এফজেড বিশেষ


@ জেএমএস: কি সম্পর্কে? জেড~এন(0,1)?
নিম্নমানের

2
@ ব্রায়ানজিড: আমি এটি বিশ্বাস করি না। দিন একটি অবিচ্ছিন্নভাবে কঠোর একঘেয়ে ফাংশন হতে (তাই, একটি ভাল সংজ্ঞায়িত বিপরীতমুখী থাকার -1) যা আপনার শর্তগুলি সন্তুষ্ট করে। তারপরে, এটি অবশ্যই হবেΦα(তোমার দর্শন লগ করা)=পি((জেড)তোমার দর্শন লগ করা)=পি(জেড-1(তোমার দর্শন লগ করা))=Φ(-1(তোমার দর্শন লগ করা)) এবং তাই -1(তোমার দর্শন লগ করা)=Φ-1(Φα(তোমার দর্শন লগ করা))। সুতরাং বিপরীতটি মোটামুটি সুস্পষ্টভাবে চিহ্নিত করা হয়েছে, তবে তা নয়নিজেই। এটি সম্পর্কে আমার আগের মন্তব্যে আমি বোঝাতে চাইছিস্পষ্টভাবে পাওয়া যাচ্ছে ।
কার্ডিনাল

@ ব্রায়ানজিড - @ কার্ডিনাল যা বলেছে :) আমি এর জন্য একটি বিশেষ মামলার কথা ভাবতেও পারিনি এফজেডযেখানে আপনি একটি বদ্ধ ফর্ম পাবেন (অবশ্যই এটির একটি নেই)।
জেএমএস 23

@JMS: জেড~ইউ[0,1]একটি ইতিবাচক উদাহরণ হতে হবে।
কার্ডিনাল

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