#### Abstract

The study attempts to determine the impact of government policies of import of gold in India on the domestic price of gold during 2013 using Autoregressive Integrated Moving Average (ARIMA) intervention model. 2013 was an amazing year for Indian gold market where the price had reached its zenith. In April 2013, to curb a record trade deficit, India imposed an import duty of 10 percent on gold and tied imports for domestic consumption to exports, creating scarce supply of the yellow metal and boosting premiums to curtail the Current Account Deficit (CAD). The objective of the paper is to model the impact of this intervention by the government on the domestic price of Indian gold. Suitable ARIMA model is fit on the preintervention period and thereafter the effects of the interventions are analysed. The results indicate that ARIMA is the most suitable model during preintervention period. Intervention analysis reveals that there is significant decrease in domestic price of gold by 56% from 2013. The model may be used by policymakers to analyse the future of gold before framing regulations and policies.

#### 1. Introduction

Gold has always been an integral part of the socioeconomic ethos of Indian household. Demand for gold is either as bullion for investment or jewellery for gift. India is the second largest importer of gold after China. This is mainly because India limits the use of recycled gold. Import of gold has added to the burden of Indian economy. Import of gold leads to Current Account Deficit as majority of the gold imported is not being monetised. In 2014, Indian household has held more than 22000 tonnes of gold and the gold imported was worth $34 billion. If the import of gold has been lower than $34 billion, the trade deficit would have been narrower than $138 billion [1]. India has formulated various policies since 1947 with different perspectives. 2013 was an important year in the Indian gold market as we witnessed extreme variations in the price of gold. This can be attributed to the changes of import duties which were increased from 2 to 10 percent to impose restrictions in the gold market [2]. The objective of the paper is to model the impact of a government policy introduced in 2013 [3] to reduce its import on the domestic price of gold. ARIMA intervention method is compared with Holt linear trends methods.

#### 2. Review of Literature

The concept of intervention analysis was introduced by Box and Tiao [4]. In their article, the model is developed which can find the impact of an intervention and noise on a response variable using difference equation model. The article also discusses the properties of maximum likelihood estimator of parameters and two applications of the model. Chung et al. [5] in their study have modelled the financial crisis in China’s manufacturing industry that occurred in 2008. They modelled the impact of the recession that had influenced the entire world abruptly. They claim the model to be more precise at explaining and analysing. Ray et al. [6] investigated the effect of introduction of Bt Cotton variety in the year 2002 on the yield of cotton for Gujarat and Maharashtra. This study also concludes that intervention ARIMA model is superior to conventional ARIMA models. Acquah et al. [7] applied intervention ARIMA model to analyse the impact of two different interventions of which one is a step function and the other a pulse function. The objective of the study is to quantify the effects of the review of the timber export levy downwards and the regulation of trade of a wood on the trade policies of wood products in Ghana. Bianchi et al. [8] in their research to find an improved forecasting model for telemarketing centers found that ARIMA intervention model performs better than the simple models such as additive and multiplicative versions of Holt-Winters exponentially weighted moving average models.

#### 3. Materials and Methods

##### 3.1. Materials

The daily traded value of gold in Multicommodity Exchange of India (MCX) Limited from March 2011 to March 2016, consisting of 1409 observations obtained from the official website of MCX [9], is used for analysis in the study. The whole set of data is divided into two pre- and postintervention periods. This division is based on Zivot-Andrews Unit Root test. 692 observations till 29th June 2013 are considered preintervention period and the remaining 717 observations as postintervention observations.

##### 3.2. Methodology

###### 3.2.1. ARIMA Model

The ARIMA model (8) of a time series refers to the model which is stationary at level with autoregressive terms and moving average terms. It may be denoted by where – is the noise term and is the backshift operator.

ARIMA model is developed in three steps: model identification, parameter estimation, and diagnostic checking. The correlogram of a stationary series will give us a rough picture of the lags of the model. Parameters are tentatively selected at the identification stage and are estimated by method of iterative least squares. The best model among a few selected tentative models is decided on the basis of the least values of Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC). AIC and BIC values are computed as follows:where and are the parameters of ARIMA model, is the standard deviation of the model, and is the number of observations used in estimating the model.

Diagnostic checking does the tests of adequacy of the selected model. This is done by checking if the error terms of the selected model are white noise using Ljung-Box test which has null hypothesis that the noise terms are random and the test statistic where is the number of observations used in estimating the model and is the maximum number of lags.

###### 3.2.2. ARIMA Model with Intervention

The ARIMA model with intervention suggested by Box and Tiao [4] is widely used to analyse the impact of external events on time dependent variables. The methodology of intervention analysis is as follows:(i)Develop a model for the time series before intervention.(ii)Add one or more dummy variables that represent the timing of the intervention.(iii)Reestimate the model, including the new dummy variables for the entire series.(iv)Interpret the coefficients of the dummy variables as measures of the effect of the intervention.The intervention model has two components: the noise component which represents the preintervention period and the intervention component that incorporates the effect of the intervention in the model. It may be represented as where is the response series, is the indicator variable coded as a binomial variable with values 0 and 1 to identify the intervention period, is the slope parameter, is the impact parameter, is the autoregressive parameter, is the moving average parameter, and is the backshift operator.

###### 3.2.3. Holt’s Linear Trend Method

This method [10] is adopted to check if a better model can be developed by forecasting with trend. The forecast equation and the two smoothing equations developed are as follows:where denotes an estimate of the level of the series at time , denotes an estimate of the trend (slope) of the series at time , is the smoothing parameter for the level, , and is the smoothing parameter for the trend,

###### 3.2.4. Comparing the Model Performance

The forecasting performance of this model can be assessed by the Mean Absolute Percent Error (MAPE) value which is supposed to be minimum. MAPE is defined as where is the actual value and is the forecast value.

#### 4. Result and Discussion

The data used for the study has 1409 observations; the descriptive statistics of the data are depicted in Figure 1.

The series has an average of Rs. 765541.40 with standard deviation Rs. 588200.00 and skewness coefficient 1.78 which indicates a drastic fall of price after the first half of 2013. However, the assumption of normality is rejected (Jarque-Bera = 2231.61, significance value = 0.00 < 0.05). To find the exact point of break for the time series to divide the series into two, Zivot-Andrews Unit Root test is used. The result of the test indicates 29th June 2013 is a point of break.

The hypothesised intervention period, 29th June 2013, is considered in the study and it is represented in the graph (Figure 2).

The preintervention period consists of 692 observations and the postintervention period has 717 observations. The descriptive statistics of the price during the two periods are given in Table 1.

It can be observed that the mean price of gold was 1070801 till June 2013 but has depreciated by 56% to 470925.67 since June 2013. To find the impact of the changes, the -test for means can be used but it will fail to describe the trend of the data; hence, a statistically rigorous method of Intervention Model Analysis is implemented.

As the first step, the most suitable ARIMA model is fitted for the preintervention period. Prior to this, Unit Root Test is carried out to find stationarity of the series. Figure 3 shows the time series as it is and the first differenced series, which indicated stationarity at first difference. This is confirmed by ADF and KPSS test.

The result of Unit Root Test by Augmented Dickey Fuller (ADF) test and KPSS test is provided in Table 2.

The series attains stationarity at first difference. Now, the correlograms for the autocorrelation function (ACF) and partial autocorrelation function (PACF) are observed to determine the parameters and of ARIMA model. The correlogram is given in Figure 4.

The correlogram suggests the value of AR model as 1 and MA also as 1; however, few other possible values are also tested and one model is accepted based on the least values of AIC and BIC. The result of these are provided in Table 3.

ARIMA is considered as the best model since it corresponds to the least values of AIC and BIC. The model developed is represented asThe MAPE of this model is 103.4 and RMSE is 585.3. The -test for significance of the coefficients also accepts the coefficients as significant.

##### 4.1. ARIMA Intervention Model

The complete series is used for developing intervention ARIMA model. The developed model is provided in Table 4.

The developed model is considered for forecast.

##### 4.2. Diagnostic Checking

Residual analysis of the model is carried out; the result is displayed in Figures 5 and 6.

The adequacy of the model is validated using Ljung-Box test for randomness and Shapiro-Wilk’s test for normality of the residuals of the model. The value of Ljung-Box test for randomness is 0.5493, accepting the null hypothesis of randomness. Similarly, the null hypothesis of normality of the residuals by Jarque-Bera test is also accepted ( value = 0.002).

##### 4.3. Holt’s Linear Trend Method

This method is adopted to check if a better model can be developed by forecasting with trend. The forecast equation and the two smoothing equations developed are presented in Table 5.

The MAPE values are compared to find a better model. The MAPE value of ARIMA intervention model is 103.4 while that of Holt’s linear trend model is 195.67. Hence, it can be concluded that the ARIMA intervention model is better than Holt’s linear trend method.

#### 5. Conclusion

The study identified the significance and impact of the government policy on import of gold which was put forth in May 2013 by employing ARIMA intervention analysis. The findings suggest that the impact of the policy can be best portrayed by this method. The value of impact is found to be Rs. 167550.7 with a slope Rs. 5.78. This method is better than the -test for difference of means before and after the interventions. The model is compared with Holt’s linear trend method and it was found that ARIMA intervention is better. This study also informs policy makers of the impact of decisions and suggests the average value and slope of the impact.

#### Competing Interests

The authors declare that they have no competing interests.