Estimación de parámetros de densidad de Weibull (error: “... en el valor inicial '' vmmin no es finito”)

Question

Estoy tratando de estimar la forma y la escala de un conjunto de datos.
He utilizado dos maneras diferentes y por tanto me dieron un mensaje de error:

En primer lugar, he intentado por momentos utilizando el paquete de survey:

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

Me dieron el mensaje de error:

tiempos de supervivencia no válidos para esta distribución

En segundo lugar, he intentado utilizar la función fitdistr():

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

Tengo un mensaje de error:

Error en optim (x = c (2.2, 2.1,1.9 ....:
valor inicial en 'vmmin' no es finito

Lo que está mal con lo que estoy haciendo?

Answer 1

Una búsqueda de google "fitdistr Weibull de error" muestra esta pregunta exacta se discutió hace un año en la lista de R-ayuda de correo: http://r.789695.n4.nabble.com/Problems-with-fitdistr-td1334772.html

Algunos puntos de enlace:

• ceros en sus datos serán causar problemas

• uso pelwei () la función de paquete lmom

Answer 2

Yo tenía un problema similar cuando se utiliza fitdistr () con una distribución beta. En ese caso, ambos los y ceros en los datos producidos este error.

Además I encontró que cuando el argumento límite se utiliza en la llamada fitdistr (), se produce un error diferente: por ejemplo, (donde x es un vector de muestras que contienen un 1,0 o 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'`