EDA: Magnitudes vergelijken voor en na 2010
Maak een ECDF van aardbevingsmagnitudes van 1980 tot en met 2009. Zet in dezelfde plot een ECDF van de magnitudes van aardbevingen van 2010 tot en met medio 2017. De tijden van de aardbevingen, als decimale jaren, staan in de NumPy-array time en de magnitudes in de NumPy-array mags.
Deze oefening maakt deel uit van de cursus
Casestudies in statistisch denken
Oefeninstructies
- Gebruik booleaanse indexering om de magnitudes van alle aardbevingen van vóór 2010 te selecteren en sla het resultaat op in
mags_pre. Maak op dezelfde manier eennumpy-arraymags_postmet alle magnitudes van aardbevingen in en na 2010. - Gebruik
plt.plot()met een argument*dcst.ecdf(____)om ECDF's te maken voor de magnitudes van vóór en na 2010. Vergeet niet argumenten op te geven voor de parametersmarkerenlinestyle. - Klik op 'Antwoord verzenden' om de plot te bekijken.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Get magnitudes before and after 2010
mags_pre = ____[____ < ____]
mags_post = ____[____ >= ____]
# Generate ECDFs
# Label axes and show plot
_ = plt.xlabel('magnitude')
_ = plt.ylabel('ECDF')
plt.legend(('1980 though 2009', '2010 through mid-2017'), loc='upper left')
plt.show()