问题标签 [fitdistrplus]

我正在用传统分布拟合数据。我正在使用 qqcomp 在 qqplot 中绘制我的拟合分布。我还想在这个图中添加数据 qqplot。

我已经尝试通过 + 获取消息以通常的 ggplot() 方式添加它们：

q+p 错误：不知道如何将 p 添加到绘图中

比方说：

我想将图 q 和 p 合二为一。我知道带有 plotstyle= "ggplot" 的 qqcomp 创建了一个 ggplot 对象，所以我想有一种方法可以将这两个图结合起来。谢谢您的帮助。

我正在尝试为 Gumbel distribtuion 运行以下命令：

出现以下错误：

我的数据头看起来像：

我对数据进行了对数转换，并尝试拟合多个分布，包括威布尔和伽马分布。我在 fitdistrplus 库中使用了 fitdist() 函数。我想使用分布的 PDF 在给定范围内使用 Weibull 分布找到期望值。对于伽马分布，很容易编写公式，因为有一个名为 dlgamma 的函数，它为我提供了对数伽马分布的 PDF。我想知道 R 中是否有类似的函数用于 log-weibull 分布？

首先，我是一个 R 新手。我正在尝试将密度图应用于我的数据中的各个组。使用 fitdistrplus，我为我的所有数据创建了一个分布密度图。

在 ggplot 中使用 facet_grid，我为每个数据分组创建了一个直方图网格。

我正在寻找的是结合两个结果，在其中我得到分组网格中每个直方图的密度图。感谢您的帮助。

样本数据：

我有一些关于特定感染流行率的数据，为每个国家/地区提供了 6 个不同年龄组。我正在尝试找到一个合适的分布，它可能适合使用 fitdistrplus 对 prev 进行建模。

流行变量的直方图显示它似乎不是单峰的，但 logit 或 log 转换表明存在一些偏斜。现在，由于起始值或我不知道的某事，我无法找到适合数据的分布。有人可以建议吗？

这是我的数据

如果这不起作用，我可以考虑使用样条曲线吗？或者可能是一些 GAM 模型，因为我有研究级别的数据，并且对于每个研究，我都有 6 个组的 prev。

当我尝试一些凝视值或将选项添加lower=c(0,0)到 fitdist 时，我遇到了一些错误the function mle failed to estimate the parameters, with the error code 100

我正在尝试使用 GLMM 分析一些反应时间数据。找到最适合我的数据的分布。我将 fitdist() 用于伽马和对数正态分布。结果表明对数正态更适合我的数据。然而，最近我读到逆高斯分布可能更适合反应时间数据。

我使用 nigFitStart 来获取起始值：

这给了我这个：

所以我尝试使用 fitdist 的启动参数：

但我收到以下错误：

我还使用了 ig_fit{goft} 并得到了以下结果：

所以，这次我将这两个参数用于 fitdist 中的 start 参数，但仍然得到完全相同的错误：

有人提到将参数名称从 mu 和 lambda 更改为 mean 和 shape 已经解决了他们的问题，但我尝试了它仍然得到同样的错误。

知道如何解决这个问题吗？或者您能否建议一种将逆高斯拟合到我的数据的替代方法？

谢谢你

我正在使用 distrplus 包中的 Fitdistcens 函数拟合指数分布。我的数据集包含完整的数据和暂停数据。但是，“exp”在下面给我一条错误消息：

这是我的代码：

有人可以帮我解决这个错误代码吗？我可以分享我的数据集。

我需要使用 fitdistrplus 但我无法区分它何时是离散的和何时是连续分布。假设我有

我们如何“告诉” fitdistrplus 它是什么分布，之后我们如何计算它的参数？

我正在尝试将某个截断的分布拟合到数据集。例如，对于对数正态分布，我将截断分布的密度函数定义为：

其中low是截断点。

我的数据是以下向量

为了适应分布，我写：

但是会出现以下错误：

我不知道发生了什么。有人可以帮助我吗？谢谢！

我想在一个图中将不同分布的拟合与我的数据进行比较。包中的qqcomp功能fitdistrplus几乎完全符合我的要求。然而，我唯一的问题是它主要是使用基本 R 情节编写的，而我所有其他情节都是用ggplot2. 我基本上只是想自定义qqcomp情节，使其看起来像是在ggplot2.

从文档（https://www.rdocumentation.org/packages/fitdistrplus/versions/1.0-14/topics/graphcomp）我知道这完全可以通过设置来实现plotstyle="ggplot"。但是，如果我这样做，则情节上不会出现任何点，即使它在没有 plotstyle 参数的情况下也能完美运行。这是一个可视化我的问题的小例子：

这给出了以下输出：

虽然这样：

给出以下输出：

为什么积分不显示？我在这里做错了什么还是这是一个错误？

编辑：

所以我还没有弄清楚为什么这不起作用，但有一个非常简单的解决方法。函数调用qqcomp(model.list, plotstyle="ggplot")仍然返回一个 ggplot 对象，其中包括用于制作绘图的数据。使用这些数据，可以轻松编写自己的绘图函数，完全按照自己的意愿行事。它不是很优雅，但是直到有人发现它为什么不能按预期工作，我才会使用这种方法。