r语言写内曼最优分配_请教：用R 做分布拟合，为什么不是最优分布？

最新推荐文章于 2022-06-10 10:06:44 发布
最新推荐文章于 2022-06-10 10:06:44 发布
阅读量298 收藏 2
点赞数
文章标签： r语言写内曼最优分配
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_34832150/article/details/112900947
版权
博主在使用R语言对损失率数据进行beta分布拟合时，发现得到的参数估计与他人使用STATA得到的结果不同，且R的K-S检验P值显著，表明R的拟合可能不如STATA。作者寻求帮助，分享了R代码和数据，以找出问题所在。
摘要由CSDN通过智能技术生成

请教各位，我用R对别人论文的损失率数据，用beta 分布拟合，参数估计结果是3.97和3.85，而原文用STATA估计的结果是6.36和5.7。做K-S检验，我的P值是0.01835，而原文的P值是0.3657，说明似乎原文的估计结果更为合理。

不知道问题出在哪里，怎么R估计出来的好像不是最优拟合参数呢？

PS：我用R中其他函数的运行结果也类似。

R程序如下：

res =fitdistr(rate, "beta", start=list(shape1 =4, shape2 =5))

ks.test(rate,"pbeta",3.97,3.85)

ks.test(rate,"pbeta",6.36,5.7)

数据为(附件为csv格式)：

rate

0.599771943        0.511062689        0.577314431        0.466203107        0.576348136        0.540878457        0.543063376        0.56423037        0.516957536        0.616257004        0.647675622        0.539647152        0.582352102        0.569858168        0.658330823        0.378698909        0.627298259        0.584641644        0.718358437        0.659514051        0.563963013        0.545025851        0.576463468    

最低0.47元/天 解锁文章
缪加
关注
根据一组数据拟合出最佳分布函数, arena, Python, Matlab, R语言
天天向上的专栏
11-22 3万+
不论是 matlab 还是 R 语言，根据一组数据拟合分布函数时，只能拟合出特定分布的参数值，并不能给出最适合这组数据的分布函数。例如，matlab 有 normfit 函数来计算一组数据拟合成正态分布时的均值和方差，wblfit 函数来给出拟合成威布尔分布的参数值，还有 gamfit, binofit 函数等。R 语言中也有 fitdistr函数。 它们都不能自动给出最佳拟合分布，个人觉得可...
R语言采用优化方法拟合曲线并计算AIC,BIC,LRT
qq_50492541的博客
11-28 1万+
文章目录前言一、RMSD(均方根误差)优化1.RMSD(均方根误差)优化代码2.运行结果二、梯度下降(Gradient Descent)1.梯度下降(Gradient Descent)优化代码2.运行结果 前言 采用冥函数为模型： Y=x^(-1) 采用优化标准和方法： RMSD(均方根误差)，梯度下降(Gradient Descent) 对比BIC和AIC 一、RMSD(均方根误差)优化 1.RMSD(均方根误差)优化代码 #冥函数RMSD powdiscrep <- function(pa.
使用R拟合分布
weixin_30530339的博客
06-29 5023
使用R拟合分布 几个常用的概率函数介绍 这里，参考R语言实战，以及[Fitting Distribution with R]的附录。 一.认识各种分布的形态 1.1 连续型随机变量的分布 首先，我们来回顾一遍各类分布函数的表达式，及其关系。 先逐一介绍与标准正态分布相关的一些分布：正态分布，卡方分布，t− 分布，F−分布，Wishart 分布。 先上个图，一睹为快。 以上几个分布之...
r语言写内曼最优分配_R语言：TOPSIS综合评价法进行多属性最优方案选择
weixin_42299888的博客
12-23 1622
TOPSIS综合评价法是进行多个对象选择排序的方法，其通过计算各对像与最优方案、最劣方案的相对距离进行排序，它用于解决多对象多属性决策问题。本文主要是介绍在R语言中如何用TOPSIS综合评价法进行多对象选择决策。TOPSIS综合评价法一般地，TOPSIS综合评价法主要包含两个步骤：计算权重和计算相对接近度。如需详细了解TOPSIS综合评价法的原理和方法，请自行百度，网上有许多非常详尽的原理说明和案...
R语言B（beta）分布
qq_27390023的博客
06-09 8245
贝塔分布（Beta Distribution) 是一个作为伯努利分布和二项式分布的共轭先验分布的密度函数，在机器学习和数理统计学中有重要应用。在概率论中，贝塔分布，也称Β分布，是指一组定义在(0,1) 区间的连续概率分布。 2.概率密度分布函数dbeta 3.累积概率函数pbeta 4.qbeta函数(pnorm的反函数) 5.beta分布参数估算 6 从平均值和方差计算beta分布参数：alpha,belta......
r语言写内曼最优分配_R中最优化函数optim
weixin_30794755的博客
12-23 2677
最优化函数optim目标函数：$$f(x_1,x_2)=(1-x_1)^2+100(x_2-x_1^2)^2$$该函数全局最小值在($x_1=1,x_2=1$)时取到。下面这种写法是因为有多个自变量函数，传入一个参数x，每个自变量用向量x的分量来表示，从而定义出目标函数。obj x1 x2 y return(y)}x1梯度：$-400*x_1*(x_2-x_1^2)-2*(1-x_1)$x2梯度：...
r语言写内曼最优分配_R语言解决最优化问题-线性优化(LP)问题
weixin_35942339的博客
12-23 901
原文：R语言解决最优化问题-线性优化(LP)问题​tecdat.cn​线性优化简介优化是一种为所有可能的解决方案找到给定问题的最佳解决方案的技术。优化使用严格的数学模型来找出给定问题的最有效解决方案。要从优化问题开始，首先确定目标非常重要。目标是绩效的量化衡量。例如：最大化利润，最小化时间，最小化成本，最大化销售。优化问题可分为两组线性规划(LP)：它也被称为线性优化，在这个问题中，目标是在数学模...
r语言分层抽样内曼分配
最新发布
10-10
R语言中可以使用"strata"函数进行分层抽样，内曼分配（Neyman allocation）是一种常用的分层抽样方法之一。 在R中，可以使用以下步骤进行分层抽样内曼分配： 1. 首先，确定需要分层的变量，并使用"strata"函数创建...
R语言-统计分布和模拟
weixin_30821731的博客
07-16 1804
R语言中统计分布和模拟 前言   很多应用都需要随机数。像interlink connection，密码系统、视频游戏、人工智能、优化、问题的初始条件，金融等都需要生成随机数。但实际上目前我们并没有“真正”的随机数生成器，尽管有一些伪随机数生成器也是非常有效的。 目录  1. 概率统计分布概述  2. 随机函数模拟介绍  3. 密度函数模拟介绍  4. 分布函数模拟介绍  5. 分位数函数模拟介绍...
如何在 R 中使用 Gamma 分布（附示例）
Mrrunsen的博客
06-03 3905
在统计学中，伽马分布通常用于模拟与等待时间相关的概率。我们可以使用以下函数来处理 R 中的 gamma 分布：dgamma(x, shape, rate) – 找到具有特定形状和速率参数的 gamma 分布的密度函数值。 pgamma(q, shape, rate) – 找到具有特定形状和速率参数的伽马分布的累积密度函数值。 qgamma(p, shape, rate) – 找到具有特定形状和速率参数的 gamma 分布的逆累积密度函数值。 rgamma(n, shape, rate) – 生成 n 个随机
R：使用R内置的mtcars数据框生成散点图及拟合曲线
追
08-27 2万+
如果以前没有绑定 mtcars 数据框，则绑定之： > attach(mtcars) 显示一下mtcars数据框的内容： > mtcars                      mpg cyl  disp  hp drat    wt  qsec vs am gear carb Mazda RX4           21.0   6 160.0 110 3.90 2.620 1
简单介绍一下R中的几种统计分布及常用模型
weixin_33769125的博客
11-23 1819
        统计学上分布有很多，在R中基本都有描述。因能力有限，我们就挑选几个常用的、比较重要的简单介绍一下每种分布的定义，公式，以及在Ｒ中的展示。 统计分布每一种分布有四个函数：d――density（密度函数），p――分布函数，q――分位数函数，r――随机数函数。比如，正态分布的这四个函数为dnorm，pnorm，qnorm，rnorm。下面我们列出各分布后缀，前面加前缀d、p、q或r就构...
应用统计学与R语言实现学习笔记（七）——拟合优度检验
胖胖雕的GIS博客
05-10 2万+
Chapter 7 Goodness of Fit本篇是第七章，内容是拟合优度检验。
R fitting R语言数据拟合总结
zefi279175732的博客
07-17 2万+
本来想翻译过来，工作了比较大先附上原文，后续再来翻译，有兴趣的可以相互交流 QQ 279175732 TKS Fit models to data This page provides tips and recommendations for fitting linear and nonlinear models to data. Updated and revised
R语言曲线拟合函数（绘图）
热门推荐
ch的专栏
09-09 9万+
本文转载自:http://blog.sina.com.cn/s/blog_8eee7fb60101g25j.html 曲线拟合：（线性回归方法：lm） 1、x排序 2、求线性回归方程并赋予一个新变量     z=lm(y~x+I(x^2)+...) 3、plot(x,y)    #做y对x的散点图 4、lines(x,fitted(z))    #添加拟合值对x的散点图并连线
r语言写内曼最优分配_R语言之内存管理
weixin_39621427的博客
01-03 855
来源 | R友舍在处理大型数据过程中，R语言的内存管理就显得十分重要，以下介绍几种常用的处理方法。1，设置软件的内存memory.size(2048) #设置内存大小 memory.size(NA) #查看当前设置下最大内存#or memory.limit() memory.size(F) #查看当前已使用的内存 #or library(pryr) mem_used() ...
常见概率分布及在R中的应用
eyejavaeye的博客
02-06 5455
  R提供工具来计算累计分布函数p(cummulative distribution function CDF),概率密度函数d和分位数函数q，另外在各种概率分布前加r表示产生随机序列 （R这种直接在分布前面加前缀的语法太难读了，pt() 误以为还是一个函数，实际上的含义是p(t())，为什么不写成这个格式呢？ 不过t()返回什么好...）   常见概率分布   离散型 1.二项分布Binomia...
R语言Γ(gamma)分布
qq_27390023的博客
06-10 1万+
Density, distribution function, quantile function and random generation for the Gamma distribution with parameters and . dgamma(x, shape, rate = 1, scale = 1/rate, log = FALSE) pgamma(q, shape, rate = 1, scale = 1/rate, lower.tail = TRUE, log.p = F
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值