基于R软件的统计模拟

ncp是什么  时间:2021-01-11  阅读:()

奚潭(南京财经大学统计系2006级)主要内容1.
统计模拟的基本概念2.
赶火车问题3.
R软件的统计模拟功能4.
应用R软件模拟验证大数定律5.
应用R软件模拟验证中心极限定理一、统计模拟的基本概念(一)统计模拟的定义统计模拟即是计算机统计模拟,它实质上是计算机建模,而这里的计算机模型就是计算机方法、统计模型(如程序、流程图、算法等),它是架于计算机理论和实际问题之间的桥梁.
它与统计建模的关系如下图.

实际问题统计、逻辑模型计算机模拟(程序、算法)统计、计算机解实际解一、统计模拟的基本概念(二)统计模拟方法一般地,统计模拟分类如下:若按状态变量的变化性质分为连续随机模拟和离散随机模拟.
而按变量是否随时间变化又可分为动态随机模拟和静态随机模拟.
常用的统计模拟方法主要有以下几种:1.
蒙特卡罗法2.
系统模拟方法3.
其它方法:包括Bootstrap(自助法)、MCMC(马氏链蒙特卡罗法)等.
一、统计模拟的基本概念(三)统计模拟的一般步骤二、赶火车问题0.
10.
20.
7概率13:1013:0513:00火车离站时刻一列列车从A站开往B站,某人每天赶往B站上车.
他已经了解到火车从A站到B站的运行时间是服从均值为30min,标准差为2min的正态随机变量.
火车大约下午13:00离开A站,此人大约13:30到达B站.
火车离开A站的时刻及概率如表1所示,此人到达B站的时刻及概率如表2所示.
问此人能赶上火车的概率有多大表1:火车离开A站的时刻及概率表2:某人到达B站的时刻及概率0.
10.
20.
40.
3概率13:3413:3213:3013:28人到站时刻二、赶火车问题——问题的分析——这个问题用概率论的方法求解十分困难,它涉及此人到达时刻、火车离开站的时刻、火车运行时间几个随机变量,而且火车运行时间是服从正态分布的随机变量,没有有效的解析方法来进行概率计算.
在这种情况下可以用计算机模拟的方法来解决.

:火车从A站出发的时刻;:火车从A站到B站的运行时间;:某人到达B站的时刻;:随机变量服从正态分布的均值;:随机变量服从正态分布的标准差;二、赶火车问题进行计算机统计模拟的基础是抽象现实系统的数学模型为了便于建模,对模型中使用的变量作出如下假定:此人能及时赶上火车的充分必要条件为:所以此人能赶上火车的概率模型为:二、赶火车问题为了分析简化,假定13时为时刻t=0,则变量、的分布律为:0.
10.
20.
710500.
10.
20.
40.
334323028二、赶火车问题R软件求解的总算法:关系式成立产生随机数验证模型成立次数k=k+1否是计算估计结果k/n成立次数不变试验次数是否达到n次是否编写R程序①借助区间(0,1)分布产生的随机数,对变量、概率分布进行统计模拟;②根据变量、、概率分布及模拟程序、命令产生n个随机分布数;③使用随机产生的n组随机数验证模型中的关系表达式是否成立;④计算n次模拟实验中,使得关系表达式成立的次数k;⑤当时,以作为此人能赶上火车的概率p的近似估计;进入演示windows(7,3)prb=replicate(100,括号内程序重复100次x=sample(c(0,5,10),1,prob=c(0.
7,0.
2,0.
1))y=sample(c(28,30,32,34),1,prob=c(0.
3,0.
4,0.
2,0.
1))plot(0:40,rep(1,41),type="n",xlab="time",ylabaxes=FALSE)axis(1,0:40)r=rnorm(1,30,2)points(x,1,pch=15)i=0while(i=y)points(y,1,pch=19)Sys.
sleep(0.
1)points(y,1,pch=19)title(ifelse(x+ry})mean(prb)进入模拟三、R软件的统计模拟功能1、R软件优秀的随机数模拟功能n,min,maxunif()unifomn,df,ncpt()Student'stn,lambdapois()Poissonn,mean,sdnorm()normaln,df1,df2,ncpF()Fn,rateexp()exponentialn,df,ncprchisq()chi-squaredn,size,probrbinom()binomial参数设置产生随机数序列命令分布生产某概率分布的随机数是实现统计模拟的前提条件,而使用R命令可以生成以下常用分布的随机数:三、R软件的统计模拟功能2、优良的编程环境和编程语言R所拥有的好的兼容性、拓展性和强大的内置函数有利于统计模拟的实现.

3、高效率的向量运算功能使用R拥有的向量运算功能可以大大减少程序运行的时间,提高程序运行的效率.

下面以求解Pi的程序为例加以说明未采用R向量运算功能的程序为:mc1下面用R软件分别执行两个程序,看看有什么差异程序1.
.
.
.
.
.
程序2三、R软件的统计模拟功能四、应用R软件模拟验证大数定律1、验证的大数定律有:(1)伯努利大数定理——设是次独立重复试验中事件发生的次数.
是事件在每次试验中发生的概率,则对于任意正数>0,有(2)辛钦定理:设随机变量相互独立,服从同一分布,且具有数学期望则对于任意正数,有四、应用R软件模拟验证大数定律2、在R软件实现的算法思想:由大数定律可知,当,样本的均值趋向与理论分布的期望,因此利用样本容量逐渐增大这一趋势来模拟这一趋势,在这种趋势下,样本的均值与理论分布期望的误差应该呈现出越来越小的趋势,同时,根据上述思想,分别对五种常用分布下的大数定律进行验证.

华为云年中聚惠618活动,新用户专区,云服务器低至88元/年,3年仅580.98元,热销抢购中,最后2天!

华为云怎么样?华为云用在线的方式将华为30多年在ICT基础设施领域的技术积累和产品解决方案开放给客户,致力于提供稳定可靠、安全可信、可持续创新的云服务,做智能世界的“黑土地”,推进实现“用得起、用得好、用得放心”的普惠AI。华为云作为底座,为华为全栈全场景AI战略提供强大的算力平台和更易用的开发平台。本次年终聚惠618活动相当给力,1核2G内存1m云耀云服务器仅88元/年起,送主机安全基础版套餐,...

HostKvm($4.25/月)俄罗斯/香港高防VPS

HostKvm又上新了,这次上架了2个线路产品:俄罗斯和香港高防VPS,其中俄罗斯经测试电信CN2线路,而香港高防VPS提供30Gbps攻击防御。HostKvm是一家成立于2013年的国外主机服务商,主要提供基于KVM架构的VPS主机,可选数据中心包括日本、新加坡、韩国、美国、中国香港等多个地区机房,均为国内直连或优化线路,延迟较低,适合建站或者远程办公等。俄罗斯VPSCPU:1core内存:2G...

Gcore(75折)迈阿密E5-2623v4 CPU独立服务器

部落分享过多次G-core(gcorelabs)的产品及评测信息,以VPS主机为主,距离上一次分享商家的独立服务器还在2年多前,本月初商家针对迈阿密机房限定E5-2623v4 CPU的独立服务器推出75折优惠码,活动将在9月30日到期,这里再分享下。G-core(gcorelabs)是一家总部位于卢森堡的国外主机商,主要提供基于KVM架构的VPS主机和独立服务器租用等,数据中心包括俄罗斯、美国、日...

ncp是什么为你推荐
美国虚拟主机美国云主机与美国虚拟主机有什么区别海外主机为什么国外的主机和国内的有这么大的差别?免费国外空间免费国外空间网站域名一个网站要几个域名域名服务商域名服务商所属区域怎么填写免费vps服务器如何免费搭建自己的vps服务器免费网站空间哪里有永久免费的网站空间?虚拟主机是什么什么是虚拟主机1g虚拟主机想买个1G虚拟主机,不限流量的,但不知道哪个建站网站靠谱,求推荐!天津虚拟主机天津哪个是新网互联代理呢,我打算购买邮局?
安徽虚拟主机 花生壳免费域名 如何注册中文域名 西安服务器 重庆服务器托管 仿牌空间 dreamhost la域名 20g硬盘 圣诞节促销 南昌服务器托管 qingyun 架设服务器 qq对话框 电信虚拟主机 33456 双线机房 百度云空间 域名转入 umax 更多