コピュラ (統計学)のソースを表示

'''コピュラ''' ({{lang-en-short|copula}}) とは、[[統計学]]において多変数の[[累積分布関数]]とその[[周辺分布関数]]の関係を示す関数のことである。

[[確率変数]]の相関を表す指標として代表的なものに[[相関係数]]があるが、相関係数が 1 個の数値であるのに対してコピュラは関数であることから、確率変数の間のきわめて多様な依存関係を表すことができる。なお、名称はラテン語で相異なる物同士の「つなぎ」や「結び付き」を意味する名詞 copula（{{lang-en-short|couple}} の語源）に由来する。この単語は元々[[音楽]]や[[言語学]]で使われていたが、統計学の用語として用いたのは、1959 年にスクラー (Abe Sklar) がパリ大学統計学会誌 (the Statistical Institute of the University of Paris) で発表したのが最初である<ref>Roger B. Nelsen (1999), ''An Introduction to Copulas''. ISBN 0-387-98623-5.</ref>。

== 定義 ==
''n'' 次元単位立方体 [0, 1]<sup>''n''</sup> から[[単位区間]] [0, 1] への関数 ''C'': [0, 1]<sup>''n''</sup> → [0, 1] が次の性質をもつとき、''C'' を ''n'' 次元コピュラ（または ''n'' コピュラ）という。

* <math>\mathbf u \in [0,1]^n</math> のうち少なくとも 1 つの要素が 0 であるとき、すなわち '''u''' = (''u''<sub>1</sub>, ..., ''u''<sub>''i''-1</sub>, 0, ''u''<sub>''i''+1</sub>, ..., ''u<sub>n</sub>''); ''i'' = 1, 2, ..., ''n'' であるとき ''C''('''u''') = 0

* <math>\mathbf u \in [0,1]^n</math> が 1 つの要素を除いてすべて 1 であるとき、すなわち '''u''' = (1, ..., 1, ''u<sub>i</sub>'', 1, ..., 1); ''i'' = 1, 2, ..., ''n'' であるとき ''C''('''u''') = ''u<sub>i</sub>''

* ''C''('''u''') は ''n''-increasing である、すなわち ''n'' 次元単位立方体内の任意の直方体 <math>B = \times_{i = 1}^{n}[x_i, y_i] \subseteq [0,1]^n</math> について

: <math>V_{C} \left( B \right) = \Delta_{x_n}^{y_n} \Delta_{x_{n - 1}}^{y_{n - 1}} \cdots \Delta_{x_1}^{y_1} C(\mathbf t) \geq 0</math>

ここで <math>\Delta_{x_i}^{y_i} C(\mathbf t) = C(t_1, \ldots, t_{i - 1}, y_i, t_{i + 1}, \ldots, t_n) - C(t_1, \ldots, t_{i - 1}, x_i, t_{i + 1}, \ldots, t_n)</math> である。

== スクラーの定理 ==
'''スクラーの定理'''は 1959 年にスクラーが示したもので、コピュラに関する基本的な定理である。定理は次のとおり。

: ''n'' 次元分布関数 ''H'' が1次元周辺分布関数 ''F''<sub>1</sub>, ''F''<sub>2</sub>, ..., ''F<sub>n</sub>'' をもつとき、''n'' 次元コピュラ ''C'' が存在して以下が成り立つ。

:: ''H''(''x''<sub>1</sub>, ''x''<sub>2</sub>, ..., ''x<sub>n</sub>'') = ''C''(''F''<sub>1</sub>(''x''<sub>1</sub>), ''F''<sub>2</sub>(''x''<sub>2</sub>), ..., ''F<sub>n</sub>''(''x<sub>n</sub>''))

: 周辺分布関数 ''F''<sub>1</sub>, ''F''<sub>2</sub>, ..., ''F<sub>n</sub>'' が連続であるとき、コピュラ ''C'' は一意に定まる。

== フレシェ-ヘフディング境界 ==
=== フレシェ-ヘフディング上界 ===
次の式で与えられる ''M'' は '''フレシェ-ヘフディング上界''' (Fréchet-Hoeffding upper bound) と呼ばれる。

:<math>M(u_1, \ldots, u_n) = \min_{1 \leq j \leq n} u_j</math>

任意のコピュラ ''C'' および任意の <math>(u_1, \ldots, u_n) \in [0, 1]^n</math> に対して <math>M(u_1, \ldots, u_n) \geq C(u_1, \ldots, u_n)</math> であることから、''M'' はコピュラの中で最大のものである。

=== フレシェ-ヘフディング下界 ===
次の式で与えられる ''W'' は '''フレシェ-ヘフディング下界''' (Fréchet-Hoeffding lower bound) と呼ばれる。

:<math>W(u_1, \ldots, u_n) = \max \left\{ \sum_{i = 1}^n {u_i} + 1 - n, 0 \right\}</math>

任意のコピュラ ''C'' および任意の <math> (u_1, \ldots, u_n) \in [0, 1]^n </math> に対して <math> W(u_1, \ldots, u_n) \leq C(u_1, \ldots, u_n) </math> が成り立つ。ただし、''W'' は 2 次元以外の場合にはコピュラではない。

== 代表的なコピュラ ==
以下に代表的なコピュラを示す。

=== 積コピュラ ===
Π(''u'', ''v'') = ''uv'' は'''積コピュラ'''と呼ばれる。

確率変数 ''X'' と ''Y'' がそれぞれ確率分布関数 ''F'' および ''G'' に従い、また ''X'' と ''Y'' の結合分布関数を ''H'' とする。このときスクラーの定理によって ''H''(''x'', ''y'') = ''C''(''F''(''x''), ''G''(''y'')) をみたすコピュラ ''C'' が存在することとなるが、''X'' と ''Y'' が互いに独立であることと、''C'' = Π であることとは同値となる。この意味で積コピュラを'''独立コピュラ'''と呼ぶこともある。

=== アルキメデスコピュラ ===
ある関数 φ を用いて次のように表せるコピュラを'''アルキメデスコピュラ''' (Archimediean copula) という。

: ''C''(''u'', ''v'') = φ<sup>-1</sup>(φ(''u'') + φ(''v''))

関数 φ は [0, 1] から <math>[0, \infty]</math> への連続な単調減少関数であって φ(1) = 0 を満たすものである。この φ を '''ジェネレーター''' (generator) という。

上記の ''W'' や Π もアルキメデスコピュラである（ジェネレーターはそれぞれ 1 - ''t'' および - ln ''t''）。なお、''M'' はアルキメデスコピュラではない。

==== クレイトンコピュラ ====
: <math>C(u, v) = \left[ \max \left( u^{-\theta} + v^{-\theta} - 1, 0 \right) \right]^{-1/\theta}</math> &nbsp;&nbsp;&nbsp;&nbsp; (<math>-1 \leq \theta < 0</math> または <math>0 < \theta</math>)

で表されるアルキメデスコピュラは'''クレイトン''' (Clayton) '''コピュラ'''と呼ばれる。このコピュラのジェネレーターは <math>\varphi(t) = \frac{1}{\theta} \left( t^{- \theta} - 1 \right)</math> である。

θ = 0 のときは積コピュラに、θ = -1 のときはフレシェ-ヘフディング下界になる。

==== グンベルコピュラ ====
: <math>C(u, v) = \exp \left( - \left[ (-\ln u)^\theta + (-\ln v)^\theta \right]^{1/\theta} \right)</math> &nbsp;&nbsp;&nbsp;&nbsp; (<math>1 \leq \theta</math>)

で表されるアルキメデスコピュラは'''グンベル''' (Gumbel) '''コピュラ'''または'''ガンベルコピュラ'''と呼ばれる。このコピュラのジェネレーターは φ(''t'') = (- ln ''t'')<sup>θ</sup> である。

θ = 1 のときは積コピュラとなる。

==== フランクコピュラ ====
:<math>C(u, v) = - \dfrac{1}{\theta} \ln \left(1 + \frac{(e^{-\theta u} - 1) (e^{-\theta v} - 1)}{e^{-\theta} - 1}\right)</math> &nbsp;&nbsp;&nbsp;&nbsp; (<math>\theta \neq 0</math>)

で表されるアルキメデスコピュラは'''フランク''' (Frank) '''コピュラ'''と呼ばれる。このコピュラのジェネレーターは <math>\varphi(t)= \ln \left( \frac{e^{\theta t} -1}{e^{\theta} -1}\right)</math> である。

θ = 0 のときは積コピュラとなる。

=== 正規コピュラ ===
'''ガウス・コピュラ'''、'''ガウス型コピュラ''' ({{lang-en-short|Gaussian copula}}) ともいう。2000年に[[:en:David X. Li]]が発表し、以後[[金融工学]]にて[[Collateralized Debt Obligation|債務担保証券]] (CDO) のリスク評価等に広く使われた。

確率変数 ''X'', ''Y'' に対して相関行列 Σ を持つ 2 変数正規分布関数を Φ<sub>2</sub>(''x'', ''y''; Σ) で表し、1 変数標準正規分布関数を Φ<sub>1</sub>(''x'') で表すものとする。このときスクラーの定理によって

: Φ<sub>2</sub>(''x'', ''y''; Σ) = ''C''(Φ<sub>1</sub>(''x''), Φ<sub>1</sub>(''y''))

を満たすコピュラ ''C'' が存在する。このコピュラを'''正規コピュラ'''という。

同様に [[t分布|''t'' 分布]]から作られるコピュラを ''t'' '''コピュラ'''という。

== コピュラの応用 ==
コピュラの実務面への応用例としては、[[Collateralized Debt Obligation|CDO]] の価格評価やリスク評価が挙げられる。CDO は複数の債務をまとめて[[証券化]]したものであるので、複数の債務がどのような確率でデフォルトを起こすかが問題となる。平常時においてはデフォルト確率の相関が低い債務であっても、景気悪化時には連鎖倒産などで相関が高まるといったことが考えられるため、1 個の相関係数では十分に価格やリスクを表せないことから、コピュラが用いられる。

コピュラはこのように、発生率の低い（すなわち分布関数の値が 0 または 1 に近い値を取る）部分で相関が高まるような場合（これを'''テール依存性'''という）での応用がしばしば考えられる。

===批判===
[[金融工学]]の世界では『悪魔の関数』また『粉飾の関数』と揶揄される。[[債務担保証券|CDO]]（Collateralized Debt Obligation）の格付けやリスク評価に関して特に正規コピュラ（ガウス・コピュラ）が広く使われた。正規コピュラを用いる実務上の利点は、2者間の相関を調べる際に、過去に生じたリスク事象の統計分析等によらずに [[クレジット・デフォルト・スワップ|CDS]] の値動きだけに着目すれば良いことにあった。これによってリスク評価実務は大幅に簡略化され、[[証券化]]市場は劇的な拡大を見た。ところが、正規コピュラには予測不能性が織り込まれておらず、現実世界と乖離しうることから、これを無批判に応用したことが 21世紀初頭の金融恐慌（[[リーマン・ショック]]）とその後の証券化市場衰退を招いた原因の一つだとする説がある<ref>{{Citation | last = Salmon | first = Felix | year = 2009 | authorlink=:en:Felix Salmon | title = Recipe for Disaster: The Formula That Killed Wall Street | journal = Wired | date = 2009-02-23 | url = http://www.wired.com/techbiz/it/magazine/17-03/wp_quant|accessdate=2010-07-07}}</ref>。

== 脚注 ==
{{reflist}}

{{統計学}}
{{Normdaten}}
{{DEFAULTSORT:こひゆら}}
[[Category:統計学]]
[[Category:数学に関する記事]]
[[Category:金融工学]]