Augmented Dicky-Fuller
I den här övningen kör du augmented Dicky-Fuller-testet på jordbävningarnas tidsserie för att testa stationaritet. Du plottade den här tidsserien i förra övningen. Den verkade kunna vara stationär – men jordbävningar är mycket förödande. Om du vill göra förutsägelser om dem är det viktigt att vara säker.
Kom ihåg att om tidsserien inte är stationär betyder det att antalet jordbävningar per år har en trend och förändras över tid. Det vore oroande om trenden går uppåt, eftersom det innebär mer skada. Men det vore lika oroande om trenden går nedåt – det kan tyda på att jordens inre håller på att förändras, vilket i sin tur kan få stora konsekvenser för oss.
DataFrame:en med jordbävningsdata har lästs in åt dig som earthquake.
Den här övningen är en del av kursen
ARIMA-modeller i Python
Övningsinstruktioner
- Importera augmented Dicky-Fuller-funktionen
adfuller()frånstatsmodels. - Kör funktionen
adfuller()på kolumnen'earthquakes_per_year'i DataFrame:enearthquakeoch tilldela resultatet tillresult. - Skriv ut teststatistikan, p-värdet och de kritiska värdena.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# Import augmented dicky-fuller test function
from statsmodels.tsa.stattools import ____
# Run test
result = ____
# Print test statistic
print(____)
# Print p-value
print(____)
# Print critical values
print(____)