संभवतः अनिर्णायक परिणामों के साथ एक परीक्षण का एक अन्य उदाहरण अनुपात के लिए द्विपद परीक्षण है जब केवल अनुपात, नमूना आकार नहीं, उपलब्ध है। यह पूरी तरह से अवास्तविक नहीं है - हम अक्सर "73% लोग इस बात से सहमत हैं ..." फॉर्म के दावों को खराब देखते हैं या सुनते हैं, जहां और कहीं नहीं है।
एच0: π= 0.5एच1: π≠ 0.5α = 0.05
पी = ५ %1195 %α = 0.05
पी = ४ ९ %
पी = ५० %एच0
पी = ० %पी = ५० %पी = ५ %पी = ० % । पी = % के समान थोड़ा महत्वपूर्ण नमूना,जो महत्वपूर्ण होगा; के लिए पीपी = १०० %पी = 16 %पीआर ( एक्स)≤ 3 ) ≈ .००,२२१ < 0.025पी = १ = %पीआर ( एक्स)≤ 1 ) 109 0.109 > 0.025पी = 16 %पी = 18 %पीआर ( एक्स)≤ 2 ) ≈ 0.0327 > 0.025पी = 19 %पीआर ( एक्स)≤ 3 ) ≈ 0.0106 < 0.025
पी = 24 %पी = १३ %α = 0.05: रेखा के नीचे के अंक असंदिग्ध रूप से महत्वपूर्ण हैं, लेकिन इसके ऊपर के लोग अनिर्णायक हैं। पी-वैल्यू का पैटर्न ऐसा है कि परिणामों के लिए स्पष्ट रूप से महत्वपूर्ण होने के लिए मनाया प्रतिशत पर एकल कम और ऊपरी सीमाएं नहीं हैं।
आर कोड
# need rounding function that rounds 5 up
round2 = function(x, n) {
posneg = sign(x)
z = abs(x)*10^n
z = z + 0.5
z = trunc(z)
z = z/10^n
z*posneg
}
# make a results data frame for various trials and successes
results <- data.frame(successes = rep(0:100, 100),
trials = rep(1:100, each=101))
results <- subset(results, successes <= trials)
results$percentage <- round2(100*results$successes/results$trials, 0)
results$pvalue <- mapply(function(x,y) {
binom.test(x, y, p=0.5, alternative="two.sided")$p.value}, results$successes, results$trials)
# make a data frame for rounded percentages and identify which are unambiguously sig at alpha=0.05
leastsig <- sapply(0:100, function(n){
max(subset(results, percentage==n, select=pvalue))})
percentages <- data.frame(percentage=0:100, leastsig)
percentages$significant <- percentages$leastsig
subset(percentages, significant==TRUE)
# some interesting cases
subset(results, percentage==13) # inconclusive at alpha=0.05
subset(results, percentage==24) # unambiguously sig at alpha=0.05
# plot graph of greatest p-values, results below red line are unambiguously significant at alpha=0.05
plot(percentages$percentage, percentages$leastsig, panel.first = abline(v=seq(0,100,by=5), col='grey'),
pch=19, col="blue", xlab="Rounded percentage", ylab="Least significant two-sided p-value", xaxt="n")
axis(1, at = seq(0, 100, by = 10))
abline(h=0.05, col="red")
( इस StackOverflow प्रश्न से राउंडिंग कोड छीन लिया गया है ।)