问题标签 [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 回答
1337 浏览

r - mle 估计参数失败,错误代码 7

我正在尝试估计 Weibull-Gamma 分布参数,但遇到以下错误:

“函数 mle 无法估计参数,错误代码 7”

我该怎么办?

Weibull-Gamma 分布

密度函数

累积分布函数

危险功能

生存功能

估计

0 投票
1 回答
1370 浏览

r - 如何使用 fitdist 函数(负二项式)?

我有以下数据集

我想找到按set和set_up分组的负二项式函数的参数k。

fitdist(data = aa$mass, distr = "nbinom", method = "mle")$estimate[[1]]给出 k 参数的值。我想估计每组 set_up 和 set 的 k。

这是它的 dplyr 代码

我得到一个输出,但每行重复的值相同。这个估计值[[1]] 与所有大量数据都被汇集(而不是分组)相同。关于如何解决这个问题的任何建议?

0 投票
1 回答
137 浏览

r - 修正 Weibull 错误 - 函数估计参数失败,错误代码 100

我正在尝试估计 Almalki 和 Yuan 的修正 Weibull 分布 (NMW) 参数,但遇到以下错误:

AIC 的值只是非常负。出了点问题。我知道在文献中 AIC 可能是负数,但我相信估计或函数中发生了一些错误。错误可以在估计、fitdist 或类似的东西中吗?来人帮帮我?

文章

累积功能

密度函数

危险功能

生存功能

估计

图片

样本:

[ https://i.stack.imgur.com/U0KwD.png][1]图片

0 投票
1 回答
1057 浏览

r - 将双模态/二高斯分布拟合到 R 中的数据集

我有一些数据集看起来像是由两个正态分布的叠加组成,形成一个双峰图。我想估计这些数据集分布的最佳拟合参数。通常我会使用 fitdistrplus 包,但我找不到一个分布函数来提供它的算法。

有人可以指出我的方向或建议我自己如何做吗?

0 投票
0 回答
181 浏览

r - 左截断 Gamma 拟合

我正在尝试在 R 中拟合左截断伽马分布。我有超过 25,000 美元的大额索赔。低于 25,000 美元的索赔是未知的。我正在尝试通过 R 中的 MLE 程序来拟合伽马参数。有人可以为此建议代码吗?我目前的代码如下:

有人可以解释一下如何为此编写代码吗?我已经尝试使用 fitdist 函数,但我不确定截断的 gamma 函数的名称以及它需要什么开始。

失败的代码:

谢谢!

0 投票
2 回答
2481 浏览

r - FitDist 错误:'数据必须是长度大于 1 的数字向量'

我正在尝试使用以下代码拟合分布:

我收到以下错误:

fitdist(x, distr = "gamma", method = "mle") 中的错误:数据必须是长度大于 1 的数值向量

X 是一个数值变量。绘制时看起来像这样。1

为什么我会收到此错误。非常感谢任何提示。

0 投票
1 回答
133 浏览

r - 转换为数据框:模型拟合数据框的全因子子集的摘要

我有一个带有变量和两个因素的数据框

我将 Weibull 分布拟合到全因子子集并收集拟合参数,如下所示:

我本质上想要作为最终输出的是一个数据框,其中包括scale,scale_sdfactor_2列。但是我目前的输出很难转换成这样的数据框。

任何帮助深表感谢。

0 投票
0 回答
169 浏览

r - 试图在 r 中为 GLMM 选择合适的分布

我有一个相对较小的数据集(160 个观察值),其中很大一部分我的响应变量为零或非常小(例如,131/160 值为 0;范围 0-1.56)。我最初计划使用 GLMM 并选择合适的分布,但一直未能找到合适的分布。我尝试过 gamma、log-normal、Weibull 和 beta。

我想尝试(并且曾向我推荐)泊松和负二项式,以及零膨胀泊松,但是当我尝试使用fitdist(in fitdistr) 查看这些拟合时,我收到以下错误消息:

优化中的简单错误(par = vstart,fn = fnobj,fix.arg = fix.arg,obs = 数据,gr = 梯度,ddistnam = ddistname,hessian = TRUE,方法 = meth,lower = lower,upper = Upper,.. .): 'vmmin' 中的初始值不是有限的> fitdist(variable_scaled, "pois") 中的错误:函数 mle 无法估计参数,错误代码为 100

我尝试使用以下方法缩放数据以避免数据中的值恰好为零:

关于导致错误代码的任何建议?

太感谢了!

0 投票
1 回答
45 浏览

r - 使用 fitdistplus 拟合分布导致 charToDate(x) 出现错误

我在我的数据框中拟合列毒物分布但是,它一直警告这个错误:“charToDate(x) 中的错误:字符串不是标准的明确格式”

这是我的数据,我想将它拟合到 tri_1 列。即使我更改了日期的类型,仍然会导致错误。

这是我的代码: 估计 <- df %>% fitdist(data= Tri_1,distr = "pois")

它不断警告此错误:“charToDate(x) 中的错误:字符串不是标准的明确格式”

0 投票
1 回答
438 浏览

r - 将 Weibull 累积分布拟合到 R 中的质量传递数据

我有一些碎岩材料的粒度质量传递累积数据,我想使用 R 拟合 Weibull 分布。我已经使用 WEIBULL.DIST() 函数在 Excel 中使用累积开关设置为 TRUE 来完成此操作。
然后我使用 excel SOLVER 使用 RMSE 推导出 alpha 和 beta 参数以获得最佳拟合。我想在 R 中重现结果。

(见附件电子表格在这里

粒子数据和累积质量通过百分比是以下向量

这是我想拟合 Weibull 结果的图:

...根据电子表格计算直径值向量的函数

我使用 Solver 在 Excel 中确定的 Weibull alpha 和 beta 的最佳值分别为 1.41 和 3.31 所以我的问题是如何在 R 中重现此分析(不一定是 Solver 部分)但将 Weibull 拟合到该数据集?