Scientific Research

An Academic Publisher

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

Share and Cite:

*Open Journal of Statistics*,

**9**, 299-307. doi: 10.4236/ojs.2019.92021.

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 a_{t}-ARCH(q) [5] , the regression equation can be established:

${a}_{t}={\alpha}_{0}+{\alpha}_{1}{a}_{t-1}^{2}+\cdots +{\alpha}_{q}{a}_{t-q}^{2}$ (1)

The null hypothesis and alternative hypothesis to be tested are:

$\begin{array}{l}{H}_{0}:{\alpha}_{1}={\alpha}_{2}=\cdots ={\alpha}_{q}=0\\ {H}_{1}:\exists {\alpha}_{i}\ne 0\text{\hspace{0.17em}}\left(1\le i\le q\right)\end{array}$ (2)

Test statistic

$LM=n{R}^{2}~{\chi}^{2}\left(q\right)$ (3)

where, n is the number of sample data; R^{2} is the determinant of the auxiliary regression formula. When significance level
$\alpha $ and degree of freedom q are known, if
$LM>{\chi}_{\alpha}^{2}\left(q\right)$ , 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:

$\mathrm{log}\left({h}_{t}\right)={\alpha}_{0}+{\displaystyle \underset{j=1}{\overset{p}{\sum}}{\theta}_{j}}\mathrm{log}\left({h}_{t-j}\right)+{\displaystyle \underset{i=1}{\overset{q}{\sum}}\left[{\alpha}_{i}\left|\frac{{\epsilon}_{t-i}}{\sqrt{{h}_{t-i}}}\right|+{\phi}_{i}\frac{{\epsilon}_{t-i}}{\sqrt{{h}_{t-i}}}\right]}$ (4)

where the expression of ${h}_{t}$ is:

${h}_{t}={\alpha}_{0}+{\alpha}_{1}{\epsilon}_{t-1}^{2}+\cdots +{\alpha}_{q}{\epsilon}_{t-q}^{2}$ (5)

If $\phi \ne 0$ , the information function is asymmetric. When $\phi $ 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=\frac{\frac{1}{n}{\displaystyle \underset{i=1}{\overset{n}{\sum}}{\left({x}_{i}-\stackrel{\xaf}{x}\right)}^{\text{3}}}}{{\left(\frac{1}{n}{\displaystyle \underset{i=1}{\overset{n}{\sum}}{\left({x}_{i}-\stackrel{\xaf}{x}\right)}^{\text{2}}}\right)}^{\frac{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=\frac{\frac{1}{n}{\displaystyle \underset{i=1}{\overset{n}{\sum}}{\left({x}_{i}-\stackrel{\xaf}{x}\right)}^{\text{4}}}}{{\left(\frac{1}{n}{\displaystyle \underset{i=1}{\overset{n}{\sum}}{\left({x}_{i}-\stackrel{\xaf}{x}\right)}^{\text{2}}}\right)}^{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=\sqrt{\frac{1}{n-1}{\displaystyle \underset{i=1}{\overset{n}{\sum}}{\left({x}_{i}-\stackrel{\xaf}{x}\right)}^{\text{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 X_{t}, 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 X_{t}.

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 R_{t} 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 R_{t}, 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 R_{t}, 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 R_{t} 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 R_{t}.

Table 2. Unit root test.

Table 3. R_{t} sequence statistics.

and its P value is 0, indicating that the sequence R_{t} does not obey the normal distribution.

3.3. ARCH Effect Test

In order to study the relationship within the return sequence R_{t}, 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 R_{t} series, as known in Table 4.

Its parameters are estimated as follows:

${R}_{t}=-0.\text{6364}0\text{9}{R}_{t-1}-0.\text{339977}{R}_{t-2}+{a}_{t}$ (9)

The sequence autocorrelation LM (Lagrange Multiplier) test for the random perturbation term a_{t} of the equation R_{t} 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 R_{t}.

Table 5. LM test for the sequence R_{t}.

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 a_{t} sequence of the ARMA model has an autocorrelation, and it has the high order ARCH effect, so the return rate sequence R_{t} 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 R_{t}, 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 R_{t}.

According to Table 6, the rate of return sequence R_{t} 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: $\begin{array}{c}\mathrm{log}\left({h}_{t}\right)=9.414443+0.450156\mathrm{log}\left({h}_{t-1}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-0.524599\frac{\left|{a}_{t-1}\right|}{\sqrt{{h}_{t-1}}}-0.775541\frac{{a}_{t-1}}{\sqrt{{h}_{t-1}}}\end{array}$ (11)

According to the variance equation, the parameters ${\alpha}_{0}=9.41443$ and ${\theta}_{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 ${\phi}_{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 ${\alpha}_{1}+{\phi}_{1}=-1$ . 30014, it shows the impact when there is good news, when ${\alpha}_{1}-{\phi}_{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.

Funds

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.

[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] | 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] | Zheng, L. (2006) Research on Monetary Policy Transmission Effect in Real Estate Market Based on VAR Model. Hunan University Press, Changsha. |

[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] |
Engle, R.F. (1982) Autoregressive Conditional Heteroscedaticity with Estimates of the Variance of U.K. Inflation. Econometrica, 50, 987-1007.
https://doi.org/10.2307/1912773 |

[6] |
Nelson, D.B. (1991) Conditional Heteroskedasticity in Asset Returns: A New Approach. Econometrica, 59, 347-370. https://doi.org/10.2307/2938260 |

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

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

Copyright © 2020 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.