নাল অনুমানের অধীনে উপদ্রব পরামিতি অজানা থাকলে বার্নার্ডের পরীক্ষা ব্যবহার করা হয়।
তবে ভদ্রমহিলা টেস্টিং পরীক্ষায় আপনি যুক্তি দিতে পারেন যে নাল অনুমানের অধীনে উপদ্রব পরামিতিটি 0.5 এ সেট করা যেতে পারে (অজানা মহিলার কাপের সঠিকভাবে অনুমান করার 50% সম্ভাবনা রয়েছে)।
তারপরে নাল অনুমানের অধীনে সঠিক অনুমানের সংখ্যাটি দ্বিপদী বিতরণে পরিণত হয়: প্রতিটি কাপের 50% সম্ভাবনা সহ 8 কাপ অনুমান করা।
অন্যান্য ক্ষেত্রে আপনার নাল অনুমানের জন্য এই তুচ্ছ 50% সম্ভাবনা নাও থাকতে পারে। এবং স্থির মার্জিন ছাড়া আপনি সম্ভবত জানেন না যে সম্ভাবনাটি কী হওয়া উচিত। সেক্ষেত্রে আপনার বার্নার্ডের পরীক্ষা দরকার।
এমনকি আপনি যদি চা পরীক্ষার ভদ্রমহিলার স্বাদ গ্রহণের বিষয়ে বার্নার্ড পরীক্ষা করেন, তবে এটি 50% হয়ে উঠবে (যদি ফলাফলটি সমস্ত সঠিক অনুমান হয়) যেহেতু সর্বাধিক পি-মান সহ উপদ্রব পরামিতি 0.5 হয় এবং তুচ্ছ দ্বিপদী পরীক্ষার ফলস্বরূপ ( এটি আসলে চারটি দুধের প্রথম কাপের জন্য দুটি দ্বি-দ্বি পরীক্ষার সংমিশ্রণ এবং চার টি প্রথম কাপের জন্য একটি)।
> library(Barnard)
> barnard.test(4,0,0,4)
Barnard's Unconditional Test
Treatment I Treatment II
Outcome I 4 0
Outcome II 0 4
Null hypothesis: Treatments have no effect on the outcomes
Score statistic = -2.82843
Nuisance parameter = 0.5 (One sided), 0.5 (Two sided)
P-value = 0.00390625 (One sided), 0.0078125 (Two sided)
> dbinom(8,8,0.5)
[1] 0.00390625
> dbinom(4,4,0.5)^2
[1] 0.00390625
নীচে এটি কীভাবে আরও জটিল ফলাফলের দিকে যাবে (যদি সমস্ত অনুমানগুলি সঠিক না হয় উদাহরণস্বরূপ 2 2 বনাম 4), তবে কী এবং চরম কী নয় তার গণনা কিছুটা আরও কঠিন হয়ে যায়
(নোট পাশাপাশি পরীক্ষা করে দেখুন যে বার্নার্ডের পরীক্ষার ক্ষেত্রে 4-2 ফলাফলের ক্ষেত্রে একটি উপদ্রব পরামিতি পি = 0.686 যা আপনি যুক্তি দিতে পারেন তা সঠিক নয়, 'চায়ে প্রথমে' উত্তর দেওয়ার 50% সম্ভাবনার পি-মান 0.08203125 হবে)। ওয়াল্ডের পরিসংখ্যানের ভিত্তিতে একটি আলাদা অঞ্চল বিবেচনা করার সময় এটি আরও ছোট হয়ে যায়, যদিও অঞ্চলটি নির্ধারণ করা এত সহজ নয় )
out <- rep(0,1000)
for (k in 1:1000) {
p <- k/1000
ps <- matrix(rep(0,25),5) # probability for outcome i,j
ts <- matrix(rep(0,25),5) # distance of outcome i,j (using wald statistic)
for (i in 0:4) {
for (j in 0:4) {
ps[i+1,j+1] <- dbinom(i,4,p)*dbinom(j,4,p)
pt <- (i+j)/8
p1 <- i/4
p2 <- j/4
ts[i+1,j+1] <- (p2-p1)/sqrt(pt*(1-pt)*(0.25+0.25))
}
}
cases <- ts < ts[2+1,4+1]
cases[1,1] = TRUE
cases[5,5] = TRUE
ps
out[k] <- 1-sum(ps[cases])
}
> max(out)
[1] 0.08926748
> barnard.test(4,2,0,2)
Barnard's Unconditional Test
Treatment I Treatment II
Outcome I 4 2
Outcome II 0 2
Null hypothesis: Treatments have no effect on the outcomes
Score statistic = -1.63299
Nuisance parameter = 0.686 (One sided), 0.314 (Two sided)
P-value = 0.0892675 (One sided), 0.178535 (Two sided)