---
title: 'Partial Copula: Adjusted Dependence in Statistics'
url: https://www.emergentmind.com/topics/partial-copula
type: topic
---

# Partial Copula: Adjusted Dependence in Statistics

A partial copula is the copula of conditional probability integral transforms and therefore a dependence object for two variables after removing the influence of a conditioning variable through their full conditional distributions. For continuous random variables, if $U=F_{Y|X}(Y\mid X)$ and $V=F_{Z|X}(Z\mid X)$, then the partial copula of $(Y,Z)$ given $X$ is $C_{Y,Z;X}(u,v)=P(U\le u,V\le v)$; equivalently, it is the unconditional copula of $(U,V)$ and the average of the conditional copulas over the distribution of the conditioning variable [1511.06665]. In this sense it generalizes partial correlation from linear residual association to a rank-based, fully distributional notion of residual dependence, and it is central both to conditional independence testing and to simplified vine copula modeling [1101.4607].

## 1. Formal definition and construction

Let $(X,Y,Z)$ be continuous random variables with conditional distribution functions $F_{Y|X}(y|x)$ and $F_{Z|X}(z|x)$. The conditional probability integral transforms are
$$
U := F_{Y|X}(Y|X), \qquad V := F_{Z|X}(Z|X).
$$
Then $U$ and $V$ are uniform on $[0,1]$, and each is independent of $X$. The partial copula of $(Y,Z)$ given $X$ is the joint distribution of $(U,V)$:
$$
C_{Y,Z;X}(u,v):=P(U\le u,V\le v), \qquad (u,v)\in[0,1]^2.
$$
Equivalently, if $C_{Y,Z|X=x}$ denotes the conditional copula of $(Y,Z)$ given $X=x$, then
$$
C_{Y,Z;X}(u,v)=\int C_{Y,Z|X=x}(u,v)\,dF_X(x)=\mathbb{E}\!\left[C_{Y,Z|X}(u,v)\right].
$$
Thus the partial copula is the expected conditional copula, or an unconditional summary of a $x$-varying family of conditional copulas [1511.06665].

This construction is the copula-theoretic analogue of residualization. In the Gaussian case, if $\varepsilon_X$ and $\varepsilon_Y$ are residuals from linear regression on $Z$, then
$$
U=\Phi(\varepsilon_X/\sigma_X), \qquad V=\Phi(\varepsilon_Y/\sigma_Y),
$$
and the copula of $(U,V)$ is Gaussian with correlation equal to the classical Pearson partial correlation. For a trivariate normal vector, that correlation is
$$
\rho_{XY\cdot Z}=\frac{\rho_{XY}-\rho_{XZ}\rho_{YZ}}{\sqrt{(1-\rho_{XZ}^2)(1-\rho_{YZ}^2)}}.
$$
Hence, in the Gaussian setting, the partial copula reduces exactly to the Gaussian copula indexed by the familiar partial correlation parameter [1511.06665].

The same definition underlies the “partial copula transform” used in testing conditional independence. If $Y\ind Z\mid X$, then $U\ind V$. This makes conditional independence testable by first transforming the data with conditional distribution functions and then applying an ordinary independence test to the transformed pairs [1101.4607].

## 2. Interpretation as residual dependence and core properties

The central interpretation is “net-of-$X$” dependence. Because $U\perp X$ and $V\perp X$, the partial copula measures dependence between $Y$ and $Z$ after removing the influence of $X$ at the level of the full conditional distributions, not merely through linear projections. This gives a nonlinear analogue of partial correlation and explains why the partial copula can capture asymmetry, tail behavior, and other non-Gaussian dependence features [1510.06971].

A basic structural property is invariance under strictly monotone transformations of the target variables. If $\tilde X=T(X)$ with $T$ strictly increasing, then $F_{\tilde X|Z}(\tilde X|Z)=F_{X|Z}(X|Z)$, so the partial copula is unaffected. Recent work extends this to strictly increasing transformations of $X$ and $Y$ that may be $z$-specific: if $X'=T_x(X;Z)$ and $Y'=T_y(Y;Z)$ are strictly increasing in the first argument for each $z$, then
$$
F_{X'|Z}(X'|Z)=F_{X|Z}(X|Z), \qquad F_{Y'|Z}(Y'|Z)=F_{Y|Z}(Y|Z),
$$
and the partial copula is unchanged [2603.10941].

As a copula, the partial copula satisfies the Fréchet–Hoeffding bounds
$$
W(u,v)\le C^{P}(u,v)\le M(u,v),
$$
where $W(u,v)=\max(u+v-1,0)$ and $M(u,v)=\min(u,v)$. Because it is an average of copulas, pointwise order properties can transfer under sign-stable conditional dependence. If for all $z$, $C_{X,Y|Z=z}(u,v)\ge uv$, then $C^{P}(u,v)\ge uv$; similarly, if $C_{X,Y|Z=z}(u,v)\le uv$ for all $z$, then $C^{P}(u,v)\le uv$ [2603.10941].

Several associated dependence measures are defined directly from the partial copula. Partial Spearman’s rho and partial Kendall’s tau are
$$
\rho_P = 12 \int_0^1\!\!\int_0^1 \big(C^{P}(u,v)-uv\big)\,du\,dv,
$$
and
$$
\tau_P = 4 \int_{[0,1]^2} C^{P}(u,v)\,dC^{P}(u,v)-1.
$$
A key distinction is that partial Spearman’s rho equals the expected conditional Spearman’s rho, whereas partial Kendall’s tau does not, in general, equal the expected conditional Kendall’s tau [1511.06665]. Tail dependence coefficients also aggregate cleanly: if $\lambda_U^P$ and $\lambda_L^P$ denote the upper and lower tail dependence coefficients of $C^{P}$, then $\lambda_U^P=\mathbb{E}[\lambda_U(C_{XY|Z})]$ and $\lambda_L^P=\mathbb{E}[\lambda_L(C_{XY|Z})]$ [1511.06665].

## 3. Conditional copulas, simplifying assumption, and classes where equality holds

The conditional copula $C_{XY|Z=z}$ and the partial copula $C^{P}$ coincide under the simplifying assumption, namely when the conditional copula does not depend on the realized conditioning value:
$$
C_{XY|Z=z}(u,v)\equiv C_0(u,v).
$$
In that case,
$$
C^{P}(u,v)=C_0(u,v)=C_{XY|Z}(u,v).
$$
Thus the partial copula is exact when the conditional dependence structure is conditioning-invariant [1205.4844].

This equality holds for important model classes. The multivariate Normal and Student-$t$ copulas give rise to simplified pair copula constructions, so their conditional copulas do not depend on the conditioning values. In particular, for the Student-$t$ copula,
$$
C_{1,2\mid 3{:}d}(u_1,u_2)
=
C_t\Big(F_{t,1}^{-1}(u_1;\nu+d-2),F_{t,1}^{-1}(u_2;\nu+d-2);R_{A\mid B},\nu+d-2\Big),
$$
which is independent of the realized values of the conditioning variables. In the Archimedean class for $d\ge 4$, the only simplified pair copula construction is generated by the gamma Laplace transform
$$
\varphi_{\Gamma}(s;\theta)=\big(1+\theta s\big)_{+}^{-1/\theta}, \qquad \theta\ge -\frac{1}{d-1}.
$$
Conditioning then yields conditional generators
$$
\varphi_{m}(s;\theta)
=
\bigg(1+\frac{\theta}{m\theta+1}\,s\bigg)_{+}^{-(m\theta+1)/\theta},
$$
so parameters depend on the number of conditioned variables but not on the conditioning values themselves [1205.4844].

Outside these special families, the simplifying assumption generally fails. The trivariate Frank copula provides a standard example: the conditional copula of $(U_1,U_2\mid U_3=u_3)$ is Ali–Mikhail–Haq with parameter $\theta(\alpha,u_3)=1-e^{-\alpha u_3}$, so the dependence varies with $u_3$. The partial copula is then the mixture
$$
C_{12;3}(u_1,u_2)
=
\int_0^1
\frac{u_1u_2}{1-\big(1-e^{-\alpha u_3}\big)(1-u_1)(1-u_2)}
\,dP_{U_3}(u_3),
$$
which is not equal to the conditional copula unless $U_3$ is degenerate [1205.4844].

This distinction matters conceptually. The conditional copula describes dependence at a fixed conditioning level; the partial copula is an average over such levels. A plausible implication is that the partial copula is informative as a global adjusted summary, but it can mask heterogeneity when the conditional copula varies strongly with the conditioning variable.

## 4. Higher-order partial copulas and the partial vine copula

In vine copula constructions, the notion extends from a single conditioning variable to general conditioning sets. Consider a D-vine on $U_{1:d}$ with an edge $(i,i+j)$ and conditioning set $K=\{i+1,\dots,i+j-1\}$. Define the higher-order transforms
$$
U_{i|K}:=F_{i|K}(U_i\mid U_K), \qquad U_{i+j|K}:=F_{i+j|K}(U_{i+j}\mid U_K).
$$
Each is uniform on $[0,1]$. The $(j-1)$-th order partial copula on that edge is
$$
C_{i,i+j;K}(a,b):=P(U_{i|K}\le a, U_{i+j|K}\le b), \qquad (a,b)\in[0,1]^2.
$$
For $j=2$, this is the ordinary bivariate partial copula; for $j\ge 3$, it is a higher-order partial copula [1510.06971].

A crucial nuance is that these higher-order partial probability integral transforms are uniform but need not be independent of the full conditioning vector. The key lemma states:
$$
U_{k|K}\perp U_K \iff U_{k|K}=F_{k|K}(U_k|U_K) \text{ almost surely.}
$$
Hence independence from the conditioning vector characterizes coincidence with the true conditional probability integral transform [1510.06971].

The partial vine copula (PVC) is the simplified vine copula that assigns to each edge the corresponding partial copula. In a D-vine, the first tree uses the true bivariate margins, the second tree uses first-order partial copulas, and trees $j\ge 3$ use recursively defined higher-order partial copulas. Its density factorizes as
$$
c_{1:d}(u_{1:d})
=
\prod_{(i,j)} c_{i,i+j;K}(u_{i|K},u_{i+j|K}),
$$
with $K=\{i+1,\dots,i+j-1\}$ and $c_{i,i+j;K}$ the density of the relevant partial copula [1510.06971].

The PVC is treewise optimal but not globally optimal in general. Sequential minimization of the Kullback–Leibler divergence over trees yields the PVC at each tree conditional on the previous ones, which explains why stepwise estimators converge to it. However, when the simplifying assumption fails, the PVC need not minimize the global Kullback–Leibler divergence over all simplified vine copulas; the best global approximation may instead be a vine pseudo-copula, which allows non-uniform pseudo-observations in upper trees [1510.06971].

Later analytic work sharpened the approximation picture. In dimension three, the family of simplified copulas is dense in the uniform metric $d_\infty$, but nowhere dense with respect to $D_1$, weak conditional convergence, total variation, and Kullback–Leibler divergence. Within this framework, the partial vine copula is never the optimal simplified copula approximation under $d_\infty$ for a non-simplified copula, and the mapping that assigns a copula its partial vine copula is discontinuous with respect to $d_\infty$ at non-simplified copulas [2006.06542]. This suggests that the PVC is best understood as a practically important sequential target rather than as a universal projection operator.

## 5. Estimation, testing, and inference

A standard estimation strategy is two-step. First estimate the conditional distribution functions $F_{X|Z}$ and $F_{Y|Z}$; then compute pseudo-observations
$$
\hat U_i=\hat F_{X|Z}(X_i|Z_i), \qquad \hat V_i=\hat F_{Y|Z}(Y_i|Z_i),
$$
and estimate the copula of $(\hat U_i,\hat V_i)$ either nonparametrically, semiparametrically, or parametrically. Suggested first-stage methods include kernel or local-likelihood regression for conditional distribution functions, as well as parametric or semiparametric models such as GLMs, quantile regression, location–scale models, and ARMA–GARCH in time series [1511.06665].

For conditional independence testing, the core procedure is to transform the sample using estimated conditional distribution functions and then apply an ordinary independence test to the transformed pairs. Kernel estimators of the form
$$
\hat F_{Y|X}(y \mid x)
=
\frac{\sum_{i=1}^n K\!\left(\frac{X_i - x}{h_n}\right) \mathbf{1}\{Y_i \le y\}}
     {\sum_{i=1}^n K\!\left(\frac{X_i - x}{h_n}\right)},
\qquad
\hat F_{Z|X}(z \mid x)
=
\frac{\sum_{i=1}^n K\!\left(\frac{X_i - x}{h_n}\right) \mathbf{1}\{Z_i \le z\}}
     {\sum_{i=1}^n K\!\left(\frac{X_i - x}{h_n}\right)}
$$
are one option for Euclidean covariates. The resulting independence statistics can be covariance, Kendall’s tau, Hoeffding’s $\Delta$, Bergsma’s $\kappa$, or $\tau^*$. Under mild conditions, the estimation effect vanishes asymptotically:
$$
T_n(\hat U,\hat V)-T_n(U,V)=o_p(1).
$$
Thus the asymptotic null and alternative behavior is the same as if the true transformed variables were observed [1101.4607].

For vine copulas, estimation practice is typically sequential. Parametric and nonparametric simplified vine estimators are usually fitted tree-by-tree. Under mild regularity conditions, when the simplifying assumption fails, stepwise estimators converge in probability to the PVC, while joint maximum likelihood converges to pseudo-true parameters minimizing the global Kullback–Leibler divergence. In general, these probability limits differ. This discrepancy is both computational and statistical: the combinatorial search for joint optimization grows as $K^{d(d-1)/2}$ for $K$ candidate families, versus $K\cdot d(d-1)/2$ in a sequential scheme [1510.06971].

A recurrent practical difficulty is the first-stage estimation of conditional distributions. Nonparametric estimation suffers from the curse of dimensionality for high-dimensional conditioning variables, and feasible nonparametric procedures typically require low-dimensional conditioning sets or structural assumptions such as additivity or single-index structure [1511.06665].

## 6. Limitations, misconceptions, examples, and recent developments

A fundamental limitation is that conditional independence and partial independence are not equivalent. If $X\ind Y\mid Z$, then $C^{P}=\Pi$, the independence copula. The converse fails: $C^{P}=\Pi$ does not imply $X\ind Y\mid Z$. The standard counterexample is the exchangeable trivariate Farlie–Gumbel–Morgenstern copula, for which the conditional copula varies with $z$ while the partial copula equals $\Pi$ [1511.06665]. The same phenomenon appears in recent covariate-adjusted formulations: opposite-signed conditional dependence can cancel in expectation, leaving an independence partial copula although conditional dependence is nonzero for most conditioning values [2603.10941].

This non-equivalence also separates the partial copula from partial correlation. In the example with $Z\sim N(0,1)$ and conditionally independent variables
$$
Y_1=-1+Z^2+\varepsilon_1, \qquad Y_2=-1+Z^2+\varepsilon_2,
$$
with independent $\varepsilon_i$ of variance $\sigma^2$, one has $Y_1\ind Y_2\mid Z$ and therefore $C^{P}=\Pi$, yet the Pearson partial correlation converges to $1$ as $\sigma^2\to 0$ [1511.06665]. This is one of the clearest demonstrations that the partial copula is not merely a rank-transformed variant of linear partial correlation.

The partial copula also lacks several closure and optimality properties one might expect. Even if the joint copula of $(X,Y,Z)$ is Archimedean, the partial copula need not be Archimedean; the Frank example shows associativity can fail. Likewise, if conditional copulas belong to a parametric family $C^F(\cdot;\theta(z))$, there need not exist a fixed $\theta$ such that the partial copula remains in the same family [1511.06665]. In $L^2$, the partial copula is generally not the projector onto the conditional copula process; the relevant $L^2$ projection is $\mathbb{E}[C_{XY|Z}(U,V|Z)\mid U,V]$, which differs in general from $\mathbb{E}[C_{XY|Z}(U,V|Z)]$ [1511.06665].

In multivariate vine settings, further interpretive caution is required. Partial independence on an upper tree edge does not imply conditional independence, nor vice versa, for higher-order partial copulas. The paper on the partial vine copula states that conditional independence and higher-order partial independence are generally unrelated for $j\ge 2$ [1510.06971]. Moreover, examples show that replacing true conditional copulas by partial copulas can induce spurious dependence in upper trees; in a five-dimensional FGM construction, true conditional copulas in the third and fourth trees were independence copulas, but the PVC induced non-trivial upper-tree copulas. This can reduce total Kullback–Leibler divergence, illustrating why the PVC may work well empirically even when the simplifying assumption fails [1510.06971].

Recent work has repositioned the partial copula as a covariate-adjusted dependence representation with direct implications for sign recovery. Because partial Spearman’s rho equals the expected conditional Spearman’s rho, if conditional dependence has constant sign across $z$, then the sign of the partial rank measure reflects that adjusted sign. The same work argues that this can clarify Simpson’s paradox settings: marginal association may have a different sign from the adjusted association, whereas partial-copula-based rank measures align with the adjusted sign when the usual adjustment assumptions hold and conditional dependence does not switch sign [2603.10941]. A plausible implication is that partial copulas are especially informative when the aim is global adjusted dependence rather than local, $z$-specific heterogeneity.

Across these developments, the partial copula has acquired a dual status. It is both a rigorous dependence measure—defined by conditional probability integral transforms and endowed with invariance and averaging properties—and a modeling device that targets tractable approximations in simplified vine constructions. Its usefulness is greatest when interpreted as an adjusted, global summary of conditional dependence rather than as a substitute for the full family of conditional copulas.

Source: https://www.emergentmind.com/topics/partial-copula