语言模拟蒲丰问题_R语言小数定律的保险业应用:泊松分布模拟索赔次数

c3bf34a26d76be9e870b6e8e56957559.png

原文链接:

拓端数据科技 / Welcome to tecdat​tecdat.cn
6b2568466b2dd27927e4f7e3bda8bfb4.png

在保险业中,由于分散投资,通常会在合法的大型投资组合中提及大数定律。在一定时期内,损失“可预测”。当然,在标准的统计假设下,即有限的期望值和独立性。由于在保险业中,灾难通常很少发生,而且代价非常高昂,精算师可能有兴趣对少量事件的发生进行建模。背后的定理有时也被称为小数定律

泊松分布

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

表示一个计数随机变量,然后它是服从泊松分布,如果有

De Moivre从二项式分布的近似值获得了该分布。回想一下,二项式分布是精算科学中的标准分布,例如,用来模拟

被保险人死亡人数 。如果单个死亡概率相同,例如

,并且如果死亡是独立事件,则


而如果

,然后

再次,这是一个渐近定理,当我们有很多观察值时(

)成立,它也成立,而且出现的可能性应该非常小(因为

),这就是为什么要使用术语“ 小数”的原因。SiméonPoisson对数学近似值不感兴趣:他的主要观点是针对他正在处理的数据获得具有良好拟合优度的分布。

小数定律

与Poisson分布有关的主要定理的启发式如下:

表示iid随机变量采用值

(一般情况下,一个分量可以是时间,另一分量可以是感兴趣的上部区域,其中某些随机过程是可能)。让

。如果

作为假设

(或

更具体地假设),则

表示事件的(随机变量表征)计数

,则

可以通过带有参数的泊松分布来近似


启发式方法是,如果考虑大量观察值,并且计算给定(小)区域中有多少观察值,则此类观察值的数量就是泊松分布。


  1. n=1000

  2. polygon(c(u,rev(u)),c(v,rev(-v)),col="yellow",border=NA)

  3. I=(X^2+Y^2)<1

  4. points(X[I],Y[I],cex=.6,pch=19,col="red")

15ea596e4502b3eebe1841ca08a6ec48.png

如果我们进行一些模拟


  1. > n=1000

  2. > ns=100000

  3. > N=rep(NA,ns)

  4. >

  5. +

  6. +

  7. +

  8. +

  9. +

  10. >

  11. > mean(N)

  12. [1] 31.41257

泊松分布的参数是黄色圆盘的面积,即正方形的面积,即


  1. > lines(0:60-.5,dpois(0:60,lambda),type="b",col="red")

77839b9e5fd2f0ca1d9cdca42b224f42.png

为了获得与保险模型有关的解释,让我们

在再保险合同中表示上层,即

某些可扣除额

。让我们

来表示个人损失。然后,可以使用泊松分布对到达该上层的索赔的数量进行建模。更准确地说,如果自付额

变得非常大(和

),我们将获得极值理论中的阈值点以上模型:如果

有一个泊松分布,并在有条件的

是独立同分布的广义帕累托随机变量,然后

具有广义的极值分布。因此,超出模型(针对罕见事件)与泊松过程密切相关。

泊松过程

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

泊松分布和索赔发生

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

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

c816063de931e3b0473e8ad774b1c7be.png

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

8d4c64e4031b3818ad5248b6ac162dcc.png

泊松分布和回归期

返回期是由Emil Gumbel在水文学中介绍的,用于链接概率和持续时间。十年事件的发生概率为1/10。那么10是发生之前的平均等待时间。这并不意味着该事件不会在10年之前发生,或者必须在10年之前发生。考虑一个返回期

(以年为单位),则每年不出现的概率为

多年未发生的概率为

。通常用下表来总结此属性,

27021d52d9d98315f3b3e9d039fafde1.png

上表中的对角线非常有趣。似乎在某种程度上趋向极限值(此处为63.2%)。在n年内观察到的事件数量具有二项式分布,其概率为

,将收敛到参数为1的泊松分布。那么

,没有灾难的概率为,等于0.632。

稀有概率与泊松分布

计算稀有事件的概率时,泊松分布不断出现。例如,在50年的时间里,至少有一次在核电厂发生事故的可能性。假设在反应堆中发生事故的年概率

很小,例如0.05%。进一步假设反应堆在时间上相互独立。在50年内发生超过80个反应堆的事件的概率是

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

另一方面


  1. >

  2. >

  3. [1] 0.1812733

  4. >

  5. [1] 0.1812692

这是具有参数为的泊松分布

时为零 的概率 。我们在这里清楚地看到近似在风险管理中的应用。

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


  1. >

  2. [1] 0.4262466

参考文献
1.R语言泊松Poisson回归模型分析案例
2.R语言进行数值模拟:模拟泊松回归模型
3.r语言泊松回归分析
4.R语言对布丰投针(蒲丰投针)实验进行模拟和动态可视化
5.用R语言模拟混合制排队随机服务排队系统
6.GARCH(1,1),MA以及历史模拟法的VaR比较
7.R语言做复杂金融产品的几何布朗运动的模拟
8.R语言进行数值模拟:模拟泊松回归模型
9.R语言对巨灾风险下的再保险合同定价研究案例:广义线性模型和帕累托分布Pareto distributions

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/371827.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

THINKPHP

路径 /index.php/home/...一般路径应用或者U方法转载于:https://www.cnblogs.com/lidepeng/p/6180631.html

JavaScript下的进制转换

JavaScript下的进制转换 //十进制转其他进制 var num 99; console.log(十进制: , num); console.log(八进制:, (num).toString(8)) console.log(十六进制:, (num).toString(16)) console.log(三十二进制:, (num).toString(32))//其他转十进制 var x 110; console.log(二进制&…

Spring Security第2部分–密码加密,自定义404和403错误页面

这是Spring安全站的第二部分。 在这篇文章中&#xff0c;我将向您展示如何使用MD5加密密码以及自定义403和404状态代码错误页面。 如果您尚未阅读第1部分&#xff0c;请单击 此处 。 因为我们在这里继续第1部分项目。 下载已完成的项目&#xff1a; http : //www.mediafire.com…

浅谈 PHP 与手机 APP 开发(API 接口开发)

本文内容转载自:http://www.thinkphp.cn/topic/5023.html 这个帖子写给不太了解PHP与API开发的人一、先简单回答两个问题&#xff1a;1、PHP 可以开发客户端&#xff1f;答&#xff1a;不可以&#xff0c;因为PHP是脚本语言&#xff0c;是负责完成 B/S架构 或 C/S架构 的S部分&…

获取人口_「微科普」14亿人口数据是如何得到的?

中国经济交出了2019年终答卷GDP总量近百万亿元人均GDP突破1万美元……小伙伴们在关心经济发展的同时也非常关注人口数据14亿人口的话题嗖的一下就上了热搜大家想不想知道14亿人口的数据是怎么得到的&#xff1f;我们今天就来科普一下如何获取人口总量&#xff1f;通常情况下&am…

8.动态规划(1)——字符串的编辑距离

动态规划的算法题往往都是各大公司笔试题的常客。在不少算法类的微信公众号中&#xff0c;关于“动态规划”的文章屡见不鲜&#xff0c;都在试图用最浅显易懂的文字来描述讲解动态规划&#xff0c;甚至有的用漫画来解释&#xff0c;认真读每一篇公众号推送的文章实际上都能读得…

更改Java包名称如何改变我的系统架构

即使只是少量更改角度&#xff0c;也可能对您如何使用系统产生深远影响。 假设您正在用Java编写Web应用程序。 在系统中&#xff0c;您处理订单&#xff0c;客户和产品。 作为Web应用程序&#xff0c;您的类包括诸如Controller&#xff0c;PersonRepository&#xff0c;Custome…

静态属性_Java面试题—内部类和静态内部类的区别

内部类和静态内部类的区别内部类&#xff1a;1、内部类中的变量和方法不能声明为静态的。2、内部类实例化&#xff1a;B是A的内部类&#xff0c;实例化B&#xff1a;A.B b new A().new B()。3、内部类可以引用外部类的静态或者非静态属性及方法。静态内部类&#xff1a;1、静态…

储存与更新 access_token

做微信的项目&#xff0c;一开始就是 access_token 的申请&#xff0c;微信文档上写的比较清楚&#xff1a; 1、为了保密appsecrect&#xff0c;第三方需要一个access_token获取和刷新的中控服务器。而其他业务逻辑服务器所使用的access_token均来自于该中控服务器&#xff0c;…

Eclipse安装以及JDK环境变量配置

首先是下载Eclipse&#xff1b;点击链接打开Eclipse官网eclipse官网点击DownLoad Packages&#xff0c;注意是点击“DownLoad Packages”点击你需要的版本开始下载&#xff08;一般是64bit Eclipse IDE&#xff09;等待几秒钟&#xff0c;开始下载这样Eclipse已经下载好了&…

完整的Web应用程序Tomcat JSF Primefaces JPA Hibernate –第1部分

我们创建了这篇文章&#xff0c;将展示如何使用以下工具创建完整的Web应用程序&#xff1a;Tomcat7&#xff0c;带有Primefaces的JSF2&#xff08;Facelets和Libraries&#xff09;&#xff08;具有AutoComplete&#xff09;&#xff0c;JPA / Hibernate&#xff08;具有NxN关系…

mysql主从架构升级_实战项目——mysql主从架构的实现

一主一从1.1 环境准备&#xff1a;centos系统服务器2台、 一台用户做Mysql主服务器&#xff0c; 一台用于做Mysql从服务器&#xff0c; 配置好yum源、 防火墙关闭、 各节点时钟服务同步、 各节点之间可以通过主机名互相通信1.2 准备步骤&#xff1a;1)iptables -F && s…

FastReport.Net使用:[30]对话框使用

使用对话框需要知道的地方 1.按钮的DialogResult属性。 假如DialogResult属性值为OK的按钮被点击&#xff0c;报表将会展现后面的对话框或者报表页&#xff1b;如果属性值为None&#xff0c;则停留在当前窗体&#xff1b;如果为其他值&#xff0c;则直接退出报表打印&#xff0…

模拟聊天室显示语句保持最新显示

<!DOCTYPE html> <html lang"en"> <head> <meta charset"UTF-8"> <title>模拟聊天室显示语句保持最新显示</title> <style> *{ border-collapse: collapse; } .dialog_box{ width:400px; height: 600px; margin…

改善Java EE生产支持技能的8种方法

参与Java EE生产支持的每个人都知道这项工作可能很困难。 7/24传呼机支持&#xff0c;定期处理的多个事件和错误修复&#xff0c;来自客户和管理团队的压力&#xff0c;要求它们尽快解决生产问题并防止再次发生。 在日常工作中&#xff0c;您还必须照顾由多个IT交付团队驱动的多…

plsql连接mysql_安装了mysql和pl/sql,怎么配置让pl/sql能联接mysql数据库

64位环境下&#xff0c;使用PL/SQL Developer连接Oracle&#xff1a;•1. 下载32位Oracle InstantClient&#xff0c;并展开到某目录&#xff0c;例如C:\instantclient-basic-nt-11.2.0.2.0&#xff1b;•2. 将系统的tnsnames.ora拷贝到该目录下&#xff1b;•3. 在PLSQL Devel…

varnish基础

varnish概念 初步认识 首先来跟我学习&#xff0c;v~a~r~n~i~s~h~~ &#xff0c;学会了没有~ 当然还有很重要的一个概念&#xff0c;它是高性能缓存服务器&#xff0c;举个例子。 好比我们要去买东西&#xff0c;所有的我们需要的东西是在超市厂家生产出来的&#xff0c;我们需…

引入Spring集成

在本文中&#xff0c;我们介绍Spring Integration 。 如果您以前没有使用过Spring Integration&#xff0c;那么可能会帮助您复习Gregor Hohpe的Enterprise Integration Patterns 。 我还将推荐Josh Long 撰写的这篇出色的介绍性文章 。 上下文设置 简而言之&#xff0c; 企业…

PAT 1024. 科学计数法 (20)

科学计数法是科学家用来表示很大或很小的数字的一种方便的方法&#xff0c;其满足正则表达式[-][1-9]"."[0-9]E[-][0-9]&#xff0c;即数字的整数部分只有1位&#xff0c;小数部分至少有1位&#xff0c;该数字及其指数部分的正负号即使对正数也必定明确给出。 现以科…

Mac上Hive环境搭建

本文介绍在Mac上搭建Hive环境。 建议首先配置好Hadoop&#xff0c;搭建与配置可以参考我之前的博文Mac Hadoop的安装与配置。 当然你也可以选择使用Docker搭建环境&#xff0c;本文不作介绍。 安装 对于MacOs&#xff0c;推荐使用HomeBrew安装hive&#xff0c;一步到位。 $ bre…