R语言Copula函数股市相关性建模:模拟Random Walk(随机游走)

在引入copula时,大家普遍认为copula很有趣,因为它们允许分别对边缘分布和相依结构进行建模。

 

由Kaizong Ye,Liao Bao撰写

给定一些边缘分布函数和一个copula,那么我们可以生成一个多元分布函数,其中的边缘是前面指定的。

考虑一个二元对数正态分布

我们可以从边缘分布开始。

×

当边缘分布(marginal probability distribution)不同的随机变量(random variable),互相之间并不独立的时候,此时对于联合分布的建模会变得十分困难。此时,在已知多个已知边缘分布的随机变量下,Copula函数则是一个非常好的工具来对其相关性进行建模。

什么是Copula函数?

copula这个单词来自于拉丁语,意思是“连接”。最早是由Sklar在1959年提出的,即Sklar定理:

以二元为例,若 [公式] 是一个具有连续边缘分布的 [公式] 与 [公式] 的二元联合分布函数,那么存在唯一的Copula函数 [公式] ,使得 [公式] 。反之,如果 [公式] 是一个copula函数,而 [公式] 和 [公式] 是两个任意的概率分布函数,那么由上式定义的 [公式] 函数一定是一个联合分布函数,且对应的边缘分布刚好就是 [公式] 和 [公式] 。

Sklars theorem : Any multivariate joint distribution can be written in terms of univariate marginal distribution functions and a copula which describes the dependence structure between the two variable.

Sklar认为,对于N个随机变量的联合分布,可以将其分解为这N个变量各自的边缘分布和一个Copula函数,从而将变量的随机性和耦合性分离开来。其中,随机变量各自的随机性由边缘分布进行描述,随机变量之间的耦合特性由Copula函数进行描述。换句话说,一个联合分布关于相关性的性质,完全由其Copula函数决定。

如果已知 [公式] , [公式] 和 [公式] ,则Copula函数可以表达为:

[公式]

这里 [公式] 代表 [公式] 的反函数,或者叫CDF的逆变换、逆累积分布函数。

Copula理论的数学表达

假设 [公式] 是 [公式] 个随机变量,它们各自的边缘分布分别为 [公式] ,它们的联合分布为 [公式] ,则存在一个将边缘分布和联合分布“连接”起来的函数 [公式] ,使得:

[公式]

而根据边缘分布的CDF的逆变换,即 [公式] ,则可以得到Copula函数的表达形式:

[公式]

补充知识:概率积分变换(Probability Integral transform)

在概率论中,概率积分变换(也称为均匀的普适性)是指将任意给定连续分布的随机变量的数据值转换成具有标准均匀分布的随机变量的结果。

简单解释一下概率积分变换[1]:如果 [公式] 和 [公式] 都是随机变量(Random Variables,RV),而 [公式] 、 [公式] 分别是二者的累计概率分布函数,即 [公式]

那么:[公式] 和 [公式] 都服从均匀分布, [公式]

简要的数学证明如下:

[公式]

换句话说,任何边际分布的CDF值都均匀分布在区间[0,1]上。如果从任意分布中随机抽取,那么抽取该分布的最大值(U=1)的概率与抽取可能的最小值(U=0)或中值(U= 5)的概率相同。而copula实际上是它所建模的随机变量的CDFs的联合分布。

基于这些边缘分布,并考虑从该伪随机样本获得的copula参数的最大似然估计值,从数值上讲,我们得到

但是,由于相依关系是边缘分布的函数,因此我们没有对相依关系进行单独处理。如果考虑全局优化问题,则结果会有所不同。可以得出密度

差别不大,但估计量并不相同。从统计的角度来看,我们几乎无法分别处理边缘和相依结构。我们应该记住的另一点是,边际分布可能会错误指定。例如,如果我们假设指数分布,

高斯copula的参数估计


热门课程

R语言数据分析挖掘必知必会

面对扑面而来的数据浪潮,包含Google、Facebook等国际企业,都已采用R语言进行数据分析

探索课程

由于我们错误地指定了边缘分布,因此我们无法获得统一的边缘。如果我们使用上述代码生成大小为500的样本,

如果边缘分布被很好地设定时,我们可以清楚地看到相依结构依赖于边缘分布,

copula模拟股市中相关随机游走

接下来我们用copula函数模拟股市中的相关随机游走


R语言实现 Copula 算法建模相依性案例分析报告

阅读文章



随时关注您喜欢的主题


FB

HO

FB

HO

标准偏差相对于均值而言非常大,接近于零;因此,在某些情况下,我们很有可能获得不稳定的结果。


可下载资源

关于作者

Kaizong Ye拓端研究室(TRL)的研究员。

本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。

​非常感谢您阅读本文,如需帮助请联系我们!

 
QQ在线咨询
售前咨询热线
15121130882
售后咨询热线
0571-63341498