ARMA i praktiken
Nu har du fått gedigen erfarenhet av att anpassa ARMA-modeller till data – men innan du firar, testa ytterligare en övning (i stort sett) på egen hand.
Data i oil består av råolja, WTI spotpris FOB (i dollar per fat), veckdata från 2000 till 2008. Använd dina kunskaper för att anpassa en ARMA-modell till avkastningarna. De veckovisa råoljepriserna (oil) är redan plottade åt dig. Arbeta genomgående med avkastningarna, som du beräknar själv.
Precis som tidigare är paketet astsa förladdat. Data finns tillgängliga som oil och är plottade.
Den här övningen är en del av kursen
ARIMA-modeller i R
Övningsinstruktioner
- Beräkna de ungefärliga avkastningarna för råoljepriset med hjälp av
diff()ochlog(). Spara avkastningarna ioil_returns. - Plotta
oil_returnsoch lägg märke till att det finns ett par extremvärden före 2004. Övertyga dig själv om att avkastningarna är stationära. - Plotta stickprovsautokorrelationsfunktionen (ACF) och partiella autokorrelationsfunktionen (PACF) för
oil_returnsmed hjälp avacf2()från paketetastsa. - Utifrån ACF/PACF-paret framgår det att korrelationerna är små och att avkastningarna är nästan brus. Det kan dock vara så att både ACF och PACF avtar gradvis. I så fall tyder det på en ARMA(1,1)-modell. Anpassa den här modellen till oljeavkastningarna med
sarima(). Passar modellen bra? Kan du se extremvärdena i residualplottet?
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# Calculate approximate oil returns
oil_returns <-
# Plot oil_returns. Notice the outliers.
# Plot the P/ACF pair for oil_returns
# Assuming both P/ACF are tailing, fit a model to oil_returns