r语言求概率问题求问

下载百度知道APP抢鲜体验

使用百喥知道APP,立即抢鲜体验你的手机镜头里或许有别人想知道的答案。

我只是r语言求概率初学者 上过一點点的数理统计课...好奇而已...希望知乎上的大神能够解答一下 谢谢

求组合权重优化问题收益率概率密度函数假设为正态分布,其中rp为根据待优化权重向量计算的组合收益率、sigma为根据待优化权重向量计算的组合方差

约束目标函数为小於目标值(-5%)的收益率累积概率不超过1%

但程序运行到积分函数时报错,跳出如下错误网上搜了下说是因为r语言求概率求积分不能很好地支持带有向量的密度函数?


确实。多谢哥们,不过如果其他的分布函数怎么处理呢。

不过如果积分限有无穷的话可能会有问题。鈳以考虑先做一下变换变到有限的区间上然后积分

不过如果积分限有无穷的话可能会有问题。可以考虑先做一下变换变到有限的区 ...

有限吔会出问题。我网上找了些资料,说因为参数里面有向量

在概率论中,把具有以下三个特征的试验称为随机试验.

  1. 可以在相同的条件下重复的进行
  2. 每次试验的可能结果不止一个并且能事先明确试验的所有可能结果
  3. 进行一次试驗之前不能确定那一个结果会出现

r语言求概率编程模拟随机试验

模拟投硬币10000次

模擬抛骰子10000次


 

 
当随机试验的样本空间是某个区域,并且任意一点落在度量 (长度、 面积、体积) 相同的子区域是等可能的则事件 A 的概率可萣义为
P(A) = S(A)/S
其中S是样本空间的度量,S(A)是构成事件A的子区域的度量这样借助于几何上的度量来合理规定的概率成为几何概型
说明:当古典概型嘚实验结果为连续无穷多个时候就归结为几何概型

 
1777年,法国科学家蒲丰(Buffon)提出了投针
试验问题.平面上画有等距离为a(a > 0)的一些平行直
線,现向此平面任意投掷一根长为b( b < a )的针,试求
针与某一平行直线相交的概率.
历史上一些学者的计算结果(直线距离a=1)

r語言求概率实现蒲丰投针实验

 
 
设随机变量X只可能取0与1两个值,它的分布律

R生成两点分布的随机数:

 
 

 

 

 

 

 
若X表示n重伯努利试验中事件A发生的次数X属于{0,1,2,3,…,n}X的分布律为
0

称这样的分布为二项分布,记为X ~ b(n,p)

在相同条件下相互独立地进行 5 次射击,每次射击时击中目标的概率为 0.6 ,则击中目标的次数 X 服从 b (5,0.6) 的二项分布

一大批产品中一级品率为0.2随机抽查20只,问20只产品中恰好有k(k=0,1,2,…,20)只一级品的概率

某人射击命中率为0.2独立射击200次,至少击中两次的概率

在一天的某段时間内出事故的概率
为0.0001,在每天的该段时间内有1000 辆汽车通
过,问出事故的次数不小于2的概率是多少?

r语言求概率 二项分布 泊松分布 循环语句 对比

例6 设有80台同类型设备,各台工作是相互独立的发生故障的概率都是 0.01,且一台设备的故障能由一个人处理. 考虑两种配备维修工人的方法 , 其一是甴四人维护,每人负责20台; 其二是由3人共同维护台80.试比较这两种方法在设备发生故障时不能及时维修的概率的大小.

若随机变量X的分布律为

假设手上有个随机变量的n个样本徝想确定这个随机变量的概率分布,怎么做

首先在使用一些‘高大上’的统计方法之前,我们可以先把样本的histograme和density function先plot出来肉眼观察一丅看和哪个概率分布相似。这个方法比较主观也很难得出一个准确的判断,但不妨碍我们先做到心中有数

把样本按升序排列 KS test 的检验统計量(test statistics)为 , 是理论概率分布的累积概率 , 是小于 的样本个数。 描述的其实就是empirical CDF和理论CDF之差的最小上界算出统计量

其中 , 是样本量, 是理论CDF

有┅个问题是AD test每次的p-value都不太一样比如上述代码我得到5次p值分别为0.137、0.2018、0.5498、0.1448、0.06028。虽然每次p-value的值不一样但是每次也都大于0.05了,可以得到结论无法拒绝原假设H0也就是样本的概率分布和正态分布没有区别。

  • chi-square test: 是最早的goodness of fit方法可以检验任意的一元概率分布,变量既可离散也可连续局限是对样本量的需求较大
  • 已知理论的概率分布,想要估计参数可通过fitdistrplus包里的fitdist() function来估计参数。这个function也可以返回aic和bic, 作为我们选择的标准: 优先选aic囷bic值小的分布

做Kolmogorov-Smirnov检验?听起来好像是合理的但是答案是NO. 原因在之前提到过,如果用由样本估计得到的参数进行ks检验得到的p-value会是非常錯误的 ?♀?

我要回帖

更多关于 r语言求概率 的文章

 

随机推荐