拟合 Weibull 函数和参数估计

     2023-02-19     254

关键词:

【中文标题】拟合 Weibull 函数和参数估计【英文标题】:Fitting Weibull function and parameter estimates 【发布时间】:2021-02-25 08:38:53 【问题描述】:

我正在尝试将 R 中的函数拟合到以下数据:

y<-c(80.32000,  55.78819,  89.23141,  30.75780, 418.26000, 254.30000,  200,316.42667, 406.83435, 364.00304, 218.27867, 153.88019, 235.49971, 148.77052,
273.20171, 123.54065, 157.75650, 120.88961, 134.64092, 177.44000, 123.62948,87.03000,  63.22455, 132.62000, 120.95000, 129.60000, 116.32000,  60.49000,  66.59000)

x<-c(0,  0,  0,  0,  1,  1,  1,  3,  3,  3,  3,  3,  3,  7,  7,  7,  7,  7,  7, 11, 11, 11, 11, 16, 16, 16, 16, 16, 16)

我已经尝试了几次拟合,但都没有真正奏效……我想 Weibull 函数最适合。我已经搜索了互联网以找到解决此问题的方法,并尝试按照此处的建议调整代码:https://groups.google.com/g/r-help-archive/c/rym6b1K54-4?pli=1

nls(y~127*dweibull(x,shape,scale), start=c(shape=3,scale=100))

但我收到以下错误:

numericDeriv(form[[3L]], names(ind), env) 中的错误:评估模型时产生的缺失值或无穷大。 在 dweibull(x, shape, scale) 中:产生了 NaN。

【问题讨论】:

NaN 的出现可能是因为参数值超出范围:例如,优化器可能会尝试使用负值来表示形状或比例。对它们设置界限可能会有所帮助,或者编写一个更复杂的目标函数来处理超出范围的值(例如,将它们推回边界,或将它们反映在边界中,等等) 在我的例子中,拟合导致形状接近 1,这样一些 x 值为零时,您可能会产生 0^0。 【参考方案1】:

数据高度分散。此外,它们在 $0

例如我尝试了下面绘制的函数。但由于分数如此分散,这个选择是有争议的。

【讨论】:

非常感谢您的回答和对另一个功能的思考!我想,在从左侧到 y=-30 到 0 附近并在 0 处增加(如您给定的函数),不可能更容易地达到 x=0 作为高原状线,对吧?无论如何,我认为这种方法对我的目的来说已经足够了!

使用 R 中的矩法拟合 Weibull

】使用R中的矩法拟合Weibull【英文标题】:FittingWeibullusingmethodofmomentsinR【发布时间】:2021-06-0218:08:58【问题描述】:我正在尝试使用矩量法将Weibull分布拟合到我在RStudio中的数据。我不知道适合Weibull或Pareto等发行版所需的必要命... 查看详情

拟合 3 参数 Weibull 分布

】拟合3参数Weibull分布【英文标题】:Fittinga3parameterWeibulldistribution【发布时间】:2012-08-0219:26:51【问题描述】:我一直在R中进行一些数据分析,我试图弄清楚如何将我的数据拟合到3参数Weibull分布。我找到了如何使用2参数Weibull... 查看详情

R 中的 Weibull 参数估计,同时考虑 X(时间)和 Y(累积观察)

】R中的Weibull参数估计,同时考虑X(时间)和Y(累积观察)【英文标题】:WeibullparameterestimationinR,consideringbothX(time)andY(cumulativeobservation)【发布时间】:2021-06-0606:36:09【问题描述】:早上好,我正在尝试从以下数据集中计算R中的... 查看详情

使用 Apache Commons Math 的 Weibull 参数估计

】使用ApacheCommonsMath的Weibull参数估计【英文标题】:WeibullParameterestimationusingApacheCommonsMath【发布时间】:2018-05-2613:13:44【问题描述】:我需要使用MLE从样本数据中估计Weibull分布的形状和尺度参数。我检查了ApachecommonsMath的WeibullDi... 查看详情

使用 scipy.stats 将 Weibull 分布拟合到数据是不是表现不佳?

】使用scipy.stats将Weibull分布拟合到数据是不是表现不佳?【英文标题】:DoesfittingWeibulldistributiontodatausingscipy.statsperformpoor?使用scipy.stats将Weibull分布拟合到数据是否表现不佳?【发布时间】:2021-01-0819:33:46【问题描述】:我正在... 查看详情

用 MLE 和 Newton-Raphson 估计 Weibull

】用MLE和Newton-Raphson估计Weibull【英文标题】:EstimatingWeibullwithMLEandNewton-Raphson【发布时间】:2021-06-2015:54:54【问题描述】:我一直在尝试使用牛顿法估计二参数Weibull分布。在阅读有关使用Newton-Raphson算法的一些信息时,我发现理... 查看详情

估计 Weibull 密度参数(错误:“...'vmmin' 中的初始值不是有限的”)

】估计Weibull密度参数(错误:“...\\\'vmmin\\\'中的初始值不是有限的”)【英文标题】:EstimatingWeibulldensityparameters(error:"...initialvaluein\'vmmin\'isnotfinite")估计Weibull密度参数(错误:“...\'vmmin\'中的初始值不是有限的”)【... 查看详情

R中3参数weibull的拐点?

】R中3参数weibull的拐点?【英文标题】:inflectionpointfor3-parameterweibullinR?【发布时间】:2021-09-2919:52:30【问题描述】:R中的drc包包含自启动函数,用于将3参数Weibull模型拟合到数据。该软件包提供了3参数weibull(https://cran.r-project.org... 查看详情

Scipy Weibull 参数置信区间

】ScipyWeibull参数置信区间【英文标题】:ScipyWeibullparameterconfidenceintervals【发布时间】:2019-06-0510:12:25【问题描述】:我一直在使用Matlab使用[paramhat,paramci]=wblfit(data,alpha)将数据拟合到Weibull分布。这给出了Weibull分布的形状和尺度... 查看详情

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

】将Weibull累积分布拟合到R中的质量传递数据【英文标题】:FitaWeibullcumulativedistributiontomasspassingdatainR【发布时间】:2019-09-0613:14:22【问题描述】:我有一些碎岩材料的粒度质量传递累积数据,我想使用R拟合Weibull分布。我已经设... 查看详情

用 ggplot2 拟合带有 Weibull 曲线的散点图

】用ggplot2拟合带有Weibull曲线的散点图【英文标题】:fitascatterplotwithWeibullcurvewithggplot2【发布时间】:2021-12-1111:03:18【问题描述】:我不确定这是否是重复的问题。但我真的希望能从这里得到帮助。我想绘制如下图所示的图表,... 查看详情

从 survreg 解释 Weibull 参数

】从survreg解释Weibull参数【英文标题】:InterpretingWeibullparametersfromsurvreg【发布时间】:2012-02-2504:45:55【问题描述】:我正在尝试使用从R中的survreg估计的参数生成逆Weibull分布。我的意思是,对于给定的概率(这将是在MSExcel中实... 查看详情

尝试 MLE 拟合 Weibull 分布时 scipy.optimize.minimize 中的 RuntimeWarning

】尝试MLE拟合Weibull分布时scipy.optimize.minimize中的RuntimeWarning【英文标题】:RuntimeWarninginscipy.optimize.minimizewhentryingtoMLEfitWeibulldistribution【发布时间】:2022-01-1419:32:06【问题描述】:我有以下函数to_minimize,它应该等于Weibull分布的数... 查看详情

将分布拟合到 R 中的给定频率值

...为某些分布的密度函数的数据点。问:假设这些频点来自Weibull分布T,我如何将最佳Weibull密度函数拟合到这些点从而推断分布@987654325@参数从哪里来的?sam 查看详情

Weibull 不超过图(轴“x”)

】Weibull不超过图(轴“x”)【英文标题】:WeibullNon-exceedanceplot(axis\'x\')【发布时间】:2021-06-1616:18:15【问题描述】:我正在研究Weibull分布。我已经估计了我的数据的参数,并且得到了概率图。如以下代码所示,在概率图上,我... 查看详情

如何用matlab进行参数以及非参数概率密度估计

...参数估计的函数。 参考技术B试试cftool中的混合高斯叠加拟合 参考技术C用Ksdensity 查看详情

SciPy LeastSq 拟合优度估计器

】SciPyLeastSq拟合优度估计器【英文标题】:SciPyLeastSqGoodnessofFitEstimator【发布时间】:2011-11-2703:38:20【问题描述】:我有一个使用SciPy的leastsq函数拟合的数据表面。我想对leastsq返回后的合身质量进行一些估计。我希望这将作为函... 查看详情

6.逻辑归回

...辑回归损失函数为交叉熵)  2.自述一下什么是过拟合和欠拟合?欠拟合 因为对于给定数据集,欠拟合的成因大多是模型不够复杂、拟合函数的能力不够。 为此可以增加迭代次数继续训练、尝试换用其他算法、增... 查看详情