parametri di densità Estimating Weibull (errore: “... valore iniziale in 'vmmin' non è finita”)

Question

Sto cercando di valutare la forma e la portata di un insieme di dati.
Ho usato due modi diversi e per entrambi ho ricevuto un messaggio di errore:

Per prima cosa, ho provato da momenti che utilizzano il pacchetto survey:

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

ho ricevuto il messaggio di errore:

  

tempi di sopravvivenza non validi per questa distribuzione

In secondo luogo, ho provato ad utilizzare la funzione fitdistr():

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

ho ricevuto un messaggio di errore:

  

Errore in Optim (x = c (2.2, 2.1,1.9 ....:
  valore iniziale in 'vmmin' non è finita

Cosa c'è di sbagliato con quello che sto facendo?

Answer 1

Una ricerca su google "fitdistr Weibull errore" mostra questa domanda esatta è stato discusso un anno fa sul R-help mailing list: http://r.789695.n4.nabble.com/Problems-with-fitdistr-td1334772.html

Alcuni punti da quel link:

• zeri nei dati causeranno problemi

• uso pelwei () funzione dal pacchetto lmom

Answer 2

Ho avuto un problema simile quando si utilizza fitdistr () con una distribuzione Beta. In questo caso entrambi i quelli e zero nei dati prodotto questo errore.

Inoltre ho scoperto che quando l'argomento limite viene utilizzato nella chiamata fitdistr (), viene prodotto un errore diverso: es (dove x è un vettore di campioni contenenti 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'`