问题标签 [fitdistrplus]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
28 浏览

r - pp plot 非线性的原因

我使用 library(fitdistrplus) 包来拟合数据,如下所示:

以下是绘图结果: 在此处输入图像描述

我的问题是:x1数据实际上是用rlnorm生成的,但是拟合后,pp图不是那么完美,如何解释?

谢谢。光明

0 投票
2 回答
109 浏览

r - 如何通过 fitdistrplus 包选择分布参数?

bootdist()之后总是用来获得更准确的参数fitdis()

如上例,summary(res) 提供几个参数,res$fitpart$estimate提供一个参数,应该选择哪一个?

0 投票
0 回答
82 浏览

r - 从年龄分组的人口数据中截断负二项分布

我有两个按年龄分类的人群的数据,每个人群都有不同的分类。

人口 1 中的年龄区间:18-24、25-29、30-34、35-45、46-60、61+

人口分布 1

人口 2 中的年龄区间:15-19、20-24、25-29、30-34 ... 85-89、90+

人口分布2

我想从这些分箱数据中推断出一个连续分布,以便更直接地比较这两个群体。我尝试拟合未截断的负二项式分布,但它低估了较低的 bin:

nbinom 适合弹出 1

所以,现在我想尝试截断的负二项分布。我做了以下事情:

但我收到以下错误:

关于如何处理/解决这个问题的任何建议?

这是pop 1的数据:

0 投票
0 回答
64 浏览

r - 如何在 R 中建立一个简单的伽马分布模型?

我目前正在尝试建立一个具有一个预测变量的模型。我有保险数据:损失金额及其保单限额。

我能够绘制带有损失量的直方图并拟合数据的伽马分布,但我不知道限制因素在哪里。我建立的唯一具有限制和损失的模型是一个简单的线性 y = mx + b 模型。

这是我第一次尝试建模,如果我缺少一些基本或基本的东西,请告诉我。

0 投票
2 回答
139 浏览

r - 用 distrplus 拟合 truncnorm 分布会导致致命错误(R 崩溃)

这是在没有任何有用错误消息的情况下使 R 崩溃的代码。有任何想法吗?我在这里要疯了。谢谢,

0 投票
1 回答
351 浏览

r - 使用 fitdist 在 R 中拟合 Pearson Type III

我已经为我的数据拟合了许多分布,但我在将 Pearson III 型分布拟合到数据时遇到了困难。我还使用 plotdist 函数来查找起始值或初始值,以便使用迭代方法拟合分布。从 plotdist 获得的图表明,该图非常适合给定起始值的数据但是 fitdist 函数不起作用并给出错误代码 100。我还研究了有关拟合 log Pearson 类型的堆栈溢出的问题和答案III 分发并应用了代码,但我再次面临运行 fitdist 函数并再次获得错误代码 100 的问题。数据可从以下链接下载。

https://ptagovsa-my.sharepoint.com/:x:/g/personal/kkhan_tga_gov_sa/EfzCE5h0jexCkVw0Ak2S2_MBWf3WUywMd1izw41r0EsLeQ?e=EiqWDc

数据可在

0 投票
2 回答
72 浏览

r - 如何返回在函数中创建的对象并忽略错误/不适用的对象?

我已经编辑了我的问题

目标

我只想保留那些成功创建的对象并忽略那些抛出错误的对象。

例子

请注意,这只是一个可重现的示例。我的原始数据集不同。

以下函数采用mtcars数据集的任何变量,拟合三个理论分布,然后返回拟合优度统计数据:

例如:

问题

一些理论分布无法成功拟合变量,并fitdist引发错误:

此错误发生在f2尝试拟合nbinom连续变量时mpg。但normgamma成功契合。

我想返回gofstat成功拟合的分布并忽略那些抛出错误的分布。

预期产出

即使f2在函数中指定了,如果它抛出错误,我仍然想要以下输出:

我试过的

显然,我可以f2从函数中删除。但这意味着重复每个变量的所有代码。那是很多代码!所以,我还是想用这个功能。

我希望能够将该函数用于任何变量。使用mtcars$mpg时,函数会失败nbinom,但使用mtcars$vs时,函数会失败gamma。在任何情况下,我都想跳过引发错误的拟合并报告gofstat有效的拟合。

我可以用purrr::possibly它悄悄地返回一个合适的结果或抛出错误而不停止在错误上。但我不知道如何仅在gofstat.

0 投票
1 回答
154 浏览

r - 寻找移位指数分布的起始值的问题

我正在尝试将移位的指数分布拟合到我的数据中,但 fitdist 函数给出的误差为 100,并且无法估计起始值。我还使用 plotdist 函数来查找起始值或初始值以拟合分布,并且在迭代过程之后,我获得了参数 rate = 0.155 shift = 0.00001 的以下图,甚至我也在 fitdist 中使用了这些值。 在此处输入图像描述 我也使用 mledist 函数来计算分布参数的起始值,但它也不起作用。我也使用了 fitdist 函数,它给出了以下错误:

fitdist(x, "sexp", start = list(rate = 0.155, shift = 1e-05)) 中的错误:函数 mle 无法估计参数,错误代码为 100

代码如下:

数据可从以下链接下载:

https://ptagovsa-my.sharepoint.com/:x:/g/personal/kkhan_tga_gov_sa/EfzCE5h0jexCkVw0Ak2S2_MBWf3WUywMd1izw41r0EsLeQ?e=EiqWDc

任何人都可以在这方面帮助我吗?

0 投票
1 回答
339 浏览

r - 当参数已知时如何使用 fitdist(帕累托分布)

我正在对一些数据进行帕累托分布拟合,并且已经估计了数据的最大似然估计。现在我需要从中创建一个 fitdist(fitdistrplus 库)对象,但我不知道该怎么做。我需要一个 fitdist 对象,因为我想使用 denscomp 等函数创建 qq、密度等图。有人可以帮忙吗?

我首先计算 MLE 的原因是因为 fitdist 没有正确执行此操作 - 即使我将正确的 MLE 作为起始值(见下文),估计值总是会膨胀到无穷大。如果早期手动给 fitdist 我的参数的选项是不可能的,fitdist 中是否有一种优化方法可以正确估计帕累托参数?

我无权发布原始数据,但这是使用 MLE 估计原始数据的伽马分布/帕累托分布的模拟。

0 投票
1 回答
112 浏览

r - Highcart R 中的 qqcomp() 和 ppcomp()

我正在尝试使用 R 中的 Highchart 包制作一些图。我的图是 qqcomp() 和 ppcomp(),我想使用 highchart 复制它。我可以使用 ggplot2 来做到这一点,但我认为它们在 Highchart 中看起来会更好。

这是我的代码: