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

r - 'fitdistrplus' 包的问题,​​t-distribution

我正在尝试将 t 分布拟合到我的数据中,但无法这样做。我的第一次尝试是

有 41 个警告,都说产生了 NaN。我不知道如何,似乎涉及对数。所以我稍微调整了我的数据,使所有数据都 > 0,但我仍然遇到同样的问题(虽然少了 9 个警告......)。sstdFit() 也有同样的问题,会产生 NaN。

因此,我尝试使用我在 stackoverflow 和 CrossValidated 上看到的 fitdist:

然后我得到

这是什么意思?我尝试查看文档,但什么也没告诉我。我只是想可能适合 t 分布,这太令人沮丧了:P

谢谢!

0 投票
1 回答
3095 浏览

r - 使用 ggplot2 制作 fitdist 图

fitdist我用包中的函数拟合了正态分布fitdistrplus。使用denscomp, qqcomp,我们可以分别绘制, , , 和cdfcomp,如下所示。ppcomphistogram against fitted density functionstheoretical quantiles against empirical onesthe empirical cumulative distribution against fitted distribution functionstheoretical probabilities against empirical ones

在此处输入图像描述

在此处输入图像描述

在此处输入图像描述

在此处输入图像描述

我对制作这些fitdist情节非常感兴趣ggplot2。MWE如下:

在此处输入图像描述

我怎样才能开始制作这些fitdist情节ggplot2

0 投票
2 回答
5663 浏览

r - 测试伽马分布的向量

我对某个向量有疑问。我想知道它是否是伽马分布的,以及(如果是的话)参数(形状、速率)是什么。我的向量有 400 个条目,但让我们举个例子

我读了一些关于 fitdistr() 的东西。但我不太明白它的实际作用!我用我的真实(长)向量尝试了以下操作:

输出是什么意思?这些是我的拟合参数吗?我对它们进行了测试,但 KS-Test 给了我一个否定的结果:

那么您能否告诉我如何确定我的向量是否是伽马分布的以及参数是什么?

0 投票
2 回答
815 浏览

r - 使用 fitdist 函数(fitdistrplus 包)估计比例和形状参数的缩放

如标题中所述,我fitdist在 R(fitdistrplus包)中的函数存在缩放问题。

请看下面的代码:

我有三个数据向量,它们存储在一个数据框中。现在我想使用该fitdist函数分别估计每列数据(和)的比例和形状参数data1,最后分别将它们存储在和中。data2data3fit_store_scalefit_store_shape

这里的问题是,fitdist如果没有适当的比例因子,该函数将无法工作data1data2并且data3需要不同的因子。我正在寻找一种解决方案来自动为每列数据确定最佳比例因子,从而让fitdist函数最终工作。

0 投票
0 回答
1500 浏览

r - R - 用 fitdistrplus 拟合混合分布

我已经定义了 LogNormal 和 Generalized Pareto (actuar) 的混合,如下所示:

我正在尝试使用 fistdistrplus 将其拟合到一些数据中:

但是我收到错误:

如果我以通常的 fitdistrplus 方式扩展拟合函数,则错误仍然存​​在:

我应该如何将我的参数传递给拟合函数?

谢谢

西蒙

0 投票
0 回答
85 浏览

r - 在 R 中的 tryCatch 命令中使用 fitdist 时遇到的问题

我想了解示例代码中 tryCatch 的行为:

fitdist命令在命令之外执行时tryCatch会给出与nbinom拟合相关的参数。但是,当我在 try catch 命令中执行相同操作时,我收到以下警告,

有人可以帮我理解这个观察吗?

编辑:我创建了这个测试用例作为我在运行循环以适应跨大数据集的多个列的分布时遇到的问题的样本。我需要使用 tryCatch 来迭代多个分布和分布拟合方法,如 mle、mme 等我面临的主要问题是 fitdist 在独立执行但不在 tryCatch 命令中时表现如预期。

0 投票
0 回答
293 浏览

r - R fitdistrplus:绘图问题 - 'by' 参数太小

我想了解绘制 fitdistrplus 对象的问题:

我的分布拟合总结如下:

我显然看不出合适,因为我在 2 个参数估计的标准误差中看到了“NA”。plot 命令出现以下错误,我想了解一下:

请帮忙。

0 投票
3 回答
2906 浏览

r - 用 fitdistrplus 拟合 Gumbel 分布

我正在尝试从这个答案中重现代码,但是这样做有问题。我正在使用包中的 gumbel 发行版VGAMfitdistrplus. 这样做时会出现问题:

好像location而且scale不是 *gumbel 的论点。

dgumbel, pgumbel,rgumbelqgumbel由 正确VGAM。然而,该包还提供了一个名为 的函数gumbel,具有不同的语法。这可能会导致问题吗?

编辑:是的,它确实引起了问题:使用包FAdist代替工作得很好。

0 投票
1 回答
1218 浏览

r - 将变量保存在R中的矩阵中

是否可以将变量存储在 R..

我必须计算 10*324 fitdist,该函数的每个输出都必须存储在具有上述大小的矩阵中?这在R中可能吗?

如果是这样,那怎么办,我完全迷路了?

所以我尝试创建一个简单的示例

但我收到了这个错误

我不确定我明白这意味着什么......

0 投票
1 回答
582 浏览

r - Plot multiple fitdist objects in same the plot with different colors?

I have a list of fitdist objects, which I have stored using this piece of code.

Is there a way to plot all the normal distributions extracted from fittest with a different color to show how the data is distributed?

Or in general how to visualize multiple normal distributions?