Tukey pięć podsumowań liczb w języku Python

Question

Nie mogłem znaleźć tej funkcji w żadnym ze standardowych pakietów, więc napisałem tę poniżej. Jednak zanim rzuci się go w stronę Cheeseshop, czy ktoś wie o już opublikowanej wersji? Alternatywnie, proszę zasugerować wszelkie ulepszenia. Dzięki.

def fivenum(v): 
    """Returns Tukey's five number summary (minimum, lower-hinge, median, upper-hinge, maximum) for the input vector, a list or array of numbers based on 1.5 times the interquartile distance""" 
    import numpy as np 
    from scipy.stats import scoreatpercentile 
    try: 
     np.sum(v) 
    except TypeError: 
     print('Error: you must provide a list or array of only numbers') 
    q1 = scoreatpercentile(v,25) 
    q3 = scoreatpercentile(v,75) 
    iqd = q3-q1 
    md = np.median(v) 
    whisker = 1.5*iqd 
    return np.min(v), md-whisker, md, md+whisker, np.max(v), 

Answer 1

Chciałbym pozbyć się tych dwóch rzeczy:

import numpy as np 
from scipy.stats import scoreatpercentile 

Powinieneś imporcie na poziomie modułu. Oznacza to, że użytkownicy będą świadomi brakujących zależności od razu po zaimportowaniu modułu, a nie po wywołaniu funkcji.

try: 
    sum(v) 
except TypeError: 
    print('Error: you must provide a list or array of only numbers') 

kilka problemów z tym:

1. Nie typ kontroli w Pythonie. Podaj, jaka jest funkcja.
2. Skąd wiesz, że dzwoniący to zobaczą? Mogą nie działać na konsoli, a nawet jeśli są, mogą nie chcieć, aby komunikat o błędzie zakłócał ich wyniki.
3. Nie wpisuj opcji Python.

Jeżeli chcesz podnieść jakiś wyjątek dla nieprawidłowych danych (nie rodzaj kontroli), albo pozwolić istniejący wyjątek propagują lub owinąć go w swoim rodzaju wyjątku.

Answer 2

Jeśli ktoś potrzebuje wersji, która współpracuje z danymi NaN, oto moja modyfikacja. Nie chciałem zmienić oryginalnej odpowiedzi na plakat, aby uniknąć nieporozumień.

import numpy as np 
from scipy.stats import scoreatpercentile 
from scipy.stats import nanmedian 

def fivenum(v): 
    """Returns Tukey's five number summary (minimum, lower-hinge, median, upper-hinge, maximum) for the input vector, a list or array of numbers based on 1.5 times the interquartile distance""" 
    try: 
     np.sum(v) 
    except TypeError: 
     print('Error: you must provide a list or array of only numbers') 
    q1 = scoreatpercentile(v[~np.isnan(v)],25) 
    q3 = scoreatpercentile(v[~np.isnan(v)],75) 
    iqd = q3-q1 
    md = nanmedian(v) 
    whisker = 1.5*iqd 
    return np.nanmin(v), md-whisker, md, md+whisker, np.nanmax(v), 

Answer 3

Jestem nowy języka, jednak powrót obliczono nieprawidłowe: powinna być max (min (V), Q1-wąsy) w dolnej granicy i min (max (v) q3 + wiskerów) dla górnej granicy. To jak to się robi w R (funkcja summary()), a to, co pojawia się na boxplots w matplotlib.pyplot i R.

Answer 4

pandasSeries i DataFrame mają describe sposób, który jest podobny do R „s summary:

In [3]: import numpy as np 

In [4]: import pandas as pd 

In [5]: s = pd.Series(np.random.rand(100)) 

In [6]: s.describe() 
Out[6]: 
count 100.000000 
mean  0.540376 
std  0.296250 
min  0.002514 
25%  0.268722 
50%  0.593436 
75%  0.831067 
max  0.991971 

NAN są obsługiwane prawidłowo.

Answer 5

Spróbuj tego:

import numpy as np 
import numpy.random 
from statstools import run 
from scipy.stats import scoreatpercentile 

data=np.random.randn(5) 

return (min(data), md-whisker, md, md+whisker, max(data)) 

Answer 6

Minimal, ale dostaje zadanie. :)

import numpy as np 
[round(np.percentile(results[:,4], i), 1) for i in [1, 2, 5, 10, 25, 50]]