এটি লগলিস্টিভিলেন্সের জন্য কেবলমাত্র একটি সুবিধা, আরও কিছু নয়।
আমি পণ্য বনাম অঙ্কের সুবিধার বলতে চাইছেন: , অঙ্কের এমন differentialtion বা ইন্টিগ্রেশন হিসাবে অনেক শ্রদ্ধা, সঙ্গে মোকাবেলা করতে সহজ। এটি কেবল ঘনিষ্ঠ পরিবারগুলির জন্য সুবিধা নয়, আমি বলার চেষ্টা করছি।ln(∏ixi)=∑ilnxi
আপনি যখন কোনও এলোমেলো নমুনা নিয়ে কাজ করেন তখন সম্ভাবনাগুলি ফর্মের হয়: , সুতরাং লগ-লেস্টিভিলিটি পরিবর্তে এই পণ্যটিকে যোগফলের মধ্যে বিভক্ত করে দেয়, যা হেরফের এবং বিশ্লেষণ করা সহজ। এটি সাহায্য করে যে আমাদের যত্ন নেওয়া সমস্তই সর্বাধিকের বিন্দু, সর্বাধিকের মানটি গুরুত্বপূর্ণ নয়, তবে আমরা লোগারিদমের মতো কোনও একঘেয়ে রূপান্তর প্রয়োগ করতে পারি।L=∏ipi
বক্রতা স্বজ্ঞাত উপর। এটি মূলত লগলিস্টিভেনসিটির দ্বিতীয় ডেরাইভেটিভ হিসাবে শেষ পর্যন্ত একই জিনিস।
আপডেট: এটি আমি বক্ররেখাকে বোঝাতে চাইছিলাম। আপনার যদি ফাংশন থাকে , তবে এটির বক্রতাটি ওল্ফ্রামে ( দেখুন (14) হবে:
κ = f ″ ( x )y=f(x)
κ = চ''( এক্স )( 1 + চ'( এক্স )2)3 / 2
লগ সম্ভাবনা দ্বিতীয় ব্যুৎপন্ন:
এ = ( এলএন)চ( এক্স ) )''= চ''( এক্স )চ( এক্স )- ( চ'( এক্স )চ( এক্স ))2
সর্বাধিকের বিন্দুতে প্রথম ডেরাইভেটিভ স্পষ্টতই শূন্য, সুতরাং আমরা পেয়েছি:
অতএব, আমার প্রশ্নটি যে সম্ভাবনার বক্রতা এবং লগলিঙ্কিলিটিটির দ্বিতীয় ডেরাইভেটিভ একই জিনিস, সাজানো।
κআমি একটি এক্স= চ''( এক্সআমি একটি এক্স) = একটি চ( এক্সআমি একটি এক্স)
অন্যদিকে, যদি সম্ভাবনার প্রথম ব্যয়টি কেবলমাত্র সর্বাধিকের বিন্দুতে নয়, তবে সম্ভাবনা কার্যটি সমতল হয় তবে আমরা পাই:
κ≈f′′(x)≈Af(x)
এখন ফ্ল্যাট সম্ভাবনা আমাদের পক্ষে ভাল জিনিস নয়, কারণ এটি সর্বাধিক সংখ্যাসূচকভাবে অনুসন্ধান করা সহজ করে তোলে এবং সর্বাধিক সম্ভাবনা এটির আশেপাশের অন্যান্য পয়েন্টগুলির চেয়ে ভাল নয়, যেমন প্যারামিটারের অনুমানের ত্রুটি বেশি।
এবং আবারও, আমাদের এখনও বক্রতা এবং দ্বিতীয় ডেরিভেটিভ সম্পর্ক রয়েছে। তাহলে ফিশার সম্ভাবনা ফাংশনটির বক্রতাটি কেন দেখেনি? আমি মনে করি এটি সুবিধার একই কারণেই। পণ্যটির পরিবর্তে অঙ্কের কারণে লগলিফিকিলিটি ম্যানিপুলেট করা সহজ। সুতরাং, লগলিস্টিটির দ্বিতীয় ডেরাইভেটিভ বিশ্লেষণ করে তিনি সম্ভাবনার বক্রতা অধ্যয়ন করতে পারতেন। যদিও বক্রতা জন্য খুবই সহজ সমীকরণ সৌন্দর্য κmax=f′′(xmax) , বাস্তবে আপনি পণ্যটির দ্বিতীয় ডেরাইভেটিভ নিচ্ছেন, যা দ্বিতীয় ডেরিভেটিভসের যোগফলের চেয়ে মেসেজ।
আপডেট 2:
এখানে একটি বিক্ষোভ। আমি একটি (সম্পূর্ণরূপে তৈরি) সম্ভাবনা ফাংশন আঁকি, এর ক) বক্রতা এবং খ) এর লগের ২ য় ডেরাইভেটিভ। বাম দিকে আপনি সংকীর্ণ সম্ভাবনা দেখতে পাবেন এবং ডানদিকে এটি প্রশস্ত। আপনি দেখতে পাচ্ছেন কীভাবে সর্বাধিক সম্ভাবনার বিন্দুতে ক) এবং খ) একত্রিত হওয়া উচিত, যেমন তাদের উচিত। আরও গুরুত্বপূর্ণ যদিও, আপনি এর লগ-সম্ভাবনার ২ য় ডেরাইভেটিভ পরীক্ষা করে সম্ভাবনা ফাংশনের প্রস্থ (বা ফ্ল্যাটনেস) অধ্যয়ন করতে পারেন। যেমনটি আমি আগে লিখেছি যেমন বিশ্লেষণ করা তার চেয়ে প্রযুক্তিগতভাবে সহজতর।
আশ্চর্যজনকভাবে নয় যে লগলিস্টিলিটিসের গভীরতর 2 টি ডেরাইভেটিভ তার সর্বোচ্চটির চারপাশে চাটুকার হওয়ার সম্ভাবনা ফাংশনকে ইঙ্গিত দেয়, এটির জন্য এটি পছন্দসই নয় কারণ এটি বড় পরামিতি অনুমানের ত্রুটির কারণ হয়ে দাঁড়ায়।
আপনি প্লটগুলির পুনঃ উত্পাদন করতে চান ক্ষেত্রে ম্যাটল্যাব কোড:
f=@(x,a)a.^2./(a.^2+x.^2);
c = @(x,a)(-2*a.^2.*(a.^2-3*x.^2)./(a.^2+x.^2).^3/(4*a.^4.*x.^2/(a.^2+x.^2).^4+1).^(3/2));
ll2d = @(x,a)(2*(x.^2-a.^2)./(a.^2+x.^2).^2);
h = 0.1;
x=-10:h:10;
% narrow peak
figure
subplot(1,2,1)
a = 1;
y = f(x,a);
plot(x,y,'LineWidth',2)
%dy = diff(y)/h;
hold on
%plot(x(2:end),dy)
plot(x,c(x,a),'LineWidth',2)
plot(x,ll2d(x,a),'LineWidth',2)
title 'Narrow Likelihood'
ylim([-2 1])
% wide peak
subplot(1,2,2)
a=2;
y = f(x,a);
plot(x,y,'LineWidth',2)
%dy = diff(y)/h;
hold on
%plot(x(2:end),dy)
plot(x,c(x,a),'LineWidth',2)
plot(x,ll2d(x,a),'LineWidth',2)
title 'Wide Likelihood'
legend('likelihood','curvature','2nd derivative LogL','location','best')
ylim([-2 1])
আপডেট 3:
উপরের কোডে আমি বক্ররেখার সমীকরণে কিছু স্বেচ্ছাচারী বেল আকারের ফাংশন প্লাগ করেছিলাম, তারপরে এর লগের দ্বিতীয় ডেরাইভেটিভ গণনা করেছি। আমি কোনও কিছুর পুনরায় স্কেল করিনি, আমি আগে উল্লিখিত সমতাটি দেখানোর জন্য সমীকরণ থেকে মানগুলি সোজা are
এখানে বিশ্ববিদ্যালয়ে থাকাকালীন ফিশার সম্ভবত প্রকাশের সম্ভাব্যতার প্রথম প্রবন্ধটি প্রকাশ করেছেন, "ফিটিং ফ্রিকোয়েন্সি কার্ভস সম্পর্কিত একটি পরম মাপদণ্ড", গণিতের মেসেঞ্জার, 41: 155-160 (1912)
যেহেতু আমি সবাইকে জোর দিয়ে বলছিলাম যে তিনি এন্ট্রপি এবং অন্যান্য অভিনব বিষয়গুলির জন্য লগ সম্ভাবনার কোনও "গভীর" সংযোগের কথা উল্লেখ করেন না, তিনি এখনও তার তথ্যের মানদণ্ডটি সরবরাহ করেন না। তিনি কেবল সমীকরণ রাখে p.54 তারপর সম্ভাব্যতা পূর্ণবিস্তার সম্পর্কে কথা বলতে আয়। আমার মতে, এটি দেখায় যে তিনি যৌথ সম্ভাবনাগুলি নিজেরাই বিশ্লেষণের একটি সুবিধাজনক পদ্ধতি হিসাবে লোগারিদমকে ব্যবহার করছেন। এটি ক্রমাগত বক্ররেখার ফিটিংয়ে বিশেষত কার্যকর, যার জন্য তিনি p.55: লগ পি = এ একটি সুস্পষ্ট সূত্র দেন
logP′=∑n1logp
logP=∫∞−∞logfdx
P
কাগজটি পড়ার সময় একটি বিষয় লক্ষণীয় যে তিনি কেবলমাত্র সর্বাধিক সম্ভাবনা অনুমানের কাজ দিয়ে শুরু করেছিলেন এবং পরবর্তী 10 বছরে আরও কাজ করেছিলেন, সুতরাং এমএলই শব্দটি এখনও তৈরি হয়নি, যতদূর আমি জানি।