बरनार्ड के परीक्षण का उपयोग तब किया जाता है जब उपद्रव पैरामीटर शून्य परिकल्पना के तहत अज्ञात होता है।
हालांकि महिला चखने के परीक्षण में आप यह तर्क दे सकते हैं कि उपद्रव पैरामीटर को शून्य परिकल्पना के तहत 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 बनाम 4), तो जो है और जो अति नहीं है उसकी गिनती थोड़ी अधिक कठिन हो जाती है
(ध्यान दें कि बरनार्ड के परीक्षण का उपयोग करता है, 4-2 परिणाम के मामले में एक उपद्रव पैरामीटर p = 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)