Estimation des paramètres de densité de Weibull (erreur: "... la valeur initiale dans 'vmmin' n'est pas finie")

Question

J'essaie d'estimer la forme et l'échelle d'un ensemble de données.
J'ai utilisé deux façons différentes et pour les deux j'ai reçu un message d'erreur:

Tout d'abord, j'ai essayé par moments en utilisant le package survey:

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

J'ai reçu le message d'erreur:

invalid survival times for this distribution

Deuxièmement, j'ai essayé d'utiliser fitdistr() fonction:

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

j'ai reçu un message d'erreur:

Error in optim(x=c(2.2, 2.1,1.9....:
initial value in 'vmmin' is not finite

Qu'est-ce qui ne va pas avec ce que je fais?

Answer 1

Une recherche google "fitdistr Weibull erreur" montre cette question exacte a été discuté il y a un an sur la liste de diffusion R-help: http://r.789695.n4.nabble.com/Problems-with-fitdistr-td1334772.html

Quelques points de ce lien:

• zéros dans votre les données risquent de créer la fonction de package LMOM

• pelwei d'utilisation()

Answer 2

J'ai eu un problème similaire lors de l'utilisation de fitdistr() avec une distribution bêta. Dans ce cas, les deux et les zéros dans les données ont produit cette erreur.

En outre, je trouve que, lorsque l'argument de la limite est utilisée dans l'appel fitdistr(), une autre erreur se produit: par exemple (où x est un vecteur d'échantillons contenant un 1,0 ou 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'`