Global uppvärmning
Nu när du har erfarenhet av att anpassa en ARIMA-modell till simulerade data är nästa steg att tillämpa dina kunskaper på verkliga data.
Data i globtemp (från astsa) visar årliga globala temperaturavvikelser fram till 2015. I den här övningen använder du etablerade metoder för att anpassa en ARIMA-modell till data. En graf över data visar beteende som liknar en slumpvandring, vilket tyder på att du bör arbeta med de differensierade data. De differensierade data diff(globtemp) är också plottade.
Efter att ha plottat urvalets ACF och PACF för de differensierade data diff(globtemp) kan du konstatera att antingen
- ACF och PACF avtar båda gradvis, vilket tyder på en ARIMA(1,1,1)-modell.
- ACF skärs av vid lag 2 och PACF avtar gradvis, vilket tyder på en ARIMA(0,1,2)-modell.
- ACF avtar gradvis och PACF skärs av vid lag 3, vilket tyder på en ARIMA(3,1,0)-modell. Även om den här modellen passar relativt bra är den sämst av de tre (du kan kontrollera det) eftersom den använder för många parametrar för så små autokorrelationer.
När du har anpassat de första två modellerna jämför du AIC och BIC för att välja den bästa modellen.
Den här övningen är en del av kursen
ARIMA-modeller i R
Övningsinstruktioner
- Plotta urvalets ACF och PACF för de differensierade data,
diff(globtemp), för att se att två modeller verkar rimliga: en ARIMA(1,1,1) och en ARIMA(0,1,2). - Använd
sarima()för att anpassa en ARIMA(1,1,1)-modell tillglobtemp. Är alla parametrar signifikanta? - Använd ett nytt anrop till
sarima()för att anpassa en ARIMA(0,1,2)-modell tillglobtemp. Är alla parametrar signifikanta? Vilken modell är bättre?
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# Plot the sample P/ACF pair of the differenced data
# Fit an ARIMA(1,1,1) model to globtemp
# Fit an ARIMA(0,1,2) model to globtemp. Which model is better?