ÎncepețiÎncepe gratuit

Vizualizarea eșantioanelor bootstrap

În acest exercițiu, vei genera eșantioane bootstrap din setul de date privind precipitațiile anuale înregistrate la Stația Meteorologică Sheffield din Marea Britanie, în perioada 1883–2015. Datele sunt stocate în array-ul NumPy rainfall, în milimetri (mm). Reprezentând grafic eșantioanele bootstrap cu ajutorul unui ECDF, vei înțelege mai bine cum permite eșantionarea bootstrap descrierea probabilistică a datelor.

Acest exercițiu face parte din cursul

Gândire statistică în Python (Partea 2)

Vezi cursul

Instrucțiuni pentru exercițiu

  • Scrie o buclă for pentru a genera 50 de eșantioane bootstrap din datele despre precipitații și pentru a reprezenta grafic ECDF-ul acestora.
    • Folosește np.random.choice() pentru a genera un eșantion bootstrap din array-ul NumPy rainfall. Asigură-te că dimensiunea (size) array-ului reeșantionat este len(rainfall).
    • Folosește funcția ecdf() pe care ai scris-o în cursul anterior pentru a genera valorile x și y ale ECDF-ului eșantionului bootstrap bs_sample.
    • Reprezintă grafic valorile ECDF. Specifică color='gray' (pentru puncte gri) și alpha=0.1 (pentru a le face semi-transparente, deoarece suprapunem foarte multe) pe lângă argumentele cheie marker='.' și linestyle='none'.
  • Folosește ecdf() pentru a genera valorile x și y ale ECDF-ului datelor originale privind precipitațiile, disponibile în array-ul rainfall.
  • Reprezintă grafic valorile ECDF ale datelor originale.
  • Apasă Trimite răspunsul pentru a vizualiza eșantioanele!

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

for _ in range(50):
    # Generate bootstrap sample: bs_sample
    bs_sample = ____(____, size=____)

    # Compute and plot ECDF from bootstrap sample
    x, y = ____
    _ = plt.plot(____, ____, ____='.', ____='none',
                 ____='gray', ____=0.1)

# Compute and plot ECDF from original data
x, y = ____
_ = plt.plot(____, ____, ____='.')

# Make margins and label axes
plt.margins(0.02)
_ = plt.xlabel('yearly rainfall (mm)')
_ = plt.ylabel('ECDF')

# Show the plot
plt.show()
Editează și rulează codul