Tutorials › Multivariate Statistics (Graduate) › Copulas and Dependence Beyond Marginals

Multivariate Laws and Transformations · Tutorial 68 of 975

Copulas and Dependence Beyond Marginals

Learn to represent a bivariate distribution using its continuous margins and copula, and see how different dependence structures can produce the same marginal distributions.

Beginner 10 min read

What You'll Learn

  • Define a bivariate copula and identify its uniform marginal distributions.
  • State and prove Sklar’s representation for continuous margins.
  • Construct independence, countermonotone, and mixture copulas.
  • Compare joint probabilities and covariances while holding margins fixed.
  • Derive the joint density from a copula density under smoothness conditions.
  • Recognize why marginal distributions alone do not determine a joint law.

Why Marginal Distributions Do Not Determine a Joint Law

Suppose two clinical measurements each have an exponential distribution. Those marginal distributions describe each measurement separately, but they do not say whether large values tend to occur together, whether one tends to be large when the other is small, or whether the measurements are independent. A copula records this missing part of the joint distribution: its dependence structure, separated from its margins.

This separation is useful when a marginal model is well understood but dependence remains uncertain. It also makes comparisons precise: we can hold the margins fixed and change the copula, then see how joint probabilities or other summaries change.

Definition: A bivariate copula is a joint cumulative distribution function \(C\) on \([0,1]^2\) whose two marginal distributions are uniform on \([0,1]\). In particular, \(C(u,1)=u\), \(C(1,v)=v\), and \(C(u,0)=C(0,v)=0\). For every \(0\leq u_1\leq u_2\leq 1\) and \(0\leq v_1\leq v_2\leq 1\), it also satisfies $$ C(u_2,v_2)-C(u_1,v_2)-C(u_2,v_1)+C(u_1,v_1)\geq 0. $$ This last condition says that every rectangle has nonnegative probability.

The two most familiar copulas are the independence copula \(C(u,v)=uv\) and the comonotone copula \(C(u,v)=\min(u,v)\). The former describes independent uniform variables; the latter describes two equal uniform variables. They have the same uniform margins, but very different joint behavior.

Sklar’s Representation

Let \(X\) and \(Y\) have joint cumulative distribution function \(H\), and marginal cumulative distribution functions \(F\) and \(G\). Sklar’s theorem identifies a copula that connects these distributions.

Theorem (Sklar): For any bivariate joint distribution \(H\) with marginal distributions \(F\) and \(G\), there is a copula \(C\) such that $$ H(x,y)=C(F(x),G(y)). $$ If both margins are continuous, this copula is unique. Conversely, if \(C\) is a copula and \(F,G\) are univariate distribution functions, then \(C(F(x),G(y))\) is a joint distribution function with margins \(F\) and \(G\).

Here is a proof of the representation and uniqueness when the margins are continuous. This is the case used in the examples below. Define \(U=F(X)\) and \(V=G(Y)\). The probability integral transform says that \(U\) and \(V\) are each uniform on \([0,1]\). Their joint distribution function $$ C(u,v)=P(U\leq u,\ V\leq v) $$ is therefore a copula.

For continuous \(F\), the events \(\{X\leq x\}\) and \(\{F(X)\leq F(x)\}\) agree up to probability zero. Monotonicity gives one inclusion. Any discrepancy in the other direction can only occur where \(F\) is flat, and a flat interval has probability zero under the distribution \(F\). The same reasoning applies to \(G\). Consequently, $$ H(x,y)=P(F(X)\leq F(x),\ G(Y)\leq G(y))=C(F(x),G(y)). $$ This proves existence for continuous margins.

For uniqueness, let \(\widetilde C\) be another copula representing \(H\). For any \(u,v\) strictly between zero and one, continuity of \(F\) and \(G\) ensures that there are \(x,y\) with \(F(x)=u\) and \(G(y)=v\). Thus $$ \widetilde C(u,v)=H(x,y)=C(u,v). $$ Copulas are continuous: changing \(u\) by an amount \(\delta\) can change the probability by at most \(|\delta|\), because the first margin is uniform; the corresponding bound for changing \(v\) is \(|\delta|\). The equality therefore extends to the boundary of the unit square. This proves uniqueness. Conversely, composing a copula with \(F\) and \(G\) gives a joint distribution function: the rectangle increments are nonnegative because \(F\) and \(G\) are nondecreasing and \(C\) has nonnegative rectangle probabilities. Setting either argument to its marginal endpoint gives \(F\) and \(G\) as the resulting margins. \(\square\)

Worked Examples: Same Margins, Different Dependence

Worked Example: Independent Exponential Measurements

Let \(U\) and \(V\) be independent uniform random variables on \([0,1]\), and define \(X=-\log(1-U)\) and \(Y=-\log(1-V)\). For \(x\geq0\), $$ P(X\leq x)=P(U\leq 1-e^{-x})=1-e^{-x}, $$ and the same calculation applies to \(Y\). Both variables therefore have the unit-rate exponential distribution.

Because \(U\) and \(V\) are independent, the copula is \(C(u,v)=uv\). At the threshold \(x=y=\log 2\), each marginal probability is \(1-e^{-\log 2}=1/2\), so $$ P(X\leq \log 2,\ Y\leq \log 2)=C(1/2,1/2)=1/4. $$ The joint density is \(e^{-x-y}\) for \(x,y\geq0\), confirming independence.

Worked Example: Perfectly Opposing Ranks With the Same Margins

Let \(U\) be uniform on \([0,1]\), and set \(X=-\log(1-U)\), \(Y=-\log U\). Each variable is again unit-rate exponential: for \(y\geq0\), for example, \(P(Y\leq y)=P(U\geq e^{-y})=1-e^{-y}\).

The uniform variables underlying the margins are \(U\) and \(1-U\). Their copula is $$ C(u,v)=P(U\leq u,\ 1-U\leq v)=\max(u+v-1,0). $$ At \(x=y=\log 2\), the two events require \(U\leq1/2\) and \(U\geq1/2\), respectively. Their intersection has probability zero. The joint probability is therefore \(0\), rather than the \(1/4\) found for independent exponential measurements. The margins are unchanged; the copula is different.

Worked Example: A Mixture of Independence and Comonotonicity

Let \(U\) and \(W\) be independent uniform random variables, and let \(Z\) be independent of both with \(P(Z=1)=\lambda\), where \(0\leq\lambda\leq1\). Set \(V=U\) when \(Z=1\), and \(V=W\) when \(Z=0\). Both \(U\) and \(V\) are uniform. Conditioning on \(Z\) gives the copula $$ C_\lambda(u,v)=\lambda\min(u,v)+(1-\lambda)uv. $$

Use the exponential transformations \(X=-\log(1-U)\) and \(Y=-\log(1-V)\). Both margins are unit-rate exponential for every \(\lambda\). With \(\lambda=0.6\), the joint probability that both measurements are at most \(\log 2\) is $$ C_{0.6}(1/2,1/2)=0.6(1/2)+0.4(1/4)=0.4. $$ This lies between the independent value \(0.25\) and the comonotone value \(0.5\).

The covariance also changes with \(\lambda\). Conditional on \(Z=1\), \(X=Y\), so \(E[XY\mid Z=1]=E[X^2]=2\); conditional on \(Z=0\), independence gives \(E[XY\mid Z=0]=1\). Hence \(E[XY]=0.6(2)+0.4(1)=1.6\). Since \(E[X]=E[Y]=1\), \(\operatorname{Cov}(X,Y)=1.6-1=0.6\). Each unit-rate exponential margin has variance \(1\), so the correlation is also \(0.6\).

When a Copula Has a Density

If the joint distribution and margins are smooth, the copula also gives a convenient density formula. Suppose \(F\) and \(G\) are continuously differentiable with densities \(f\) and \(g\), and suppose the copula has a continuous density \(c(u,v)\), so that its mixed derivative is \(c\). Then the joint density is $$ h(x,y)=c(F(x),G(y))f(x)g(y). $$ To verify this, write \(H(x,y)=C(F(x),G(y))\). Differentiating first with respect to \(x\) and then with respect to \(y\), the chain rule gives $$ \frac{\partial^2 H(x,y)}{\partial y\,\partial x} =\frac{\partial^2 C(F(x),G(y))}{\partial v\,\partial u}f(x)g(y) =c(F(x),G(y))f(x)g(y). $$ Under the stated smoothness conditions, this mixed derivative is the joint density. For independence, \(c(u,v)=1\), and the formula reduces to \(h(x,y)=f(x)g(y)\). Not every copula has a density: the comonotone copula, for example, concentrates probability on the diagonal and is singular with respect to area.

Interpretation and a Common Pitfall

A copula describes how marginal ranks occur together. For continuous margins, \(F(X)\) and \(G(Y)\) are uniform, so a point such as \(C(0.9,0.9)\) is the probability that both variables fall at or below their respective 90th percentiles. This interpretation does not require the original measurements to share units or even to have similar marginal shapes.

A common mistake is to infer a joint distribution from the two marginal distributions alone. For example, knowing that two measurements are each unit-rate exponential does not justify multiplying their marginal probabilities. That multiplication gives the joint probability only if the variables are independent, which is the additional assumption \(C(u,v)=uv\). The independent, countermonotone, and mixture examples all have the same exponential margins but different probabilities of joint events.

Another useful implication is invariance under strictly increasing transformations: increasing a measurement’s scale or changing its units does not alter its rank, and hence does not alter its copula. A decreasing transformation reverses that variable’s ranks instead. This is why separating margins from dependence can help when comparing dependence across different measurement scales.

Check Your Understanding

Use the definitions and examples above to check your understanding.

  1. What two properties distinguish a bivariate copula from an arbitrary distribution function on the unit square?
  2. For continuous margins, how does Sklar’s representation express the joint distribution in terms of the margins and a copula?
  3. For two unit-rate exponential variables with the mixture copula \(C_\lambda\), what is the probability that both are at most \(\log 2\), in terms of \(\lambda\)?
  4. Why does multiplying two marginal probabilities require an independence assumption?
  5. Can every copula be represented by a density with respect to area? Explain using the comonotone example.
Clinical Biostats