我正在对一些数据拟合Pareto分布，并且已经估计了这些数据的最大似然估计。现在我需要从它创建一个fitdist (fitdistrplus库)对象，但我不确定该如何做。我需要一个fitdist对象，因为我想创建qq，密度等图的函数，如密集。有人能帮帮忙吗？我首先计算MLEs的原因是因为fitdist没有正确地这样做-即使我给出了正确的MLEs作为起始值，估计值也总是放大到无穷大(见下文)。如果之前手动给fi
浏览 78提问于2020-10-20得票数 0
回答已采纳
我正在尝试使用fdistrplus包中的函数fdist来拟合分布到一些数据。ppareto = function(x,alpha,theta){之后，它说它需要一个顺序
浏览 4提问于2014-06-16得票数 1
为了简要介绍背景知识，我习惯于描述火灾大小的分布，该分布假定遵循对数正态分布(小火灾多，大火灾少)。对于我的特定应用程序，我只对落入特定大小范围(> min，< max)的fires感兴趣。因此，我正在尝试将对数正态分布拟合到两端都已被审查的数据集。本质上，我希望找到对数正态分布(µ和sigma)的参数，在审查之前最适合完全分布。考虑到我只看了分布<e
浏览 2提问于2013-06-06得票数 4
回答已采纳
我试图使t-分布符合我的数据，但无法做到这一点。我的第一次尝试是有41个警告，都说NaNs是生产出来的。我不知道怎么会涉及对数。因此，我对数据做了一些调整，使所有数据都大于0，但我仍然存在相同的问题(虽然警告减少了9个.)。sstdFit()也存在同样的问题，会产生NaNs。因此，我尝试使用我在堆栈溢出和CrossValidated上看到的适配者：
fitdist(my
浏览 1提问于2015-05-27得票数 3
回答已采纳
1回答
我在包rplcon()中使用poweRlaw函数生成一些随机变量。现在，我想知道哪些已知的发行版最适合数据。洛格诺姆？经验？伽马？权力法？所以我在包fitdist()中使用函数fitdistrplusfit.gammadl <- fitdist(data, "gamma", lower
浏览 5提问于2016-05-11得票数 3
1回答
我需要计算Fréchet分布的参数。
生成以下错误
Error in computing default starting values.Error in start.arg.default(obs, distname) :
浏览 57提问于2019-11-23得票数 2
2回答
当我使用绘图时，如何使下列数据更加贴合clear all; y = [23 23 23 -7.59 23 22.82 22.40 13.54
浏览 4提问于2014-11-12得票数 0
回答已采纳
我想要定义我自己的发行版，用它来拟合我的每月降水数据，从现在起称为“月”。我正在使用“lmomco”函数来帮助我定义发行版，但无法使它工作。例如，我定义了广义极值(gev)分布，如下所示： pgev<-cdfgev
qgev<-quagev由于"fitdistrplus“需要参数"start”(由所需分布的初始参数
浏览 8提问于2015-04-27得票数 7
回答已采纳
2回答
我一直在努力将一个发行版安装到我在R中的数据样本中，我已经研究过使用fitdistr函数和fitdistr函数，但我似乎在这两个方面都遇到了问题。一个快速的背景；我的代码的输出应该是最适合的发行版(从发行版的列表中)到提供的数据，带有参数。这需要在没有人类交互的情况下发生，所以比较图表不是一种选择。我在想，我可以将每个分布与数据相匹配，从x-平方检验中提取p值，并找到p值最高<em
浏览 3提问于2015-04-24得票数 5
回答已采纳
2回答
library(VGAM) start = list(lower=1,lower = lower, upper = upper, ...): function cannot be evaluated at initial parameters> with the error code 100
浏览 10提问于2022-02-05得票数 0
回答已采纳
1回答
我有一个数据集，其中观察来自高度不同的群体。每个组可能有一个大不相同的分布，所以我尝试使用fitdistrplus中的fitdist找到最佳分布，然后使用gamlss包中的gamlssML找到最佳参数。
我的问题是在这一步之后转换数据。对于某些分布，如Box-Cox t，我可以找到使用BCT系数对数据进行规范化的方程，但对于许多
浏览 15提问于2022-09-15得票数 0
回答已采纳
我正在将波高数据集(m)拟合为广义的帕累托分布。我已经选择了0.15m的阈值。我正在从Matlab过渡到R，所以我想比较两个程序的结果。我知道这两个程序不会提供完全相同的答案，但我希望得到的结果会比较接近。 当然，对于这两个拟合，我使用了相同的数据和相同的阈值。两种方法都是通过极大似然估计来估计参数的。在Matlab中，fitdist函数不估计阈值。假设它是已知<
浏览 54提问于2020-09-16得票数 0
然后利用中心极限定理求出正态分布的参数，并绘制出正态分布的PDF。因此，直方图和PDF大致上应该是“相似的”(并且随着n的增长变得更加“相似”)。我选择了Pareto发行版，使用这个Python代码，import matplotlib.pyplot as pltfrom math import sqrt
random_value
浏览 0提问于2018-06-17得票数 4
回答已采纳
1回答
我正在将SAS代码转换为R，在SAS单变量过程中使用直方图和中点有一个使用对数正态分布的特征。结果是一个包含以下变量的表，
OBSPCT -直方图区间变量值的百分比MIDPT -直方图间隔中点在应用分布时，SAS中有一个选项可以考虑zeta、theta和sigma参数的</e
浏览 1提问于2015-09-09得票数 0
2回答
library(fitdistrplus)serving <- groundbeef$serving> res$fitpart$estimate4.00825257 0.05441911
bootdist()通常用于在fitdis()之后获得更精确的参数如上例所示，summary(res)提供了多个参数，res$fitpart
浏览 5提问于2020-05-21得票数 2
回答已采纳
scipy.stats.entropy计算连续随机变量的微分熵。用哪种估计方法和公式准确地计算微分熵？(即norm分布的微分熵与beta分布的微分熵)
下面是它的github代码。微分熵是p.d.f的负积分和。乘以日志p.d.f，但我在任何地方都看不到这个或日志。
浏览 7提问于2020-08-04得票数 2
所以在问这个问题之前，我对堆栈溢出和google做了大量的调查。我有一个我正在做的模拟，并且需要能够生成0到1.0之间的随机输入数据，这些数据遵循一定的统计分布。到目前为止，我得到了一个正态分布和一个均匀的实分布，但我仍然坚持Pareto分布。
前两种是boost/随机/，但pareto仅作为原始分布可用(即不能在变量生成器中使用)。有谁知道产
浏览 2提问于2015-03-19得票数 1
回答已采纳
我试图用R中已知的参数对R分布进行采样，我在另一个软件中使用了参数位置、β、alpha1和alpha2。我如何在R中生成样本？我一直在尝试使用函数实数：：rburr(shape1，shape2，scale)。在这里，shape1、shape2和scale似乎与已知的参数β、alpha1、alpha2相匹配。如何介绍位置参数？在R中是否有替代的参数
浏览 0提问于2018-12-14得票数 2
回答已采纳
我想在单个图中比较不同分布与我的数据的拟合程度。fitdistrplus包中的qqcomp函数几乎完成了我想要做的事情。然而，我唯一的问题是，它主要是使用基本R图编写的，而我所有的其他图都是用ggplot2编写的。我基本上只想自定义qqcomp图，使其看起来像是在ggplot2中创建的。然而，如果我这样做，绘图上没有显示任何点，即使它在没有plotstyle参数的情况下</
浏览 47提问于2019-12-30得票数 3
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
