Simultaneous false discovery rate control in location families

तकनीकी सारांश: लोकेशन फैमिलीज़ में एक साथ फॉल्स डिस्कवरी रेट (FDR) नियंत्रण

समस्या का विवरण
मल्टीपल हाइपोथीसिस टेस्टिंग में, फॉल्स डिस्कवरी रेट (FDR) को पारंपरिक रूप से केवल नुल वैल्यू (जैसे, $\theta = 0$) पर परिभाषित और नियंत्रित किया जाता है। हालांकि, लोकेशन फैमिलीज़ से जुड़े व्यावहारिक अनुप्रयोगों में, शोधकर्ताओं को न केवल नुल के लिए, बल्कि उन पैरामीटर मानों के लिए भी कम FDR की आवश्यकता होती है जिन्हें "व्यावहारिक रूप से महत्वहीन" (जैसे, छोटे प्रभाव आकार) माना जाता है। मानक बेन्जामिनी-होचबर्ग (BH) प्रक्रिया नुल के लिए FDR को नियंत्रित करती है, लेकिन यह स्पष्ट रूप से संपूर्ण "FDR कर्व" पर नियंत्रण की गारंटी नहीं देती है—अर्थात, पैरामीटर $\theta$ के एक फलन (function) के रूप में FDR पर। यह शोध पत्र लोकेशन फैमिलीज़ के भीतर पैरामीटर मानों के एक निरंतरता (continuum) में एक साथ FDR नियंत्रण की आवश्यकता को संबोधित करता है।

कार्यप्रणाली
लेखक एक उपयोगकर्ता-निर्दिष्ट स्तर के नीचे FDR कर्व को नियंत्रित करने के लिए मानक BH प्रक्रिया के सामान्यीकरण का प्रस्ताव करते हैं। सेटअप में स्वतंत्र डेटा $X_i \sim P_{\theta_i}$ शामिल है जो एक लोकेशन फैमिली से आता है जहाँ $P_\theta(X \le x) = F(x - \theta)$ है। लक्ष्य एक-तरफा परिकल्पनाओं $H_{i,\theta}: \theta_i \ge \theta$ का परीक्षण करना है।

1. FDR कर्व की परिभाषा: FDR कर्व $FDP_S(\theta)$ के एक्सपेक्टेशन (expectation) को एक $FDP$कर्व के रूप में परिभाषित किया गया है, जो एक दिए गए थ्रेशोल्ड$\theta$ के लिए खारिज की गई परिकल्पनाओं के बीच गलत खोजों (false discoveries) के अनुपात को मापता है।
2. सामान्यीकृत BH प्रक्रिया: एक स्केलर नॉमिनल लेवल $q$ के बजाय, यह विधि एक उपयोगकर्ता-निर्दिष्ट गैर-बढ़ते (non-increasing) FDR कर्व $q(\theta)$ का उपयोग करती है।
   • सामान्यीकृत P-मान (Normalized P-values): प्रत्येक हाइपोथीसिस $i$ के लिए, एक सामान्यीकृत p-मान की गणना की जाती है:
     $$\bar{P}_i = \sup_{\theta \in \mathbb{R}} \frac{P_{i,\theta}}{q(\theta)}$$
     जहाँ $P_{i,\theta} = F(X_i - \theta)$ $\theta$ के लिए परीक्षण हेतु p-मान है।
   • अस्वीकृति नियम (Rejection Rule): यह प्रक्रिया इन सामान्यीकृत p-मानों $\{\bar{P}_1, \dots, \bar{P}_m\}$ पर लक्ष्य स्तर 1 के साथ मानक BH एल्गोरिदम लागू करती है। अस्वीकृति सेट $S$ सबसे बड़ा उपसमुच्चय (subset) है कि $i \in S$ यदि और केवल यदि $\bar{P}_i \le |S|/m$।

मुख्य सैद्धांतिक परिणाम
यह शोध पत्र स्थापित करता है कि सामान्यीकृत BH प्रक्रिया FDR कर्व को एक साथ नियंत्रित करती है। विशेष रूप से, मान लीजिए $q^*(\theta)$ उपयोगकर्ता-निर्दिष्ट कर्व $q(\theta)$ का एक रूपांतरित संस्करण है, जिसे इस प्रकार परिभाषित किया गया है:
$$q^*(\theta) = \inf_{\theta' \in \mathbb{R}} \sup_{a \le x - \theta' \le b} q(\theta') \cdot \frac{F(x - \theta)}{F(x - \theta')}$$
जहाँ $a = F^{-1}(q(\theta')/m)$ और $b = F^{-1}(q(\theta'))$ है।

• थ्योरम 3.1: यह प्रक्रिया FDR कर्व को एक साथ नियंत्रित करती है। विशेष रूप से, $\sup_{\theta \in \mathbb{R}} \frac{FDR(\theta)}{q^*(\theta)} \le 1$। फलस्वरूप, किसी भी $\theta$ के लिए $FDR(\theta) \le q^*(\theta)$ लगभग निश्चित रूप से (almost surely) सत्य है।
• "फ्री लंच" घटना (The "Free Lunch" Phenomenon): एक महत्वपूर्ण परिणाम यह है कि मानक BH प्रक्रिया (जो एक स्टेप फंक्शन $q_{BH}(\theta) = q \cdot \mathbb{1}_{\theta \ge 0} + \mathbb{1}_{\theta < 0}$ का उपयोग करती है) स्वाभाविक रूप से पूरे FDR कर्व पर एक साथ नियंत्रण प्रदान करती है। वास्तविक नियंत्रण स्तर $q^*_{BH}(\theta)$ गैर-शून्य $\theta \neq 0$ के लिए नॉमिनल स्टेप फंक्शन से सख्ती से नीचे है, जबकि $\theta = 0$ पर इसे छूता है।
• टचिंग प्रॉपर्टी (Touching Property): "मोनोटोन रेश्यो प्रॉपर्टी" के तहत (जो सामान्य परिवारों जैसे कि गॉसियन द्वारा संतुष्ट की जाती है), प्रपोजिशन 3.1 सिद्ध करता है कि किसी भी घटते, राइट-कंटीन्यूअस $q$ के लिए, कम से कम एक बिंदु $\theta^*$ मौजूद है जहाँ वास्तविक नियंत्रण कर्व $q^*(\theta^*)$ लक्ष्य कर्व $q(\theta^*)$ को छूता है। इसका तात्पर्य है कि नियंत्रण हर जगह नॉमिनल स्पेसिफिकेशन से बेहतर नहीं हो सकता।
• निचली सीमाएँ (Lower Bounds): प्रपोजिशन च्या 3.2 वर्स्ट-केस FDR कर्व के लिए एक निचली सीमा प्रदान करता है, जो दिखाता है कि विश्लेषण सटीक है, विशेष रूप से जब $q^*(\theta)$ छोटा होता है।

विस्तार और संख्यात्मक उदाहरण
इस ढांचे को उन सेटिंग्स में विस्तारित किया गया है जहाँ हाइपोथीसिस अलग-अलग लोकेशन फैमिलीज़ (अलग-अलग $F_i$) से आते हैं। इस सामान्य मामले में, $q^*$ की परिभाषा विशिष्ट वितरणों पर सुप्रेमम लेने के लिए संशोधित की जाती है।

• प्रपोजिशन 4.1: पीसवाइज कांस्टेंट टारगेट कर्व्स के लिए, ऐसे बाधाओं (constraints) का एक उपसमुच्चय मौजूद है जिससे प्रेरित $q^*$ कर्व सभी निर्धारित बाधाओं को नियंत्रित करता है और कम से कम एक बार लक्ष्य को छूता है।
• ब्रेस्ट कैंसर डेटासेट: लेखक इस पद्धति को स्टोरी और तिबशिरानी (2003) ब्रेस्ट कैंसर डेटासेट पर लागू करते हैं। वे प्रभाव-आकार (effect sizes) बनाम सिग्नल-टू-नॉइज़ रेशियो (जहाँ एक सामान्य वितरण लागू होता है) के आधार पर परीक्षणों की तुलना करते हैं।
  • प्रभाव-आकार परिदृश्य में, सभी बाधाओं को एक साथ उपयोग करने से एकल बाधा की तुलना में कम अस्वीकृतियाँ हुईं, जो यह सुझाव देता है कि $q^* \le q$ स्थिति को संतुष्ट करने वाले न्यूनतम बाधा सेट का चयन करना इष्टतम है।
  • सिग्नल-टू-नॉइज़ परिदृश्य (सामान्य $F$) में, $q^*$ कर्व ने सभी निर्दिष्ट स्थानों पर लक्ष्य को छुआ, और अस्वीकृति सेट विभिन्न बाधा विकल्पों के बीच स्थिर रहे।

महत्ता और दावे
शोध पत्र एक "फ्री लंच" थ्योरम का दावा करता है: मानक BH प्रक्रिया, जो केवल नुल कंट्रोल के लिए डिज़ाइन की गई है, स्वाभाविक रूप से बिना किसी अतिरिक्त कम्प्यूटेशनल जटिलता या डेटा के गैर-नुल मानों के लिए पूरे FDR कर्व पर मजबूत एक साथ नियंत्रण प्रदान करती है। सामान्यीकृत प्रक्रिया शोधकर्ताओं को व्यावहारिक रूप से महत्वहीन पैरामीटर मानों पर FDR नियंत्रण को स्पष्ट रूप से लक्षित करने की अनुमति देती है, जो डिस्कवरी के लिए एक लचीला उपकरण प्रदान करती है। लेखक नोट करते हैं कि यह दृष्टिकोण चयनात्मक कंपाउंड निर्णय समस्याओं (Gao, Hu, and Zhao 2024 का संदर्भ देते हुए) के व्यापक ढांचे का एक विशिष्ट उदाहरण है, और अनुमानित प्रायर्स (priors) के आधार पर लक्षित कर्व $q$ चुनने के लिए एम्पेरिकल बायेसियन दृष्टिकोणों वाले संभावित भविष्य के दिशाओं का सुझाव देते हैं, हालांकि वे ऐसे कार्यान्वयन में "डबल-डिपिंग" (double-dipping) के मुद्दे के प्रति आगाह करते हैं।