fitting distributions with r

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

Question

Estoy tratando de estimar la forma y la escala de un conjunto de datos.
Usé dos formas diferentes y para ambos recibí un mensaje de error:

Primero, lo intenté por momentos usando el paquete de survey :

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

Recibí el mensaje de error:

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

En segundo lugar, traté de usar la función fitdistr() :

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

Recibí un mensaje de error:

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

¿Qué está mal con lo que estoy haciendo?

Answer 1

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

Algunos puntos de ese enlace:

• ceros en sus datos causará problemas

• use la función pelwei () del paquete lmom

Answer 2

Tuve un problema similar al usar fitdistr () con una distribución Beta. En ese caso, unos y ceros en los datos produjeron este error.

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