Open Journal of Statistics
Vol.09 No.02(2019), Article ID:92120,9 pages

Empirical Analysis of Commercial Housing Sales Based on EARCH(1,1) Model

Shichang Shen, Chao Feng

School of Mathematics and Statistics, Qinghai Nationalities University, Xining, China

Copyright © 2019 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

Received: April 4, 2019; Accepted: April 25, 2019; Published: April 28, 2019


Since the 1980s, China’s commercial housing market has shown an unprecedented rapid development, and the commercial houses still has a high price. This paper studies the sales rate of commercial housing sales to find an appropriate model, and it analyzes the volatility of the commercial housing market to describe the sustainable development of the commercial housing market. By selecting month data of China’s commercial housing sales from January 2006 to October 2018, this paper uses EViews7.2 and the ARMA Model as the tool in order to establish EARCH(1,1) through the method of quantitative analysis. It is found that the yield of commercial housing sales has obvious cluster, asymmetry and leverage effect, and the impact of adverse news on the commercial housing market is more significant than the impact of favorable news.


Sales Volume of Commercial Housing, ARMA Model, EARCH Model, Leverage Effect

1. Introduction

In recent years, China’s development rate has been very rapid, and the people’s living standards have also been significantly improved. China’s commercial housing market has ushered in an unprecedented period of development especially in the context of the establishment of a comprehensive well-off society. Many domestic scholars have conducted a lot of researches for the future development prospects of the commercial housing market. Wei Junjia and Zhang Chi [1] (2014) explored the commercial housing in Nanning through the method of cointegration analysis and stationarity analysis. Shi Huiling [2] (2018) used the cross-sectional data of commercial housing sales and GDP in 30 provinces to explore the relationship between commercial housing prices and economic development. Zheng Lan [3] (2006) used the VAR model to discuss the transmission effect of monetary policy in the commercial housing market. Tian Tian and Dong Weijuan [4] (2009) studied the impact of macroeconomics on the demand for commercial housing. None of the above documents have studied the volatility of commercial housing market. This paper establishes commodity housing sales rate of return sequence and also makes a series of quantitative analysis of the residual of the model in order to EARCH model. The development status of China’s commercial housing sales market is analyzed, and relevant suggestions are put forward as a reference.

2. Theoretical Basis

2.1. ARCH Effect Test

If the random perturbation term of the subject model at-ARCH(q) [5] , the regression equation can be established:

a t = α 0 + α 1 a t 1 2 + + α q a t q 2 (1)

The null hypothesis and alternative hypothesis to be tested are:

H 0 : α 1 = α 2 = = α q = 0 H 1 : α i 0 ( 1 i q ) (2)

Test statistic

L M = n R 2 ~ χ 2 ( q ) (3)

where, n is the number of sample data; R2 is the determinant of the auxiliary regression formula. When significance level α and degree of freedom q are known, if L M > χ α 2 ( q ) , the null hypothesis is not accepted, indicating that there is ARCH effect in this sequence. Otherwise, there is no ARCH effect.

2.2. EARCH Model

The EARCH model is also called exponential GARCH model, which was proposed by Nelson [6] in 1991. The conditional variance expression of the model is:

log ( h t ) = α 0 + j = 1 p θ j log ( h t j ) + i = 1 q [ α i | ε t i h t i | + φ i ε t i h t i ] (4)

where the expression of h t is:

h t = α 0 + α 1 ε t 1 2 + + α q ε t q 2 (5)

If φ 0 , the information function is asymmetric. When φ is less than 0, leverage is significant.

2.3. Skewness

The skewness measures the symmetry of the data. Normal distribution has a skewness of 0, the right is positive, and the left is negative. The calculation formula of the skewness is as follows:

S = 1 n i = 1 n ( x i x ¯ ) 3 ( 1 n i = 1 n ( x i x ¯ ) 2 ) 3 2 (6)

2.4. Kurtosis

The kurtosis measures the flatness of the data distribution. The large data distribution at the tail has a large kurtosis value. The normal distribution has a kurtosis value of 3. Its formula is as follows:

K = 1 n i = 1 n ( x i x ¯ ) 4 ( 1 n i = 1 n ( x i x ¯ ) 2 ) 2 (7)

2.5. Standard Deviation

The standard deviation can reflect the degree of dispersion of the data set, and the calculation formula is as follows:

S = 1 n 1 i = 1 n ( x i x ¯ ) 2 (8)

3. Empirical Analysis of ARCH Model

This paper takes the data of China’s commercial housing sales from January 2006 to October 2018 as the sample time series, and conducts data analysis with the help of EViews7.2 [7] [8] . In order to study the ARCH model of commodity housing sales, the ARMA model is used to make an empirical analysis of China’s commodity housing market by determining the return rate sequence of commodity housing sales and the required sequence must be stable. The data comes from China statistical yearbook (2018).

3.1. Data Preprocessing

It can be seen from Figure 1 (commercial housing sales) that the sales volume of commercial housing presents an exponential trend on the whole. The annual December is much larger than the annual January, showing seasonal influence. And the unit root test is further used (Table 1).

According to the results in Table 1, the ADF (Augmented Dickey-Fuller test statistic) value is 0.358366. At the significance level of 1%, the sequence accepts the original false, and the P value (Probable value) is 0.9805. It is believed that the unit root exists in the sequence Xt, indicating that the commercial housing sales volume sequence is a non-stationary time series.

In order to study the sequence of commercial housing sales, relevant data are

Figure 1. Time series diagrams of Xt.

Table 1. Unit root test.

processed as follows: two first-order and one seasonal difference were performed on the data, respectively eliminating the exponential trend and seasonal influence of the sequence, so that the new yield sequence obtained does not have the correlation trend and periodic influence, and at this time, the yield sequence is a stationary sequence.

It can be seen from Figure 2 that the sequence Rt is stable, and the continuity and clustering of the sequence can be observed.

According to the results in Table 2, the ADF value is −8.295501. At the significance level of 1%, the sequence does not accept the null hypothesis, and the P value is 0. It is believed that there is no unit root in the sequence Rt, indicating that the rate of return at this time is a stable time series.

3.2. Normality Test

Calculate each statistic for the new sequence, as known in Table 3.

Table 3 lists the results of statistics of the return rate of commercial housing sales Rt, in which the skewness is −0.488687, less than 0, indicating the existence of left-deviation phenomenon of the return rate sequence. The kurtosis value was 15.75282, while the kurtosis value of the normal distribution was 3, which was much higher than its range value, indicating that the sequence Rt had a distribution of “sharp peak and fat tail”. The value of the J-B statistic is 954.2733,

Figure 2. Time series diagrams of Rt.

Table 2. Unit root test.

Table 3. Rt sequence statistics.

and its P value is 0, indicating that the sequence Rt does not obey the normal distribution.

3.3. ARCH Effect Test

In order to study the relationship within the return sequence Rt, ARMA (auto-regressive moving average) model is taken as the main model, ARMA(2,0) was finally determined as the subject model by analyzing the autocorrelation and partial autocorrelation graphs of the yield Rt series, as known in Table 4.

Its parameters are estimated as follows:

R t = 0. 6364 0 9 R t 1 0. 339977 R t 2 + a t (9)

The sequence autocorrelation LM (Lagrange Multiplier) test for the random perturbation term at of the equation Rt obtained a P value of 0.0034, indicating the existence of autocorrelation in the sequence. Next, the author performs the autocorrelation LM test with order 7 on the sequence, as known in Table 5.

Table 4. Parameter estimation of sequence Rt.

Table 5. LM test for the sequence Rt.

It can be seen from Table 5 that the p-value of the seventh-order statistic is 0.0017, less than 0.01, indicating that the random perturbation term at sequence of the ARMA model has an autocorrelation, and it has the high order ARCH effect, so the return rate sequence Rt can build the GARCH model.

3.4. Model Recognition

GARCH(1,1) model should be selected to fit the high-order ARCH effect in the yield sequence Rt, and other types of ARCH effect models should be fitted. Table 6 lists GARCH(1,1), GARCH-M(1,1), TARCH(1,1), EARCH(1,1), logarithmic likelihood function values, AIC values and SC values of the return sequence Rt.

According to Table 6, the rate of return sequence Rt is more suitable for EARCH(1,1) model under AIC (Akaike Information Criterion) and SC (Schwarz Criterion) criteria, and its value reaches the minimum, which is 19.68798 and 19.81525 respectively. Meanwhile, the value of the logarithmic likelihood function is the largest, and the value is −1352.471. Therefore, the fitting model EARCH(1,1) is finally determined.

It can be seen from Table 7 that all P values of the coefficients of EARCH(1,1) model are 0, indicating good fitting effect:

Mean value equation: R t = 0.266744 R t 1 0.312926 R t 2 + a t (10)

Variance equation: log ( h t ) = 9.414443 + 0.450156 log ( h t 1 ) 0.524599 | a t 1 | h t 1 0.775541 a t 1 h t 1 (11)

According to the variance equation, the parameters α 0 = 9.41443 and θ 1 = 0.450156 are both greater than 0, indicating that the fluctuation of income is positive and in line with the actual situation of commercial housing sales today. the asymmetric coefficient φ 1 = 0.775541 is less than 0, indicating that the fluctuation phenomenon of the rate of return on sales of commercial housing in China is asymmetric, when α 1 + φ 1 = 1 . 30014, it shows the impact when there is good news, when α 1 φ 1 = 0.250942 , it shows the impact of bad news.

3.5. Adaptability Test

The ARCH effect test is performed by the sequence in the variance equation to

Table 6. Model identification.

Table 7. Parameter estimation of EARCH model.

find out whether this effect has been eliminated. Figure 3 is an autocorrelation and partial autocorrelation test of the standard deviation residual sequence of the EARCH(1,1) model of the sales rate of commercial housing in China.

It can be seen from Figure 3 that the P value is 0.017 when proceeding to the fifth step, and the autocorrelation hypothesis that the significance level is 0.01 standard deviation residual sequence is rejected, indicating that the residual sequence has no ARCH effect. It may be reasonable to establish an EARCH(1,1) model. Next, ARCH effect test can be performed on the residual sequence.

It can be seen from Table 8 that the P value of LM statistic is 0.7131, greater than 0.01, which further indicates that the ARCH effect no longer exists in the residual sequence. Therefore, fitting EARCH(1,1) model is appropriate.

4. Conclusion

This paper mainly studies the characteristics of the return rate of commercial housing sales and finds that the return rate sequence has ARCH effect, and EARCH(1,1) model is still effective. From 2006 to 2018, China’s commercial housing sales have witnessed unprecedented rapid development. What’s more,

Figure 3. Correlation diagram of standard deviation residual sequence.

Table 8. ARCH effect test results of sequence.

the fluctuation of the return rate of commercial housing sales is positive, and the leverage effect coefficient of model is negative, indicating that the fluctuation of China’s commercial housing sales is asymmetric and has significant leverage effect. At the same time, by analyzing the variance equation of the model, it can be seen that the rate of return on sales of commercial housing has a more significant response to adverse news. In other words, adverse factors in the market have a greater impact on commercial housing market. With the increase of favorable news, the commercial housing market will eventually level off. Through the above empirical analysis, several suggestions can be made for the development of the commercial housing market. The government can regulate and control the commercial housing market; it can increase good news that is conducive to market development. Buyers can determine the purchase intention according to the price increase. The government should adopt policy measures in the commodity housing market to steadily increase sales and consolidate the development of the commercial housing market.


This work is supported by the National Natural Science Foundation of China (No. 11561056) and Natural Science Foundation of Qinghai (No. 2016-ZJ-914).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Shen, S.C. and Feng, C. (2019) Empirical Analysis of Commercial Housing Sales Based on EARCH(1,1) Model. Open Journal of Statistics, 9, 299-307.


  1. 1. Wei, J.J. and Zhang, C. (2014) Empirical Analysis on the Influencing Factors of Commercial Housing Sales in Nanning City. Economic and Social Development, 12, 16-19.

  2. 2. Shi, H.L. (2018) Empirical Analysis of Real Estate Price and Economy—Based on Sectional Data Analysis of Commodity Housing Sales and Gross Domestic Product in 30 Provinces. Era Economy and Trade, 7, 57-58.

  3. 3. Zheng, L. (2006) Research on Monetary Policy Transmission Effect in Real Estate Market Based on VAR Model. Hunan University Press, Changsha.

  4. 4. Tian, T. and Dong, W.J. (2009) Research on the Impact of Macroeconomic Changes on the Demand for Commercial Housing. Modern Economic Information, 9, 54-55.

  5. 5. Engle, R.F. (1982) Autoregressive Conditional Heteroscedaticity with Estimates of the Variance of U.K. Inflation. Econometrica, 50, 987-1007.

  6. 6. Nelson, D.B. (1991) Conditional Heteroskedasticity in Asset Returns: A New Approach. Econometrica, 59, 347-370.

  7. 7. Gao, T.M. (2009) Econometric Analysis Methods and Modeling—EViews Applications and Examples. 2nd Edition, Tsinghua University Press, Beijing.

  8. 8. Yi, D.H. (2008) Data Analysis and Application of Eviews. China Renmin University Press, Beijing.