Forecasting with AutoRegressive (AR) Model in R

Now that we know how to estimate the AR model using ARIMA, we can create a simple forecast based on the model.

Step 1: Fit the model

The first step is to fit the model as ARIMA(1, 0, 0). We have already seen this in the previous lesson.

Step 2: Create Forecast

We can now use the predict() function to create a forecast using the fitted AR model. It takes as its inputs, the model object that we created in step 1, and an additional parameter n.ahead which establishes the forecast horizon, that is, how many steps (periods) in the future we want to create the forecast. In our example, we will provide n.ahead=20, which will create forecast for next 20 steps which corresponds to 20 days for our daily data.

The object generated by the predict() command contains two time series: 1) $pred which contains the forecasted values and 2) $se which contains the standard error for the forecast. We will use the $pred time series to plot the forecast and the $se time series to add confidence intervals to our plot.

For simplicity sake, let’s also extract the two series in their own respective variables.

Step 3: Plot the Forecast

We can now use the plot.ts() function to first plot the original data and then add points for the forecasted values using the points() function as shown below:

Notice that while creating the initial plot, we’ve specified scale limits for x and y axis in order to provision for the forecast values.

Step 4: Add Confidence Intervals to Forecast

We can add a 95% confidence interval to our forecast using the standard error values.

You have successfully created your first forecast.

Series Navigation

‹ Estimating AutoRegressive (AR) Model in RMoving Average (MA) Model in R ›