2011-01-11 77 views
5

Ich versuche, die Form und den Maßstab eines Datensatzes zu schätzen.
habe ich zwei verschiedene Arten und für beide bekam ich eine Fehlermeldung:Schätzung der Weibull-Dichte-Parameter (Fehler: "... Anfangswert in 'vmmin' ist nicht endlich")

Zuerst habe ich versucht, durch Momente mit dem survey Paket:

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

ich die Fehlermeldung bekam:

invalid survival times for this distribution

Zweitens versuchte ich mit fitdistr() Funktion:

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

Ich habe eine Fehlermeldung:

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

Was mit dem, was falsch ist, ich tue?

+2

versuchen Sie vielleicht auf http://www.crossvalidated.com. –

Antwort

6

Eine Google-Suche "fitdistr Weibull Fehler" genau diese Frage zeigt, wurde vor einem Jahr auf der R-Hilfe-Mailingliste diskutiert: http://r.789695.n4.nabble.com/Problems-with-fitdistr-td1334772.html

Einige Punkte aus diesem Link:

  • Nullen in Ihrem Daten werden Probleme verursachen

  • Verwendung pelwei() Funktion von Paket lmom

+0

danke !! die Nullen ist das Problem, ich verstehe nicht, warum – eliavs

+1

anscheinend, weil es ein Protokoll dauert, nicht sicher, obwohl. –

+0

prasad danke amillion die pelwei funktion funktioniert wunderbar !!! – eliavs

0

Ich hatte ein ähnliches Problem bei der Verwendung von fitdistr() mit einer Beta-Verteilung. In diesem Fall erzeugten sowohl als auch und Nullen in den Daten diesen Fehler.

Zusätzlich fand ich, dass, wenn die Grenze Argument in dem fitdistr() Aufruf verwendet wird, wird ein anderer Fehler erzeugt hat: z (wobei x ein Vektor der Proben mit einem 1,0 oder 0 enthalten):

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

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