查看“︁耦合 (概率)”︁的源代码

{{otheruses|耦合}}
{{NoteTA|G1=Math}}
'''耦合'''，或稱'''关联结构'''（{{lang-en|Copula}}），為处理统计中[[随机变量]][[相关性]]问题的一种方法，由一组随机变量的[[邊際分布]]来确定它们的[[联合分布]]。通过关联结构来确定一个联合分布的方法是基于如下思想，一个简单转换可以通过分别将每个边缘分布都转换为[[平均分布]]的转换组成。这样，一个关联结构（dependence structure）就可以表达为一个基于上述所得平均分布之上的联合分布，而关联结构（copula）即是边缘均匀随机变量之上的一个联合分布。在实际应用中，上述的转换可能被设置为每个边缘变量的初始化步骤，或者上述转换的参数可能根据具体'''关联结构'''的对应参数设置。

按照所表达的关联关系的不同，关联结构被分为很多不同类别。典型情况下，一个种类的关联结构有多个参数用来表达不同的关联强度和关联类型。下面将大概描述一些有代表性的关联结构。关联结构的一个典型应用是，通过选择某一种类的关联结构来定义某一适合特定样本数据分布的联合分布，当然关联结构也可以来自于任何相应的给定联合分布。

== 基本思想 ==
考察两个随机变量<math>X</math>, <math>Y</math>，分别具有连续[[累积分布函数]]<math>F_X</math>和<math>F_Y</math>。通过分别在两个随机变量上应用[[概率积分转换]]，得到<math>X'=F_X(X)</math> 和<math>Y'=F_Y(Y)</math>。因此<math>X'</math>和<math>Y'</math>都是具有连续[[均匀分布]]的变量，相关性通常取决于''<math>X</math>''和''<math>Y</math>''是否是相关（自然，如果''<math>X</math>''和''<math>Y</math>''是不相关的，那么''<math>X'</math>''和''<math>Y'</math>''也是不相关的）。因为这个转换是可逆的，可以定义''<math>X</math>''和''<math>Y</math>''之间的相关性等于''<math>X'</math>''和''<math>Y'</math>''之间的相关性。因为''<math>X'</math>''和''<math>Y'</math>''是均匀分布的随机变量，所以问题被简化为定义一个在两个均匀分布之上的[[二项分布]]，这就是关联结构。所以，这一基本思想就是，通过把边缘变量转化为均匀分布变量而不再需要考察很多不同的边缘分布以简化问题，然后再把相关性定义为一个在均匀分布之上的联合分布。

== 定义 ==

一个 '''关联结构'''是一个定义在<math>n</math>维[[单位立方体]]<math>[0,1]^n</math>上的多元[[联合分布]]，其每个[[边缘分布]]都是在<math>[0,1]</math>区间上的[[均匀分布]]。

特别的，<math>C:[0,1]^n\to [0,1]</math>是一个''n''维关联结构，有

: <math> C\left(\mathbf u\right)=0</math>当<math>\mathbf u\in [0,1]^n</math>有至少一个分量为<math>0;</math>

: <math> C\left(\mathbf u\right)=u_i</math>当<math>\mathbf u\in [0,1]^n</math>所有分量为<math>1</math>除了第''i''个分量等于<math>u_i;</math>

: <math> C\left(\mathbf u\right) </math>是[[n维递增]]的，也即，有每个[[hyperrectangle]] <math>B=\times_{i=1}^{n}[x_i,y_i]\subseteq [0,1]^n;</math> 

: <math> V_{C}\left( B\right):=\sum_{\mathbf z\in \times_{i=1}^{n}\{x_i,y_i\}} (-1)^{N(\mathbf z)} C(\mathbf z)\ge 0;</math>

其中<math>N(\mathbf z)=\operatorname{card}\{k\mid z_k=x_k\}</math>。<math> V_{C}\left( B\right)</math>所谓的<math>B</math>的''C''-体积（volume）。

== Sklar定理 ==

由Sklar提出的这条定理<ref name="Sklar 1959">{{cite journal | last=Sklar | first=A. | title=Fonctions de répartition à n dimensions et leurs marges | journal=Publ. Inst. Statist. Univ. Paris | year=1959 | volume=8 | pages=229–231}}</ref>是大多数关联结构的应用的基础。Sklar定理指出，一个给定的<math>p</math>个变量的联合分布函数<math>H</math>，<math>F_1, F_2,F_3,\ldots,F_p</math>为其边缘分布函数，必存在这样一个关联结构<math>C</math>使<math>H = C(F_1, F_2,F_3,\ldots,F_p)</math>

以[[二项分布]]为例，Sklar定理应用如下。对任一二项分布函数<math>H(x,y)</math>，令<math>F(x)=H(x,\infty)</math> 而<math>G(y)=H(\infty,y)</math>为其单变量边缘概率分布函数。那么存在关联结构''<math>C</math>''以使

:<math>H(x,y)=C(F(x),G(y))\,</math>

（此处已知分布''<math>C</math>''和它的累积分布函数）。此外，如果边缘分布<math>F(x)</math> 和<math>G(y)</math>连续，那么关联结构函数''<math>C</math>''是唯一的。否则，关联结构''<math>C</math>''在边缘分布的值域上是唯一确定的。

== 弗雷歇–霍夫丁（Fréchet–Hoeffding）关联结构边界 ==
[[File:Copule ord.svg|thumb|right|Graphs of the Fréchet–Hoeffding copula limits and of the independence copula (in the middle).]]

'''最小（反单调）关联结构：'''是所有关联结构的下边界。仅在二项分布中，变量间表现为完全负相关。

:<math> W(u,v) = \max(0,u+v-1).\,</math>

对''n''-元关联结构，下边界为
:<math> W(u_1,\ldots,u_n) := \max\left\{1-n+\sum\limits_{i=1}^n {u_i} , 0 \right\} \leq C(u_1,\ldots,u_n). </math>

'''最大 （共单调 ） 关联结构：'''是所有关联结构的上边界。其在二项分布中，变量间表现为完全正相关：

:<math> M(u,v) = \min(u,v).\,</math>

对''n''-元关联结构，上边界为
:<math>C(u_1,\ldots,u_n)\le \min_{j \in \{1,\ldots,n\}} u_j =: M(u_1,\ldots,u_n).</math>

'''结论：'''对所有关联结构''C''（''u'', ''v''）,
:<math> W(u,v) \le C(u,v) \le M(u,v).</math>

对于多元关联的情况为
:<math> W(u_1,\ldots,u_n) \le C(u_1,\ldots,u_n) \le M(u_1,\ldots,u_n).</math>

== 关联结构种类 ==
=== 正态关联结构 ===
[[File:Copula_gaussian.svg|thumb|right|Cumulative distribution and probability density functions of Gaussian copula with <math>\rho=0.4</math>]]

在金融建模中常用到的一个关联结构是正态关联结构，正态关联结构是根据Sklar定理由二元[[正态分布]]构成。设<math>\Phi_\rho</math>是标准二元正态累积分布函数，[[相关系数]]为''ρ''，则正态关联结构函数为

:<math> C_\rho(u,v) = \Phi_\rho \left(\Phi^{-1}(u), \Phi^{-1}(v) \right)</math>

其中，<math>u, v \in [0,1]</math>而<math>\Phi</math>表示标准正态[[累积分布函数]]。

对''C''[[微分]]得出关联结构的密度函数：

:<math> c_\rho(u,v) = \frac{\varphi_{X,Y, \rho} (\Phi^{-1}(u), \Phi^{-1}(v))}
{\varphi(\Phi^{-1}(u))\varphi(\Phi^{-1}(v))}</math>

其中

:<math> \varphi_{X,Y, \rho}(x,y) = \frac{1}{2 \pi\sqrt{1-\rho^2}} \exp \left ( -\frac{1}{2(1-\rho^2)}  \left [{x^2+y^2} -2\rho xy  \right ] \right ) </math>

是[[相关#皮尔逊积差系数|皮尔逊矩相关系数]]为<math>\rho</math>标准二元正态分布的概率密度函数，其标准正态密度为<math>\varphi</math>。

== 参考资料 ==
{{reflist}}

{{Authority control}}
[[Category:獨立 (機率論)]]
[[Category:多變量統計]]
[[Category:保险学]]