Analiza danych – współczynnik urodzeń
Teraz wykorzystasz swoje nowe umiejętności, aby starannie dopasować model SARIMA do szeregu czasowego birth z pakietu astsa. Dane przedstawiają miesięczną liczbę żywych urodzeń (skorygowaną) w tysiącach dla Stanów Zjednoczonych w latach 1948–1979 i obejmują wyż demograficzny po II wojnie światowej.
Szereg birth jest wykreślony w konsoli R. Zwróć uwagę na długoterminowy trend (błądzenie losowe) oraz sezonowy składnik danych.
To ćwiczenie jest częścią kursu
Modele ARIMA w R
Instrukcje do ćwiczenia
- Użyj funkcji
diff(), aby zróżnicować dane (d_birth). Następnie użyjacf2(), aby wyświetlić próbkową ACF i PACF tych danych do opóźnienia 60. Zauważ sezonową persistencję. - Użyj kolejnego wywołania
diff(), aby wyznaczyć sezonową różnicę danych. Zapisz wynik jakodd_birth. Wywołaj ponownieacf2(), aby wyświetlić ACF i PACF tych danych – również do opóźnienia 60. Stwierdź, że model SARIMA(0,1,1)x(0,1,1)12 wydaje się rozsądnym wyborem. - Dopasuj model SARIMA(0,1,1)x(0,1,1)12. Co się dzieje?
- Dodaj dodatkowy parametr AR (niesezonowy,
p = 1), aby uwzględnić dodatkową korelację. Czy model dobrze dopasowuje dane?
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# Plot P/ACF to lag 60 of differenced data
d_birth <- diff(birth)
# Plot P/ACF to lag 60 of seasonal differenced data
dd_birth <- diff(d_birth, lag = 12)
# Fit SARIMA(0,1,1)x(0,1,1)_12. What happens?
# Add AR term and conclude