বিকল্প অনুমানের অধীনে সঠিক সংখ্যার অনুমানের বন্টন একটি অ-কেন্দ্রীয় হাইপারজোমেট্রিক বিতরণকে অনুসরণ করে , যা প্রতিকূলতার অনুপাতের দিক দিয়ে পরামিতি করা হয়, অর্থাৎ ভদ্রমহিলা "চা প্রথমে" অনুমান করার সময় যে প্রতিকূলতা কত বেশি যখন আসলে দুধ আগে যুক্ত করা হয়েছিল (বা অন্যদিকে চারপাশে) আসলে তার বিপরীতে প্রথমে চা যুক্ত করা হয়েছিল। যদি বিজোড় অনুপাত 1 হয় তবে আমরা কেন্দ্রীয় হাইপারজমেট্রিক বিতরণ পাই।
দেখা যাক এটি কাজ করে কিনা। আমি MCMCpack
প্যাকেজটি ব্যবহার করে উদাহরণের প্রয়োজনে আর ব্যবহার করব , যা dnoncenhypergeom()
একটি (নন-সেন্ট্রাল) হাইপারজমেট্রিক বিতরণের ঘনত্বের কম্পিউটিংয়ের জন্য কাজ করে । এটা তোলে আর্গুমেন্ট নেই x
অনুমান সঠিক সংখ্যার জন্য (সতর্ক থাকুন: এই দুই শর্ত এক অধীনে অনুমান সঠিক সংখ্যা, উদাহরণস্বরূপ, যখন চা সত্যিই প্রথম যোগ করা হয়েছিল), আর্গুমেন্ট n1
, n2
এবং m1
চার মার্জিন তিন জন্য, এবং psi
জন্য সত্য প্রতিকূল অনুপাত। x
আসল প্রতিকূল অনুপাত 1: যখন আসল ঘনত্বটি 0 থেকে 4 এর সমান (4 টির সমান মার্জিন সহ) গণনা করা যাক
install.packages("MCMCpack")
library(MCMCpack)
sapply(0:4, function(x) dnoncenhypergeom(x, n1=4, n2=4, m1=4, psi=1))
এই ফলন:
[1] 0.01428571 0.22857143 0.51428571 0.22857143 0.01428571
সুতরাং, সেখানে একটি 1.43% সম্ভাবনা রয়েছে যে ভদ্রমহিলা 8 টি সঠিক অনুমান করবেন (অর্থাত্, তিনি সমস্ত 4 কাপ সঠিকভাবে অনুমান করেছিলেন যেখানে চা আগে যুক্ত হয়েছিল এবং তাই তিনি সমস্ত 4 কাপ সঠিকভাবে অনুমান করেছিলেন যেখানে দুধ আগে যুক্ত হয়েছিল) নাল অনুমানের অধীনে। এটি ফিশার নাল অনুমানকে প্রত্যাখ্যান করার জন্য যথেষ্ট পরিমাণে প্রমাণ হিসাবে বিবেচনা করেছিলেন This
প্রশ্নের মধ্যে বর্ণিত সম্ভাব্যতাগুলি বিজোড় অনুপাত গণনা করতে ব্যবহার করা যেতে পারে, (.90/(1−.90))/(.10/(1−.10))=81 (অর্থাত, odds(guessA|trueA)/odds(guessA|trueB))। মহিলাটি এখন 8 টি কাপ সঠিকভাবে অনুমান করবে (এখন, সে প্রথমে 4 টি কাপ অনুমান করবে যেখানে প্রথমে চা যুক্ত হয়েছিল এবং তাই 4 কাপও যেখানে দুধ আগে যুক্ত করা হয়েছিল) এখন কী সম্ভাবনা রয়েছে?
dnoncenhypergeom(4, n1=4, n2=4, m1=4, psi=81)
এই ফলন:
[1] 0.8312221
সুতরাং শক্তি তখন প্রায় 83%।