---
title: Maximum A Posteriori Ratio Test (MAPRT)
url: https://www.emergentmind.com/topics/maximum-a-posteriori-ratio-test-maprt
type: topic
---

# Maximum A Posteriori Ratio Test (MAPRT)

Maximum A Posteriori Ratio Test (MAPRT) denotes a Bayesian decision rule that compares posterior support between competing hypotheses or models and selects the alternative whose posterior odds, posterior density ratio, or analogous posterior quantity exceeds a prescribed threshold. In the canonical binary form,
\[
\Lambda_{\mathrm{MAP}}(y)=\frac{p(y\mid H_1)\,P(H_1)}{p(y\mid H_0)\,P(H_0)}
\underset{H_0}{\overset{H_1}{\gtrless}} \eta,
\]
with \(\eta=1\) for standard MAP under equal priors and costs. The same Bayesian logic reappears in model-order selection through posterior odds between \(M_a\) and \(M_b\), in adaptive radar through EM-estimated posterior responsibilities, and in inverse problems through posterior mass or posterior-density comparisons over constrained parameter sets [2002.07673] [1902.05809] [2503.02214] [2402.00686] [1909.09402].

## 1. Decision-theoretic formulation

For binary hypotheses, MAPRT is a posterior-odds test. Using Bayes’ rule, the detector multiplies the likelihood ratio by the prior odds and compares the result with a threshold. In this form, MAPRT differs from the pure likelihood ratio test, which compares \(p(y\mid H_1)/p(y\mid H_0)\) alone, and from maximum likelihood, which ignores priors entirely. In the Gaussian case, the log-MAP ratio is the difference of two Gaussian log-likelihoods plus the log prior odds, so the decision boundary is quadratic when the covariances differ and linear when they coincide [2002.07673].

The same structure extends beyond simple binary testing. In model-order selection for HAADF STEM, one compares posterior odds
\[
PO_{a,b}=\frac{p(M_a\mid y)}{p(M_b\mid y)}
       =\frac{p(y\mid M_a)\,p(M_a)}{p(y\mid M_b)\,p(M_b)},
\]
and with a uniform prior over model orders the decision reduces to comparing marginal likelihoods. The selected model is
\[
M^*=\arg\max_k p(M_k\mid y).
\]
This makes MAPRT a natural bridge between hypothesis testing and Bayesian evidence maximization [1902.05809].

A distinct but related formulation arises in statistical inverse problems. There, the paper’s primary object is a MAP mass test that compares posterior probabilities of half-spaces,
\[
\Pi(X_K\mid Y=y)\quad \text{versus}\quad \Pi(X_H\mid Y=y),
\]
where \(X_H=\{x:\langle \phi,x\rangle\le 0\}\) and \(X_K=\{x:\langle \phi,x\rangle>0\}\). A posterior-density MAPRT can also be defined through
\[
T(y):=\frac{\sup_{u\in X_K}\pi(u\mid y)}{\sup_{u\in X_H}\pi(u\mid y)},
\]
but this is explicitly distinguished from both Bayes factors and the posterior-mass MAP test. A recurrent misconception is to identify all Bayesian tests with Bayes factors; the inverse-problem formulation shows that posterior set comparison, posterior density comparison, and model-evidence comparison are different constructions [2402.00686].

## 2. Gaussian discriminants in networked sensing

In networked sensing, MAPRT is instantiated on a discrete-time linear time-invariant network
\[
\mathbf{x}[k+1]=G\,\mathbf{x}[k]+\Pi\,\mathbf{w}[k],
\]
with measurements
\[
\mathbf{y}_{\mathcal J}[k]=C\,\mathbf{x}[k]+\mathbf{v}[k],\qquad
\mathbf{v}[k]\sim\mathcal N(0,\sigma_v^2 I),
\]
and stacked observation vector \(\mathbf Y_{\mathcal J}\). Under simple Gaussian input hypotheses
\[
H_i:\ \mathbf w[k]\overset{\text{i.i.d.}}{\sim}\mathcal N(\mu_i,\Sigma_i),\qquad i\in\{1,2\},
\]
the measurement mean and covariance are
\[
\overline{\mu}_i=\mathcal F(\mathbf 1_N\otimes \mu_i),\qquad
\overline{\Sigma}_i=\mathcal O\Sigma_0\mathcal O^\top+\mathcal F(I_N\otimes \Sigma_i)\mathcal F^\top+\sigma_v^2 I.
\]
The matrices \(\mathcal O\) and \(\mathcal F\) encode, respectively, observability and impulse-response structure, so the network filtering induced by \((G,\Pi,C)\) is embedded directly into the MAPRT statistic [2002.07673].

For the mean-shift model, \(\Sigma_1=\Sigma_2\equiv \Sigma_c\) and \(\mu_1\neq \mu_2\), the MAP detector reduces to a linear discriminant with sufficient statistic
\[
s=2\,\mu_\Delta^\top \overline{\Sigma}_c^{-1}\mathbf Y_{\mathcal J},
\qquad \mu_\Delta=\mu_2-\mu_1,
\]
and a threshold shifted by \(\gamma=\ln(\pi_1/\pi_2)\). For the covariance-shift model, \(\mu_1=\mu_2\equiv \mu_c\) and \(\Sigma_1\neq \Sigma_2\), the paper uses a one-dimensional linear discriminant \(y=b^\top \mathbf Y_{\mathcal J}\) chosen to maximize the \(I\)-divergence. The decision rule becomes
\[
\ln\!\left(\frac{d_1}{d_2}\right)-2\gamma
\underset{\widehat H=H_1}{\overset{\widehat H=H_2}{\gtrless}}
(y-b^\top \mu_c)^2\left(\frac{1}{d_2}-\frac{1}{d_1}\right),
\qquad d_i=b^\top \Sigma_i b,
\]
so the statistic is the centered energy \((y-b^\top \mu_c)^2\), with \(b\) selected through a generalized eigenproblem or, equivalently, as a principal eigenvector of \(\Sigma_2^{-1}\Sigma_1\) [2002.07673].

The central asymptotic characterization is expressed through the transfer function
\[
T(z)=C\,(zI-G)^{-1}\Pi.
\]
For mean shifts, the asymptotic signal-to-noise ratio is
\[
\eta^2
=
N\,\widetilde{\mu}_\Delta^\top
\Big([L^\top L+\sigma_v^2 I]^{-1}L^\top L\Big)\widetilde{\mu}_\Delta,
\qquad
L=T(1)\Sigma_c^{1/2},
\]
with \(\widetilde{\mu}_\Delta=\Sigma_c^{-1/2}(\mu_2-\mu_1)\). For covariance shifts, the key quantity is
\[
R=1+\sigma_v^{-2}\,\|T(z)\Sigma_1^{1/2}\|_\infty^2,\qquad
\tau=\frac{\ln R}{R-1}.
\]
The paper proves that the corresponding error probabilities decrease monotonically with \(\eta\) or \(R\), respectively. This yields a network-theoretic interpretation of detection performance in terms of input-output gain rather than only graph distance [2002.07673].

Sensor-placement consequences are formulated through node cutsets and downstream partitions. In the noiseless case, measuring the cutset is always at least as good as measuring the downstream partition. In the noisy case, however, the subnetwork transfer
\[
T_s(z)=(zI-G_{pp})^{-1}G_{pc}
\]
controls whether cutset or downstream sensors perform better. If the gain from cutset to partition is at most one, cutset sensing is superior; if it exceeds one, downstream sensing can outperform cutset sensing because the downstream subnetwork can amplify the signal relative to noise. The numerical examples make this reversal explicit: in a 10-node synthetic network with \(\sigma_v^2=1.2\), partition nodes achieved lower error probabilities than the cutset node when the “row sums \(>1\)” condition held, while in a 50-node random network the cutset uniformly outperformed every other 3-node subset when \(\|[G_{pp}\ G_{pc}]\|_\infty<1/\sqrt 7\) [2002.07673].

## 3. Distributed hypothesis testing over graphs

In distributed detection over a pairwise Markov random field, MAPRT appears as a per-node posterior-odds test derived from max-product beliefs. The posterior factorization is
\[
p(\mathbf x\mid \mathbf Y)\propto
\prod_{n\in\mathcal V}\phi_n(x_n)\,
\prod_{(i,j)\in\mathcal E}\psi_{ij}(x_i,x_j),
\]
with \(\phi_k(x_k)=p(\mathbf y_k\mid x_k)\) and \(\psi_{kj}(x_k,x_j)=e^{J_{kj}x_kx_j}\). At node \(j\), the log-ratio of max-product beliefs is
\[
\lambda_j^{(l)}=
\ln\frac{b_j^{(l)}(x_j=+1)}{b_j^{(l)}(x_j=-1)}
\gtrless \tau_j,
\]
which is the nodewise MAPRT in log form. The local observation contribution is
\[
\gamma_j=
\ln\frac{\phi_j(+1)}{\phi_j(-1)}
=
\ln\frac{p(\mathbf y_j\mid x_j=+1)}{p(\mathbf y_j\mid x_j=-1)},
\]
and for coherent Gaussian detection
\[
\gamma_j=\mathbf s_j^\top \mathbf y_j-\frac{1}{2}E_j,\qquad E_j=\|\mathbf s_j\|^2.
\]
The paper’s key result is that the max-product decision variable is a linear combination of local log-likelihood ratios:
\[
\lambda_j^{(l)}=\sum_{i\in\mathcal M_j^{(l)}} w_{ji}^{(l)}\,\gamma_i,
\qquad
\lambda_j=\sum_{i\in\mathcal V} w_{ji}\,\gamma_i
\]
after enough iterations [1909.09402].

This linearity is established by showing that the maximizers inside the message updates are affine functions of neighboring states, which yields linear message-LLR recursions. A one-hop approximation gives
\[
\lambda_j\approx \gamma_j+\sum_{k\in\mathcal N_j}J_{kj}\,\gamma_k.
\]
Thus, in the binary-Gaussian setting studied, max-product behaves as linear data fusion rather than as an intrinsically nonlinear detector. The paper reaches the same conclusion for a linearized sum-product recursion, where the coefficients become \(c_{jk}=\tanh(J_{kj}/2)\), and reports very similar behavior of max-product and sum-product under the distributed hypothesis test [1909.09402].

Performance analysis follows because, conditioned on a fixed network state, each \(\gamma_i\) is Gaussian or approximately Gaussian, so any linear fusion \(\lambda_j\) is Gaussian as well. This permits explicit false-alarm and detection probabilities
\[
P_{\mathrm f}^{(j)}=\Pr\{\lambda_j>\tau_j\mid x_j=-1\},\qquad
P_{\mathrm d}^{(j)}=\Pr\{\lambda_j>\tau_j\mid x_j=+1\},
\]
with thresholds selected either by Bayesian tuning or, as in the paper, by Neyman–Pearson constraints. A convergence condition is given for the sum-product linearization,
\[
|c_{jk}|<\frac{1}{\max_n|\mathcal N_n|-1},
\]
which ensures contraction. Simulations on a five-node network with \(K=100\) samples per node and a training window \(T=2500\) show that the optimal centralized linear fusion (“linOpt”) achieves the best \(P_D\), the decentralized one-hop optimization (“linProp” and “linPropB”) is close to it, and max-product and sum-product exhibit similar performance while remaining sensitive to the coupling-learning factor \(\zeta\) [1909.09402].

## 4. Evidence-based MAPRT in HAADF STEM

In HAADF STEM, MAPRT is used for atom-column detection and model-order selection. Competing hypotheses \(M_k\) represent images containing \(k\) projected atom columns. With a uniform prior over \(k\), the posterior satisfies
\[
p(M_k\mid y)=\frac{p(y\mid M_k)\,p(M_k)}{\sum_j p(y\mid M_j)\,p(M_j)},
\qquad
M^*=\arg\max_k p(M_k\mid y),
\]
so selection reduces to maximizing the marginal likelihood \(p(y\mid M_k)\). Pairwise MAPRT between two model orders is therefore a posterior-odds test, with threshold \(\tau=1\) for pure MAP selection [1902.05809].

The image model represents \(N=k\) columns by Gaussian peaks on a constant background:
\[
f_{kl}(\theta)=
\zeta+\sum_{n=1}^N
\eta_n\exp\!\left\{
-\frac{(x_k-\beta_{x,n})^2+(y_l-\beta_{y,n})^2}{2\rho_n^2}
\right\},
\]
with reduced-parameter variants for equal widths or equal widths and equal heights. The observation noise is physically Poisson, but the implementation adopts a normal approximation to Poisson with plug-in variance \(w_{kl}\), which yields
\[
p(w\mid \theta,N)=\frac{\exp(-\chi^2(\theta)/2)}
{\prod_{k,l}(2\pi w_{kl})^{1/2}},
\qquad
\chi^2(\theta)=\sum_{k,l}\frac{[w_{kl}-f_{kl}(\theta)]^2}{w_{kl}}.
\]
Maximum likelihood estimation therefore becomes weighted least squares, and the marginal likelihood is evaluated with a Laplace approximation around the ML estimate \(\hat \theta\) [1902.05809].

For the general unequal-width model, the practical posterior expression is
\[
p(N\mid w)\propto
N!(4\pi)^{2N}e^{-\chi^2_{\min}/2}
[\det(\nabla\nabla\chi^2)]^{-1/2}
\times
\Big[
(\beta_{x,\max}-\beta_{x,\min})
(\beta_{y,\max}-\beta_{y,\min})
(\eta_{\max}-\eta_{\min})
(\rho_{\max}-\rho_{\min})
\Big]^N,
\]
with analogous forms for the equal-width and equal-width/equal-height cases. The factors \(N!\), the Hessian determinant, and the prior-volume terms are the characteristic finite-sample corrections of this MAPRT. In “information criterion” form,
\[
-2\log p(N\mid w)
=
-2\log \hat L
+\log\det(\nabla\nabla\chi^2)
-2\log(N!)
-M\log(4\pi)
-2\log p(\theta\mid N)
+\mathrm{cst},
\]
and the paper shows that as \(W=K\times L\to\infty\), the leading terms reduce to \(-2\log \hat L+M\log W\), so the MAP rule is asymptotically equivalent to BIC while retaining curvature, labeling, and prior-volume penalties in finite samples [1902.05809].

Algorithmically, the method performs forward selection from \(N_0\) to \(N_{\max}\), adds peaks incrementally, tests many initial positions for each added peak, minimizes \(\chi^2\) under prior-range constraints, computes the Hessian determinant, and then evaluates \(p(N\mid w)\). In simulations of 1000 images of size \(12.5\ \text{\AA}\times 12.5\ \text{\AA}\) with pixel size \(0.25\ \text{\AA}\), 1–5 randomly positioned Au atoms, dose uniformly \(5\times 10^3\)–\(10^5\ e^{-}/\text{\AA}^2\), and a Gaussian model with equal widths and equal heights, the MAP rule detected the correct number most frequently at \(\mathrm{ICNR}=3.0\); AIC, GIC, and HQC tended to overfit, whereas BIC tended to underfit. At \(\mathrm{ICNR}=4.0\) and \(5.0\), all criteria improved and the differences narrowed [1902.05809].

A major by-product of this work is the integrated CNR,
\[
\mathrm{ICNR}=
\frac{2\pi \eta \rho^2}
{\sqrt{2\pi \eta \rho^2+\pi(3\rho)^2\zeta}},
\]
defined from the integrated signal and integrated background over a footprint of radius \(3\rho\), containing \(99.46\%\) of the Gaussian volume. The paper reports that ICNR correlates with atom detectability better than conventional SNR and CNR, and that columns with ICNR less than around \(5\) become challenging, with detection rate rapidly dropping below this value. This distinguishes the MAPRT from a purely formal model-selection device: it also supports quantitative detectability analysis in low-dose imaging [1902.05809].

## 5. EM-based adaptive radar detection

In adaptive radar, MAPRT is instantiated through a hybrid ML-Bayesian detector for a cell under test \(x\in\mathbb C^N\) and secondary target-free data \(D=\{z_k\}_{k=1}^K\). The hypotheses are
\[
H_0:\ x=n,\quad z_k=n_k,\qquad
H_1:\ x=\alpha s+n,\quad z_k=n_k,
\]
where \(n,n_k\sim\mathcal{CN}(0,\Sigma)\), \(s\) is a known steering vector, \(\alpha\in\mathbb C\) is an unknown target amplitude, and \(\Sigma\) is an unknown Hermitian positive definite covariance matrix. Defining
\[
S=\sum_{k=1}^K z_k z_k^H,
\]
the joint likelihoods are
\[
p(x,D\mid H_0,\Sigma)=
\left[\frac{1}{\pi^N\det\Sigma}\right]^{K+1}
\exp\{-\operatorname{Tr}[\Sigma^{-1}(x x^H+S)]\},
\]
\[
p(x,D\mid H_1,\alpha,\Sigma)=
\left[\frac{1}{\pi^N\det\Sigma}\right]^{K+1}
\exp\{-\operatorname{Tr}[\Sigma^{-1}((x-\alpha s)(x-\alpha s)^H+S)]\}.
\]
A binary latent variable \(c\in\{0,1\}\) indexes the hypothesis, with priors \(p_t=P(c=t)\) estimated inside the EM procedure rather than fixed externally [2503.02214].

The MAPRT statistic is the posterior ratio
\[
\Lambda_{\mathrm{MAP}}(x)=
\frac{P(H_1\mid x,D)}{P(H_0\mid x,D)}
=
\frac{\hat p_1 f_1(Z;\hat\alpha,\hat\Sigma)}
{\hat p_0 f_0(Z;\hat\Sigma)}
\underset{H_0}{\overset{H_1}{\gtrless}} \eta.
\]
Because EM computes the posterior responsibilities
\[
q_t(Z)=P(c=t\mid Z;\text{current parameters}),
\]
the same detector can be written as
\[
\Lambda_{\mathrm{MAP}}(x)=\frac{q_1(Z)}{q_0(Z)}.
\]
Its log form is
\[
\ell(x)=\log(\hat p_1/\hat p_0)
+x^H\hat\Sigma^{-1}x
-(x-\hat\alpha s)^H\hat\Sigma^{-1}(x-\hat\alpha s),
\]
or equivalently
\[
2\Re\{\hat\alpha\, s^H\hat\Sigma^{-1}x\}
-|\hat\alpha|^2 s^H\hat\Sigma^{-1}s
+\log(\hat p_1/\hat p_0).
\]
This makes the radar detector a MAPRT whose posterior terms are estimated jointly with nuisance parameters [2503.02214].

The EM iteration alternates between an E-step, which computes \(q_0(Z)\) and \(q_1(Z)\), and an M-step, which updates
\[
\hat p_t\leftarrow q_t(Z),
\]
\[
\hat\Sigma\leftarrow \frac{1}{K+1}
\Big[q_0\,x x^H+q_1\,(x-\hat\alpha s)(x-\hat\alpha s)^H+S\Big],
\]
\[
\hat\alpha\leftarrow
\frac{s^H\hat\Sigma^{-1}x}{s^H\hat\Sigma^{-1}s}.
\]
The detector is named EM-BML-D in the paper and is explicitly described as a “generalized MAP” rule because it compares posterior probabilities with ML-estimated nuisance parameters. Initialization with \(\hat p_0^{(0)}=\hat p_1^{(0)}=1/2\), \(\hat\Sigma^{(0)}=\frac1K\sum_{k=1}^K z_k z_k^H\), and the generalized matched-filter estimate for \(\hat\alpha^{(0)}\) yields a CFAR statistic with respect to the unknown \(\Sigma\) [2503.02214].

The CFAR proof is based on whitening:
\[
y=\Sigma^{-1/2}x,\qquad
\tilde s=\Sigma^{-1/2}s,\qquad
\tilde S=\Sigma^{-1/2}S\Sigma^{-1/2}\sim \mathcal W_{\mathbb C}(K,I_N).
\]
Under \(H_0\), the distribution of \(q_1/q_0\) depends only on whitened invariants and is therefore independent of \(\Sigma\). By induction over EM iterations, this invariance propagates to the final statistic, so a threshold \(\eta\) can be calibrated for a desired \(P_{FA}\) without knowing the interference covariance. The paper contrasts this detector with the GLRT and Kelly’s detector, emphasizing that MAPRT incorporates hypothesis priors and fuses information from both hypotheses through the responsibilities, whereas the GLRT remains purely ML-based [2503.02214].

Numerically, for \(N=8\), \(K=16\), \(\mathrm{CNR}=30\ \mathrm{dB}\), \(\rho=0.9\), and \(\mathrm{SCNR}=15\ \mathrm{dB}\), the proposed MAPRT achieved \(P_D>0.8\), versus approximately \(0.6\) for GLRT and approximately \(0.5\) for AMF. The paper also reports CFAR preservation across \(\mathrm{CNR}\) from \(30\) to \(110\ \mathrm{dB}\) and correlation \(\rho\) from \(0.5\) to \(0.9\), improved performance on MIT Lincoln Laboratory Phase One data, and robustness under mismatch similar to AMF. Empirically, about five EM iterations suffice, with average execution times for \(N=8,K=16\) over 1000 trials of approximately \(2.29\times 10^{-5}\ \mathrm{s}\) for GLRT, \(1.93\times 10^{-5}\ \mathrm{s}\) for AMF, \(3.81\times 10^{-4}\ \mathrm{s}\) for EM-BML-D5, and \(5.37\times 10^{-4}\ \mathrm{s}\) for EM-BML-D7 [2503.02214].

## 6. MAPRT and MAP testing in statistical inverse problems

For linear-Gaussian inverse problems on separable Hilbert spaces,
\[
Y=T u^\dagger+\sigma Z,
\]
with Gaussian prior \(U\sim N(m_0,C_0)\), the posterior is Gaussian \(N(m,C)\) with
\[
m=
m_0
+
C_0^{1/2}
\big(C_0^{1/2}T^*TC_0^{1/2}+\sigma^2 I\big)^{-1}
C_0^{1/2}T^*(y-Tm_0),
\]
\[
C=
\sigma^2 C_0^{1/2}
\big(C_0^{1/2}T^*TC_0^{1/2}+\sigma^2 I\big)^{-1}
C_0^{1/2}.
\]
The tested feature is \(g(u)=\langle \phi,u\rangle\), with hypotheses
\[
H:\ \langle \phi,u^\dagger\rangle\le 0,
\qquad
K:\ \langle \phi,u^\dagger\rangle>0.
\]
Since
\[
g\mid y\sim N(m_g,\sigma_g^2),
\qquad
m_g=\langle \phi,m\rangle,
\qquad
\sigma_g^2=\langle \phi,C\phi\rangle,
\]
the paper’s MAP mass test rejects \(H\) iff \(\Pi(\langle \phi,U\rangle>0\mid Y=y)>1/2\), equivalently iff \(m_g>0\) [2402.00686].

This mass test admits a linear-statistic representation. Defining
\[
\Phi_{\mathrm{MAP}}
=
T C_0^{1/2}
\big(C_0^{1/2}T^*TC_0^{1/2}+\sigma^2 I\big)^{-1}
C_0^{1/2}\phi,
\]
the decision rule becomes
\[
\langle Y,\Phi_{\mathrm{MAP}}\rangle_Y
>
\langle m_0, T^*\Phi_{\mathrm{MAP}}-\phi\rangle_X.
\]
The paper proves an exact Gaussian power formula,
\[
P_{u^\dagger}[\mathrm{MAP}(Y)=1]
=
Q\!\left(
\frac{
\langle u^\dagger,T^*\Phi_{\mathrm{MAP}}\rangle_X
-
\langle m_0,T^*\Phi_{\mathrm{MAP}}-\phi\rangle_X
}{
\sigma\|\Phi_{\mathrm{MAP}}\|_Y
}
\right),
\]
which yields both level and power once \(m_0\) or the threshold is calibrated. It also shows that, without a priori restrictions, the level is typically \(1\); under commuting assumptions and when \(\phi\) is an eigenvector of \(C_0^{1/2}T^*TC_0^{1/2}\), one can obtain exact level \(\alpha\), and the resulting test coincides with the uniformly most powerful unregularized test [2402.00686].

A posterior-density MAPRT can also be defined in this setting. Let
\[
T(y)=\frac{\sup_{u\in X_K}\pi(u\mid y)}{\sup_{u\in X_H}\pi(u\mid y)}.
\]
Because the posterior is Gaussian and the constrained maximizers are projections of \(m\) onto the half-space boundaries in the \(C^{-1}\)-norm, the ratio has the closed form
\[
T(y)=
\begin{cases}
\exp\!\big(m_g^2/(2\sigma_g^2)\big), & m_g\ge 0,\\[4pt]
\exp\!\big(-m_g^2/(2\sigma_g^2)\big), & m_g<0.
\end{cases}
\]
Hence rejection of \(H\) at threshold \(\tau>1\) is equivalent to
\[
m_g\ge 0
\quad\text{and}\quad
\frac{|m_g|}{\sigma_g}\ge \sqrt{2\log \tau}.
\]
The paper stresses that this posterior-density MAPRT is distinct from the MAP mass test that forms the main object of the frequentist analysis [2402.00686].

The regularization viewpoint is central. Under Assumption 2.1, the MAP mass test is a regularized test of the form
\[
\Psi_{\Phi,c}(Y)=\mathbf 1\{\langle Y,\Phi\rangle>c\},
\]
with \(\Phi=\Phi_{\mathrm{MAP}}\) and
\[
c(\Phi)=
\sigma\|\Phi\|_Y Q^{-1}(1-\alpha)+\|T^*\Phi-\phi\|_{V'}.
\]
Moreover, \(\Phi_{\mathrm{MAP}}\) is the Tikhonov–Phillips regularized solution of \(T^*\Phi=\phi\), or equivalently the solution of
\[
(TT^*+\sigma^2 V C_0^{-1}V^*)\Phi_{\mathrm{MAP}}=T\phi.
\]
Under the spectral source condition \(u^\dagger=(T^*T)^{\nu/2}w\), \(\|w\|_X\le \rho\), and priors \(C_0=\gamma^2(T^*T)^\mu\) with \(\mu\ge 1\), the paper derives explicit lower bounds on power and shows that for scalings \(\gamma=\gamma_0\sigma^\omega\), \(\omega\in(-\mu,1)\), the power tends to \(1\) as \(\sigma\to 0\). Numerical experiments in deconvolution, numerical differentiation, and the backward heat equation show conservative level control and strong power even in severely ill-posed settings, where unregularized testing becomes numerically infeasible [2402.00686].

Source: https://www.emergentmind.com/topics/maximum-a-posteriori-ratio-test-maprt