2016-07-10 17 views
0

Also alle Daten, die ich tatsächlich habe, ist ein Bild eines Histogramms, aus dem ich Höhen und Bin-Breite, den Median und einen Sigma-Fehler bekommen kann.Eine Wahrscheinlichkeitsdichtefunktion finden, die ein Histogramm in Python reproduziert

Das Histogramm ist verzerrt, daher sind das 16. und 84. Quantil nicht symmetrisch. Ich habe festgestellt, dass der Median und die Fehler mit einer verzerrten Gauss-Funktion repliziert werden können, aber das resultierende Histogramm aus meiner gefundenen PDF ist schwer zu vergleichen, egal wie viel ich mit den Bin-Zahlen und den Bin-Breiten spiele.

Ich verstehe, dass ich das Histogramm möglicherweise nicht genau neu erstellen kann, aber ich werde sehr glücklich mit etwas sein, das nahe genug ist.

Meine beste Idee ist es, mögliche Parameter des Gaußschen Schräglaufs zu durchlaufen, ein Histogramm zu erstellen, den Unterschied irgendwie zu quantifizieren (wie den Höhenunterschied an allen Punkten) und den besten zu finden. Ich denke, dass das ein sehr langer Prozess sein könnte und ich bin mir sicher, dass es etwas in scipy gibt, das das schneller macht. Bitte beziehen Sie sich auf etwas Nützliches, wenn möglich.

+0

Wenn Sie eine Verteiler gesichtet haben, die passen könnten Sie das Histogramm eine normalisieren könnte Dichte und dann eine Anpassung der kleinsten Quadrate der Parameter an dieses Histogramm. Dies ist jedoch alles andere als optimal - durch die Verwendung nur des Histogramms haben Sie bereits viele Informationen verloren. – cel

Antwort