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

时间:2011-01-11 10:26:47

标签: r statistics weibull

我正在尝试估算数据集的形状和比例 我使用了两种不同的方式,并且两者都收到了错误消息:

首先,我尝试使用survey包:

survreg(Surv(all.ws)~1, dist="weibull")

我收到了错误消息:

  

此分发的无效生存时间

其次,我尝试使用fitdistr()函数:

fitdistr(all.ws, densfun=dweibull, start=list(scale=1, shape=2))    

我收到错误消息:

  

优化错误(x = c(2.2,2.1,1.9 ....:
  'vmmin'中的初始值不是有限的

我在做什么有什么问题?

2 个答案:

答案 0 :(得分:6)

谷歌搜索“fitdistr Weibull错误”显示一年前在R-help邮件列表中讨论了这个确切的问题:http://r.789695.n4.nabble.com/Problems-with-fitdistr-td1334772.html

该链接的一些观点:

  • 数据中的零会导致问题

  • 使用包lmom

  • 中的pelwei()函数

答案 1 :(得分:0)

使用带有Beta版本的fitdistr()时遇到了类似的问题。在这种情况下,数据中的 1 和0都会产生此错误。

此外,我发现在fitdistr()调用中使用limit参数时,会产生不同的错误: 例如(其中x是包含1.0或0的样本的矢量):

fitdistr(x, "beta", list(shape1 = 1, shape2 = 0.2),lower=0.001)

   `Error in stats::optim: L-BFGS-B needs finite values of 'fn'`