2025-11-25T18:34:18.645261

Simultaneous Frequentist Calibration of Confidence Regions for Multiple Functionals in Constrained Inverse Problems

Batlle, Patil, Stanley et al.
Many scientific analyses require simultaneous comparison of multiple functionals of an unknown signal at once, calling for multidimensional confidence regions with guaranteed simultaneous frequentist under structural constraints (e.g., non-negativity, shape, or physics-based). This paper unifies and extends many previous optimization-based approaches to constrained confidence region construction in linear inverse problems through the lens of statistical test inversion. We begin by reviewing the historical development of optimization-based confidence intervals for the single-functional setting, from "strict bounds" to the Burrus conjecture and its recent refutation via the aforementioned test inversion framework. We then extend this framework to the multiple-functional setting. This framework can be used to: (i) improve the calibration constants of previous methods, yielding smaller confidence regions that still preserve frequentist coverage, (ii) obtain tractable multidimensional confidence regions that need not be hyper-rectangles to better capture functional dependence structure, and (iii) generalize beyond Gaussian error distributions to generic log-concave error distributions. We provide theory establishing nominal simultaneous coverage of our methods and show quantitative volume improvements relative to prior approaches using numerical experiments.
academic

المعايرة المتزامنة للتكرار للمناطق الثقة لعدة دوال في مسائل معكوسة مقيدة

المعلومات الأساسية

  • معرّف الورقة: 2510.11708
  • العنوان: المعايرة المتزامنة للتكرار للمناطق الثقة لعدة دوال في مسائل معكوسة مقيدة
  • المؤلفون: Pau Batlle, Pratik Patil, Michael Stanley, Javier Ruiz Lupon, Houman Owhadi, Mikael Kuusela
  • التصنيف: math.ST stat.TH
  • تاريخ النشر: 13 أكتوبر 2025
  • رابط الورقة: https://arxiv.org/abs/2510.11708

الملخص

تتطلب العديد من التحليلات العلمية المقارنة المتزامنة لعدة دوال للإشارة المجهولة، مما يستلزم بناء مناطق ثقة متعددة الأبعاد مع ضمانات تغطية تكرارية متزامنة تحت القيود الهيكلية (مثل عدم السلبية والقيود المستندة إلى الفيزياء). تقدم هذه الورقة إطار عمل موحد وموسع لطرق بناء مناطق الثقة القائمة على التحسين في مسائل معكوسة خطية مقيدة من خلال منظور اختبار الانعكاس الإحصائي. تستعرض الورقة أولاً التطور التاريخي لفترات الثقة القائمة على التحسين في إعداد الدالة الواحدة، من "الحدود الصارمة" إلى حدسية بوروس وتفنيدها الأخير من خلال إطار اختبار الانعكاس. ثم تمتد الإطار إلى إعداد متعدد الدوال. يمكن استخدام الإطار لـ: (i) تحسين ثوابت المعايرة للطرق السابقة، مما ينتج عنه مناطق ثقة أصغر مع الحفاظ على التغطية التكرارية؛ (ii) الحصول على مناطق ثقة متعددة الأبعاد قابلة للتعامل معها والتي لا تكون بالضرورة فائقة المستطيلة، لالتقاط هيكل الاعتماد بين الدوال بشكل أفضل؛ (iii) التعميم من توزيعات الأخطاء الغاوسية إلى توزيعات الأخطاء اللوغاريتمية المقعرة بشكل عام.

خلفية البحث والدافع

تعريف المشكلة

تدرس هذه الورقة مشكلة بناء مناطق ثقة متزامنة لعدة دوال في المسائل المعكوسة الخطية. نعتبر المسألة المعكوسة الخطية: y=Kx+εy = Kx^* + \varepsilon حيث yRny \in \mathbb{R}^n هي القيم المرصودة، وxRpx^* \in \mathbb{R}^p هي المعاملات المجهولة، وKRn×pK \in \mathbb{R}^{n \times p} هي المشغل الأمامي المعروف، وεRn\varepsilon \in \mathbb{R}^n هي الضوضاء العشوائية.

دافع البحث

  1. الحاجة العلمية: تتطلب العديد من التحليلات العلمية الاستدلال المتزامن على عدة دوال خطية HxHx^* للإشارة المجهولة، بدلاً من تقدير المعامل عالي الأبعاد الكامل xx^*
  2. المعلومات المقيدة: عادة ما يرضي المعامل الحقيقي xx^* قيوداً مستندة إلى المعرفة الفيزيائية السابقة (مثل عدم السلبية x0x^* \geq 0)
  3. التغطية المتزامنة: يجب ضمان التغطية التكرارية المتزامنة لجميع الدوال، وليس التغطية الهامشية فقط

قيود الطرق الموجودة

  1. المحافظة: طرق الحدود الصارمة المتزامنة (SSB) التقليدية محافظة جداً، حيث تبني أولاً مجموعة ثقة لـ xx^* ثم تعينها إلى فضاء الدوال
  2. قيود المستطيل: عادة ما تنتج الطرق الموجودة مناطق ثقة فائقة المستطيلة، وغير قادرة على التقاط هيكل الاعتماد بين الدوال
  3. مشاكل المعايرة: تفتقر الطرق التاريخية مثل حدسية بوروس إلى ضمانات نظرية صارمة

المساهمات الأساسية

  1. إطار عمل موحد: توحيد طرق بناء مناطق الثقة المقيدة أحادية الدالة ومتعددة الدوال من خلال منظور اختبار الانعكاس
  2. اختراقات نظرية:
    • إثبات التحدب لدوال الكميات المقابلة لإحصائيات الاختبار λu2\lambda_u^2 و λ1\lambda^1
    • تحديد موقع الحل الأمثل لمشاكل تحسين الكميات
    • إنشاء علاقات الهيمنة العشوائية بين إحصائيات الاختبار
  3. خوارزميات عملية:
    • توفير ثوابت معايرة مثلى لمسائل القيود غير السالبة
    • تطوير طريقة اختزال TFM للمسائل عالية الأبعاد
    • اقتراح تقنيات فصل فضاء الصفوف/فضاء الصفر
  4. تحسن الأداء: تقليل حجم مناطق الثقة بشكل كبير مقارنة بالطرق الكلاسيكية مع الحفاظ على معدل التغطية الاسمي

شرح الطريقة

تعريف المهمة

بالنظر إلى المصفوفة HRk×pH \in \mathbb{R}^{k \times p}، الهدف هو بناء مجموعة ثقة تكرارية 1α1-\alpha محدودة العينة Rα(y)RkR_\alpha(y) \subseteq \mathbb{R}^k للمتجه المجهول HxRkHx^* \in \mathbb{R}^k، بحيث: PyPx(HxRα(y))1αP_{y \sim P_x}(Hx \in R_\alpha(y)) \geq 1-\alpha لجميع xXx \in \mathcal{X} (مجموعة القيود).

إطار اختبار الانعكاس

مراجعة حالة الدالة الواحدة

لكل μR\mu \in \mathbb{R}، نعتبر اختبار الفرضية: H0:xΦμXمقابلH1:xXΦμH_0: x^* \in \Phi_\mu \cap \mathcal{X} \quad \text{مقابل} \quad H_1: x^* \in \mathcal{X} \setminus \Phi_\mu حيث Φμ={xRp:hTx=μ}\Phi_\mu = \{x \in \mathbb{R}^p: h^T x = \mu\}.

التمديد متعدد الدوال

لـ μRk\mu \in \mathbb{R}^k، نعرّف Φμ={xRp:Hx=μ}\Phi_\mu = \{x \in \mathbb{R}^p: Hx = \mu\}، يصبح اختبار الفرضية: H0:xΦμXمقابلH1:xXΦμH_0: x^* \in \Phi_\mu \cap \mathcal{X} \quad \text{مقابل} \quad H_1: x^* \in \mathcal{X} \setminus \Phi_\mu

ثلاث إحصائيات اختبار

تحلل الورقة ثلاث إحصائيات اختبار:

  1. الحد الثاني المقيد λc2(μ,y)\lambda_c^2(\mu, y): λc2(μ,y)=minHx=μ,AxbKxy22minAxbKxy22\lambda_c^2(\mu, y) = \min_{Hx=\mu, Ax \leq b} \|Kx - y\|_2^2 - \min_{Ax \leq b} \|Kx - y\|_2^2
  2. الحد الثاني غير المقيد λu2(μ,y)\lambda_u^2(\mu, y): λu2(μ,y)=minHx=μ,AxbKxy22minxRpKxy22\lambda_u^2(\mu, y) = \min_{Hx=\mu, Ax \leq b} \|Kx - y\|_2^2 - \min_{x \in \mathbb{R}^p} \|Kx - y\|_2^2
  3. الحد الواحد λ1(μ,y)\lambda^1(\mu, y): λ1(μ,y)=minHx=μ,AxbKxy22\lambda^1(\mu, y) = \min_{Hx=\mu, Ax \leq b} \|Kx - y\|_2^2

مشكلة المعايرة

لكل إحصائية اختبار، يجب تحديد عتبة لضمان معدل تغطية 1α1-\alpha:

  • العتبة المقسمة: d(μ)=supHx=μ,AxbQx,1αd^*(\mu) = \sup_{Hx=\mu, Ax \leq b} Q_{x,1-\alpha}
  • العتبة العامة: D=supAxbQx,1αD^* = \sup_{Ax \leq b} Q_{x,1-\alpha}

حيث Qx,1αQ_{x,1-\alpha} هي الكمية (1α)(1-\alpha) لـ Zx=λ(Hx,Kx+ε)Z_x = \lambda(Hx, Kx + \varepsilon).

النتائج النظرية الرئيسية

نظرية التحدب

النظرية 5.4: لأي 0<α<10 < \alpha < 1 ثابت، دالة الكمية Qu2(x)Q_u^2(x) هي دالة محدبة في xx.

نظرية التعظيم

النظرية 5.6 (القيود الخطية): تحت القيود الخطية AxbAx^* \leq b، supxPQu2(x)=maxi=1:mQu2(pi)\sup_{x \in P} Q_u^2(x) = \max_{i=1:m} Q_u^2(p_i) حيث {pi}i=1m\{p_i\}_{i=1}^m هي مجموعة النقاط الطرفية للمتعدد الوجوه PP.

النظرية 5.7 (قيود المخروط): تحت قيود المخروط xCx^* \in C، supxCQu2(x)=Qu2(0)\sup_{x \in C} Q_u^2(x) = Q_u^2(0)

الإعداد التجريبي

الأمثلة الرقمية

نعتبر مشكلة القيود غير السالبة: y=Kx+ε,εN(0,I),x0y = Kx + \varepsilon, \quad \varepsilon \sim N(0, I), \quad x^* \geq 0 حيث: K=(211011),H=(110011)K = \begin{pmatrix} 2 & 1 & 1 \\ 0 & 1 & 1 \end{pmatrix}, \quad H = \begin{pmatrix} 1 & -1 & 0 \\ 0 & 1 & -1 \end{pmatrix}

الطرق المقارنة

  1. SSB_x: حدود الصندوق الموصوفة بـ xx للحدود الصارمة المتزامنة
  2. SSB_μ: الوصف الموصوف بـ μ\mu للحدود الصارمة المتزامنة
  3. QuantileZero_x/μ: النسخة المحسنة باستخدام الثابت الأمثل
  4. Bonferroni: فترات المنتج المصححة بـ Bonferroni
  5. طريقة Split: تقنية فصل فضاء الصفوف/فضاء الصفر

مؤشرات التقييم

  • معدل التغطية التجريبي: التحقق من خلال إعادة العينات N=105N = 10^5
  • منطقة المنطقة: حساب مساحة منطقة الثقة باستخدام التكامل بالإحداثيات القطبية

نتائج التجارب

النتائج الرئيسية

مناطق الثقة للقيمة المرصودة الثابتة

بالنسبة لـ y=(0,0)y = (0,0) و y=(20,10)y = (20,10)، تنتج الطرق الموصوفة بـ μ\mu مجموعات محدبة تحتويها بدقة حدود الصندوق الموصوفة بـ xx، مما يقلل بشكل كبير من مساحة المنطقة.

توزيع التغطية والمساحة

  • x=(0,0,0)x^* = (0,0,0): تحقق طريقة QuantileZero_μ معدل تغطية دقيق تقريباً بنسبة 68%، بأصغر متوسط مساحة
  • x=(5,5,5)x^* = (5,5,5): جميع الطرق تحقق التغطية، لكن الطريقة الموصوفة بـ μ\mu لا تزال لها ميزة كبيرة في المساحة

التحسينات الكمية

  1. تحسن ثابت المعايرة: لمستويات الثقة 68% و 95%، الثوابت المثلى هي 1.644 و 5.139 على التوالي، مع تحسن كبير مقارنة بـ 2.279 و 5.991 لتوزيع χ22\chi^2_2
  2. تقليل المساحة: الوصف الموصوف بـ μ\mu يقلل متوسط المساحة بحوالي 30-50% مقارنة بحدود الصندوق الموصوفة بـ xx

الأعمال ذات الصلة

التطور التاريخي

  1. Burrus (1964): اقترح لأول مرة طريقة التحسين للفترات الثقة المقيدة
  2. Rust & O'Leary (1986): طوروا خوارزميات عملية
  3. Stark (1992): اقترح طريقة الحدود الصارمة
  4. Tenorio وآخرون (2007): طوروا تقنية اختزال TFM

التقدم الحديث

  • Batlle وآخرون (2023): تفنيد حدسية بوروس من خلال إطار اختبار الانعكاس
  • أدبيات الاستدلال المقيد: الارتباط بنظرية توزيع χ2\chi^2-bar
  • التنبؤ المطابق: الاختلافات في الأهداف والافتراضات

الخلاصة والمناقشة

الاستنتاجات الرئيسية

  1. المساهمات النظرية: إنشاء إطار عمل موحد لاختبار الانعكاس لمناطق الثقة المقيدة متعددة الدوال
  2. المزايا الحسابية: توفير خوارزميات قابلة للتوسع للمسائل عالية الأبعاد
  3. تحسن الأداء: تقليل حجم مناطق الثقة بشكل كبير مقارنة بالطرق الكلاسيكية

القيود

  1. إحصائية λc2\lambda_c^2: تفتقر دالة الكمية إلى التحدب، وتبقى مشكلة التعظيم مفتوحة
  2. التعقيد الحسابي: قد يكون البحث عن النقاط الطرفية صعباً في الحالات عالية الأبعاد
  3. العتبات المقسمة: حساب الدالة الكاملة d(μ)d^*(\mu) عادة ما يكون تحدياً

الاتجاهات المستقبلية

  1. التمديد غير الغاوسي: التمديد إلى توزيعات الأخطاء اللوغاريتمية المقعرة بشكل عام
  2. معايرة λc2\lambda_c^2: تطوير خوارزميات معايرة لإحصائية الحد الثاني المقيد
  3. النظرية المقاربة: دراسة الخصائص في العينات الكبيرة
  4. مجالات التطبيق: التمديد إلى القيود على الشكل والمسائل الإحصائية الأخرى

التقييم المتعمق

المزايا

  1. الصرامة النظرية: توفير إطار عمل رياضي كامل يتضمن إثباتات التحدب ونتائج الأمثلية
  2. القيمة العملية: تطوير خوارزميات قابلة للتوسع لحل المسائل الفعلية عالية الأبعاد
  3. المنظور الموحد: توحيد الطرق المتفرقة تاريخياً تحت إطار اختبار الانعكاس
  4. التحسينات الملحوظة: تقليل كبير في مناطق الثقة مع الحفاظ على الضمانات النظرية

أوجه القصور

  1. الفجوات النظرية: لم يتم بعد إنشاء النظرية الكاملة لإحصائية λc2\lambda_c^2
  2. القيود الحسابية: التعقيد الحسابي في بعض الحالات عالية الأبعاد
  3. نطاق التجارب: التجارب الرقمية نسبياً بسيطة، تفتقد التطبيقات الفعلية المعقدة

التأثير

  1. المساهمة الأكاديمية: توفير أساس نظري جديد لقياس عدم اليقين في المسائل المعكوسة المقيدة
  2. التطبيقات العملية: آفاق تطبيق واسعة في العلوم الفيزيائية والهندسة والمجالات الأخرى التي تتطلب استدلالاً مقيداً
  3. الأهمية المنهجية: قد يلهم إطار اختبار الانعكاس حل المسائل الإحصائية الأخرى

السيناريوهات المعمول بها

  • الاستدلال المتزامن متعدد الدوال في المسائل المعكوسة الخطية
  • تقدير المعاملات مع القيود الفيزيائية
  • الحسابات العلمية التي تتطلب ضمانات تكرارية صارمة
  • قياس عدم اليقين في مسائل التحسين المقيدة عالية الأبعاد

المراجع

تستشهد الورقة بـ 47 مرجعاً ذا صلة، تغطي الأعمال المهمة في مجالات الاستدلال المقيد والمسائل المعكوسة ونظرية التحسين والإحصاء، مما يوفر أساساً نظرياً متيناً للبحث.