问题标签 [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 投票
0 回答
45 浏览

r - 使用 fitdistr 拟合 R 中的帕累托分布

我目前使用 R 4.0.4 版本,并且正在调试一些用 R 3.6.1 原生编写的代码。更具体地说,我有点卡在以下几点:代码的作者使用函数 fitdistrplus::fitdist 将帕累托分布拟合到一些数据,并将形状和比例参数指定为“开始”输入(如帕累托II 型分布)。当我尝试运行同一行时,我收到以下错误

'start' 必须指定作为 distr 参数的名称

深入研究 fitidt 函数,我意识到开始参数应该是(对于 Pareto)以下参数:'x、'lambda'、'a'和'log',这表明默认分布可能是 Pareto 类型 iv。有人可以给我一个见解并解释为什么代码的作者没有遇到同样的错误(我已经和他仔细检查过)。

非常感谢!

0 投票
1 回答
61 浏览

r - 如何使用 fitdist 的结果在 R 中创建概率分布?

我使用 fitdistrplus 包中的 fitdist 将(伽玛)分布拟合到我的数据中: fitg <- fitdist(mdt, "gamma")

结果是描述拟合的参数列表。我想知道是否有一种方法可以使用该结果从该分布中创建累积分布函数和随机样本生成器。

例如,如果拟合 fitdist 的分布对应于均值为 0 和 sd 1 的正态分布,我怎样才能轻松地重新创建 pnorm(..,0,1) 和 rnorm(..,0,1)?

我知道我可以手动执行此操作,但对我来说,拥有一个“自动”执行此功能会更容易,因为我必须为许多不同的数据集执行此操作,这些数据集将配备不同类型的分布。

非常感谢你的帮助!

0 投票
1 回答
32 浏览

r - tryCatch() 不会抑制错误消息

我想创建一个不打印任何错误消息的函数。

假设我有以下数据:

然后以下给出错误消息:

现在我想创建一个返回 NULL 的函数。我用 tryCatch() 试过这个。问题是 fit_fn() 仍然返回错误“计算默认起始值​​时出错”:

这样做的方法是什么?此处仅应打印 NULL:

reprex 包于 2021-11-02 创建(v2.0.1)

0 投票
0 回答
43 浏览

r - 在 R 中使用“fitdist()”函数时出错

我想使用 R 中 fitdistrplus 包中的 fitdist() 函数将 Weibull 分布拟合到我的数据中。我的数据分成箱,并且需要在每个箱中进行拟合。但是,我有一些只有一个数据值的 bin,而 fitdist() 函数不喜欢这样。我想知道是否有一种方法可以“告诉” fitdist() 不适合该垃圾箱 - 我正在考虑类似于为避免负值所做的事情,即 fitdist(x[x > 0], "weibull" )。

这是我的数据示例:

我拟合 Weibull 分布的方式如下:

0 投票
1 回答
10 浏览

r - 使用矩估计器方法在 R 中拟合 Log Pearson Type III 的问题

我是 R 的新手,我很难弄清楚如何使用方法 = mme 将 Log Pearson Type III 拟合到数据中。

我得到的错误是这样的:

使用method = mle时没有问题,但我们的研究需要使用矩量法进行参数估计。

提前感谢那些可以提供帮助的人。

0 投票
2 回答
52 浏览

r - R中截断帕累托的形状估计

length(u_NI$k_u) = 637594

我收到了这个错误:

问题出在数据集过大还是起始参数上?

可重复的例子:

这是行不通的,我不明白为什么。

这里似乎有问题:

argument 'lower' must be positive