Tutorials › Statistical Learning (Graduate) › Consistency of Response-Adaptive Partitions

Randomized Forests and Consistency · Tutorial 855 of 989

Consistency of Response-Adaptive Partitions

Learn when response-guided partitioning preserves regression consistency, and how sample splitting or complexity control limits selection effects.

Advanced 10 min read

What You'll Learn

  • Define response-adaptive partitions and distinguish them from honest partitions.
  • Explain how selecting a partition using responses can make training leaf means optimistic.
  • Apply a finite-class consistency condition with minimum cell probability and controlled partition complexity.
  • Derive the conditional error bound for leaf means estimated on an independent sample.
  • Identify why localization and leaf occupancy alone do not remove selection effects.

Why Response Adaptation Needs Its Own Condition

A partitioning method may use the observed responses to decide where to split. That can make its leaves well adapted to the training sample, but it also means the partition and the leaf averages are statistically linked. A split chosen because it separates unusually high and low responses can look predictive even when the apparent pattern is noise.

Consequently, conditions that control cell diameter and occupancy do not, by themselves, control the effect of selecting boundaries using the same responses later averaged within those cells. Two broad remedies are to estimate leaf means with responses not used to choose the partition, or to restrict the collection of partitions enough that empirical averages behave uniformly over all candidates.

Definition: A response-adaptive partition is a partition \(\hat{\Pi}_n\) of the predictor space whose choice can depend on the observed responses \(Y_1,\ldots,Y_n\), as well as on the observed predictors and auxiliary randomness. The associated regression estimate at \(x\) is the average of responses whose predictors lie in the selected cell containing \(x\). A partition is honest for that estimate if the responses used to choose the partition are separate from the responses used to calculate its leaf averages.

Here \(X\) is the predictor, \(Y\) is the response, and \(m(x)=E[Y\mid X=x]\) is the regression function. For a measurable cell \(A\) with \(P(X\in A)>0\), write \(\mu_A=E[Y\mid X\in A]\). The population cell predictor for a partition \(\Pi\) is \(m_\Pi(x)=\mu_A\) when \(x\in A\in\Pi\). Its approximation error is \(E[(m_\Pi(X)-m(X))^2]\). This error describes the coarseness of the partition, not the extra error from estimating the cell means with data.

How Selection Creates Optimism

If a partition is fixed before the responses are observed, a leaf average is an ordinary sample average within a fixed region. With response adaptation, the selected region may instead be one where the observed responses happen to be unusually large or small. The resulting leaf average can differ systematically from its population counterpart. A small training error is therefore not enough to establish prediction consistency.

Worked Example: A Search Over Noise Splits

Let \(n=6\), and suppose the observed responses are balanced: three are \(+1\) and three are \(-1\). For each of \(p=32\) candidate predictors, suppose the six observed predictor values are independent Rademacher variables, independent of the responses. A candidate split sends observations with predictor value \(+1\) to one leaf and those with value \(-1\) to the other.

For any one candidate, there are \(2^6=64\) equally likely predictor-value vectors. Exactly two vectors—the response vector and its negative—put all three positive responses in one leaf and all three negative responses in the other. Thus the probability that a candidate perfectly separates the responses is \(2/64=1/32\). Because the candidate predictors are independent, the probability that at least one of the 32 candidates does so is

$$ 1-\left(1-\frac{1}{32}\right)^{32} =1-\left(\frac{31}{32}\right)^{32} \approx 0.638. $$

On this event, choosing a perfectly separating split gives training error zero: each leaf average equals its response value. But the candidate predictor and response are independent for a new observation. The prediction is \(+1\) or \(-1\), while the new response is an independent \(+1\) or \(-1\). Its test squared error has expectation

$$ E[(Y_{\mathrm{new}}-\hat m(X_{\mathrm{new}}))^2] =E[Y_{\mathrm{new}}^2]+E[\hat m(X_{\mathrm{new}})^2] -2E[Y_{\mathrm{new}}\hat m(X_{\mathrm{new}})] =1+1-0=2. $$

The zero predictor has test squared error \(E[Y_{\mathrm{new}}^2]=1\). Thus, on a substantial-probability selection event, a split with perfect training fit has worse test error than the zero predictor. The issue is not that the leaves are too large; it is that the same responses were searched over and then averaged.

A Sufficient Condition: Uniform Control Over Candidates

A direct way to allow response-guided selection is to prove that empirical leaf averages are close to population leaf means simultaneously over the candidate cells. Then the selected partition inherits the same control, even though its selection used the responses. The following finite-class condition is sufficient. Its minimum cell-probability assumption is deliberately explicit: very rare cells cannot have reliably estimated means without additional structure.

Theorem: For each \(n\), let \(\mathcal P_n\) be a fixed finite collection of measurable partitions, and let \(J_n\) be the total number of cells across those partitions. Suppose \((X_i,Y_i)_{i=1}^n\) are independent and identically distributed, \(|Y_i|\le B\) almost surely for a finite constant \(B\), and every candidate cell \(A\) has \(P(X\in A)\ge q_n>0\). Let \(\hat\Pi_n\in\mathcal P_n\) be any measurable partition selected using the full sample, and estimate each selected leaf mean using that same sample. If \[ \frac{\log(2J_n)}{nq_n}\longrightarrow 0 \quad\text{and}\quad E[(m_{\hat\Pi_n}(X)-m(X))^2]\xrightarrow{p}0, \] then the integrated prediction risk \[ E[(\hat m_n(X_{\mathrm{new}})-m(X_{\mathrm{new}}))^2\mid (X_i,Y_i)_{i=1}^n] \xrightarrow{p}0, \] where \(X_{\mathrm{new}}\) is an independent draw from the predictor distribution. The expectation in the approximation condition is with respect to a fresh \(X\), conditional on the selected partition.

The approximation condition is important: controlling estimation error cannot make a persistently poor partition approximate \(m\). The theorem instead isolates the selection issue. A larger candidate collection raises the uniform estimation cost through \(\log J_n\); smaller cell probabilities raise it through \(1/q_n\).

Proof. For a candidate cell \(A\), set \(p_A=P(X\in A)\) and \(\mu_A=E[Y\mid X\in A]\). Let \(P_n\) denote an empirical average over the sample. Bernstein's inequality, applied to \(\mathbf 1_{\{X\in A\}}\) and \(Y\mathbf 1_{\{X\in A\}}\), gives deviations of orders \(\sqrt{p_A t/n}+t/n\) and \(B\sqrt{p_A t/n}+Bt/n\), respectively. More precisely, with probability at least \(1-\delta\), simultaneously for every candidate cell,

$$ \begin{aligned} |P_n\mathbf 1_A-p_A| &\le \sqrt{\frac{2p_A t}{n}}+\frac{2t}{3n},\\ |P_n(Y\mathbf 1_A)-E[Y\mathbf 1_A]| &\le B\sqrt{\frac{2p_A t}{n}}+\frac{2Bt}{3n}, \end{aligned} \qquad t=\log\left(\frac{4J_n}{\delta}\right). $$

Indeed, each inequality fails with probability at most \(2e^{-t}\); summing these failure probabilities over the two quantities for \(J_n\) cells gives at most \(4J_ne^{-t}=\delta\). Since \(\log(2J_n)/(nq_n)\to0\), for fixed \(\delta\) we have \(t/(nq_n)\to0\). Thus, eventually, the first deviation bound is at most \(p_A/2\) for every cell, so \(P_n\mathbf 1_A\ge p_A/2>0\).

The empirical leaf average on \(A\) is \(P_n(Y\mathbf 1_A)/P_n\mathbf 1_A\). Since \(|\mu_A|\le B\) and \(E[Y\mathbf 1_A]=p_A\mu_A\), its difference from \(\mu_A\) is bounded by

$$ \left|\frac{P_n(Y\mathbf 1_A)}{P_n\mathbf 1_A}-\mu_A\right| \le \frac{|P_n(Y\mathbf 1_A)-E[Y\mathbf 1_A]| +B|P_n\mathbf 1_A-p_A|}{P_n\mathbf 1_A}. $$

Substituting the simultaneous deviations and \(P_n\mathbf 1_A\ge p_A/2\) shows that, for a universal constant \(C\), the absolute difference is at most \(CB\{\sqrt{t/(nq_n)}+t/(nq_n)\}\) on every candidate cell. This bound holds simultaneously before the data-dependent choice is made, so it also holds on every cell of \(\hat\Pi_n\). Call the bound \(r_n\). For a fresh \(X\), decompose the prediction error into the empirical leaf-mean error and the population approximation error. The inequality \((a+b)^2\le 2a^2+2b^2\) gives

$$ E[(\hat m_n(X)-m(X))^2\mid\text{sample}] \le 2r_n^2+2E[(m_{\hat\Pi_n}(X)-m(X))^2]. $$

For fixed \(\delta\), \(r_n\to0\), and the second term tends to zero in probability by assumption. The deviation event has probability at least \(1-\delta\); since \(\delta>0\) is arbitrary, the risk converges to zero in probability. \(\square\)

Worked Example: Checking a Finite-Class Condition

Suppose a fitting procedure can choose among five candidate partitions, each with ten cells. Then \(J_n=50\). Assume every candidate cell has probability at least \(q_n=0.02\), responses are bounded by \(B\), and the sample size is \(n=10{,}000\). For a target failure probability \(\delta=0.05\), the proof uses

$$ t=\log\left(\frac{4J_n}{\delta}\right) =\log(4000)\approx 8.294, \qquad \frac{t}{nq_n}=\frac{8.294}{200}\approx 0.0415. $$

The cell-probability requirement corresponds to an expected count of at least \(nq_n=200\) observations in each candidate cell. The uniform mean-error scale in the proof is \(CB(\sqrt{0.0415}+0.0415)\), which tends to zero as the sample grows with \(J_n\) controlled and \(nq_n\) increasing. To conclude consistency, the population approximation error of the selected partition must also tend to zero; the calculation of the estimation scale alone does not establish that condition.

A Sufficient Condition: Honest Estimation

Honesty avoids direct reuse of responses by separating the partition-selection sample from the leaf-estimation sample. Conditional on the selection sample, the chosen partition is fixed, and the estimation responses are independent of its selection. The next bound makes the remaining occupancy requirement explicit.

Result: Suppose a partition is chosen using data independent of an estimation sample of size \(m\). Let \(A\) be one of its cells, with \(p=P(X\in A)>0\) and \(|Y|\le B\). Estimate \(\mu_A=E[Y\mid X\in A]\) by the average of estimation responses in \(A\), using zero if the cell is empty. Conditional on the chosen partition, the mean squared error of this cell estimate is at most \[ B^2\left(e^{-mp}+\frac{2}{(m+1)p}\right). \]

Derivation. Conditional on the partition, the number \(N\) of estimation observations in \(A\) has a \(\operatorname{Binomial}(m,p)\) distribution. Given \(N=k>0\), the observations in the cell are independent draws from the conditional distribution given \(X\in A\). Their average is unbiased for \(\mu_A\) and has variance at most \(B^2/k\), since \(\operatorname{Var}(Y\mid X\in A)\le B^2\). When \(N=0\), the squared error of the zero default is \(\mu_A^2\le B^2\). Therefore the mean squared error is at most

$$ B^2P(N=0)+B^2E\left[\frac{\mathbf 1_{\{N>0\}}}{N}\right]. $$

Now \(P(N=0)=(1-p)^m\le e^{-mp}\). For integers \(k\ge1\), \(1/k\le2/(k+1)\). Also, direct expansion of the binomial probabilities gives

$$ E\left[\frac{1}{N+1}\right] =\frac{1-(1-p)^{m+1}}{(m+1)p} \le \frac{1}{(m+1)p}. $$

It follows that \(E[\mathbf 1_{\{N>0\}}/N]\le2/((m+1)p)\), proving the bound. \(\square\)

Worked Example: An Honest Response-Guided Split

Suppose a structure sample selects a split using responses, and an independent estimation sample of size \(m=20\) is used to calculate its leaf averages. Conditional on the structure sample, consider a selected leaf \(A\) with \(p=P(X\in A)=0.5\) and \(Y\mid X\in A\sim\operatorname{Bernoulli}(0.2)\). If the estimation sample happens to contain \(N=10\) observations in that leaf, the leaf average is unbiased for \(0.2\), and its conditional variance is

$$ \frac{0.2(1-0.2)}{10}=\frac{0.16}{10}=0.016. $$

The probability that this leaf receives no estimation observations is \(P(N=0)=(1-0.5)^{20}=2^{-20}\), about \(0.000000954\). Thus the response-guided choice of the split does not bias the independent-sample leaf average, but the estimate still depends on adequate estimation-sample occupancy. To obtain overall consistency, the selected cells must also approximate the regression function as the sample grows.

What the Conditions Do—and Do Not—Guarantee

Honesty and complexity control address different parts of the problem. Honesty removes the direct link between the responses used to select a partition and those used to estimate its means. Uniform control allows that link, but requires all candidate-cell means to be sufficiently well estimated at once. Neither remedy alone guarantees that the selected partition has small approximation error, and neither can compensate for cells that remain too sparsely populated.

The finite-class theorem is one concrete complexity restriction, not a claim that every response-adaptive partition is consistent. For a large or continuous candidate collection, a corresponding argument needs a suitable uniform concentration result for that collection. Merely counting leaves in the final fitted tree does not count all the alternative partitions searched to obtain it.

Common Mistake

  • Wrong reasoning: “The leaf diameters shrink and the fitted leaves contain many observations, so the response-adaptive estimator is consistent.”
  • Correction: Those facts can control approximation and occupancy, but do not by themselves control selection based on the same responses. Also establish honesty or a uniform estimation bound over the candidate partitions, and require the selected partition's population approximation error to vanish.
Key takeaway: Response-guided boundaries create a selection problem because the same responses can determine both where leaves are placed and what their averages are. Honesty separates those roles; complexity control bounds the selection effect uniformly. Consistency still requires both reliable leaf estimation and vanishing population approximation error.

Check Your Understanding

Use the definitions and conditions above to answer these questions.

  1. Why does a small training error fail to establish consistency when the partition was selected by searching over response-guided splits?
  2. In the finite-class theorem, what roles do \(J_n\) and \(q_n\) play in controlling the leaf-mean estimation error?
  3. Why does honesty make a selected partition fixed for the purpose of analyzing the estimation-sample leaf mean?
  4. Which two distinct sources of error must vanish for the finite-class theorem to yield consistency?
  5. What can go wrong if a leaf has a probability that is too small relative to the estimation-sample size?
Clinical Biostats