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

সীমাবদ্ধ বিপরীত সমস্যায় একাধিক কার্যকরের জন্য আস্থা অঞ্চলের সমসাময়িক ফ্রিকোয়েন্টিস্ট ক্যালিব্রেশন

মৌলিক তথ্য

  • পত্র ID: 2510.11708
  • শিরোনাম: সীমাবদ্ধ বিপরীত সমস্যায় একাধিক কার্যকরের জন্য আস্থা অঞ্চলের সমসাময়িক ফ্রিকোয়েন্টিস্ট ক্যালিব্রেশন
  • লেখক: Pau Batlle, Pratik Patil, Michael Stanley, Javier Ruiz Lupon, Houman Owhadi, Mikael Kuusela
  • শ্রেণীবিভাগ: math.ST stat.TH
  • প্রকাশনা সময়: ২০২৫ সালের ১৩ অক্টোবর
  • পত্র লিঙ্ক: 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} দেওয়া, লক্ষ্য অজানা ভেক্টর HxRkHx^* \in \mathbb{R}^k এর জন্য সীমিত নমুনা 1α1-\alpha ফ্রিকোয়েন্টিস্ট আস্থা সেট Rα(y)RkR_\alpha(y) \subseteq \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} হল Zx=λ(Hx,Kx+ε)Z_x = \lambda(Hx, Kx + \varepsilon) এর (1α)(1-\alpha) কোয়ান্টাইল।

মূল তাত্ত্বিক ফলাফল

উত্তলতা উপপাদ্য

উপপাদ্য 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, χ22\chi^2_2 বিতরণের 2.279 এবং 5.991 এর তুলনায় উল্লেখযোগ্য উন্নতি
  2. এলাকা হ্রাস: μ\mu-বর্ণনা xx-বর্ণনা সীমানা বাক্সের তুলনায় গড় এলাকা প্রায় 30-50% হ্রাস করে

সম্পর্কিত কাজ

ঐতিহাসিক বিকাশ

  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টি সম্পর্কিত সংদর্ভ উদ্ধৃত করে, গবেষণার জন্য একটি দৃঢ় তাত্ত্বিক ভিত্তি প্রদান করে।