Ich habe einige Datenanalysen in R durchgeführt und versuche herauszufinden, wie ich meine Daten an eine 3-Parameter-Weibull-Verteilung anpassen kann. Ich habe herausgefunden, wie man es mit einem 2-Parameter-Weibull macht, bin aber bei der Suche nach einem 3-Parameter zu kurz gekommen.
So passe ich die Daten mit der Funktion fitdistr () aus dem MASS-Paket an:
%Vor% x[[6]] ist eine Teilmenge meiner Daten und y speichert das Ergebnis der Anpassung.
Zuerst sollten Sie sich das FAdist-Paket ansehen. Das ist jedoch nicht so schwer von rweibull3 nach rweibull :
und ähnlich von dweibull3 bis dweibull
also haben wir das
%Vor%Bearbeiten: Wie im Kommentar erwähnt, erscheinen verschiedene Warnungen, wenn Sie versuchen, die Verteilung auf diese Weise anzupassen.
%Vor% Für mich zuerst war es nur NaNs produced , und das ist nicht das erste Mal, wenn ich es sehe, also dachte ich, dass es nicht so aussagekräftig ist, da die Schätzungen gut waren. Nach einigem Suchen schien es ein recht beliebtes Problem zu sein und ich konnte weder Ursache noch Lösung finden. Eine Alternative könnte stats4 package und mle() function verwenden, aber es scheint auch einige Probleme zu haben. Aber ich kann Ihnen anbieten, eine modifizierte Version des Codes von danielmedic zu verwenden, die ich ein paar Mal überprüft habe:
Tags und Links r model-fitting distribution weibull