Zacznij terazZacznij za darmo

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

Zobacz kurs

Instrukcje do ćwiczenia

  • Użyj funkcji diff(), aby zróżnicować dane (d_birth). Następnie użyj acf2(), 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 jako dd_birth. Wywołaj ponownie acf2(), 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

Edytuj i uruchom kod