** Using the a simple autoregressive model I calculate that 4th quarter GDP was only 14,900 billion dollars. The model also forecasts that GDP growth in the first quarter of 2011 will be at a sluggish 1.1% compared to the 30 year quarterly average of 1.47%. **At the time of this post the government is still compiling the 4th quarter GDP so the information is not public and the 1st quarter of 2011 has only barely begun. I use seasonally adjusted quarterly GDP data starting from 1947 from the St. Louis Federal Reserve Bank’s FRED database. Then set by using the Box-Jenkins methodology to fit an ARIMA model to this time series that will be the basis for the GDP forecast in this post.

**What is ARIMA?**

The ARIMA model consist of three parts: 1) The autoregressive part consist of a linear regression that establishes how past values of GDP are related to future values, 2) The “Integrated” part refers to how many times we have to take a difference to get a stationary series (see post on stationary series), and 3) the moving average part consist of how past forecast errors are related to future values of GDP. An ARIMA model will be developed using the Box-Jenkin’s methodology that will take into account past values and forecast errors to predict future GDP levels.

**What is Box-Jenkins Methodology**

I will use the Box-Jenkins methodology to fit an Autoregressive Integrated Moving Average (ARIMA) model to previous GDP data to support the forecast of future values of GDP. This methodology aids in identifying a forecast model, estimating its parameters, checking the model’s performance, and finally using it to forecast. All of these steps will be illustrated below as I develop a simple model to forecast GDP.

**BOX-JENKINS METHODOLOGY **

**1) MODEL IDENTIFICATION**

The first part of the identification phase is to ensure that the model is stationary (see previous post on stationary time series) and if not then taking the first difference until it is. The plot of GDP on time above shows that gdp is not a stationary process so first differencing will be needed to ensure this assumption of the ARIMA model is satisfied. The number of times first differencing is used will determine the order of integration in the ARIMA model.

The second difference of GDP turns out to be a stationary process. The Augmented Dickey Fuller test confirms that we can reject the null hypothesis that the second difference of GDP has a unit root (non-stationary) at the 1% significance level. The next step is to choose the ARMA model that best fits this data using a correlelogram.

The first table above is called a correlelogram of the second difference in GPD and the second table is used to identify which kind of model we should be using AR, MA, or mixture model ARIMA. Given the fact that there can be multiple models that are likely candidates I will use the (Schwartz Information Criterion) SIC statistics to choose the best one. The model with the smallest SIC value will be the better fit. Two model candidates, ranked by order of complexity, are below:

**ARIMA(1,3,8,12,13; 2; 1,3,8)**

Eviews code: ddgdp c ddgdp(-1) ddgdp(-3) ddgdp(-8) ddgdp(-12) ddgdp(-13) ddgdp(-18) u(-1) u(-3) u(-8)

**ARIMA(1,3; 2 ; 1,3)**

Eviews code: ddgdp c ddgdp(-1) ddgdp(-3) u(-1) u(-3)

The first step is to set up a forecast error which will be the basis for the Moving Average term. A simple regression on the second difference of GDP to its lag value will serve as this error term.

**ESTIMATING THE PARAMETERS OF THE MODEL**

The table above shows the regression on the smaller model ARIMA (1,3,;2,1,3) this model has a Schwarz criterion of 10.42, while the expanded model of ARIMA(1,3,8,12,13;2,1,3,8,) has a Schwarz criterion of 10.46. I choose the simple model because it has the smallest Schwarz criterion value.

**CHECKING THE MODELS PERFORMANCE**

The errors of the model appear to be a white noise process, mean 0 and constant variance. The Augmented Dickey Fuller test rejects the hypothesis that the model’s error term is non-stationary. The diagnostic check of the models performance validates unbiased estimates and the SIC criteria above used to select the simple model ensures maximum efficiency is estimation.

**FORECASTING GDP WITH THE MODEL SELECTED**

In order to use the model that has been selected a mathematical formulation of the ARIMA (1,3,;3,1,3) must be outlines. After the model has been constructed using algebra the parameters estimated in this model above will be useful in building our forecasting model. Breaking down the lags and first differences yields the final forecasting model that will be used to forecast the GDP numbers.

**Using the Autoregressive model above I calculate the 4th quarter GDP was only 14,900 billion dollars resulting in a 1% increase that quarter. The model also predicts that GDP growth in the first quarter of 2011 will grow at 1.1%. Since 1980 the average quarterly growth rate for GDP has been 1.47%, so the forecasted growth will be sluggish in the next two quarters compared to historical averages.**

thanks for the expose on arima. However, i dont comprehend the coding aspect of the gdp data used. i thought ARIMA(1,3; 2; 1,3) should be coded as:

ddgdp c ar(1) ar(3) ma(1) ma(3).

secondly, in your statement “The first step is to set up a forecast error which will be the basis for the Moving Average term. A simple regression on the second difference of GDP to its lag value will serve as this error term.” did you mean ddgp c ddgdp(-1)?

May God bless you. I am from Ethiopia. I am a student .I found this discussion informative and resourceful. Keep up , please.