小数定律的保险业应用:泊松分布模拟索赔次数

在保险业中,由于分散投资,通常会在合法的大型投资组合中提及大数定律。在一定时期内,损失“可预测”。

当然,在标准的统计假设下,即有限的期望值和独立性。由于在保险业中,灾难通常很少发生,而且代价非常高昂,精算师可能有兴趣对少量事件的发生进行建模。背后的定理有时也被称为小数定律

泊松分布

所谓的泊松分布(请参阅http://en.wikipedia.org/…)由SiméonPoisson于1837年进行了介绍。亚伯拉罕·德·莫伊夫(Abraham De Moivre)于1711年在De Mensura Sortis seu对其进行了定义

让  https://latex.codecogs.com/gif.latex?N​表示一个计数随机变量,然后它是服从泊松分布,如果有  https://latex.codecogs.com/gif.latex?\lambda\in(0,\infty)​这样

https://latex.codecogs.com/gif.latex?\mathbb {P}(N = k)= e ^ {-\ lambda} \ frac {\ lambda ^ k} {k!},\ forall%20k \ in \ mathbb {N}

De Moivre从二项式分布的近似值获得了该分布。回想一下,二项式分布是精算科学中的标准分布,例如,用来模拟https://latex.codecogs.com/gif.latex?n​被保险人死亡人数  。如果单个死亡概率相同,例如https://latex.codecogs.com/gif.latex?p​,并且如果死亡是独立事件,则

https://latex.codecogs.com/gif.latex?\mathbb{P}(N=k)=\binom{n}{k}p^k(1-p)^{nk},\forall%20k \ in \ {0,1,\ cdots,n \}


而如果  https://latex.codecogs.com/gif.latex?n\rightarrow\infty​和  https://latex.codecogs.com/gif.latex?np\rightarrow%20\lambda​,然后

https://latex.codecogs.com/gif.latex?\mathbb{P}(N=k)\rightarrow%20e^{-\lambda}\frac{\lambda^k}{k!}

​再次,这是一个渐近定理,当我们有很多观察值时(https://latex.codecogs.com/gif.latex?n\rightarrow\infty​)成立,它也成立,而且出现的可能性应该非常小(因为https://latex.codecogs.com/gif.latex?p\sim\lambda/n​),这就是为什么要使用术语“ 小数”的原因。SiméonPoisson对数学近似值不感兴趣:他的主要观点是针对他正在处理的数据获得具有良好拟合优度的分布。


This image has an empty alt attribute; its file name is image.png

专栏

精算科学

关于结合数学、统计方法以及程序语言对经济活动来做风险分析、评估的见解。

探索专栏

小数定律

与Poisson分布有关的主要定理的启发式如下:  https://latex.codecogs.com/gif.latex?X_1,%20\cdots,X_n​表示iid随机变量采用值https://latex.codecogs.com/gif.latex?%20\mathbb{R}^d​(一般情况下,一个分量可以是时间,另一分量可以是感兴趣的上部区域,其中某些随机过程是可能)。让https://latex.codecogs.com/gif.latex?\mathcal{A}_n\subset\mathbb{R}^d​。如果  https://latex.codecogs.com/gif.latex?\mathbb{P}(X_i%20\in%20\mathcal{A}_n)\rightarrow%200​作为假设https://latex.codecogs.com/gif.latex?n\rightarrow\infty​(或https://latex.codecogs.com/gif.latex?\mathbb{P}(X_i%20\in%20\mathcal{A}_n)=O(n^{-1})​更具体地假设),则  https://latex.codecogs.com/gif.latex?N​表示事件的(随机变量表征)计数https://latex.codecogs.com/gif.latex?\ {X_i%20 \ in%20 \ mathcal {A} _n \}​,则  https://latex.codecogs.com/gif.latex?N​可以通过带有参数的泊松分布来近似https://latex.codecogs.com/gif.latex?\lambda%20=n%20\times%20\mathbb%20P(X_i%20\in%20\mathcal{A}_n)​。
启发式方法是,如果考虑大量观察值,并且计算给定(小)区域中有多少观察值,则此类观察值的数量就是泊松分布。

n=1000






polygon(c(u,rev(u)),c(v,rev(-v)),col="yellow",border=NA)
I=(X^2+Y^2)<1
points(X[I],Y[I],cex=.6,pch=19,col="red")

为了获得与保险模型有关的解释,让我们  https://latex.codecogs.com/gif.latex?\mathcal{A}​在再保险合同中表示上层,即https://latex.codecogs.com/gif.latex?\mathcal {A} = \ {x%3Ed \}​某些可扣除额https://latex.codecogs.com/gif.latex?d​。让我们  https://latex.codecogs.com/gif.latex?X_i​来表示个人损失。然后,可以使用泊松分布对到达该上层的索赔的数量进行建模。更准确地说,如果自付额  https://latex.codecogs.com/gif.latex?d​变得非常大(和https://latex.codecogs.com/gif.latex?\mathbb{P}(X_i%20\in%20\mathcal{A})\rightarrow%200​),我们将获得极值理论中的阈值点以上模型:如果  https://latex.codecogs.com/gif.latex?N​有一个泊松分布,并在有条件的https://latex.codecogs.com/gif.latex?N​,  https://latex.codecogs.com/gif.latex?X_1,\cdots,X_N​是独立同分布的广义帕累托随机变量,然后 https://latex.codecogs.com/gif.latex?\ max \ {X_1,\ cdots,X_N \}​具有广义的极值分布。因此,超出模型(针对罕见事件)与泊松过程密切相关。

泊松过程

如上所述,当事件以某种方式随机且独立地随时间发生时,就会出现泊松分布。然后很自然地研究两次事件之间的时间(或在保险范围内两次索赔)。

泊松分布和索赔发生

既不是SiméonPoisson也不是De Moivre,而是Ladislaus Von Bortkiewicz首先提到了Poisson分布是小数定律。1898年,他研究了1875年至1894年间被马踢倒杀死的士兵的人数,其中有200个兵团。

他确实获得了以下分布(此处,泊松分布的参数为0.61,即每年的平均死亡人数)


在很多情况下,泊松分布都非常适合。例如,如果我们考虑1850年后在佛罗里达州的飓风数量,

泊松分布和回归期

返回期是由Emil Gumbel在水文学中介绍的,用于链接概率和持续时间。十年事件的发生概率为1/10。那么10是发生之前的平均等待时间。这并不意味着该事件不会在10年之前发生,或者必须在10年之前发生。考虑一个返回期  https://latex.codecogs.com/gif.latex?T​(以年为单位),则每年不出现的概率为https://latex.codecogs.com/gif.latex?1-(1/T)​。

https://latex.codecogs.com/gif.latex?n​多年未发生的概率为  https://latex.codecogs.com/gif.latex?1-[1-(1/T)]^n​。通常用下表来总结此属性,

上表中的对角线非常有趣。似乎在某种程度上趋向极限值(此处为63.2%)。在n年内观察到的事件数量具有二项式分布,其概率为https://latex.codecogs.com/gif.latex?1/T=1/n​,将收敛到参数为1的泊松分布。那么https://latex.codecogs.com/gif.latex?1-\exp(-1)​,没有灾难的概率为,等于0.632。

稀有概率与泊松分布

计算稀有事件的概率时,泊松分布不断出现。例如,在50年的时间里,至少有一次在核电厂发生事故的可能性。假设在反应堆中发生事故的年概率  https://latex.codecogs.com/gif.latex?p​很小,例如0.05%。进一步假设反应堆在时间上相互独立。在50年内发生超过80个反应堆的事件的概率是

https://latex.codecogs.com/gif.latex?\mathbb{P}(N\neq0)=1-(1-p)^{50%20\times%2080}

当然,线性近似是不正确的

https://latex.codecogs.com/gif.latex?\mathbb%20P(N\neq%200)\neq%2050\times%2080\times%20p

另一方面

https://latex.codecogs.com/gif.latex?\mathbb%20P(N\neq 0)= 1-(1-p)^ {50 \ times80%20}%20 \ sim1- \ exp \ left( -50 \ times80 \ times%20p%20 \ right)

> 
> 
[1] 0.1812733
> 
[1] 0.1812692

这是具有参数为的泊松分布https://latex.codecogs.com/gif.latex?N​时为零  的概率  。我们在这里清楚地看到近似在风险管理中的应用。https://latex.codecogs.com/gif.latex?Nhttps://latex.codecogs.com/gif.latex?\lambda=50\times80\times%20p

解决这个问题的另一种方法是基于以下思想:鉴于在对全球450座反应堆进行的45年观察中(,观察到了三起重大事故,包括“三哩岛”(1979年)和“福岛”(2011年),即两次事故之间的平均时间估计为16年。对于单个反应堆,我们可以假设事件发生之前等待的平均时间是16年的450倍,即7200年。或者,一个反应堆在一年内发生一次事件的概率是7200以上的事件之一(这是“返还期”概念背后的想法)。如果我们假设事故的到来是随机且彼此独立发生的(如上定义),则在50年内观察到的重大事故数量遵循参数为50 /(7200/80)的泊松分布。也,

https://latex.codecogs.com/gif.latex?1-\exp(-50\times%2080/7200)

> 
[1] 0.4262466


可下载资源

关于作者

Kaizong Ye拓端研究室(TRL)的研究员。在此对他对本文所作的贡献表示诚挚感谢,他在上海财经大学完成了统计学专业的硕士学位,专注人工智能领域。擅长Python.Matlab仿真、视觉处理、神经网络、数据分析。

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

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

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

关注有关新文章的微信公众号


永远不要错过任何见解。当新文章发表时,我们会通过微信公众号向您推送。

技术干货

最新洞察

This will close in 0 seconds