---
title: Real Log Canonical Threshold
url: https://www.emergentmind.com/topics/real-log-canonical-threshold
type: topic
---

# Real Log Canonical Threshold

The real log canonical threshold (RLCT) is a birational and analytic invariant measuring the singularity of real-analytic and algebraic structures, and plays a central role in both real algebraic geometry and singular learning theory. In the context of statistics and machine learning, the RLCT (also called the learning coefficient) governs the leading behavior of marginal likelihood and generalization error in singular models, allowing for refined asymptotic model comparison well beyond classical regular cases.

## 1. Analytic and Geometric Definition

Given a real-analytic function $f \colon \mathbb{R}^d \to \mathbb{R}$, the RLCT, denoted $\lambda_R(f)$, is defined via the analysis of the integrability of $|f|^{-p}$ near its zero locus, or equivalently through resolution of singularities:
\[
\lambda_R(f) = \inf_{P,\,i:\ a_i>0} \frac{b_i+1}{a_i}
\]
where for a real log resolution $\sigma: U \to \mathbb{R}^d$, locally $f\circ \sigma(u) = \eta(u)\, u_1^{a_1}\cdots u_d^{a_d}$ and $\det D\sigma(u) = \eta'(u)\, u_1^{b_1}\cdots u_d^{b_d}$ for non-vanishing analytic functions $\eta,\eta'$, $a_i, b_i \in \mathbb{N}$ for $i=1,\ldots,d$ [2411.13392].

The multiplicity $m_R(f)$ is defined as the maximal number of exponents $(b_i+1)/a_i$ achieving the minimum at any point. This pair $(\lambda_R, m_R)$ fully characterizes the leading singular behavior of sharp analytic volume or zeta integrals around the singular locus:
\[
V(\epsilon) = \int_{|f(w)| \leq \epsilon} \varphi(w)\, dw \sim C\, \epsilon^{\lambda_R}(-\ln \epsilon)^{m_R-1}
\]
as $\epsilon \to 0$ for compact $W\subset\mathbb{R}^d$, analytic $\varphi>0$ [2411.13392].

## 2. RLCT in Complex and Real Settings

The RLCT is an extension of the log canonical threshold (lct) defined in the complex algebraic context, often calculated from complex resolutions of singularities. For $f\in\mathbb{R}[x_1, \ldots, x_d]$, the complex lct $\lambda_C(f)$ is calculated via divisorial data from a complex log resolution $\rho: X \to \mathbb{C}^d$:
\[
\lambda_C(f) = \min_{1 \leq i \leq m} \frac{b_i+1}{a_i}
\]
where $b_i, a_i$ are the discrepancy and multiplicity coefficients along exceptional divisors $E_i$ [2411.13392].

When the log resolution is defined over $\mathbb{R}$ and each exceptional divisor $E_i$ meets real points, $\lambda_R(f) = \lambda_C(f)$, with the multiplicities also coinciding under mild conditions [2411.13392].

## 3. RLCT in Singular Learning Theory

The RLCT, also called the learning coefficient, is the key invariant dictating the small-sample asymptotics for both Bayesian evidence and generalization error in singular or non-identifiable statistical models. For a statistical model $p(x|w)$ and true data distribution $q(x)$, define the KL distance function:
\[
K(w) = \int q(x) \log\frac{q(x)}{p(x|w)}\, dx
\]
The RLCT $(\lambda, m)$ is then determined by the leading pole $- \lambda$ (order $m$) of the zeta function $\zeta(z) = \int_\Omega K(w)^z \varphi(w) dw$ [1906.01341, 2303.11619].

For data of size $n$, the log marginal likelihood (free energy) asymptotically satisfies (Watanabe's Main Formula II):
\[
-\log L(n) = n L_0 + \lambda \log n - (m-1)\log \log n + O_p(1)
\]
In regular cases $\lambda = m = d/2$; for singular models, $\lambda$ is typically smaller, yielding weaker Bayesian Occam penalties [1906.01341, 2303.11619].

## 4. Combinatorial Formulas: Hyperplane Arrangements

For $f(x) = \prod_{j=1}^n L_j(x)^{s_j}$, $L_j$ linear, $s_j > 0$ (not necessarily reduced), set $\mathcal{B} = \{H_j = \{L_j=0\}\}$ the arrangement, and let $L(\mathcal{B})$ be its intersection poset (all proper nontrivial intersections $W \neq \mathbb{R}^d$):
- $\mathrm{codim}\, W = d - \dim W$
- $s(W) = \sum_{j: L_j|_W \equiv 0} s_j$

Explicitly:
\[
\lambda_R(f) = \min_{W \in L(\mathcal{B})} \frac{\mathrm{codim}\, W}{s(W)}
\]
\[
m_R(f) = \max_{\text{chains } W_0 \subsetneq \ldots \subsetneq W_{d-1}} |\{i : \mathrm{codim}\, W_i/s(W_i) = \lambda_R(f)\}|
\]
If $f$ has real coefficients, $\lambda_R(f) = \lambda_C(f)$ and $m_R(f) = m_C(f)$ [2411.13392].

**Examples:**
- In $\mathbb{R}^2$, for $f=\prod_{i=1}^n L_i^{s_i}$ ($s_1 \leq \cdots \leq s_n$), $\lambda_R = \min\{1/s_n, 2/\sum_i s_i\}$.
- For $f=x\,y^2\,z^2\,(x+y+z)$ in $\mathbb{R}^3$, $\lambda_R = 1/2$, $m_R = 3$ [2411.13392].

## 5. Algorithmic and Blow-up Approaches

Computation of RLCTs for general polynomial models hinges on real-analytic resolution of singularities. For sum-of-products (sop) polynomials, dedicated blow-up algorithms—iteratively replacing singular charts via coordinate blow-ups—yield normal or locally normal crossing forms [2303.11619]. For binomials $f = w^p + w^q$ (with $g(w)=w^{s-1}$), critical poles can be computed explicitly:
\[
1/\lambda = \max\left\{ \max_i \{\min(\rho_i^l, \rho_i^r)\},\, \max_{i \neq j} \frac{\nu_i^l \nu_j^r - \rho_i^l \rho_j^r}{\nu_i^l + \nu_j^r - \rho_i^l - \rho_j^r} \right\}
\]
where $\rho_i^l = p_i/s_i$, $\rho_i^r = q_i/s_i$, $\nu_i^l = (p_i + r_i)/s_i$, etc. For higher terms ($n>2$), a linear programming simplex bound provides $1/\lambda \ge \max_{\alpha \in \Delta} \min_j \sum_i \alpha_i (a_{ij}/s_i)$ [2303.11619].

For hyperplane arrangements, the combinatorial algorithm computes all intersection flats, their codimension and $s(W)$, and determines $(\lambda_R, m_R)$ through inclusion poset chains [2411.13392].

## 6. RLCT at Non-Singular Points and Upper Bounds

At non-singular points of the true parameter manifold in statistical models, RLCT is immediately computable. Under analytic assumptions and after a variable split $(\theta, \tau)$, assumptions (feasibility, independence, vanishing order) yield
\[
\lambda_0 = \frac{d_1 - r + r m}{2m}
\]
where $d_1 =$ effective parameter dimension, $r$ counts directions with quadratic expansion, and $m$ minimal order for higher directions [2408.13030].

For global estimation in singular models, the RLCT is the minimum over non-singular points, with the above formula providing an effective upper bound, guiding practical Bayesian model comparison [2408.13030, 2303.05731].

## 7. Statistical Estimation and Model Selection

In practice, analytic RLCTs are known only in restricted cases; for broader models, Monte Carlo-based estimators exploit the relationship between the variance of log-likelihoods under tempered posteriors and the RLCT:
\[
\hat{\lambda}_V^{(1)} = t^2 \cdot \operatorname{Var}_t[\log p(D|w)]
\]
for $t \sim 1/\log n$, and averaging over multiple simulated datasets increases accuracy. This enables the singular Bayesian information criterion (sBIC) to be widely deployed as WsBIC, using $n^{-\hat{\lambda}_V}$ as the complexity penalty, which empirically outperforms WBIC/BIC in mixtures, reduced-rank, and other singular models [1906.01341].

## Table: RLCT Calculation Methods and Their Domains

| Method                        | Domain/Applicability                    | Reference         |
|-------------------------------|-----------------------------------------|-------------------|
| Real log resolution           | General analytic or algebraic $f$       | [2411.13392]      |
| Combinatorial poset (arrangements) | Hyperplane arrangements                   | [2411.13392]      |
| Blow-up algorithm (sop polynomials) | Sum-of-products, binomial polynomials    | [2303.11619]      |
| Non-singular point expansion  | Regular points in parametric models     | [2408.13030]      |
| Monte Carlo (variance-based)  | General statistical models              | [1906.01341]      |

The RLCT is thus a unifying invariant at the interface of algebraic geometry, singularity theory, and Bayesian statistics, encoding local analytic complexity and model selection criteria in both theoretical and computational practice. Calculating or estimating the RLCT, exactly or via bounds, is essential for the asymptotic characterization of evidence and generalization in nonregular models, including applications in high-dimensional arrangements, tensor decompositions, and non-identifiable latent-variable models [2411.13392, 2303.05731].

Source: https://www.emergentmind.com/topics/real-log-canonical-threshold