এখানে R
একটি ছবি উত্পন্ন করার জন্য কিছু কোড রয়েছে। এটি তাদের আদেশের বিরুদ্ধে প্লট করা 15 সিমুলেটেড পি-মানগুলি দেখায়। সুতরাং তারা একটি আরোহণ পয়েন্ট প্যাটার্ন গঠন। লাল / বেগুনি রেখার নীচের পয়েন্টগুলি 0.1 বা 0.2 স্তরে উল্লেখযোগ্য পরীক্ষার প্রতিনিধিত্ব করে। এফডিআর লাইনের নীচে কালো পয়েন্টগুলির সংখ্যা লাইনের নীচে পয়েন্টের মোট সংখ্যার দ্বারা বিভক্ত হয়।
x0 <- runif(10) #p-values of 10 true null hypotheses. They are Unif[0,1] distributed.
x1 <- rbeta(5,2,30) # 5 false hypotheses, rather small p-values
xx <- c(x1,x0)
plot(sort(xx))
a0 <- sort(xx)
for (i in 1:length(x0)){a0[a0==x0[i]] <- NA}
points(a0,col="red")
points(c(1,15), c(1/15 * 0.1 ,0.1), type="l", col="red")
points(c(1,15), c(1/15 * 0.2 ,0.2), type="l", col="purple")
আমি আশা করি এটি অর্ডার করা পি-মানগুলির বন্টনের যে আকৃতি সম্পর্কে কিছুটা অনুভূতি দিতে পারে। যে লাইনগুলি সঠিক এবং উদাহরণস্বরূপ কিছু নীতিগর্ভ আকারের বক্ররেখার নয়, তা অর্ডার বিতরণের আকারের সাথে করতে হবে। এটি সুস্পষ্টভাবে গণনা করতে হবে। আসলে, লাইনটি কেবল একটি রক্ষণশীল সমাধান।