Ich versuche, aus Daten in Python eine ECDF (und eine Vertrauensgrenze) zu erstellen. Ich kann das ECDF ziemlich einfach numpy
durch Sortieren und Verwenden erzeugen linspace
. Ich bin mir jedoch nicht ganz sicher, welche Konfidenzgrenzen angemessen sind, und es scheint keine integrierten Bibliotheken zu geben, die die Grenzen berechnen ( statsmodels
scheint nur die ECDF anzugeben).
Wenn ich eine punktweise Konfidenzgrenze von möchte, ist es angebracht, die DKW-Ungleichung zu verwenden, um meine Region mit zu berechnen
Wo ist die Anzahl der Beobachtungen in meiner Stichprobe? Wenn also mein ECDF ist, wären meine oberen und unteren Grenzen
MATLAB hat eine eingebaute ECDF -Funktion , aber ich hatte nicht viel Glück zu verstehen, wie man die Greenwood-Formel (unten angegeben) anwendet, um die Grenzen zu generieren.
Antworten:
Geben Sie in Matlabs Konsolentyp Folgendes ein:
Es öffnet den Quellcode im Editor.
Gehen Sie zu Zeile 194:
Dies ist der Beginn des Codeblocks, der die unteren und oberen (Vertrauens-) Grenzen berechnet :
[Flo, Fup]
. Der Codeblock ist 30 Zeilen lang und ziemlich einfach. Gepostet unten für Ihre Bequemlichkeit:Die Quadratwurzel von Greenwoods Formel, dh
wird in Zeile 203 implementiert als:
Kannst du es von hier nehmen? Gib mir Bescheid.
quelle