lmoment#
- scipy.stats.lmoment(sample, order=None, *, axis=0, sorted=False, standardize=True, nan_policy='propagate', keepdims=False)[Quelle]#
Berechnet L-Momente einer Stichprobe aus einer kontinuierlichen Verteilung
Die L-Momente einer Wahrscheinlichkeitsverteilung sind zusammenfassende Statistiken, die ähnliche Anwendungen wie herkömmliche Momente haben, aber sie sind anhand der Erwartungswerte von Ordnungsstatistiken definiert. Stichproben-L-Momente sind analog zu Populations-L-Momenten definiert und können als Schätzer für Populations-L-Momente dienen. Sie sind tendenziell weniger empfindlich gegenüber extremen Beobachtungen als herkömmliche Momente.
- Parameter:
- samplearray_like
Die realwertige Stichprobe, deren L-Momente gewünscht werden.
- orderarray_like, optional
Die (positive Ganzzahl-)Ordnungen der gewünschten L-Momente. Muss ein Skalar oder ein nicht-leeres 1D-Array sein. Standard ist [1, 2, 3, 4].
- axisint oder None, Standard: 0
Wenn es sich um eine ganze Zahl handelt, ist dies die Achse des Eingabearrays, entlang der die Statistik berechnet wird. Die Statistik jedes Achsen-Slices (z. B. Zeile) der Eingabe erscheint dann in einem entsprechenden Element der Ausgabe. Wenn
None, wird die Eingabe vor der Berechnung der Statistik geglättet.- sortedbool, default=False
Ob sample bereits auf der axis in aufsteigender Reihenfolge sortiert ist. Wenn False (Standard), wird sample sortiert.
- standardizebool, default=True
Ob L-Momentenverhältnisse für Ordnungen ab 3 zurückgegeben werden sollen. L-Momentenverhältnisse sind analog zu standardisierten herkömmlichen Momenten: Sie sind die nicht-standardisierten L-Momente geteilt durch das L-Moment der Ordnung 2.
- nan_policy{‘propagate’, ‘omit’, ‘raise’}
Definiert, wie Eingabe-NaNs behandelt werden.
propagate: Wenn ein NaN in der Achsen-Slice (z. B. Zeile) vorhanden ist, entlang der die Statistik berechnet wird, wird der entsprechende Eintrag der Ausgabe NaN sein.omit: NaNs werden bei der Berechnung weggelassen. Wenn im Achsen-Slice, entlang dem die Statistik berechnet wird, nicht genügend Daten verbleiben, wird der entsprechende Eintrag der Ausgabe NaN sein.raise: Wenn ein NaN vorhanden ist, wird einValueErrorausgelöst.
- keepdimsbool, Standard: False
Wenn dies auf True gesetzt ist, bleiben die reduzierten Achsen im Ergebnis als Dimensionen mit der Größe eins erhalten. Mit dieser Option wird das Ergebnis korrekt gegen das Eingabearray gestreut (broadcasted).
- Rückgabe:
- lmomentsndarray
Die Stichproben-L-Momente der Ordnung order.
Siehe auch
Hinweise
Seit SciPy 1.9 werden
np.matrix-Eingaben (für neuen Code nicht empfohlen) vor der Berechnung innp.ndarraykonvertiert. In diesem Fall ist die Ausgabe eine Skalar- odernp.ndarraymit geeigneter Form anstelle eines 2D-np.matrix. Ebenso werden, während maskierte Elemente von Masked Arrays ignoriert werden, die Ausgabe eine Skalar- odernp.ndarrayanstelle eines Masked Arrays mitmask=Falsesein.Referenzen
[1]D. Bilkova. „L-Moments and TL-Moments as an Alternative Tool of Statistical Data Analysis“. Journal of Applied Mathematics and Physics. 2014. DOI:10.4236/jamp.2014.210104
[2]J. R. M. Hosking. „L-Moments: Analysis and Estimation of Distributions Using Linear Combinations of Order Statistics“. Journal of the Royal Statistical Society. 1990. DOI:10.1111/j.2517-6161.1990.tb01775.x
[3]„L-moment“. Wikipedia. https://en.wikipedia.org/wiki/L-moment.
Beispiele
>>> import numpy as np >>> from scipy import stats >>> rng = np.random.default_rng() >>> sample = rng.exponential(size=100000) >>> stats.lmoment(sample) array([1.00124272, 0.50111437, 0.3340092 , 0.16755338])
Beachten Sie, dass die ersten vier standardisierten Populations-L-Momente der Standard-Exponentialverteilung 1, 1/2, 1/3 und 1/6 sind; die Stichproben-L-Momente liefern vernünftige Schätzungen.