Next Article in Journal
Some Results on Majorization of Matrices
Next Article in Special Issue
Some Notes for Two Generalized Trigonometric Families of Distributions
Previous Article in Journal
Existence of Mild Solutions for Hilfer Fractional Evolution Equations with Almost Sectorial Operators
Previous Article in Special Issue
An Effective Approximation Algorithm for Second-Order Singular Functional Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Financial Market Factors on House Prices: An Expected Utility Three-Asset Approach

1
Department of Economics, Shandong University of Finance and Economics, Jinan 250014, China
2
Center of Excellence in Econometrics, Faculty of Economics, Chiang Mai University, Chiang Mai 50200, Thailand
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(4), 145; https://doi.org/10.3390/axioms11040145
Submission received: 31 January 2022 / Revised: 13 March 2022 / Accepted: 15 March 2022 / Published: 22 March 2022
(This article belongs to the Collection Mathematical Analysis and Applications)

Abstract

:
This study aimed to theoretically identify the impact factors of the financial market on house prices. Developed upon the two-asset model and with the consideration of risky financial assets, our three-asset model reveals a new derivation of house prices. Compared with the two-asset model, the newly emerged term is similar to the Sharpe β ; therefore, it is a risk premium term. Based on China’s 2001–2018 panel data, theoretical derivations are examined. However, the short-term effect of this risk term on house prices is practically small. Given the nonlinear pattern, the long-term effect of the risk term is checked by repeated stochastic simulation. The results imply the following: (i) real house prices are nonlinearly affected by three financial market factors, namely, the expected financial market return, financial market volatility, and the correlation between housing and financial markets; (ii) the correlation determines the signs and the significance of the effects of the other two factors; and (iii) the naturally changed correlation causes periodic house price fluctuations. Therefore, to stabilize real house prices, it is recommended that the government control the money flow between the two markets.

1. Introduction

House prices have been sharply inflated internationally in recent decades. Most people believe that the financial market plays an important role in housing market appreciation. It seems that house prices always increase in boom periods and decrease during financial crises. However, related empirical works have not supported this view. Muellbauer and Murphy [2] pointed out that house price changes are related to financial market liberalization in the UK. Green [3] showed that the impact of stock prices on house prices is high in Northern California but low in Southern California in the US. Lee et al. [4] focused on the relationship between Australian housing and stock markets and found that stock market prices were affected by house prices before the recent GFC but affected house price after the crisis. Kakes and Van Den End [5] found that house prices are significantly affected by stock market prices in the Netherlands. Takala and Pere [6] found that house prices are Granger-caused by stock prices in Finland. Ibrahim [7] applied the same method but found that there was no significant Granger-causality between housing and financial markets in Thailand. All of these related empirical works have suggested the following: (i) the impact factors of the financial market on house prices are unclear; (ii) their influence mechanisms are unclear; and (iii) the effects of the financial market factors on house prices should be much more complicated than we believed. Therefore, related theoretical work on the effects of financial market factors on house prices is necessary.
Housing is different from most financial assets since houses not only have investment demand but also have consumption demand, usually known as the dual property [8]. Houses provide housing services at every moment for people living in the house. Due to the consumption demand, house prices are affected by several economic factors, such as household income, family wealth, and construction costs [8]. These impacts are usually theoretically derived by consumption type models, especially the housing life-cycle model [9]. In contrast, housing is a type of asset that can be selected and invested in by investors. Based on the determinants of housing investment demand, the prices of houses are affected by many financial factors, such as (mortgage) interest rates, expected financial market return, the variances (or the standard deviations) of both housing and financial market returns, and the covariance (or the correlation) between the returns of the two markets [10,11,12]. These factors are derived from financial theories, for example, the Consumption-based Capital Asset Pricing Model (C-CAPM) and Housing CAPM [10,13,14,15,16]. However, although some studies [17] have nested the economic and financial factors in their empirical analyses and found that both of these factors are significant, it remains controversial whether it is appropriate to assume housing to be a type of pure financial assets and then to analyse house prices using financial theories [18].
In economics, the housing life-cycle model [9] is popularly used in analysing housing issues since houses are clearly defined as both consumption and investment products in the model. Despite the standard housing life-cycle model allowing for the consideration of the dual property of housing, as a result of the underlying assumption of certainty, it still has no position in risk issues. The expected utility two-asset housing life-cycle model (hereinafter, two-asset model) that [1] developed had a theoretical framework based on the standard model, in which (i) housing is assumed to be risky, and thus future house prices are assumed to be uncertain, like a normal distribution; and (ii) households maximize their expected lifetime utility. The results showed that the volatility of housing capital return has significant and non-negligible negative impacts on real house prices, particularly in long-run forecasting. The missing consideration of the effect of the volatility will cause over-estimation. However, the two-asset model remains incomplete since to assume houses are the only risky asset in the market is not sufficient. The effects of financial market factors on housing variables are also worth being analysed. In addition, the time-varying [19] and regional differed [20] distributions of house prices as well as the volatilises are also possibly caused by the factors from financial market or related government policies. These issues will be discussed in the future studies.
In light of these discussions, there have been three main gaps recently in this area. First, it is necessary to have a theoretical base in which, at least, risky housing and risky financial assets are jointly considered. Second, given the consumption-based model, households’ choices of risky assets, including housing and risky financial assets, should be derived and compared with those implied in financial theories. Further, the reasons why the results are consistent or different should be discussed. Third, the theoretical derivation should reflect the pattern of how financial market factors affect house prices, which would be particularly helpful for related empirical works.
Therefore, this article is organized as follows. Section 2 is the theoretical section, which provides the details of the three-asset model, including related theoretical derivations. Section 3 is the empirical section that practically estimates the short-term effect of housing market risk premiums on house prices. In Section 4, based on reasonable assumptions, the dynamic impact of the cyclically fluctuated correlation with real house prices is simulated. Section 5 provides the conclusion and some political suggestions. References and appendices are provided at the end.

2. Theory

In this section, the expected utility three-asset housing life-cycle model (hereinafter, three-asset model) is introduced. In related work [1], housing was assumed to be risky, and the housing market risk premium appeared in the housing user cost of capital (UCC), as well as the house price derivation. However, when risk issues are considered, the assumption of merely unique risk housing assets is insufficient. At least, all other financial assets should be categorized into two groups: risk-free and risky. Therefore, the theoretical structure and some details of the three-asset model are first provided, and then the general solution of the model is derived. Under two specific utilities, the Constant Absolute Risk Aversion (CARA) and the Constant Relative Risk Aversion (CRRA) utilities, special solutions are derived. Because of the joint consideration of two risk assets, the potential households’ portfolio decisions are derived and discussed.

2.1. The Expected Utility Three-Asset Housing Life-Cycle Model

Based on the standard inter-temporal housing life-cycle model [9] and the expected utility two-asset housing life-cycle model [1], household temporal utility is determined by non-housing consumption (C) and the used volume of housing services, which is proportional to the housing stock (H) as μ [ H ( t ) , C ( t ) ] . For “now,” household’s utility μ [ H ( 0 ) , C ( 0 ) ] is certain, but all future t > 0 utilities are uncertain. Inter-temporally, given an assumed constant real discount rate (r), the present value of the lifetime utility is described for an infinite horizon discrete time (t) as 0 1 ( 1 + r ) t μ [ H ( t ) , C ( t ) ] or continuous time (t) as 0 e r t μ [ H ( t ) , C ( t ) ] d t . A rational household now will make decisions to maximize the expectation of the continuous type present value of its uncertain lifetime utility (U), given by Equation (1).
U = E 0 e r t μ [ H ( t ) , C ( t ) ] d t
Households’ expected lifetime utility is maximized with respect to budget constraint (2) and technical constraints (3)–(5), which describe changes in the three real asset stocks (risky housing, risk-free assets, and risky financial assets, respectively) over time. Generally, based on the expected returns of these three assets, households will make an investment decision on the three assets for the purpose of maximizing their expected lifetime utility. In addition, for the purpose of simplicity, the issues of taxes [21,22], transaction costs [23], liquidity [24], and restrictions on lending [25] are removed.
p ( t ) C ( t ) + p H ( t ) X ( t ) + p ( t ) S f ( t ) + p A ( t ) S A ( t ) = Y ( t ) + i ( t ) p ( t ) A f ( t ) + i d ( t ) p A ( t ) A A ( t )
H ˙ ( t ) = X ( t ) δ ( t ) H ( t )
A f ˙ ( t ) = S f ( t ) π ( t ) A f ( t )
A A ˙ ( t ) = S A ( t ) π ( t ) A A ( t )
where
p ( t ) = price level/Consumer Price Index(CPI);
C ( t ) = quantity of non-housing consumption;
p H ( t ) = nominal purchase price of dwellings;
X ( t ) = quantity of new purchases of dwellings;
S f ( t ) = quantity of savings net of real new loans volume of risk-free assets;
p A ( t ) = nominal purchase price of risky financial assets;
S A ( t ) = quantity of savings net of real new loans volume of risky assets;
Y ( t ) = nominal disposable income;
i ( t ) = market risk-free/interest rate;
A f ( t ) = quantity of net non-housing risk-free assets;
i d ( t ) = dividend yield of the financial risky asset;
A A ( t ) = quantity of net non-housing risky assets;
δ ( t ) = depreciation rate on housing;
π ( t ) = general inflation rate; and
x ˙ d x ( t ) / d t , denotes the time derivative for any variable x(t).
Budget constraint (2) implies the equivalence between nominal terms. By dividing the price level (p) for both sides, Equation (6) indicates the equivalence between real terms. It shows that the flow of net expenditures on non-housing consumption, housing, risk-free savings, and savings on risky financial assets is equal to the net earnings from labour, risk-free investments, and risky financial investments. Now, housing acts as both consumption and investment since housing stock (service) directly determines household’ instantaneous utility and, at the same time, inflates in a similar pattern to financial assets. Additionally, some have asked why there is no housing rent/yield in the budget constraint. It is because, housing rent, in principle paid by all people, is equal to the housing rental income obtained by all households, so they are cancelled on both sides of Equations (2) and (6).
C ( t ) + g H ( t ) X ( t ) + S f ( t ) + g A ( t ) S A ( t ) = R Y ( t ) + i ( t ) A f ( t ) + i d ( t ) g A ( t ) A A ( t )
where
g H ( t ) = real purchase price of dwellings;
g A ( t ) = real purchase price of the financial risky assets; and
R Y ( t ) = real disposable income.
Mathematically, to maximize lifetime utility (1) subjected to (3)–(6), a Lagrange equation is structured. From the first-order conditions (see Appendix A.1), we achieved three Equations (7)–(9).
L ( t ) H ( t ) : E [ μ H ( t ) ] = g H ( t ) E [ μ C ( t ) ] δ ( t ) E [ r H ( t ) μ C ( t ) ] e r t E d d t [ e r t μ C ( t ) ]
L ( t ) A f ( t ) : [ i ( t ) π ( t ) ] E [ μ C ( t ) ] = e r t E d d t [ e r t μ C ( t ) ]
L ( t ) A A ( t ) : E [ [ r A ( t ) + i d ( t ) π ( t ) ] μ C ( t ) ] = e r t E d d t [ e r t μ C ( t ) ]
where
μ H ( t ) = μ [ H ( t ) , C ( t ) ] / H ( t ) , the (uncertain) marginal utility of housing;
μ C ( t ) = μ [ H ( t ) , C ( t ) ] / C ( t ) , the (uncertain) marginal utility of consumption;
r H ( t ) = g H ˙ ( t ) / g H ( t ) , the (uncertain) capital return of housing; and
r A ( t ) = g A ˙ ( t ) / g A ( t ) , the (uncertain) capital return of risky financial asset.
Since r H and r A are functions of future states (e.g., t + d t ), returns on risky assets are uncertain. Similarly, the marginal utilities of both housing and non-housing consumption are uncertain since they are functions of r H and r A , respectively. Due to the same term e r t E d d t [ e r t μ C ( t ) ] in all three of these equations, shown in Appendix A.2, the relationship between risk-free and risky financial assets is derived in (10) through the combination of (8) and (9); similarly, the relationship between housing and risk-free assets is derived in (11) through the combination of (7) and (8).
E [ r A D ( t ) μ C ( t ) ] E [ μ C ( t ) ] = i ( t ) π ( t )
E [ μ H ( t ) ] E [ μ C ( t ) ] = i ( t ) π ( t ) + δ ( t ) E ( r H ( t ) ) + τ ( t ) g H ( t )
where
E [ μ H ( t ) ] = E [ μ [ H ( t ) , C ( t ) ] ] / H ( t ) , the expected marginal utility of housing;
E [ μ C ( t ) ] = E [ μ [ H ( t ) , C ( t ) ] ] / C ( t ) , the expected marginal utility of consumption;
r A D ( t ) = r A ( t ) + i d ( t ) π ( t ) , the (uncertain) total return of risky financial assets; and
τ ( t ) = C o v [ r H ( t ) , μ C ( t ) ] / E [ μ C ( t ) ] , the housing market risk premium.
Equation (10) is identical to the key derivation of C-CAPM [13], and Equation (11) is identical to the key derivation indicated in the two-asset model [1], where a housing market risk premium term, τ , appeared in the house price derivation compared to the standard derivation [9].
The left-hand side of (11) is the expected marginal rate of substitution between H and C, meaning that households must cut down on unit housing service if they want to obtain R units of non-housing goods as compensation; therefore, R is the imputed real housing rental price. The term in [ ] on the right-hand side of (11) is the widely used standard definition of the real housing user cost of capital (UCC). All of these factors are summarized by Equation (12). This equation can be also rewritten as g ( t ) = R ( t ) / U C C ( t ) , the well-known arbitrage condition.
R ( t ) = E [ M R S H , C ( t ) ] = E [ μ H ( t ) ] E [ μ C ( t ) ] = U C C ( t ) g H ( t )
where
R ( t ) = the imputed real housing rental price;
E [ M R S H , C ( t ) ] = the expected marginal rate of substitution between H and C; and
U C C ( t ) = i ( t ) π ( t ) + δ ( t ) E [ r H ( t ) ] + τ ( t ) , the user cost of capital.
Equations (10)–(12), particularly (11), are the general solutions of our key interest. However, since they are too general, we must try very hard to use them to estimate housing market risk premiums in empirical works.

2.2. Special Solutions under the Specific Utility Functions

The last section theoretically derived and described the general solution of house prices and housing market risk premiums. In principle, the solution can be adopted for all countries, which is also why a general solution is too abstract. To be able to derive intuitive results of τ for a specific region, specific utility functions and distribution of the returns are needed. For the purpose of showing the calculation process, the distribution of the risky returns is here simplistically assumed to be joint normal as an example, given by (13). Similarly, the Constant Absolute Risk Aversion (CARA) and Constant Relative Risk Aversion (CRRA) utilities [26] given by (14) and (15), respectively, are employed.
r H ( t ) r A D ( t ) N r H e ( t ) r A D e ( t ) , σ H 2 ( t ) σ A H ( t ) σ A H ( t ) σ A D 2 ( t )
μ [ H ( t ) , C ( t ) ] = e x p [ φ H ( t ) ] e x p [ φ C ( t ) ]
μ [ H ( t ) , C ( t ) ] = H ( t ) ( 1 θ ) 1 θ + C ( t ) ( 1 θ ) 1 θ
where
φ = the parameter of CARA;
θ = the parameter of CRRA;
σ A H = ρ × σ A D × σ H , the covariance between r A D and r H ; and
ρ = the correlation between r A D and r H .
Since there are two utility functions, Appendix A.3 and Appendix A.4 separately show the process of deriving housing market risk under the CARA and CRRA utilities, and the special solutions are given by (16)–(18) and (19)–(21), respectively.
τ C A R A = τ 1 , C A R A + τ 2 , C A R A
τ 1 , C A R A ( t ) = [ r A D ( t ) ( i ( t ) π ( t ) ) ] × [ ρ ( t ) σ H ( t ) / σ A D ( t ) ]
τ 2 , C A R A ( t ) = φ g H ( t ) H ( t ) ( 1 ρ 2 ( t ) ) σ H 2 ( t )
τ C R R A = τ 1 , C R R A + τ 2 , C R R A
τ 1 , C R R A ( t ) = [ r A D ( t ) ( i ( t ) π ( t ) ) ] × [ ρ ( t ) σ H ( t ) / σ A D ( t ) ]
τ 2 , C R R A ( t ) = θ C e ( t ) θ 1 E [ μ C ( t ) ] g H ( t ) H ( t ) ( 1 ρ 2 ( t ) ) σ H 2 ( t ) [ θ / C e ( t ) ] g H ( t ) H ( t ) ( 1 ρ 2 ( t ) ) σ H 2 ( t )
First, both τ C A R A and τ C R R A consist of two parts: τ 1 and τ 2 . Second, the τ 1 of the two cases is identical. When risky financial assets are regarded as the whole market, and housing is regarded as a particular asset in this market, r A D ( i π ) becomes the expected exceed return of the financial market, and ρ × σ H / σ A D (or σ A H / σ A D 2 ) becomes the systematic risk of housing in finance, also known as Sharpe β [27]. Therefore, τ 1 is the housing market risk premium affected by the financial market elements. Third, τ 2 is the housing market risk directly caused by the uncertainty of the housing market yield, which has been discussed in the two-asset model [1].
In light of these considerations, τ 1 is the main focus in this article. It is determined by three financial market factors: the expected exceed return on the risky financial asset ( r A D ( i π ) ), the standard deviation ( σ A D ), and the correlation ( ρ ). In principle, the former two are positive, and the latter is in the interval [ 1 , 1 ] . The correlation is key in determining the effects of other financial market elements. When ρ > 0 , relatively higher expected returns or lower volatility of the financial market indicates lower real house prices; when ρ < 0 , real house prices are positively determined by the expected return and negatively determined by the volatility. Simply speaking, when the two markets are strongly “substituted,” financial market flourishing causes lower house prices; when the two markets are strongly “compensated,” financial market flourishing indicates higher house prices; and when the correlation is relatively tiny, the effects of financial market factors on housing market risk premiums, as well as house prices, can be neglected.
Our τ 2 derivations (18) and (21) are similar to Equations (18) and (19) shown in the two-asset model [1]. The effects of risk aversion level ( φ and θ / C ), real housing wealth ( g H H ), and the volatility of housing capital return ( σ H 2 ) have been discussed and analysed [1]. The only difference is that term ( 1 ρ 2 ) newly appears in our τ 2 in both CARA and CRRA cases. Whether positive or negative, a strong correlation leads to ( 1 ρ 2 ) 0 , implying that the housing market risk premium caused by the uncertainty of the housing market ( τ 2 ) will be eliminated when the two markets reflect strong correlations in their market returns because, when the two markets have no interaction, the volatility of the housing capital return ( σ H 2 ) can fully capture the uncertainty of the housing market. In contrast, when the two markets are related, the correlated part ( σ H 2 ρ 2 ) will decrease from the volatility since this risk is dispersed by households’ asset portfolios.
Therefore, the special solutions indicate that financial market factors determine the housing market risk premium and thus affect real house prices. The correlation is the key in both of τ 1 and τ 2 . A tiny correlation between the two markets means tiny financial market impacts but significant τ 2 . However, despite significant correlations being able to hide τ 2 , the effects of financial market factors on house prices are significant. Hence, we can see that both de- and over-financialization are bad for housing markets. To answer the question of the appropriate level of financialization of the housing market, the risk dispersion mechanism through market correlation must be further studied.

2.3. Portfolios between Housing and Risky Financial Assets

Derivation τ 1 reflects a highly similar structure to the systematic risk derivation or the Sharpe β . Remember that both CAPM and its predecessor, the mean-variance model [28], are concentrated on portfolios among risky assets and then derive the price/return to risk premium relation. Here, we undergo a reversed process since we have the price/return to risk premium relation of housing ( τ 1 ) and want to know the portfolio solution.
As shown in Appendix A.5, conditional on the CARA utility function, the optimal ratio of risky financial wealth to housing wealth is derived and given by (22); similarly, conditional on the CRRA utility function, the derivation of the optimal ratio is also (22). Appendix A.6 provides the main body of the mean-variance model [28]. When housing and risky financial assets are the two specific assets A and B, its derivation of the optimal ratio is identical to (22), suggesting that our three-asset model has in-built households optimal portfolio decisions.
k * ( t ) = g A ( t ) A A ( t ) g H ( t ) H ( t ) = r A D e ( t ) [ i ( t ) π ( t ) ] σ H 2 ( t ) + r H T e ( t ) [ i ( t ) π ( t ) ] σ A H ( t ) r H T e ( t ) [ i ( t ) π ( t ) ] σ A D 2 ( t ) + r A D e ( t ) [ i ( t ) π ( t ) ] σ A H ( t )
where r H T e ( t ) = r H e ( t ) + R ( t ) / g ( t ) δ ( t ) , the expected total housing return.
One thing that must be mentioned is that, despite the variance of risky financial assets considering both the variances of capital and dividend yields, the variance of housing is only the variance of the capital yield without housing rental yield because the dividends of financial instruments are usually uncertain, but, practically, almost all housing rents are agreed and paid before tenants move in. In light of this difference, the housing imputed rent is assumed to be certain and thus has no risk.
All of the above results suggest that, when households make their optimal decisions to maximize their expected lifetime utility, their decisions are not only optimal in consumption–investment balance [9] but also include instantaneous optimal portfolio decisions in investing in housing and risky financial assets.

3. Empirical Section

The main purpose of this section is to check whether housing market risk premiums, especially τ 1 , are significant or not. As given in Wang et al. [1], Chinese data are employed for three reasons. First, the China Statistic Office (CSO) provides both monthly and annual house prices. Thus, the volatility of the latter (annual variance) could be calculated from the former. Second, provincial data are available, and therefore, empirical relationships can be analysed by panel data regressions, which are more robust than single equation time series regressions. Third, the housing price index provided by CSO is the price per square meter. Compared to the price per house index, it is less likely that the prices have strong regional heterogeneity.

3.1. Data and Testing

The complete data set is annual panel data including China’s 30 provinces except for Tibet, Taiwan, Hong Kong, and Macao for the period of 2001–2018. Details of the variables are listed in Table 1. The interest rate (i) is proxied by the base interest rate on RMB (CNY, the Chinese currency) loans provided by China Financial Statistical Yearbook; the market price of risky financial assets of China is approximated by the Shanghai Securities Composite Index (SSEC). Except for these two factors, all of the other variables are obtained from China Statistical Yearbook. As shown in Table 1, the annual expected returns, standard deviations, and correlations are calculated through the monthly growths/returns of house prices and SSEC. Since there are no high-quality housing stock data, as discussed in the two-asset model [1], the flow of real housing wealth ( H E / C P I ) is selected to be the approximation of g H ( t ) H ( t ) . Furthermore, for the purpose of avoiding confusion, we only focus on the CRRA case in the empirical section.
In empirical analysis, the stationarity of variables determines the panel model selection. The stationarity of variables in the panel structure is usually tested by LLC, IPS, ADF-Fisher, and Hadri LM stationary tests. The stationarity of l n r y , l n g , and i were already tested by Wang et al. [1], and the results suggested that (log) real house prices and (log) real income contain unit-roots in all panels, but the interest rate (i) is stationary in all panels. Similarly, the stationarity of τ 1 and τ 2 , C R R A is tested by all four of these stationary tests, as shown in Table 2. The results support that both τ 1 and τ 2 , C R R A are stationary for all panels.

3.2. Regression Analysis

The empirically used equations are slightly different compared to the theoretical derivations because: (i) the imposed housing rent (R) is unavailable; and (ii) the logarithm of user cost of capital (UCC) cannot be simply decomposed in the linear form of its components. The logarithm form of standard (12) is Equation (23). Many studies (e.g., [8]) have described how to derive a simple empirical equation such as Equation (24) based on the theoretical formula (23). At the same time, our three-asset derivation implies that risk terms should be added to the standard UCC. Therefore, our empirical equation should be (25). The vector X is used to show the explanatory variables of imposed housing rent. Personal disposable income is selected to capture the change in l n ( R ) . Next, in line with previous studies [8], the sum of expected real housing capital return and inflation is approximated by the last nominal house price growth rate ( r H ( t ) + π ( t ) = r h ( t 1 ) ). In addition, we estimate the households’ relative risk aversion level ( θ × H E / C O N ), which is 1.015 [29]. Then, we rely on the values of 2010, and θ is set to be 5.0 ( 1.015 × C O N / H E ).
l n ( G ) = l n ( R ) l n ( U C C ) = l n ( R ( X ) ) l n ( i π + δ r H )
l n g t = a 0 + a 1 l n ( X ) t + a 2 i t + a 3 r h t 1 + ε t
l n g t = a 0 + a 1 l n ( X ) t + a 2 i t + a 3 r h t 1 + a 4 t a u 1 t + a 5 t a u 2 t + ε t
Since the key variables l n g and l n r y are non-stationary, co-integration should be noted. The error correction type of Equation (25) has three panel relations: they are: (26)–(28). Heterogeneity is not considered in (26); (27) only considers heterogeneity in the short-term relation; and (28) implies that heterogeneity exists in all of the short-term, long-term, and error-correction processes. In addition, u i is the individual effect, which is necessary for panel data analysis.
Δ l n g t = β Δ l n r y t + η [ l n g t 1 a 1 l n r y t 1 a 2 i t 1 a 3 r h t 2 a 4 t a u 1 t 1 a 5 t a u 2 t 1 ] + u i + ε t
Δ l n g t = β i Δ l n r y t + η i [ l n g t 1 a 1 l n r y t 1 a 2 i t 1 a 3 r h t 2 a 4 t a u 1 t 1 a 5 t a u 2 t 1 ] + u i + ε t
Δ l n g t = β i Δ l n r y t + η i [ l n g t 1 a i 1 l n r y t 1 a i 2 i t 1 a i 3 r h t 2 a i 4 t a u 1 t 1 a i 5 t a u 2 t 1 ] + u i + ε t
Equations (26)–(28) are regressed by the three regressions in the non-stationary heterogeneous panel model: pooled mean-group (PMG), mean-group (MG), and dynamic fixed-effect (DFE). Due to the allowance of heterogeneity in the regression of PMG and MG, the estimated coefficients are the averages of all panels. Estimated coefficients, as well as their t-values, are given in Table 3.
Given the results of Hausman tests among MG&PMG (Prob>chi2 = 0.51), MG&DFE (Prob>chi2 = 1.00) and PMG&DFE (Prob>chi2 = 1.00), Equation (26) is better than (28) and then better than (27), implying that the spatial heterogeneity is not significant, and therefore, all provinces can be regarded as the same in co-integration, including all of the short-term, long-term, and error-correction process.
The UCC is formed based on the estimated values of the DFE shown in Table 3. The ratios among estimated coefficients of the interest rate, expected housing capital return, t a u 1 , and t a u 2 ( a 2 , a 3 , a 4 , and a 5 ) are approximate to 1 : 0.05 : 0.05 : 0.05 . The annual housing depreciation rate, a necessary component of housing UCC, is commonly assumed to be 0.01 [8]. Thus, standard UCC is organized as (29) for comparison. Our UCC is organized as (30). Since the DFE regression is better than the MG and PMG regressions, the empirical equation with log UCC is (31). The results are given in Table 4.
U C C t = i t + 0.01 0.05 × r h t 1
U C C 3 a s s e t , t = i t + 0.01 0.05 × r h t 1 + 0.05 × t a u 1 t + 0.05 × t a u 2 t
Δ l n g t = β Δ l n r y t + η [ l n g t 1 a 1 l n r y t 1 a 2 l n ( U C C ) t 1 ] + u i + ε t
The estimated values in Table 4 imply that, when housing market risk premium terms (both t a u 1 and t a u 2 ) are considered, the significance of l n ( U C C ) becomes greater ( | t | increases from 3.6 to 4.5), and R 2 is also increased. This evidence suggests that the three-asset house price derivation is more suitable in house price analysis and forecasting.
In the two-asset model [1], the empirical effect of t a u 2 was discussed, where t a u 2 plays an important role in restricting the price-return positive circle. In contrast, the mean values in Table 1 give us the impression that t a u 1 is significantly smaller than any of the interest rate (i), housing capital return ( r h ), and t a u 2 . Conditional on their similar marginal effects, the practical effect of t a u 1 seems to be tiny. However, since the correlation between the returns of the two markets practically has a strong cyclical trend rather than being stochastically distributed in the interval [ 1 , 1 ] , the effect of t a u 1 on real house prices could be apparent.

3.3. Historical Correlation between Housing and Financial Market Returns

The correlation between housing and financial market returns ( ρ ) is the key to the housing market risk premium for both τ 1 and τ 2 . If the correlation is close to zero, the impact of financial market factors on housing market risk premiums will be tiny, and our τ 2 will decrease to the two-asset derivation; in contrast, if the correlation is close to 1 or 1 , τ 1 will be significantly non-zero, and τ 2 will converge to zero.
There is a fascinating phenomenon by which, although the correlation between housing and financial market returns is usually close to zero in the long-term measurement, the correlation actually is significantly non-zero in the short-run. Furthermore, the correlation is always negative before financial/economic crises and positive during and after crises. Meen [8] measured the correlation of the UK and showed that, for example, the correlation between house capital return and FTSE return is 0.08 during the period of 1970–2012 but −0.66 during the period of 2001–2005 and 0.69 during the period of 2006–2012. Related works [8] have mentioned this phenomenon too. We also measured the correlation between US housing capital return and S&P return and obtained supporting evidence. For China, the long-term correlation was 0.04 during the period of 2000M04–2020M06, and the short-term correlations were calculated based on the sliding window method, with the window width equal to the last 12 months, as shown in Figure 1.
Figure 1 indicates that the correlation in the short run can be quite large or quite small (≈±0.7). One of the possible interpretations is that, in boom periods, some assets will be more attractive than others, causing changes in households’ optimal portfolio decisions ( k * ) and then causing reversed changes in demand for different assets, but in bust periods, the consumption–investment decision dominates households’ behaviours, and therefore, the demands of all risky assets will fall (and recover afterward). If it makes sense that the correlation is negative in the boom period, τ 1 will be negative, and τ 2 will be positive (but small) so that the housing market risk premium could be negative, causing relatively lower UCC and relatively higher house prices; when a crisis occurs, the correlation changes from negative to positive. τ 1 should be positive, and τ 2 should also be positive (but small) so that the housing market risk premium should be positive; therefore, UCC goes up, and house prices go down. In light of this outcome, if the real estate market is the origin of the financial crisis, the changed correlation will aggravate the severity of the financial crisis. In the next section, on the basis of reasonable assumptions about the correlation and other factors, we simulate the influence of the change in correlation on real house prices.

4. Simulation Analysis

In this section, we attempt to analyse the dynamic impact of the correlation between housing capital return and financial market return on real house prices under our three-asset model. Generally, housing market risk premiums are generated as the CRRA derivation (19)–(21), housing user cost of capital (UCC) is generated as the empirical Equation (30), and real house prices are generated as the long-term housing price Equation (32), which is the long-run relationship implied by the ECM Equation (31) with the empirical estimations given in Table 4. Note that, since real income will be generated as a non-stationary series, co-integration is in-built when the long-term relationship (32) is applied.
l n g t = a 0 + 0.812 × l n r y t 0.345 × l n ( U C C ) t
The exogenous shocks are contained by the real income series. The mean and standard deviation of the historical real income growth rate are 0.0776 and 0.0807, respectively. We simply assume that the growth rate of real income is distributed normally and thus generate the exogenous noises, while to simulate the cyclical correlation, the correlation is generated as a sine function ( c o r r t = 0.7 × s i n ( t × 2 π / 6.5 ) ), where π is approximate to 3.14, 0.7 is the amplitude, and 6.5 (years) is the cycle period roughly estimated by Figure 1. All of the other details are summarized in Appendix B (Table A1).
Three curves are simulated for the purposes of comparison. The first is the baseline, where exogenous noises are not considered in real income growth. Therefore, house prices increase steadily. Noises are considered in both the second and third. The only difference is that housing market risk premiums are set to be zero in the second, referring to the standard house price model, but the risk premium is considered in the third. It is important to mention that calibration simulation (such as the first curve) of our three-asset model is meaningless because the generation of the risk premium requires volatilities. However, the consideration of noises will lead to a case in which we are not sure whether the change in real house prices is caused by the consideration of the risk premium or the noises. In light of this uncertainty, the stochastic process was replicated 100 times, and we calculated the averages to remove the influence of noise.
Figure 2 shows the key results of the simulation. The first frame shows the stochastic and non-stochastic real income growth rates over 25 years (periods). The second frame shows the generated sine shape correlation. To reflect the properties of risk premium terms, the third frame shows 55 years of results. Doubtlessly, since tau1 is determined by the correlation, it also reflects a sine shape movement. As in the analysis and simulation shown by Wang et al. [1], tau2 is increasing in recent years in China. Tau, the sum of the two, also reflects a strong cyclicity, but compared with tau1, the peaks are wider, and the troughs are narrower, indicating that housing market risk premiums in most of the periods have the relatively strong ability to restrain house price growth; when the financial market and housing market are negatively correlated, the restrictive ability will Decrease and thus cause housing price bubbles. As clearly shown by the fourth frame, when the risk premium terms are involved, there will be obvious periodicity in house prices. Therefore, although housing market risk premiums have a negative impact on real house prices, the natural change is not conducive to stabilizing house price. For instance, real house prices decreased 9.6% in the two years from period 18 to 20 and fell 9.3% in one year from period 25 to 26. These simulated values support our view that the effect of financial market factors on housing market risk premiums, as well as real house prices, cannot be ignored in practice and should be of sufficient concern in both theoretical studies and government policy making.

5. Conclusions

This paper aims to analyse the effect of financial market factors on real house prices. By introducing the risky financial asset, we first develop the theoretical framework from a two-asset base to a three-asset model. Theoretical derivations show that three financial market elements—expected return, volatility and the correlation—determine housing market risk premiums and thus affect the housing user cost of capital, as well as real house prices. If correlations between the returns on housing and financial markets are positive, real house prices are negatively affected by expected returns and positively affected by the volatility of the risky financial asset; if the correlation is negative, the effects are reversed. Empirical works have supported the existence of the theoretical derivations of the housing market risk premiums ( τ 1 and τ 2 ), but the estimated marginal effects are relatively small. Statistics have also shown that, despite the correlation between the two markets being close to zero in the long-term measurement, it is significantly negative in boom periods and significantly positive in bust periods. Given reasonable assumptions, our dynamic simulations point out that financial markets cause cyclical changes in real house prices. More specifically, the results show that China’s real house prices could decrease by approximately 10% within one to two years in the future if the effects of financial market factors are still as large as they are now.
Therefore, for the purpose of stabilizing housing prices and avoiding housing market collapse, we recommend the following: (i) government should control the money flow from one to another between the housing and financial markets; and (ii) compared with the housing market, financial markets should have an absolute advantage (relatively higher expected return and lower risk).
Although the three-asset model is very complex, it is still incomplete. Housing mortgage debt must necessarily be introduced into the model next. Then, the difference between risk-free interest rates and mortgage interest rates will be allowed. The effect of monetary policy on housing prices can be more precisely analysed.

Author Contributions

Conceptualization, Y.W.; methodology, Y.W.; software, Y.W. and J.L.; validation, Y.W. and J.L.; formal analysis, Y.W.; investigation, Y.W.; resources, Y.W.; data curation, Y.W. and J.L.; writing—original draft preparation, Y.W., J.L., Z.Q. and S.S.; writing—review and editing, Y.W., J.L., Z.Q. and S.S.; visualization, Y.W.; supervision, Y.W.; project administration, Y.W.; funding acquisition, Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by (China’s) Shandong Social Science Planning Research Project (grant number: 19DTJJ02).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Acknowledgments

This work was supported by the Faculty of Economics and the Centre of Excellence in Econometrics at Chiang Mai University, the China-CASEAN High-Quality Development Research Center in Shandong University of Finance and Economics, and the education plan of the youth and creative talents in Shandong higher education institutions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Mathematical Derivations

Appendix A.1. The Expected Utility Three-Asset Housing Life-Cycle Model, Part I

According to the lifetime utility function (1) and constraints (3)–(6), the households’ decisions can be found by maximizing the Hamiltonian function:
L ( t ) = E { 0 e r t { μ [ H ( t ) , C ( t ) ] + λ ( t ) [ R Y ( t ) C ( t )
g H ( t ) H ˙ ( t ) + δ ( t ) H ( t )
[ A f ˙ ( t ) [ i ( t ) π ( t ) ] A f ( t ) ]
g A ( t ) [ A A ˙ ( t ) [ i d ( t ) π ( t ) ] A A ( t ) ] ] } d t } ,
The first-order conditions are
L ( t ) C ( t ) = E e r t μ C ( t ) e r t λ ( t ) = 0 ,
L ( t ) H ( t ) = E e r t μ H ( t ) e r t λ ( t ) g H ( t ) δ ( t ) + d d t [ e r t λ ( t ) g H ( t ) ] = 0 ,
L ( t ) A f ( t ) = E e r t λ ( t ) [ i ( t ) π ( t ) ] + d d t [ e r t λ ( t ) ] = 0 ,
L ( t ) A A ( t ) = E e r t λ ( t ) g A ( t ) [ i d ( t ) π ( t ) ] + d d t [ e r t λ ( t ) g A ( t ) ] = 0 .
Note that all variables of the are certain in time period (t) and uncertain in the future ( t + d t , t + 2 d t , t + 3 d t , ) Therefore, certain elements can be removed to the outside of expectations.
From (A2),
μ C ( t ) = λ ( t ) .
From (A3),
E [ e r t μ H ( t ) ] = E [ e r t λ ( t ) g ( t ) δ ( t ) ] E d d t [ e r t λ ( t ) g H ( t ) ]
e r t E [ μ H ( t ) ] = e r t E [ λ ( t ) ] g H ( t ) δ ( t ) E g H ˙ ( t ) e r t λ ( t ) g H ( t ) d d t [ e r t λ ( t ) ]
e r t E [ μ H ( t ) ] = e r t E [ λ ( t ) ] g H ( t ) δ ( t ) e r t E g H ˙ ( t ) λ ( t ) g H ( t ) E d d t [ e r t λ ( t ) ]
E [ μ H ( t ) ] = g H ( t ) E [ λ ( t ) ] δ ( t ) E r H ( t ) λ ( t ) e r t E d d t [ e r t λ ( t ) ]
where r H ( t ) = g H ˙ ( t ) / g H ( t ) , real capital return of housing.
From (A4),
[ i ( t ) π ( t ) ] E [ λ ( t ) ] = e r t E [ d d t [ e r t λ ( t ) ] ] .
From (A5),
e r t g A ( t ) E [ λ ( t ) [ i d ( t ) π ( t ) ] ] = E d d t [ e r t λ ( t ) g A ( t ) ]
e r t g A ( t ) E [ λ ( t ) [ i d ( t ) π ( t ) ] ] = E g A ˙ ( t ) e r t λ ( t ) g A ( t ) d d t [ e r t λ ( t ) ]
E [ [ i d ( t ) π ( t ) ] λ ( t ) ] = E r A ( t ) λ ( t ) e r t E d d t [ e r t λ ( t ) ]
E [ [ r A ( t ) + i d ( t ) π ( t ) ] λ ( t ) ] = e r t E d d t [ e r t λ ( t ) ]
where r A ( t ) = g A ˙ ( t ) / g A ( t ) is real capital return of risky financial assets.
Substitute (A6) into (A7)–(A9),
L ( t ) H ( t ) : E [ μ H ( t ) ] = g H ( t ) E [ μ C ( t ) ] δ ( t ) E [ r H ( t ) μ C ( t ) ] e r t E d d t [ e r t μ C ( t ) ]
L ( t ) A f ( t ) : [ i ( t ) π ( t ) ] E [ μ C ( t ) ] = e r t E [ d d t [ e r t μ C ( t ) ] ]
L ( t ) A A ( t ) : E [ [ r A ( t ) + i d ( t ) π ( t ) ] μ C ( t ) ] = e r t E d d t [ e r t μ C ( t ) ]
Equations (A10)–(A12) are Equations (7)–(9) in the main text.

Appendix A.2. The Expected Utility Three-Asset Housing Life-Cycle Model, Part II

In this part, conditions are Equations (A10)–(A12).
Substitute (A11) into (A12),
E [ [ r A ( t ) + i d ( t ) π ( t ) ] μ C ( t ) ] = [ i ( t ) π ( t ) ] E [ μ C ( t ) ]
E [ r A D ( t ) μ C ( t ) ] E [ μ C ( t ) ] = i ( t ) π ( t )
where r A D ( t ) = r A ( t ) + i d ( t ) π ( t ) and this is Equation (10) in main text.
Or since E [ A B ] = E [ A ] E [ B ] + C o v [ A , B ] ,
E [ r A D ( t ) ] ( i ( t ) π ( t ) ) = C o v [ r A D ( t ) , μ C ( t ) ] E [ μ C ( t ) ] .
Substitute (A11) into (A10):
E [ μ H ( t ) ] = { E [ μ C ( t ) ] δ ( t ) E [ r H ( t ) μ C ( t ) ] + [ i ( t ) π ( t ) ] E [ μ C ( t ) ] } g H ( t )
E [ μ H ( t ) ] E [ μ C ( t ) ] = i ( t ) π ( t ) + δ ( t ) E [ r H ( t ) μ C ( t ) ] E [ μ C ( t ) ] g H ( t ) .
Or since E [ A B ] = E [ A ] E [ B ] + C o v [ A , B ] ,
E [ μ H ( t ) ] E [ μ C ( t ) ] = i ( t ) π ( t ) + δ ( t ) E [ r H ( t ) ] C o v [ r H ( t ) , μ C ( t ) ] E [ μ C ( t ) ] g H ( t ) .
Equations (A14) and (A15) are Equations (11)–(12) in the main text.

Appendix A.3. The Three-Asset Expected Utility Housing Life-Cycle Model under CARA Utility

In this part, conditions are Equations (13) and (14).
Since r H and r A D are joint normally distributed, the probability density function of the standard joint normal distribution between z 1 and z 2 is used, given as:
ϕ ( z 1 , z 2 ) = 1 2 π 1 ρ 2 e x p 1 2 ( 1 ρ 2 ) [ z 1 2 2 ρ z 1 z 2 + z 2 2 ]
where ρ is the correlation between z 1 and z 2 .
Equations (3)–(6) can be organized as one equation:
C ( t ) = R Y ( t ) + i ( t ) A f ( t ) + i d ( t ) A A ( t ) g H ( t ) [ H ˙ ( t ) + δ ( t ) H ( t ) ]
[ A ˙ ( t ) + π ( t ) A ( t ) ] g A ( t ) [ A A ˙ ( t ) + π ( t ) A A ( t ) ]
or
C ( t ) R Y ( t ) + [ r H ( t ) δ ( t ) ] g H ( t ) H ( t )
              + [ i ( t ) π ( t ) ] A f ( t ) + r A D ( t ) g A ( t ) A A ( t )
This equation implies that, if the quantities of households’ housing and non-housing assets are unchanged at the moment, all of their labour income, housing net yield and non-housing financial yield will be consumed. Further, since the housing capital return ( r H ( t ) ) and the total return of risky financial assets ( r A D ( t ) ) are uncertain, households’ consumption (A17) can be divided into certain and uncertain parts as:
C ( t ) = R Y ( t ) + [ r H e ( t ) δ ( t ) ] g H ( t ) H ( t ) + g H ( t ) H ( t ) σ H ( t ) z 1
+ [ i ( t ) π ( t ) ] A f ( t ) + r A D e ( t ) g A ( t ) A A ( t ) + g A ( t ) A A ( t ) σ A D ( t ) z 2
Therefore,
E [ C ( t ) ] C e ( t ) = R Y ( t ) + [ r H e ( t ) δ ( t ) ] g H ( t ) H ( t )
+ [ i ( t ) π ( t ) ] A f ( t ) + r A D e ( t ) g A ( t ) A A ( t )
or
C ( t ) = C e ( t ) + g H ( t ) H ( t ) σ H ( t ) z 1 + g A ( t ) A A ( t ) σ A D ( t ) z 2
Remember that, to derive a specific solution for Equation (12) requires solving E [ μ C ( t ) ] , E [ r H ( t ) μ C ( t ) ] and E [ r A D ( t ) μ C ( t ) ] .
Step 1: The derivation of E [ μ C ( t ) ] .
Under the specific CARA utility function (14),
E [ μ C ( t ) ] = E [ φ e x p [ φ C ( t ) ] ] .
Substitute (A18),
E [ μ C ( t ) ] = E [ φ e x p [ φ C e ( t ) φ g H ( t ) H ( t ) σ H ( t ) z 1 φ g A ( t ) A A ( t ) σ A D ( t ) z 2 ] ]
E [ μ C ( t ) ] = φ e x p [ φ C e ( t ) ] × E [ e x p [ φ g H ( t ) H ( t ) σ H ( t ) z 1 φ g A ( t ) A A ( t ) σ A D ( t ) z 2 ] ]
Based on the standard joint normal distribution,
E [ μ C ( t ) ] = φ e x p [ φ C e ( t ) ] × + { e x p [ φ g H ( t ) H ( t ) σ H ( t ) z 1
φ g A ( t ) A A ( t ) σ A D ( t ) z 2 ] } ϕ ( z 1 , z 2 ) d z 1 d z 2
Substituting the probability density function (A16),
E [ μ C ( t ) ] = φ e x p [ φ C e ( t ) ] × + { e x p [ φ g H ( t ) H ( t ) σ H ( t ) z 1 φ g A ( t ) A A ( t )
σ A D ( t ) z 2 ] } 1 2 π 1 ρ 2 ( t ) e x p { 1 2 ( 1 ρ 2 ( t ) ) [ z 1 2 2 ρ ( t ) z 1 z 2 + z 2 2 ] } d z 1 d z 2
E [ μ C ( t ) ] = φ e x p [ φ C e ( t ) ] × + 1 2 π 1 ρ 2 ( t ) e x p { 1 2 ( 1 ρ 2 ( t ) ) [ z 1 2 2 ρ ( t ) z 1 z 2
+ z 2 2 + 2 ( 1 ρ 2 ( t ) ) φ g H ( t ) H ( t ) σ H ( t ) z 1 + 2 ( 1 ρ 2 ( t ) ) φ g A ( t ) A A ( t ) σ A D ( t ) z 2 ] } d z 1 d z 2
Because of the mathematical law that
x 2 2 ρ x y + y 2 + 2 ( a ρ b ) x + 2 ( b ρ a ) y + ( a 2 2 ρ a b + b 2 )
= ( x + a ) 2 2 ρ ( x + a ) ( y + b ) + ( y + b ) 2 .
Let, x = z 1 ; a ρ ( t ) b = ( 1 ρ 2 ( t ) ) φ g H ( t ) H ( t ) σ H ( t ) ; y = z 2 ;
b ρ ( t ) a = ( 1 ρ 2 ( t ) ) φ g A ( t ) A A ( t ) σ A D ( t ) .
Therefore, a = φ g H ( t ) H ( t ) σ H ( t ) + ρ ( t ) g A ( t ) A A ( t ) σ A D ( t ) b = φ ρ ( t ) g H ( t ) H ( t ) σ H ( t ) + g A ( t ) A A ( t ) σ A D ( t ) .
Thus, a 2 2 ρ a b + b 2 = φ 2 ( 1 ρ 2 ( t ) ) M ( t )
where M ( t ) = g H ( t ) H ( t ) σ H ( t ) 2 + 2 ρ ( t ) g H ( t ) H ( t ) σ H ( t ) g A ( t ) A A ( t ) σ A D ( t ) + g A ( t ) A A ( t ) σ A D ( t ) 2 .
Therefore, E [ μ C ( t ) ] = φ e x p [ φ C e ( t ) ] × + 1 2 π 1 ρ 2 ( t ) e x p { 1 2 ( 1 ρ 2 ( t ) )
[ ( x + a ) 2 2 ρ ( x + a ) ( y + b ) + ( y + b ) 2 ( a 2 2 ρ a b + b 2 ) ] } d z 1 d z 2
E [ μ C ( t ) ] = φ e x p [ φ C e ( t ) ] × e x p { a 2 2 ρ a b + b 2 2 ( 1 ρ 2 ( t ) ) } + 1 2 π 1 ρ 2 ( t )
e x p { 1 2 ( 1 ρ 2 ( t ) ) [ ( x + a ) 2 2 ρ ( x + a ) ( y + b ) + ( y + b ) 2 ] } d z 1 d z 2
E [ μ C ( t ) ] = φ e x p [ φ C e ( t ) ] × e x p φ 2 ( 1 ρ 2 ( t ) ) M ( t ) 2 ( 1 ρ 2 ( t ) )
× + ϕ [ ( z 1 + a ) , ( z 2 + b ) ] d z 1 d z 2
E [ μ C ( t ) ] = φ e x p [ φ C e ( t ) ] × e x p 1 2 φ 2 M ( t ) .
Step 2: The derivation of E [ r H ( t ) μ C ( t ) ] .
E [ r H ( t ) μ C ( t ) ]
= E [ [ r H e ( t ) + σ H ( t ) z 1 ] μ C ( t ) ]
= r H e ( t ) E [ μ C ( t ) ] + E { σ H ( t ) z 1 μ C ( t ) }
= r H e ( t ) E [ μ C ( t ) ] + E { σ H ( t ) z 1 φ e x p [ φ C ( t ) ] ]
= r H e ( t ) E [ μ C ( t ) ] + E { σ H ( t ) z 1 φ e x p [ φ [ C e ( t ) + g H ( t ) H ( t ) σ H ( t ) z 1
+ g A ( t ) A A ( t ) σ A D ( t ) z 2 ] ] }
= r H e ( t ) E [ μ C ( t ) ] + + { σ H ( t ) z 1 φ e x p [ φ [ C e ( t ) + g H ( t ) H ( t ) σ H ( t ) z 1
+ g A ( t ) A A ( t ) σ A D ( t ) z 2 ] ] } ϕ ( z 1 , z 2 ) d z 1 d z 2
= r H e ( t ) E [ μ C ( t ) ] + + { σ H ( t ) z 1 φ e x p [ φ C e ( t ) φ g H ( t ) H ( t ) σ H ( t ) z 1
φ g A ( t ) A A ( t ) σ A D ( t ) z 2 ] } 1 2 π 1 ρ 2 ( t ) e x p { 1 2 ( 1 ρ 2 ( t ) ) [ z 1 2 2 ρ z 1 z 2 + z 2 2 ] } d z 1 d z 2
= r H e ( t ) E [ μ C ( t ) ] + + { σ H ( t ) z 1 φ × e x p [ φ C e ( t ) ] 1 2 π 1 ρ 2 ( t ) e x p [ φ g H ( t ) H ( t ) σ H ( t ) z 1
φ g A ( t ) A A ( t ) σ A D ( t ) z 2 ] } e x p { 1 2 ( 1 ρ 2 ( t ) ) [ z 1 2 2 ρ ( t ) z 1 z 2 + z 2 2 ] } d z 1 d z 2
= r H e ( t ) E [ μ C ( t ) ] + σ H ( t ) φ × e x p [ φ C e ( t ) ] + { 1 2 π 1 ρ 2 ( t ) z 1 × e x p { 1 2 ( 1 ρ 2 ( t ) )
[ z 1 2 2 ρ ( t ) z 1 z 2 + z 2 2 + 2 ( 1 ρ 2 ( t ) ) φ g H ( t ) H ( t ) σ H ( t ) z 1 + 2 ( 1 ρ 2 ( t ) ) φ g A ( t ) A A ( t ) σ A D
( t ) z 2 ] } d z 1 d z 2 .
Because of the mathematical law that ϕ ( x , y ) d x d y = ϕ ( x , y ) d y d x , where ϕ ( x , y ) d y = ϕ x ( x ) = 1 2 π e x p 1 2 x 2 is the marginal probability of x.
(Kenney and Keeping, 1951, p. 202).
Then,
E [ r H ( t ) μ C ( t ) ] = r H e ( t ) E [ μ C ( t ) ] + φ e x p [ φ C e ( t ) ] σ H ( t ) × e x p 1 2 φ 2 M ( t ) + z 1 f z 1 ( z 1 ) d z 1 , where f z 1 ( z 1 ) = 1 2 π e x p 1 2 z 1 + φ g H ( t ) H ( t ) σ H ( t ) + ρ ( t ) φ g A ( t ) A A ( t ) σ A D ( t ) 2 .
Therefore,
E [ r H ( t ) μ C ( t ) ] = r H e ( t ) E [ μ C ( t ) ] + φ e x p [ φ C e ( t ) ] σ H ( t ) × e x p 1 2 φ 2 M ( t )
× φ g H ( t ) H ( t ) σ H ( t ) ρ ( t ) φ g A ( t ) A A ( t ) σ A D ( t )
E [ r H ( t ) μ C ( t ) ] = r H e ( t ) E [ μ C ( t ) ] + φ g H ( t ) H ( t ) σ H 2 ( t ) φ g A ( t ) A A ( t ) [ ρ ( t ) σ H ( t ) σ A D ( t ) ] E [ μ C ( t ) ] .
Step 3: The derivation of E [ r A D ( t ) μ C ( t ) ] .
Similarly,
E [ r A D ( t ) μ C ( t ) ] = r A D e ( t ) E [ μ C ( t ) ] + φ e x p [ φ C e ( t ) ] σ A D ( t ) × e x p { 1 2 φ 2 M ( t ) }
× ρ ( t ) φ g H ( t ) H ( t ) σ H ( t ) φ g A ( t ) A A ( t ) σ A D ( t )
E [ r A D ( t ) μ C ( t ) ] = r A D e ( t ) E [ μ C ( t ) ] + φ g A ( t ) A A ( t ) σ A D 2 ( t ) φ g H ( t ) H ( t ) [ ρ ( t ) σ H ( t ) σ A D ( t ) ] E [ μ C ( t ) ] .
Step 4: Calculate the UCC as well as tau ( τ C A R A ) under CARA.
Substitute (A21) into (10):
i ( t ) π ( t ) = r A D e ( t ) + φ g A ( t ) A A ( t ) σ A D 2 ( t ) φ g H ( t ) H ( t ) [ ρ ( t ) σ H ( t ) σ A D ( t ) ]
φ g A ( t ) A A ( t ) = r A D e ( t ) ( i ( t ) π ( t ) ) σ A D 2 ( t ) φ g H ( t ) H ( t ) ρ ( t ) σ H ( t ) σ A D ( t ) .
Substitute (A20) into (11):
E [ μ H ( t ) ] E [ μ C ( t ) ] = { i ( t ) π ( t ) + δ ( t ) r H e ( t ) + τ C A R A ( t ) } g H ( t )
where τ C A R A ( t ) = φ g H ( t ) H ( t ) σ H 2 ( t ) + φ g A ( t ) A A ( t ) [ ρ ( t ) σ H ( t ) σ A D ( t ) ] .
Substitute (A22) into τ C A R A ( t ) ; then we have:
τ C A R A ( t ) = τ 1 , C A R A ( t ) + τ 2 , C A R A ( t )
τ 1 , C A R A ( t ) = [ r A D ( t ) ( i ( t ) π ( t ) ) ] × [ ρ ( t ) σ H ( t ) / σ A D ( t ) ]
τ 2 , C A R A ( t ) = φ g H ( t ) H ( t ) ( 1 ρ 2 ( t ) ) σ H 2 ( t )
Equations (A23)–(A25) are Equations (16)–(18) in the main text.

Appendix A.4. The Three-Asset Expected Utility Housing Life-Cycle Model under CRRA Utility

In this part, conditions are Equations (13), (15), (A17) and (A18).
Similarly, to derive the specific solution of Equation (12) requires solving E [ μ C ( t ) ] , E [ r H ( t ) μ C ( t ) ] and E [ r A D ( t ) μ C ( t ) ] .
Step 1: The derivation of E [ μ C ( t ) ] .
E [ μ C ( t ) ] = + μ C ( t ) ϕ ( z 1 , z 2 ) d z 1 d z 2
E [ μ C ( t ) ] = + C ( t ) θ ϕ ( z 1 , z 2 ) d z 1 d z 2
where ϕ ( z 1 , z 2 ) is defined in Appendix A.3.
Let us take a second-order Taylor expansion for C ( t ) θ , where C ( t ) = C e ( t ) :
C ( t ) θ C e ( t ) θ θ C e ( t ) θ 1 [ C ( t ) C e ( t ) ]
+ 1 2 θ ( θ + 1 ) C e ( t ) θ 2 [ C ( t ) C e ( t ) ] 2 Then,
E [ μ C ( t ) ] + { C e ( t ) θ θ C e ( t ) θ 1 [ C ( t ) C e ( t ) ]
+ 1 2 θ ( θ + 1 ) C e ( t ) θ 2 [ C ( t ) C e ( t ) ] 2 } ϕ ( z 1 , z 2 ) d z 1 d z 2
E [ μ C ( t ) ] + { C e ( t ) θ θ C e ( t ) θ 1 [ g H ( t ) H ( t ) σ H ( t ) z 1 + g A ( t ) A A ( t ) σ A D ( t ) z 2 ]
+ 1 2 θ ( θ + 1 ) C e ( t ) θ 2 [ g H ( t ) H ( t ) σ H ( t ) z 1 + g A ( t ) A A ( t ) σ A D ( t ) z 2 ] 2 } ϕ ( z 1 , z 2 ) d z 1 d z 2 Since the mathematical laws
+ c ϕ ( x , y ) d x d y = c (constant);
+ x ϕ ( x , y ) d x d y = 0 (mean of the joint standard normal distribution);
+ x 2 ϕ ( x , y ) d x d y = 1 (variance of the joint standard normal distribution); and
+ x y ϕ ( x , y ) d x d y = ρ (covariance of the joint standard normal distribution);
Then,
E [ μ C ( t ) ] C e ( t ) θ + 1 2 θ ( θ + 1 ) C e ( t ) θ 2 { [ g H ( t ) H ( t ) σ H ( t ) ] 2 + [ g A ( t ) A A ( t ) σ A D ( t ) ] 2 + 2 ρ ( t ) [ g H ( t ) H ( t ) σ H ( t ) ] [ g A ( t ) A A ( t ) σ A D ( t ) ] } .
Step 2: The derivation of E [ r H ( t ) μ C ( t ) ] .
E [ r H ( t ) μ C ( t ) ] = + r H ( t ) μ C ( t ) ϕ ( z 1 , z 2 ) d z 1 d z 2
E [ r H ( t ) μ C ( t ) ] = + r H ( t ) C ( t ) θ ϕ ( z 1 , z 2 ) d z 1 d z 2
Let us take a second-order Taylor expansion for r H ( t ) C ( t ) θ , where C ( t ) = C e ( t ) and r H ( t ) = r H e ( t ) :
r H ( t ) C ( t ) θ r H e ( t ) C e ( t ) θ + r H e ( t ) ( θ ) C e ( t ) θ 1 [ C ( t ) C e ( t ) ] + C e ( t ) θ [ r H ( t ) r H e ( t ) ] + 1 2 θ ( θ + 1 ) r H e ( t ) C e ( t ) θ 2 [ C ( t ) C e ( t ) ] 2 + 0 × [ r H ( t ) r H e ( t ) ] 2 + ( θ ) C e ( t ) θ 1 [ C ( t ) C e ( t ) ] [ r H ( t ) r H e ( t ) ]
r H ( t ) C ( t ) θ r H e ( t ) C e ( t ) θ + r H e ( t ) ( θ ) C e ( t ) θ 1 [ g H ( t ) H ( t ) σ H ( t ) z 1 + g A ( t ) A A ( t ) σ A D ( t ) z 2 ] + C e ( t ) θ [ σ H ( t ) z 1 ] + 1 2 θ ( θ + 1 ) r H e ( t ) C e ( t ) θ 2 [ g H ( t ) H ( t ) σ H ( t ) z 1 + g A ( t ) A A ( t ) σ A D ( t ) z 2 ] 2 + ( θ ) C e ( t ) θ 1 [ g H ( t ) H ( t ) σ H ( t ) z 1 + g A ( t ) A A ( t ) σ A D ( t ) z 2 ] [ σ H ( t ) z 1 ]
Therefore,
E [ r H ( t ) μ C ( t ) ] r H e ( t ) C e ( t ) θ + 1 2 θ ( θ + 1 ) r H e ( t ) C e ( t ) θ 2 { [ g H ( t ) H ( t ) σ H ( t ) ] 2 + [ g A ( t ) A A ( t ) σ A D ( t ) ] 2 + 2 ρ ( t ) [ g H ( t ) H ( t ) σ H ( t ) ] [ g A ( t ) A A ( t ) σ A D ( t ) ] } + ( θ ) C e ( t ) θ 1 [ g H ( t ) H ( t ) σ H 2 ( t ) + ρ ( t ) g A ( t ) A A ( t ) σ H ( t ) σ A D ( t ) ]
Step 3: The derivation of E [ r A D ( t ) μ C ( t ) ] .
Similarly,
E [ r A D ( t ) μ C ( t ) ] r A D e ( t ) C e ( t ) θ + 1 2 θ ( θ + 1 ) r A D e ( t ) C e ( t ) θ 2 { [ g H ( t ) H ( t ) σ H ( t ) ] 2 + [ g A ( t ) A A ( t ) σ A D ( t ) ] 2 + 2 ρ ( t ) [ g H ( t ) H ( t ) σ H ( t ) ] [ g A ( t ) A A ( t ) σ A D ( t ) ] } + ( θ ) C e ( t ) θ 1 [ g A ( t ) A A ( t ) σ A D 2 ( t ) + ρ ( t ) g H ( t ) H ( t ) σ H ( t ) σ A D ( t ) ]
Step 4: Calculate the UCC, as well as tau ( τ C R R A ) under CRRA.
Substitute (A26) into (A28):
E [ r A D ( t ) μ C ( t ) ] r A D e ( t ) E [ μ C ( t ) ] θ C e ( t ) θ 1 [ g A ( t ) A A ( t ) σ A D 2 ( t ) + g H ( t ) H ( t ) ρ ( t ) σ H ( t ) σ A D ( t ) ]
Substitute (A29) into (10):
[ r A D e ( t ) ( i ( t ) π ( t ) ) ] E [ μ C ( t ) ] θ C e ( t ) θ 1 g A ( t ) A A ( t ) σ A D 2 ( t ) + g H ( t ) H ( t ) ρ ( t ) σ H ( t ) σ A D ( t )
Substitute (A26) into (A27):
E [ r H ( t ) μ C ( t ) ] r H e ( t ) E [ μ C ( t ) ] θ C e ( t ) θ 1 [ g H ( t ) H ( t ) σ H 2 ( t )
+ ρ ( t ) g A ( t ) A A ( t ) σ H ( t ) σ A D ( t ) ]
Substitute (A30) into (A31):
E [ r H ( t ) μ C ( t ) ] E [ μ C ( t ) ] r H e ( t ) [ r A D ( t ) ( i ( t ) π ( t ) ) ] × [ ρ ( t ) σ H ( t ) / σ A D ( t ) ]
θ C e ( t ) θ 1 E [ μ C ( t ) ] g H ( t ) H ( t ) ( 1 ρ 2 ( t ) ) σ H 2 ( t )
If expectations are ignored, { θ C e ( t ) θ 1 } / { E [ C ( t ) θ ] } = θ / C ( t ) , it will be the Pratt-Arrow coefficient for CRRA utility under certainty ( μ C ( t ) / μ C ( t ) = θ / C ( t ) ). Under this approximation, the housing market risk ( τ C R R A ) is much more easily understood. Note that the coefficient of the uncertain utility ( { θ C e ( t ) θ 1 } / { E [ C ( t ) θ ] } is still solvable, but since (i) national risk attitude is relatively stable, (ii) its difference from θ / C ( t ) is tiny, and (iii) it has a pattern similar to that of τ C A R A derivation, the approximation is more meaningful. Under the approximation, substitute (A32) into (11):
E [ μ H ( t ) ] E [ μ C ( t ) ] { i ( t ) π ( t ) + δ ( t ) r H e ( t ) + τ C R R A ( t ) } g H ( t )
where
τ C R R A ( t ) = τ 1 , C R R A ( t ) + τ 2 , C R R A ( t )
τ 1 , C R R A ( t ) = [ r A D ( t ) ( i ( t ) π ( t ) ) ] × [ ρ ( t ) σ H ( t ) / σ A D ( t ) ]
τ 2 , C R R A ( t ) θ C e ( t ) g H ( t ) H ( t ) ( 1 ρ 2 ( t ) ) σ H 2 ( t )
Equations (A33)–(A35) are Equations (19)–(21) in the main text.

Appendix A.5. Derivation of the Optimal Ratio of Risk Financial Assets to Housing Assets from Our Three-Asset Model

E [ μ H ( t ) ] E [ μ C ( t ) ] = M R S H , C ( t ) = R ( t )
Equation (A36) is the arbitrage condition. The marginal rate of substitution (MRS) between housing and the consumption good is equal to the imputed rental price of housing services (R), indicating that households’ utility is unchanged if their housing services are cut by one unit, but they obtain some non-housing goods as compensation (more details are provided in Wang et al., 2020).
Since Equation (11), the expected total housing return will be
r H T e ( t ) = r H e ( t ) + R ( t ) / g ( t ) δ ( t ) = i ( t ) π ( t ) E [ r H ( t ) μ C ( t ) ] E [ μ C ( t ) ] + r H e ( t )
Part 1: The three-asset model under the CARA utility:
The two conditions are (A20) and (A21). Substitute (A20) into (A37):
r H T e ( t ) [ i ( t ) π ( t ) ] σ H ( t ) = φ g H ( t ) H ( t ) σ H ( t ) + φ g A ( t ) A A ( t ) [ ρ ( t ) σ A D ( t ) ]
Substitute (A21) into (10):
r A D e ( t ) [ i ( t ) π ( t ) ] σ A D ( t ) = φ g A ( t ) A A ( t ) σ A D ( t ) + φ g H ( t ) H ( t ) [ ρ ( t ) σ H ( t ) ]
Now, let φ g H ( t ) H ( t ) and φ g A ( t ) A A ( t ) be the two unknowns of the binary Equations (A38) and (A39). Then,
k * ( t ) = g A ( t ) A A ( t ) g H ( t ) H ( t ) = r A D e ( t ) [ i ( t ) π ( t ) ] σ H 2 ( t ) + r H T e ( t ) [ i ( t ) π ( t ) ] σ A H ( t ) r H T e ( t ) [ i ( t ) π ( t ) ] σ A D 2 ( t ) + r A D e ( t ) [ i ( t ) π ( t ) ] σ A H ( t )
where σ A H ( t ) = ρ ( t ) σ H ( t ) σ A D ( t ) is the covariance of the two returns.
Part 2: The three-asset model under the CRRA utility:
The two conditions are (A29) and (A31). Substitute (A29) into (A37):
r H T e ( t ) [ i ( t ) π ( t ) ] σ H ( t ) = θ C e ( t ) θ 1 E [ μ C ( t ) ] [ g H ( t ) H ( t ) σ H ( t ) + g A ( t ) A A ( t ) ρ ( t ) σ A D ( t ) ]
Substitute (A31) into (10):
r A D e ( t ) [ i ( t ) π ( t ) ] σ A D ( t ) = θ C e ( t ) θ 1 E [ μ C ( t ) ] [ g H ( t ) H ( t ) ρ ( t ) σ H ( t ) + g A ( t ) A A ( t ) σ A D ( t ) ]
Similarly, let θ C e ( t ) θ 1 E [ μ C ( t ) ] g H ( t ) H ( t ) and θ C e ( t ) θ 1 E [ μ C ( t ) ] g A ( t ) A A ( t ) be the two unknowns of the binary Equations (A41) and (A42). Then,
k * ( t ) = g A ( t ) A A ( t ) g H ( t ) H ( t ) = r A D e ( t ) [ i ( t ) π ( t ) ] σ H 2 ( t ) + r H T e ( t ) [ i ( t ) π ( t ) ] σ A H ( t ) r H T e ( t ) [ i ( t ) π ( t ) ] σ A D 2 ( t ) + r A D e ( t ) [ i ( t ) π ( t ) ] σ A H ( t )
where σ A H ( t ) = ρ ( t ) σ H ( t ) σ A D ( t ) is the covariance of the two returns.
Equations (A40) and (A43) are the same and they are Equation (22) in the main text.

Appendix A.6. Derivation of the Optimal Ratio of Risky Financial Assets to Housing Assets from the Mean-Variance Model

The mean-variance model is the predecessor of the standard CAPM, and the risk-return relationship is clearly described by the mean-variance model when all assets are categorized into three kinds: two are risky assets, and one is risk-free assets. Simply assume that there are two risky assets A and B, and one risk-free asset, RF, in the market. The means and variances (expected returns and risks) of these three assets are:
A : E [ r A ] = r A e and V a r [ r A ] = σ A 2 B : E [ r B ] = r B e and V a r [ r B ] = σ B 2 RF : E [ r R F ] = r R F and V a r [ r R F ] = 0
Investors will invest in these three assets; thus, the shares of money invested in A, B and RF are d 1 , d 2 and d 3 . respectively The mean and variance of the portfolio are:
r T e = d 1 r A e + d 2 r B e + d 3 r R F
σ T 2 = d 1 2 σ A 2 + d 2 2 σ B 2 + 2 d 1 d 2 σ A B
(where d 1 + d 2 + d 3 = 1 ; σ A B is the covariance of r A and r B ).
Given a variety of combinations ( d i ), investors’ portfolio decisions will be more efficient when the expected total return ( r T e ) is the maximum from all plots that have the same total risk ( σ T 2 ). Thus, maximize (A44) subject to (A45), and d 1 + d 2 + d 3 = 1 .
Lagrangian:
L = d 1 r A e + d 2 r B e + ( 1 d 1 d 2 ) r R F + λ [ d 1 2 σ A 2 + d 2 2 σ B 2 + 2 d 1 d 2 σ A B σ T 2 ]
First order conditions:
L d 1 = r A e r R F + 2 λ d 1 σ A 2 + 2 λ d 2 σ A B = 0
L d 2 = r B e r R F + 2 λ d 2 σ A 2 + 2 λ d 1 σ A B = 0
Combine (A46) and (A47):
d 1 = [ r A e r R F ] σ B 2 + [ r B e r R F ] σ A B 2 λ ( σ A 2 σ B 2 σ A B 2 )
d 2 = [ r B e r R F ] σ A 2 + [ r A e r R F ] σ A B 2 λ ( σ A 2 σ B 2 σ A B 2 )
To remove λ ,
d 1 d 2 = [ r A e r R F ] σ B 2 [ r B e r R F ] σ A B [ r B e r R F ] σ A 2 [ r A e r R F ] σ A B
(A48) is the key equation of the mean-variance model, which implies that investors’ portfolio decisions will be efficient when the ratio of the shares of the money invested in risky assets A and B is subject to (A48).
When the general indicators r A e , r B e , r R F , σ A 2 , σ B 2 and σ A B are assumed to be r A D e , r H T e , i π , σ A D 2 , σ H 2 and σ A H , the optimal portfolio between housing and risky financial assets for each period of time (t) will be:
d 1 d 2 = g A A A g H H = k * = r A D e [ i π ] σ H 2 + r H T e [ i π ] σ A H r H T e [ i π ] σ A D 2 + r A D e [ i π ] σ A H
Equations (A49) is the same as Equation (22).

Appendix B. Values Used in the Simulations

Table A1. Values Used in the Simulations.
Table A1. Values Used in the Simulations.
ParameterInterpretationValueRationale
g_initialInitial real house price40.74national level in 2010
ry_initialInitial real income152.03national level in 2010
HEtoC H E / C 0.203national level in 2010
rry_consGrowth rate of real income0.0776Sample average
std_rry_consStandard deviation of rry0.0807Sample average
interest_consNominal interest rate0.053Sample average (Table 1)
inf_consInflation rate0.0248Sample average
dep_consHousing depreciation rate0.01See text
theta_consCRRA parameter5See text
ra_consSSEC growth rate0.0577sample average
std_ra_consStandard deviation of ra0.0655sample average
Further details and the complete replication of our results are available via our MATLAB simulation codes, which can be provided on request.

References

  1. Wang, Y.; Liu, J.; Tang, Y.; Sriboonchitta, S. Housing Risk and Its Influence on House Price: An Expected Utility Approach. Math. Probl. Eng. 2020, 2020, 16. [Google Scholar] [CrossRef]
  2. Muellbauer, J.; Murphy, A. Booms and Busts in the UK Housing Market. Econ. J. 1997, 107, 1701–1727. [Google Scholar] [CrossRef]
  3. Green, R. Stock prices and house prices in California: New evidence of a wealth effect? Reg. Sci. Urban Econ. 2002, 32, 775–783. [Google Scholar] [CrossRef]
  4. Lee, M.; Lee, C.; Lee, M.; Liao, C. Price linkages between Australian housing and stock markets. Int. J. Hous. Mark. Anal. 2017, 10, 305–323. [Google Scholar] [CrossRef]
  5. Kakes, J.; Van Den End, J. Do stock prices affect house prices? Evidence for the Netherlands. Appl. Econ. Lett. 2004, 11, 741–744. [Google Scholar] [CrossRef]
  6. Takala, K.; Pere, P. Testing the cointegration of house and stock prices in Finland. Finn. Econ. Pap. 1991, 4, 33–51. [Google Scholar]
  7. Ibrahim, M. House price-stock price relations in Thailand: An empirical analysis. Int. J. Hous. Mark. Anal. 2010, 3, 69–82. [Google Scholar] [CrossRef]
  8. Meen, G. Modelling Spatial Housing Markets: Theory, Analysis, and Policy, 2nd ed.; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2001; pp. 42–76. [Google Scholar]
  9. Dougherty, A.; Van Order, R. Inflation, Housing Costs, and the Consumer Price Index. Am. Econ. Rev. 1982, 72, 154–164. [Google Scholar] [CrossRef]
  10. Case, K.; Cotter, J.; Gabriel, S. Housing Risk and Return: Evidence from a Housing Asset-Pricing Model. J. Portf. Manag. 2011, 37, 89–109. [Google Scholar] [CrossRef] [Green Version]
  11. Hu, J.; Sui, Y.; Ma, F. The Measurement Method of Investor Sentiment and Its Relationship with Stock Market. Comput. Intell. Neurosci. 2021, 2021, 1–11. [Google Scholar] [CrossRef]
  12. Wang, J.; Wang, J.; Fang, W.; Niu, H. Financial time series prediction using elman recurrent random neural networks. Comput. Intell. Neurosci. 2016, 105, 214–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Lucas, R., Jr. Asset Pricing in an Exchange Economy. Econometrica 1978, 46, 1429–1445. [Google Scholar] [CrossRef]
  14. Grossman, S.; Laroque, G. Asset Pricing and Optimal Portfolio Choice in the Presence of Illiquid Durable Consumption Goods. Econometrica 1990, 58, 25–51. [Google Scholar] [CrossRef]
  15. Flavin, M.; Nakagawa, S. A Model of Housing in the Presence of Adjustment Costs: A Structural Interpretation of Habit Persistence. Am. Econ. Rev. 2008, 98, 474–495. [Google Scholar] [CrossRef] [Green Version]
  16. Piazzesi, M.; Schneider, M.; Tuzel, S. Housing, Consumption and Asset Pricing. J. Financ. Econ. 2007, 83, 531–569. [Google Scholar] [CrossRef] [Green Version]
  17. Chen, N. Asset price fluctuations in Taiwan: Evidence from stock and real estate prices 1973 to 1992. J. Asian Econ. 2001, 12, 215–232. [Google Scholar] [CrossRef]
  18. Aliev, R.; Mraiziq, D.; Huseynov, O. Expected Utility Based Decision Making under Z-Information and Its Application. Comput. Intell. Neurosci. 2015, 2015. [Google Scholar] [CrossRef] [Green Version]
  19. Jiang, Y.; Wang, Y. Price dynamics of China’s housing market and government intervention. Appl. Econ. 2021, 53, 1212–1224. [Google Scholar] [CrossRef]
  20. Zhi, T.; Li, Z.; Jiang, Z.; Wei, L.; Sornette, D. Is there a housing bubble in China? Emerg. Mark. Rev. 2019, 39, 120–132. [Google Scholar] [CrossRef] [Green Version]
  21. Poterba, J. Tax Subsidies to Owner-occupied Housing: An Asset-market Approach. Q. J. Econ. 1984, 99, 729–752. [Google Scholar] [CrossRef] [Green Version]
  22. Poterba, J. House Price Dynamics: The Role of Tax Policy and Demography. Brookings Pap. Econ. Act. 1991, 1991, 143–203. [Google Scholar] [CrossRef]
  23. Andrew, M.; Meen, G. House price appreciation, transactions and structural change in the British housing market: A macroeconomic perspective. Real Estate Econ. 2003, 31, 99–116. [Google Scholar] [CrossRef]
  24. Yamashita, T. House price appreciation, liquidity constraints, and second mortgages. J. Urban Econ. 2007, 62, 424–440. [Google Scholar] [CrossRef]
  25. Meen, G. The removal of mortgage market constraints and the implications for econometric modelling of UK house prices. Oxf. Bull. Econ. Stat. 1990, 52, 1–23. [Google Scholar] [CrossRef]
  26. Dynan, K. How Prudent Are Consumers? J. Political Econ. 1993, 101, 1104–1113. [Google Scholar] [CrossRef]
  27. Sharpe, W. Mutual Fund Performance. J. Bus. 1966, 39, 119–138. [Google Scholar] [CrossRef]
  28. Markowitz, H. Portfolio Selection. J. Financ. 1952, 7, 77–91. [Google Scholar] [CrossRef]
  29. Babcock, B.; Choi, E.; Feinerman, E. Risk and probability premiums for CARA utility functions. J. Agric. Resour. Econ. 1993, 18, 17–24. [Google Scholar]
Figure 1. Correlation between housing capital return and SSEC return.
Figure 1. Correlation between housing capital return and SSEC return.
Axioms 11 00145 g001
Figure 2. Average dynamics under 100 replications (random seed: 1–100).
Figure 2. Average dynamics under 100 replications (random seed: 1–100).
Axioms 11 00145 g002
Table 1. Data Description.
Table 1. Data Description.
NameVariableLabel/MeasurementMeanS.D.
Monthly:
M H E monthly H E 100 million yuan
M H S monthly H S 10 thousand m 2
M H P monthly H P = H E / H S × 10 , 000
M P A SSEC
m r h M H P monthly growth m r h m = l n ( M H P m ) l n ( M H P m 1 )
m r a M P A monthly growth m r a m = l n ( M P A m ) l n ( M P A m 1 )
Annual:
H E housing expenditure100 million yuan1592.12194.1
H S housing space sold10 thousand m 2 2862.42765.5
H P house price (yuan/m 2 ) = H E / H S × 10 , 000 4566.03730.9
Ydisposable incomeyuan/person18,218.810,332.2
C P I consumer price index2000 = 100123.117.72
C O N final consumption100 million yuan6136.56342.3
Greal house price G = H P / C P I 35.95327.701
l n g logarithm form l n g = l n ( G )
R Y real income R Y = Y / C P I 142.02868.680
l n r y logarithm form l n r y = l n ( R Y )
i(loan) interest rate0.0530.0056
r h nominal return r h 1 = r H + π 0.0950.093
t a u 1 τ 1 Equation (20)−0.00360.084
t a u 2 τ 2 , C R R A Equation (21)0.0350.100
r h e expected housing yield r h e t = ( 1 / 12 ) m = J a n D e c m r h t , m × 12
r a e expected financial yield r a e t = ( 1 / 12 ) m = J a n D e c m r a t , m × 12
S D H S.D. of r H σ H , t = ( 1 / 12 ) m = J a n D e c ( m r h t , m r h e t ) 2
S D A S.D. of r A D σ A D , t = ( 1 / 12 ) m = J a n D e c ( m r a t , m r a e t ) 2
c o r r correlation ρ t = σ [ H , A D ] , t / ( σ H , t σ A D , t )
Notes: (i) The time subscripts t and m are year and month, respectively; (ii) SSEC = Shanghai Securities Composite Index; (iii) “×12” in rhe and rae indicates annualization; (iv) σ[H,AD] is the covariance between mrh and mra.
Table 2. Panel Data Unit-Root/Stationary Test (p-value).
Table 2. Panel Data Unit-Root/Stationary Test (p-value).
Test tau 1 tau 2
Levin-Lin-Chu (LLC) test<0.001<0.001
Im-Pesaran-Shin (IPS) test<0.001<0.001
ADF-Fisher unit-root test (lag = 2)<0.001<0.001
Hadri LM stationary test0.6800.726
Table 3. Non-stationary Heterogeneous Panel Model, Decomposed UCC, 2002–2018.
Table 3. Non-stationary Heterogeneous Panel Model, Decomposed UCC, 2002–2018.
CoefficientEquation (26), DFEEquation (27), MGEquation (28), PMG
β (or E [ β i ] )0.772 *** (4.4)0.785 *** (3.3)1.007 *** (5.3)
η (or E [ η i ] )−0.383 *** (10.8)−0.644 *** (10.0)−0.499 *** (9.7)
a 1 (or E [ a i 1 ] )0.830 *** (25.7)0.906 *** (7.4)0.853 *** (55.5)
a 2 (or E [ a i 2 ] )−7.401 *** (3.9)−4.173 (0.6)−5.806 *** (6.8)
a 3 (or E [ a i 3 ] )0.381 *** (3.3)0.607 *** (4.0)0.168 ** (2.4)
a 4 (or E [ a i 4 ] )−0.367 ** (2.1)−0.987 * (1.7)−0.319 *** (3.3)
a 5 (or E [ a i 5 ] )−0.303 ** (2.1)−0.982 (0.5)0.077 (0.9)
constant−0.079 (0.9)−0.283 * (1.7)−0.224 *** (5.0)
Notes: |t| in parentheses; *, **, and *** are significant at the 10%, 5% and 1% levels.
Table 4. Non-stationary Heterogeneous Panel Model, ln(UCC), 2002–2018.
Table 4. Non-stationary Heterogeneous Panel Model, ln(UCC), 2002–2018.
CoefficientEquation (31), UCC by (29)Equation (31), UCC by (30)
β 0.798 *** (4.5)0.815 *** (4.7)
η −0.396 *** (11.0)−0.380 *** (10.7)
a 1 0.812 *** (24.1)0.807 *** (23.2)
a 2 −0.345 *** (3.6)−0.410 *** (4.5)
constant−0.583 *** (4.3)−0.619 *** (5.2)
R 2 (within)0.2960.315
Notes: |t| in parentheses; *** is significant at the 1% levels.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Liu, J.; Qiu, Z.; Sriboonchitta, S. The Effect of Financial Market Factors on House Prices: An Expected Utility Three-Asset Approach. Axioms 2022, 11, 145. https://doi.org/10.3390/axioms11040145

AMA Style

Wang Y, Liu J, Qiu Z, Sriboonchitta S. The Effect of Financial Market Factors on House Prices: An Expected Utility Three-Asset Approach. Axioms. 2022; 11(4):145. https://doi.org/10.3390/axioms11040145

Chicago/Turabian Style

Wang, Yehui, Jianxu Liu, Zhaolin Qiu, and Songsak Sriboonchitta. 2022. "The Effect of Financial Market Factors on House Prices: An Expected Utility Three-Asset Approach" Axioms 11, no. 4: 145. https://doi.org/10.3390/axioms11040145

APA Style

Wang, Y., Liu, J., Qiu, Z., & Sriboonchitta, S. (2022). The Effect of Financial Market Factors on House Prices: An Expected Utility Three-Asset Approach. Axioms, 11(4), 145. https://doi.org/10.3390/axioms11040145

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop