基于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软件实现的算法思想:由大数定律可知,当,样本的均值趋向与理论分布的期望,因此利用样本容量逐渐增大这一趋势来模拟这一趋势,在这种趋势下,样本的均值与理论分布期望的误差应该呈现出越来越小的趋势,同时,根据上述思想,分别对五种常用分布下的大数定律进行验证.

炭云188元/年,上海CN2 VPS/2核/384MB内存/8GB空间/800GB流量/77Mbps端口/共享IP

炭云怎么样?炭云(之前的碳云),国人商家,正规公司(哈尔滨桓林信息技术有限公司),主机之家测评介绍过多次。现在上海CN2共享IP的VPS有一款特价,上海cn2 vps,2核/384MB内存/8GB空间/800GB流量/77Mbps端口/共享IP/Hyper-v,188元/年,特别适合电信网络。有需要的可以关注一下。点击进入:炭云官方网站地址炭云vps套餐:套餐cpu内存硬盘流量/带宽ip价格购买上...

IonSwitch:$1.75/月KVM-1GB/10G SSD/1TB/爱达荷州

IonSwitch是一家2016年成立的国外VPS主机商,部落上一次分享的信息还停留在2019年,主机商提供基于KVM架构的VPS产品,数据中心之前在美国西雅图,目前是美国爱达荷州科德阿伦(美国西北部,西接华盛顿州和俄勒冈州),为新建的自营数据中心。商家针对新数据中心运行及4号独立日提供了一个5折优惠码,优惠后最低1GB内存套餐每月仅1.75美元起。下面列出部分套餐配置信息。CPU:1core内存...

Virmach款低价VPS可选可以选择多个机房,新增多款低价便宜VPS主机7.2美元起

Virmach商家我们是不是比较熟悉?速度一般,但是人家价格低,而且机房是比较多的。早年的时候有帮助一个有做外贸也许需要多个机房且便宜服务商的时候接触到这个商家,有曾经帮助够买过上百台这样的低价机器。这里需要提醒的,便宜但是速度一般,尤其是中文业务速度确实不快,如果是外贸业务,那肯定是没有问题。这几天,我们有看到Virmach推出了夏季优惠促销,VPS首年8折,最低年付仅7.2美元,多机房可选,如...

ncp是什么为你推荐
云南虚拟主机云南虚拟主机,公司网站用本地客户,云南数据港怎么样?双线虚拟主机双线虚拟主机是智能的吗买域名买域名怎么买?顶级域名顶级域名有哪些备案域名购买20元左右买个已备案的域名域名抢注怎么抢注域名?子域名查询爱站网外链查询工具高防服务器租用国内最稳定高防服务器,国内千兆硬防服务器哪里租用?游戏服务器租用租用游戏服务器一般需要多少钱传奇服务器租用谁能告诉我去哪里租传奇服务器最好,要信誉要的?
长沙虚拟主机 北京域名空间 香港bgp机房 debian7 免费ftp空间申请 华为网络硬盘 52测评网 怎样建立邮箱 流量计费 服务器托管什么意思 河南移动网 web服务器安全 四核服务器 上海电信测速网站 稳定空间 免费稳定空间 后门 存储服务器 广州服务器托管 512内存 更多