问题标签 [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.
r - pp plot 非线性的原因
我使用 library(fitdistrplus) 包来拟合数据,如下所示:
以下是绘图结果: 在此处输入图像描述
我的问题是:x1数据实际上是用rlnorm生成的,但是拟合后,pp图不是那么完美,如何解释?
谢谢。光明
r - 如何通过 fitdistrplus 包选择分布参数?
bootdist()
之后总是用来获得更准确的参数fitdis()
。
如上例,summary(res)
提供几个参数,res$fitpart$estimate
提供一个参数,应该选择哪一个?
r - 从年龄分组的人口数据中截断负二项分布
我有两个按年龄分类的人群的数据,每个人群都有不同的分类。
人口 1 中的年龄区间:18-24、25-29、30-34、35-45、46-60、61+
人口 2 中的年龄区间:15-19、20-24、25-29、30-34 ... 85-89、90+
我想从这些分箱数据中推断出一个连续分布,以便更直接地比较这两个群体。我尝试拟合未截断的负二项式分布,但它低估了较低的 bin:
所以,现在我想尝试截断的负二项分布。我做了以下事情:
但我收到以下错误:
关于如何处理/解决这个问题的任何建议?
这是pop 1的数据:
r - 如何在 R 中建立一个简单的伽马分布模型?
我目前正在尝试建立一个具有一个预测变量的模型。我有保险数据:损失金额及其保单限额。
我能够绘制带有损失量的直方图并拟合数据的伽马分布,但我不知道限制因素在哪里。我建立的唯一具有限制和损失的模型是一个简单的线性 y = mx + b 模型。
这是我第一次尝试建模,如果我缺少一些基本或基本的东西,请告诉我。
r - 用 distrplus 拟合 truncnorm 分布会导致致命错误(R 崩溃)
这是在没有任何有用错误消息的情况下使 R 崩溃的代码。有任何想法吗?我在这里要疯了。谢谢,
r - 使用 fitdist 在 R 中拟合 Pearson Type III
我已经为我的数据拟合了许多分布,但我在将 Pearson III 型分布拟合到数据时遇到了困难。我还使用 plotdist 函数来查找起始值或初始值,以便使用迭代方法拟合分布。从 plotdist 获得的图表明,该图非常适合给定起始值的数据但是 fitdist 函数不起作用并给出错误代码 100。我还研究了有关拟合 log Pearson 类型的堆栈溢出的问题和答案III 分发并应用了代码,但我再次面临运行 fitdist 函数并再次获得错误代码 100 的问题。数据可从以下链接下载。
数据可在
r - 如何返回在函数中创建的对象并忽略错误/不适用的对象?
我已经编辑了我的问题
目标
我只想保留那些成功创建的对象并忽略那些抛出错误的对象。
例子
请注意,这只是一个可重现的示例。我的原始数据集不同。
以下函数采用mtcars
数据集的任何变量,拟合三个理论分布,然后返回拟合优度统计数据:
例如:
问题
一些理论分布无法成功拟合变量,并fitdist
引发错误:
此错误发生在f2
尝试拟合nbinom
连续变量时mpg
。但norm
和gamma
成功契合。
我想返回gofstat
成功拟合的分布并忽略那些抛出错误的分布。
预期产出
即使f2
在函数中指定了,如果它抛出错误,我仍然想要以下输出:
我试过的
显然,我可以f2
从函数中删除。但这意味着重复每个变量的所有代码。那是很多代码!所以,我还是想用这个功能。
我希望能够将该函数用于任何变量。使用mtcars$mpg
时,函数会失败nbinom
,但使用mtcars$vs
时,函数会失败gamma
。在任何情况下,我都想跳过引发错误的拟合并报告gofstat
有效的拟合。
我可以用purrr::possibly
它悄悄地返回一个合适的结果或抛出错误而不停止在错误上。但我不知道如何仅在gofstat
.
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
代码如下:
数据可从以下链接下载:
任何人都可以在这方面帮助我吗?
r - 当参数已知时如何使用 fitdist(帕累托分布)
我正在对一些数据进行帕累托分布拟合,并且已经估计了数据的最大似然估计。现在我需要从中创建一个 fitdist(fitdistrplus 库)对象,但我不知道该怎么做。我需要一个 fitdist 对象,因为我想使用 denscomp 等函数创建 qq、密度等图。有人可以帮忙吗?
我首先计算 MLE 的原因是因为 fitdist 没有正确执行此操作 - 即使我将正确的 MLE 作为起始值(见下文),估计值总是会膨胀到无穷大。如果早期手动给 fitdist 我的参数的选项是不可能的,fitdist 中是否有一种优化方法可以正确估计帕累托参数?
我无权发布原始数据,但这是使用 MLE 估计原始数据的伽马分布/帕累托分布的模拟。
r - Highcart R 中的 qqcomp() 和 ppcomp()
我正在尝试使用 R 中的 Highchart 包制作一些图。我的图是 qqcomp() 和 ppcomp(),我想使用 highchart 复制它。我可以使用 ggplot2 来做到这一点,但我认为它们在 Highchart 中看起来会更好。
这是我的代码: