Abstract

House price is affected by households’ expectation of future house price trend and volatility, where the expected volatility of housing capital return, indicated by variance, is defined as the housing market risk. Theoretically, risk element cannot be directly inserted in the standard housing models because most of the models are built on the underlying assumption of certainty. Extending the life-cycle model to a two-asset expected utility case with uncertainty, we show house price is affected by housing market risk premium, which is a function of households’ risk-aversion coefficient, real housing wealth, and expected housing volatility. Empirical analysis relying on China’s 2001–2018 provincial housing panel data supports the theoretical innovation. Despite the empirical results show that the effect of housing market risk on house price is tiny, simulations suggest that the consideration of risk is quite helpful in analyzing and predicting the long-run house price equilibriums.

1. Introduction

Over the past 20 years since entering the 21st century, house prices are sharply raised in most of the developing and developed countries. This worldwide event gives people an impression that investment in housing market is relatively better on asset preservation and appreciation. However, house price does not increase always. For instance, Japan experienced housing market collapse in 1991, where the price in 2003 is only 1/4 of the price in 1990; as the indicator of Southeast Asia, Hong Kong’s house price reduced from 13000 to 5000 (HK dollar) during the Southeast Asia Financial Crisis; and in the recent 2007–2008 Global Financial Crisis, house prices in lots of countries are dramatically reduced. Housing market risk gradually comes into our view. Academically, many attentions focus on house price prediction; however, although classical housing economic models work very well in common periods, they always get lost in the depreciations.

The development on house price theories has a long history such as (i) the “mark-up” model ([1]) where house price is mainly tracked by construction cost, especially the land price, (ii) the reduced-form models ([2]) where demand and supply of housing services are described, (iii) the housing stock-flow model ([3, 4]) where the housing stock is accumulated by the flow of housing purchases, (iv) the hedonic model ([5]) where house price is determined by lots of factors, and (v) the housing life-cycle model ([6]) where house is clearly defined as both consumption and investment products. Despite the standard housing life-cycle model is popularly used in analyzing housing issue, there is no position for risk element in the theoretical framework. Lots of research studies concentrate on this issue, but they find that the fundamental problem is that the standard model is built on the underlying assumption of certainty. The assumption of certainty implies that house price has no risk (unexpected future volatility); thus, people are perfectly informed even for future prices. As the series of studies (e.g., [3, 4, 710]) suggested, the assumption of rational expectations should be rejected. They also claim that, under the assumption of certainty, the house price derivation of the standard model actually reflects the price of “nonrisk” houses. However, the price of realistic houses is risky. So, even the early studies, such as Poterba [3], recognised the importance of a risk premium in an asset pricing framework.

The issue of risk is deeply studied in financial area. Financial theory originated with the mean-variance portfolio model [11] and evolved into the capital asset pricing model (CAPM) [1214]. The CAPM shows that the return of any risky asset is equal to a riskless rate plus a risk premium for a particular risky asset. However, the empirical support for the static CAPM was weak [15, 16], which may reflect theoretical failings arising from many simplifying assumptions as well as difficulties in implementing valid tests of the model [17]. Consequently, the consumption-based capital asset pricing model (C-CAPM) was developed by Lucas [18] as a dynamic generalization of the CAPM also linking it to the real macroeconomy. In contrast to the CAPM, it can theoretically show how economic fundamentals affect returns and prices of risky assets. Nevertheless, housing was not incorporated in asset allocation models of the C-CAPM type until Grossman and Laroque [19]. Furthermore, Flavin and Nakagava [20] found that their housing model, which adds nondurable consumption in the utility function to extend Grossman and Laroque [19] beyond durable housing being the sole argument in utility, provides support for the C-CAPM. Piazzesi et al. [21] proposed a nonseparable utility model aggregating “shelter” and “consumption” via a constant elasticity of substitution index and found that the change in the expenditure share emerges as a second factor, complementing the single-factor C-CAPM, in their “housing C-CAPM” [21]. In an empirical context, Case et al. [22] developed a multifactor housing asset pricing model and found that house prices are determined not only by riskless assets and a theoretically derived risk premium but also by idiosyncratic risk, momentum, and a size effect, where Banz [23] suggested that returns on small firms were relatively high. Therefore, given on all the above studies, expected future volatility of housing capital return is commonly defined as housing (market) risk and empirically measured by variance of the housing capital return; however, their theoretical definitions of the risk premium are slightly different, but most of them estimate risk premium empirically based on standard C-CAPM definitions [24]. In addition, some studies track the determinants of housing market risk premium such as the uncertainties in interest rate, inflation, and income and then research on a specific aspect and define the decomposed risk (premium) separately such as the interest risk of housing, inflation risk of housing, and households’ income risk of housing [25, 26]. However, these separately defined risks or risk premiums are not the interest in this paper.

The theoretical framework about the continuous-type life-cycle model under uncertainty can be tracked back to Merton’s “lifetime portfolio selection under uncertainty: the continuous-time case” [27]. This framework can estimate the time-varying risks but is not typically for housing issues. Rely on this framework, some further studies (e.g., [2831]) measure and explore housing risks. However, most of these do not take to much experience from housing models, especially Dougherty and Van Order’s housing life-cycle model. It means that they still treat house as a pure financial asset and ignore its consumption property (demand). However, even if houses cannot be traded, the uncertain future prices will also affect households’ consumption demand of housing services. Recently, although the most relevant study [32] derives the results based on the discrete type of households’ lifetime utility and discrete type of lifetime budget constraints, whether the continuous-type modelling has effect or not still needs to be checked. In addition, massive of studies focus on the portfolio decision between housing and financial assets, or the interaction between housing and financial markets [33], which are not the interest in this paper.

Based on the above and some other relevant studies, a couple of points can be summarized: (1) according to the comparison between the CAPM and C-CAPM, housing market risk should be automatically derived from the model rather than subjectively defined. (2) The necessary condition of getting a risk element in the theoretical derivation is the allowance and the consideration of uncertainty. (3) Pure financial models especially the CAPM-type models are not appropriate in describing and analyzing house issues because houses not only have investment property/demand but also have consumption property/demand. In light of these, it is necessary to rebuild the standard housing life-cycle model, where housing purchase is risky and households’ lifetime utility is uncertain, and then to test whether a theoretically coherent measure of the risk premium can be derived, parallel to that used in financial economics.

Therefore, this paper is organized as follows. Section 2 is the theoretical section including the mathematical background of the standard housing life-cycle model, the results when the standard model and CAPM are combined, and the expected utility housing life-cycle model. Section 3 is the empirical section in which the effect of housing market risk will be analyzed based on China’s housing data. Simulations will be given and discussed in Section 4. Conclusions are drawn in Section 5. References and appendixes are given at the end.

2. Theory

In this section, the standard housing life-cycle model will be firstly introduced. Based on its assumptions, we will then discuss why the standard user cost equation is biased in estimating house price. Given by the standard user cost equation, the more specific house price equation can be derived through a simple link between the certain and uncertain housing capital returns corresponding with the theory of uncertainty. Although this combined complement method is simple, the consistency between models is hard to be proved. After that, a fully developed and consistent theoretical model, the two-asset expected utility housing life-cycle model, is shown. The general and special derivations of house price as well as the use cost of capital (UCC) will also be shown and discussed.

2.1. Standard Housing Life-Cycle Model

The standard intertemporal housing life-cycle model extended to housing assumes two goods, housing services and a composite consumption good (C). The model has a long history, including as noted above [6]. If, for simplicity, the flow of housing services is proportional to the demand for the housing stock (H), and given an assumed constant real discount rate (r), lifetime utility is described in an infinite horizon continuous time by the following equation:

μ[H, C] denotes the instantaneous utility of the representative household. (1) is maximized with respect to the budget constraint (2) and technical constraints (3) and (4), which describe changes to real asset stocks (housing and financial, respectively) over time. In addition, for the purpose of simplicity, the issues of tax and housing transaction cost are removed [3].where H(t) = stock demand for housing,  = real purchase price of dwellings, X(t) = new purchases of dwellings, S(t) = real savings net of real new loans, RY(t) = real household income, i(t) = market interest rate, A(t) = real net nonhousing assets, δ(t) = depreciation rate on housing, π(t) = general inflation rate, and ≡ dx(t)/dt denotes the time derivative for any variable x(t).

From the first-order conditions (see Appendix A), the marginal rate of substitution (MRS) between housing and the consumption good, μH/μC, is given by (5), where rH(t) =  is the real price growth rate:

The term in [⋯] in the right-hand side is the widely-used standard definition of the real housing user cost of capital (UCC). The left-hand side is the marginal rate of substitution between housing and consumption good, and this MRS should equal the imputed rental price of housing services (R) which means that households’ utility is unchanged if their housing services cut down unit, but they obtain some nonhousing goods as compensations. Thus, (5) can be rewritten as R(t) = UCC(t), and this is the well-known arbitrage condition. Equation (5) can also be written as follows:

The left-hand side of (6) is the real net housing return, which is the sum of housing rental yield () and housing capital return (rH) minus the housing depreciation rate, and the right-hand side is the real interest rate on liquid assets, or the risk-free rate. Given the assumption of certainty, the future house price is observable, so the housing capital return is certain. Equation (6) implies that there is no difference in the final returns whether people choose buy-to-let investment or bank saving. As we know, housing investment should be risky because its capital return is practically uncertain. So, the problem is that housing in the standard model is assumed to be risk-free. The derived house price equation (not shown) actually refers to the price of the “risk-free housing asset” rather than the price of the practical risky houses.

2.2. Combined Complement Method

Now, we assume that there are two kinds of housing assets in the market: one is risky and another is risk-free. Risk-free houses can be simply regarded as the houses play on option market where owners know the future prices of their houses. The real capital return of the risk-free houses is certain and marked as rH,free, and the real capital return of the risky houses is uncertain and marked as rH,risky.

In order to guarantee two kinds of houses have no absolute advantage compared with another one in the market, the expected utilities of the final wealth of the two kinds of houses should be indifferent, given by (7), where HW is the real housing wealth (or real housing value):

The utility function, μ[…], is concave, and the general result implies that E[rH,risky] > rH,free. In order to solve the issue and calculate a more specific result, two more conditions are needed: the distribution function of the uncertain element (rH,risky) and a specific utility function. The real capital return is assumed to be distributed normally as (8), and the CARA (constant absolute risk aversion) utility function (9) is selected to be the specific utility function, in which γ is the Pratt(1964)–Arrow(1971) parameter of CARA [21, 34]:

The derivation of equation (10) is shown in Appendix B, and the capital return from the risk-free house (if such a thing existed) should be the same as the expected capital return on the risky house minus a housing risk premium. The housing risk premium is determined by the CARA parameter, real housing wealth, and the expected housing volatility. Substituting (10) into (5), the cost of capital for risky houses is described by (11):

Although oversimplified, (11) provides the first insights and shows that it is not sufficient just to include a measure of the price variance in empirical studies because risk also depends on other factors, notably the market value of the housing stock. Despite the combined complement method is simple and easily understandable, it is very difficult to illustrate the consistency between the standard housing life-cycle model and the theory of uncertainty.

2.3. Expected Utility Housing Life-Cycle Model

In this section, housing acts as both a consumption good and a risky asset. Liquidity of a risky asset including housing will affect the asset price [35], but for the purpose of simplicity, we have already removed the costs of housing in the theoretical framework in Section 2.1, and therefore, housing now is a liquid risky asset. The liquidity issue will be studied in our further work. In contrast to the standard model, the extension is to treat housing as risky in terms of the capital return. Extending the standard life-cycle housing model to the uncertain case, compared to equation (1), households maximize their expected lifetime utility (12). Although the constraints mathematically are still (2)–(4), some related mathematical derivations are different. Typically, future house price () as well as expectation on housing capital return (rH,risky(t)) in time period (t) for the following period (t + dt) is uncertain now (see Appendix C).

From the first-order conditions (see Appendix C), the MRS between housing and the consumption good (E[μH]/E[μC]) is given bywhere E[μH(t)] = ∂E[μ[H(t), C(t)]]/∂H(t), the expected marginal utility of housing, and E[μC(t)] = ∂E[μ[H(t), C(t)]]/∂C(t), the expected marginal utility of consumption.

Equation (13) shows a generalized MRS and UCC under uncertainty. When housing capital return is certain, expected marginal utility of consumption in both numerator and denominator will be canceled, and then the right-hand side of (13) would reduce to the widely-used standard definition (5). UCC can be rewritten as (14) and (15), where the newly derived term in UCC is τ. τ is quite the same as the risk premium derivation in the C-CAPM, but equation (15) is specifically for housing. Similar to the risk premium derivation in the C-CAPM, τ is the finally generalized theoretical form of the housing market risk premium. Equation (14) points out that user cost of capital is negatively affected by expected housing capital return but positively affected by the housing market risk premium:

In some financial literatures, consumption data are used to check whether risk premium derivation is significant or not. In contrast, consumption is endogenous in our expected utility model. This makes the risk premium derivation much abstract. To be able to derive intuitive results of the housing market risk premium (τ), specific utility function and distribution of housing capital return are needed. Practically, we should identify the specific households’ utility and the distribution of the expected housing capital return of a specific housing market and then calculate the specific form of τ. However, evidences from previous studies show that the different distribution functions and the utility functions will cause different special solutions. More specifically, both the kurtosis (heavy-tail issue) and skewness (asymmetric issue) have impacts on the risk premium derivation and the distribution of housing capital return compared with normal is practically heavy-tailed and asymmetric [36, 37]. Furthermore, different utility functions have different risk-averse coefficients and different curvatures, which influence the effect of risk on price [38]. For different housing markets, the distribution of households’ expectation of housing capital return and utility function are unique. In order to capture the main line and avoid the complicated works, here, for simplicity, the distribution is assumed to be normal. Similarly, the standard CARA and CRRA (constant relative risk aversion) utilities are employed [38] which are given by equations (16) and (17). The main aim is to understand the calculation process in this paper. The effects of the specific distributions and utilities will be studied in the future.

Case I. Housing capital return is distributed normal as equation (8) and utility function is CARA utility denoted by equation (16).

Case II. Housing capital return is distributed normal as equation (8) and utility function is CRRA utility denoted by equation (17).The specific analytical effect of the housing market risk (τ) under CARA and CRRA utilities are derived in Appendixes D and E, respectively, given by equations (18) and (19). Equation (18) is completely the same as (11) because of the different scales of the risk aversion (e.g., γ = 2ϕ). By definitions, ϕ is the CARA coefficient and θ/C is the CRRA coefficient. In finance, the variance of the expected asset return is usually defined as (total) risk. Similarly, the variance term in (18) and (19) is housing market risk. In light of this, the housing market risk premium is a function of the specific risk-averse coefficient, real house price, housing stock, and housing market risk. The existence of real housing wealth is unsurprised as some microstudies found that the effect of housing risk varies over different income groups [39]. The construction of τ is also supported by the empirical founding provided by Case et al. [22]. Equation (19) is analogous to the ratio of the real relative quantities of housing services and the consumption numeraire derived by Piazzesi et al. [21] under their nonseparable CRRA utility:Compared with other studies, our housing market risk premium (τ) has more than two elements including housing market risk. In other words, it cannot simply add a variance term in regression as before because the coefficient should be variable-varying with households’ risk-averse attitude and real volume of housing wealth. More details will be discussed in the empirical section.
The user cost of capital is positively affected by τ and thus positively affected by the level of households’ risk-averse, the real housing wealth, and the housing market risk. Because real house price is negatively affected by UCC, the effects of these four elements particularly house price itself on house price are negative. House price has an in-built negative feedback adjustment through housing market risk premium. When housing market is overheated, even given by the same housing market risk, house price works as an amplifier to extend the housing market risk premium and then to cool down itself through a housing collapse. As risk is always positive but return is not, this very special in-built mechanism can help to explain (i) why house price steady grows in the boom but sharply falls in the bust and (ii) why standard housing models work very well in common periods but are terribly biased in forecasting long-run house price level.

3. Empirical Section

In this section, the main purpose is to check whether housing risk term is significant or not, and if it is significant, then estimate its marginal effect on house price. Technically, China’s housing market is selected to be the research object for three reasons: (i) China’s statistic office provides both monthly and annual house price indexes. The volatility (annual variance) of the later could be captured from the former. (ii) Provincial data is available and therefore the relation can be analyzed by panel data method, which is relatively robustness than the single-equation time-series analysis. (iii) China’s house price index is the price per unit square meter which will then be less likely affected by the housing market regional heterogeneity compared to the indexes which measure the price per house.

3.1. Data and Tests

The final data set is a panel including 31 China’s provinces except Taiwan, Hong Kong, and Macao during 18 years from 2001 to 2018. 13 variables are contained in this annual data set, and their details are listed in Table 1. Except interest rate, all variables are obtained from China Statistical Yearbook; interest rate is proxied by the base interest rate on RMB loans which is provided by China Financial Statistical Yearbook. Despite the means of the annual house price index (HP) and the estimated annual house price index (HPm) are not completely the same, their correlation is 0.9992. So, the estimated housing volatility has the ability to capture the annual fluctuation. As shown in Table 1, logarithms of the real house price and real disposable income are formed which will be applied in the empirical analysis. In addition, because the CRRA parameter (θ) is unknown and the empirical effect of term (τ) is in the same position, we assume that θ = 1.

Stationary test is necessary in deciding which kind of panel models will be selected in this specific analysis. Four variables are mainly focused; they are real house price (lg), real income (lry), interest rate (i), and housing market risk (tau). Because of the panel data set, the stationarity of these four variables are checked by LLC, IPS, ADF-Fisher, and Hadri LM stationary tests, shown in Table 2. Because LLC and Hadri LM tests require a balanced panel, Tibet is removed in tau. Notice that the null hypotheses of the former three tests are that the variable is unit-root, but the null of the last one is that the variable is stationary. After comprehensive consideration, results support that real house price (lg) and real income (lry) series contain unit-roots in all panels, but interest rate (i) and housing risk (tau) are stationary in all panels.

3.2. Empirical Analysis

Because (1) the imposed housing rent (R) is unavailable and (2) the logarithm of user cost of capital (UCC) cannot be simply decomposed in a linear form of its components, the theoretical equation and the empirical equation are slightly different. Lots of studies (e.g., [1]) describe how to derive a simple empirical equation (21) based on theoretical formula (20) which is the logarithm form of (5). Given on our slightly different UCC derivation (14), the housing risk term (tau) could be added into the formula in a similar pattern as interest rate and housing capital return, shown by equation (22). In general, a vector X is used to show the explanatory variables of imposed housing rent. As the data are available, personal disposable income (lry) is selected to capture the change in ln(R). In addition, keep in line with previous studies (e.g., [1]), households’ expectation on future housing capital returns (rH in the standard model and in the expected utility model) and inflation rate (π) are jointly approximated by last nominal house price growth rate ().

Because variables are nonstationary, equation (22) should be analyzed by error correction-type models. Because our data set is a panel, there are three possible cointegration relations, given by (23)–(25). Equation (23) implies that the short-run and long-run relations are the same in all regions; equation (24) implies that the short-run relations and error-correction speed vary across regions, but they will converge in the same long-run relation; equation (25) implies that, for all regions, each has its own error-correction approach. In addition, individual effect (ui) should also be considered in this panel data analysis.

The nonstationary heterogeneous panel model could estimate (23)–(25) because it has three analogical specifications: pooled mean-group (PMG) regression, mean-group (MG) regression, and dynamic fixed-effect (DFE) regression. The only difference is that the former two regressions estimate the average coefficient rather than specific coefficients for all regions. Notice that the average is enough for our analysis and the specific regional patterns are not our interest. Regression results are shown in Table 3.

The cointegration patterns are posterior determined by the Hausman test. Based on the tests among MG&PMG (Prob > chi2 = 0.15), MG&DFE (Prob > chi2 = 1.00), and PMG&DFE (Prob > chi2 = 1.00), the dynamic fixed-effect (DFE) estimator is inconsistent compared with the other two estimators, so it is relatively more efficient. The results of the Hausman tests suggest that the spatial heterogeneity is not significant. In light of this, all regions in China follow the same pattern in all of the short-run, the long-run, and the convergence process.

Despite the coefficients of expected housing capital return (a3) and housing market risk (a4) of the DFE regression are not significant, their values are just little bit higher than 10% significant level in both the DFE and PMG regressions. This case was also found in other papers, and there are two possible reasons. First, the dependent variable, Δlg, is quite closed to housing capital return, and then the expected housing capital return (rht−1) is not significant in ECM structure. Second, in order to empirically exam the form of UCC, the delogarithmic linear decompositions of ln(UCC) leads to this insignificant case.

The annual housing depreciation rate is commonly assumed to be 0.01 in previous papers (e.g., [1]). Given by the estimated coefficients of the DFE in Table 3, the ratio between the two estimated coefficients a3/a2 = 0.166/3.714 is 0.045 and a4/a2 = −1.595/3.714 is 0.429. Based on equation (20), rounding all these three numbers, the standard UCC is organized as equation (26). Similarly, with the consideration of housing market risk premium, the UCC (UCCrisky) is organized as equation (27). We have examined that the DFE regression is more appropriate in describing the China provincial housing data; the empirical relation with logarithmic UCC is given by equation (28):

As shown in Table 4, the consideration of housing market risk makes the significance of the short-run ln(UCC) parameter (β2) much higher and the goodness-of-fit also increases. Thus, the developed model is more appropriate to be employed in the housing issue analysis. Some once argued that the effect of housing market risk on house price is tiny. Indeed, the direct effect of tau is quite small, but we believe that this direct effect will gradually work through the price-return positive circle and completely change the long-run equilibrium of house price. However, because the volatility in price now is self-adjusted, common foresting methods do not work. Therefore, in the next section, simulations are conducted.

4. Simulation Analysis

We now proceed to analysis and interpretation of the role that the housing market risk premium plays in dynamic stochastic simulations of the expected utility model. To summarize, simulation framework consists of two equations: the long-run house price equation and the UCC equation. The long-run house price equation comes from (28) with the empirical estimations given in Table 4. The UCC equations (26) and (27) are separately applied, where the former is the comparison. τ is calculated as its CRRA form (19).

Real income is volatile and is the only source of exogenous shocks to house price. Based on RY’s statistical mean and standard deviation and simply assuming that its noises are distributed normal, the noises are generated. All of the initial values and growth rates of the variables and parameters also rely on the means given in Table 5. All details are summarized in Appendix F (Table 5).

For comparison, we conduct three curves. The first is based on (28) and (26), and the real income is not volatile in order to show the long-run tendency of house price ( series) implied in the standard housing life-cycle model. The second is still based on (28) and (26), but the noises of real income are considered. Then, the generated house price ( series) will also volatile. The third is based on (28) and (27), and the noises are still considered. The difference between and will reflect the effect of housing market risk on long-run house price equilibrium. In τ, the variance of households’ expectation is captured by the variance of last four years’ housing capital returns.

Notice that calibration simulation (such as the first curve) of the developed expected utility model is meaningless because there is no stochastic element in the calibration simulation and thus the generated τ becomes to zero. The key to the issue is that the housing market risk relies on the changes in house price inflation from both the exogenous shocks in real income and its self-adjustment. This is the reason why the popularly used CGE-type and DSGE-type models cannot reach to our target.

As an initial illustration of the properties, the four frames of Figures 1 and 2 separately plot one realization for the key variables—the variance of housing capital return, housing market risk premium, real house price, and user cost of capital—over 80 years (periods). The differences between these two figures are due to variations in the income series, which is generated by different random seeds (white noise series). Figures generally show that even if the volatility of the housing capital return is stable, housing market risk premium still has an upward trend due to the real house price inflation. The UCC with the consideration of τ will also has an upward trend, thus limiting the rise in price. In the short-run, the two dynamic house price series have similar dynamic changes because of the same exogenous real income fluctuation. This is the reason why the standard model still works very well in practical data.

The results from a single replication could occur by chance; therefore, repeat the results and average over 100 replications, shown in Figure 3. The results clearly show that because the noises are strongly eliminated, the difference in the long-run trends between standard UCC and the UCC with τ is now significant. First, the average of the real house price fits to the baseline very well, which implies that the consideration of noises is not important to the standard model in forecasting the house price long-run equilibrium. However, when risk is included, simulated house price series is significantly lower than the baseline meaning that the consideration of risk is quite necessary in analyzing house price long-run equilibriums. All these simulations, especially Figure 2, supports our view: house price has a special pattern converges to its long-run equilibrium, where it steadily grows in the booms but sharply reduces in the busts, although the undesirable market collapse is extremely painful.

5. Conclusion

This paper provides a number of sections theoretically and empirically focusing on the effect of the housing market risk on house price. A commonly held view is that house prices rarely fall and, therefore, housing is a safe investment. Historical experience appears to bear this out. Lots of studies try to interpret the effect of housing risk on price related to financial theories. However, unlike financial products, house has both consumption and investment properties (demands), and it is very hard to illustrate the consistency between models when financial theories are embedded into economic theory. Risk element cannot be directly introduced into the standard housing model, especially the housing life-cycle model because of the underlying assumption of certainty. To be more realistic, expected housing capital return is redefined to be uncertain and then households maximize the expectation of their uncertain future lifetime utility.

Theoretically, starting at the standard housing life-cycle model, the combined complement method is firstly given in order to show the effect of housing risk on price although the combination between economic and financial theories is still a problem. A fully developed expected utility housing life-cycle model is given, which is purely economic-based, and one of its specific derivations is identical to the equation derived from the combined complement method. The results imply that the effect of the housing market risk—the risk-aversion parameter times real housing value and times housing volatility—is positive to UCC and is negative to house price. Given on the specific CRRA housing market risk derivation and relying on the China housing data, empirical results support our definition of the housing market risk premium.

Based on the empirical results, simulations are given in order to check whether housing market risk could have significant impact on house price long-run equilibriums. The results clearly show that, with the inclusion of housing market risk, long-run house price is no long terribly high compared with the implication derived from the standard housing life-cycle model. Evidences support the view that house price steadily grows in the common periods but backward adjusted through housing market catastrophe.

Finally, for the purpose of keeping house price more stable and reducing the probability of housing market collapse, we suggest that (i) governments can increase the feedback effect of risk on price through policies and systems such as housing security; (ii) governments can also promote residences’ risk perception through news and media; (iii) governments can smooth price cycles by the contracyclical regulation.

Appendix

A. The standard housing life-cycle model

The lifetime utility function and constraints:

The households’ lifetime housing decision can be found from maximizing the following Hamiltonian function:

The first-order conditions are as follows:

From (A.4),

From (A.3),where ġ(t) = .

Substituting (A.5) and (A.6) into (A.7),

Then,where rH(t) = .

B. The equivalence between risky and riskless houses

Conditions:

Rewrite the generalized equivalent under the CARA utility:

Then, include (B.2):where , the standard normal probability density function for z).

Then,

Therefore,

C. The two-asset expected utility housing life-cycle model

The lifetime utility function and constraints:

The households’ lifetime housing decision can be found by maximizing the expected Hamiltonian function:

The first-order conditions are as follows:

(Notice that all variables are certain at time period (t) and uncertain in the future (t + dt, t + 2dt, t + 3dt, …). So, certain elements can be removed to the outside of expectations.)

From (C.7),

From (C.6),where ġ(t) =  = .

Substituting (C.8) and (C.9) into (C 10), we obtainwhere rH,risky(t) = .

Then,

D. The two-asset expected utility housing life-cycle model under CARA utility conditions:

Equations (C.2)–(C.4) can be organized as one equation:or

This equation implies that if people do not want to change the quantities of their housing and nonhousing assets, all of their labour income, housing net yield, and nonhousing financial yield will be consumed.

Also, because housing net yield is uncertain due to the uncertain house price inflation rate, household expected consumption is given by

Remember that, to derive specific solution of equation (C.12), we need to solve E[μC(t)] and E[rH,risky(t)μC(t)] based on equation (D.5).Step 1: the derivation of E[μC(t)]:where , the standard normal probability density function for z.Method of deriving the final equation is the same as Appendix B.Therefore,Step 2: the derivation of E[rH,risky(t)μC(t)]:Therefore,Step 3: calculate the UCC as well as tau (τ) under CARA.Combine (D.9) and (D.10):Substitute (D.11) into (D.12):So,

E. The two-asset expected utility housing life-cycle model under CRRA utility conditions:

As is the standard normal probability density function for z, there are three related useful mathematical theorems:(1)Constant: (2)Constant × mean: (3)Constant × variance:

The marginal utility of consumption under CRRA is

Based on all the above conditions and mathematical theorems, solve E[μC(t)] and E[rH,risky(t)μC(t)] in order to derive the specific derivation of equation (D.12).Step 1: the derivation of E[μC(t)]:Take a second-order Taylor expression of μC(t), where C(t) = Ce(t):Then,Step 2: the derivation of E[rH,risky(t)μC(t)]:Take a second-order Taylor expression of rH,risky(t)μC(t), where C(t) = Ce(t) and :Then, Step 3: calculate the UCC as well as tau (τ) under CRRACombine the results from Steps 1 and 2:orIf expectations are ignored, {−θCe(t)θ−1}/{E[C(t)θ]} = θ/C(t); it will be the Pratt–Arrow coefficient for CRRA utility under certainty ().

Under this approximation, substituting (E.9) into (D.12), we obtain

Notice that the coefficient of the uncertain utility ({−θCe(t)θ−1}/{E[C(t)θ]} is still solvable, but because (i) national risk attitude is relatively stable, (ii) the difference with θ/C(t) is tiny, and (iii) the similar pattern with τCARA derivation, the approximation is more meaningful.

F. Key values used in the simulations

Further details and the complete replication of our results are available via our MATLAB simulation codes, which can be provided.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

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

Acknowledgments

This research was supported by Shandong Social Science Planning Research Project (China) (grant no. 19DTJJ02).