Predicting Stock Return with Economic Constraint: Can Interquartile Range Truncate the Outliers?

We find that imposing economic constraint on stock return forecasts based on the Interquartile Range of equity premium can significantly strengthen predictive performance. Specifically, we construct a judgment mechanism that truncates the outliers in forecasts of stock return. We prove that our constraint approach can realize more accurate predictive information relative to the unconstraint approach from the perspective of statistics and economics. In addition, the new constraint approach can effectively defeat CT constraint and CDA strategy. -e three mixed models we proposed can further enhance the accuracy of prediction, especially the mixed model combined with our constraint approach. Finally, utilizing our new constraint approach can help investors obtain considerable economic gains. With the application of extension and robustness analysis, our results are robust.


Introduction
Stock is a major financial product in the financial market. e uncertainty of stock market fluctuation not only involves numerous economic behaviors but also contains many political factors. erefore, researching the prediction of stock return has always been a considerable subject for many scholars and practitioners. However, predicting stock returns is notoriously difficult. e predictability of stock returns in sample has been confirmed by many studies (Campbell and Shiller [1]; Pontiff and Schall [2]; Dai et al. [3,4]; Ferson and Harvey [5]). In addition, according to Welch and Goyal [6], in terms of predicting out-of-sample stock returns, it is hard to find a certain model to go beyond the benchmark prediction model that treats stock returns as a constant. For this reason, a large number of scholars have proposed many predictors linked to macrofundamentals, which can enhance the predictability of out-of-sample stock returns, for example, oil-related variables (Liu et al. [7]), the variance risk premium (Bollerslev et al. [8,9]), technical indicators (Neely et al. [10]), economic policy uncertainty (Brogaard and Detzel [11]), manager sentiment (Jiang et al. [12]), short interest (Rapach et al. [13]), news-implied volatility (Manela and Moreira [14]).
Numerous recent studies have verified that the singlepredictor model fails to beat the benchmark prediction model in the out-of-sample prediction of stock returns. e research of Pesaran and Timmermann [15] has shown that the uncertainty of parameters brings numerous issues to the prediction of stock returns, particularly, when the data contain a large number of noises. While more and more pieces of literature have found that economic constraints can improve the prediction of stock returns, there is no common view on the implementation of constraints. For example, to identify the parameter of the term structure model, Ang and Piazzesi [16] impose no arbitrage restrictions. Campbell and ompson [17] require the negative predictive forecasts to become nonnegative and bind the sign of the slope coefficient in the predictive model of stock returns. Pástor and Stambaugh [18,19] apply a priori information to ensure a negative correlation between expected returns and shocks of unexpected returns. Yi et al. [20] propose the cycledecomposing approach (CDA hereafter) that imposes restrictions on parameters.
In this paper, we establish a new constraint approach to truncate stock return forecasts. We mainly report two types of constraints for comparison. Firstly, Campbell and ompson [17] believe that when the stock return forecasts are negative, it should be truncated to be a nonnegative value. is is because it is unrealistic for a risk-averse investor to hold stocks when their expected returns are negative. Secondly, Yi et al. [20] propose the cycledecomposing approach to decompose economic variables into short-term bias and long-term trend.
is approach achieves the limitation of model parameters. e idea of this paper is designed to truncate prediction outliers generated by univariate predictive regression model. To realize this idea, we advance an approach that is related to the principle of boxplot. is approach needs to apply the concept of statistical quartile and uses the quartile difference (that is, the difference between the third quartile and the first quartile, also known as Interquartile Range) to deal with abnormal information. erefore, we name this approach the Interquartile Range constraint (IQR, hereafter). e operation mechanism of the IQR constraint is setting up a reasonable interval for the stock return forecasts by utilizing the principle of boxplot when real stock return forecasts beyond the upper and lower range of this interval can be considered abnormal information.
Our main motivation for constructing this constraint approach is linked to but quite different from the CT constraint. CTconstraint believes that investors will not hold stock when the expected returns are negative, but our approach finds that negative expected returns are still profitable. is is because investors are just unwilling to believe in extremely negative forecasts rather than being afraid of all negative forecasts. With this line, investors will not believe in extremely positive forecasts, too. Another motivation is that although it is unlikely to generate extreme return forecasts in the future, investors will pay special attention to the extreme changes of stock return prediction information. erefore, they are more likely to react to these pieces of extreme information. In addition, there is always a positive relationship between return and risk; that is, extreme return forecasts are often accompanied by high risk. Consequently, investors will be more cautious about their portfolio allocation, and stock prices will respond to the behavior of investor. Collectively, the accuracy of stock return forecasts can be further improved after the IQR constraint model identifies return forecast outliers and rules out them. On the other hand, our motivation is partly inspired by Pettenuzzo et al. [21]. To make the portfolio performance more reasonable, they apply the Bayesian approach to restrict the Sharpe ratio to lie an upper and lower bound. e difference is that we use box graph instead of the Bayesian approach to determine the upper and lower bound.
For empirical analysis of our constraint approach, we apply 14 macroeconomic variables referred by Welch and Goyal [6] as the predictor of the single variable prediction regression model. Following Rapach and Zhou [22] and Dai and Zhu [23], we use the data from January 1927 to December 1947 as the initial period and employ January 1948 as the starting point for the recursive prediction of the expected return of the S&P 500 index out of sample. We employ out-of-sample R 2 statistic to describe the forecasting performance of stock returns and examine the null hypothesis that MSPE b ≤ MSPE i with Clark and West [24] statistics. is test is often used to check whether R 2 OS is statistically significant. e empirical results show that only the single-predictor model of SVAR and TBL can generate positive R 2 OS value, which is consistent with the previous literature results; that is, the predictive performance of macrofundamentals fails to exceed the historical average benchmark. In addition, compared with the original model, CT constraint and the CDA strategy can improve the prediction performance of 9 and 10 variables, respectively. Fortunately, we find that our constraint approach can generate larger R 2 os for all predictors relative to the original model. We also compute the average R 2 os value of each model as a comprehensive evaluation indicator for the corresponding predictive approach. e results reveal that the original model generates the average R 2 os value at −0.502%, the CT constraint increases it to −0.272%, the CDA strategy increases it to 0.094%, and notably, the IQR constraint increases it to 0.261%. Consequently, we conclude that the IQR constraint model can effectively improve the out-of-sample prediction ability compared with the original model, CT constraint, and the CDA strategy.
Furthermore, we further evaluate the forecasting performance of the model of interest by computing certainty equivalent return (CER) from an economics perspective. When the relative risk aversion coefficient is 3, we calculate the difference for the certainty equivalent return (CER) obtained from the model of interest and the benchmark prediction model. Here, we find that the constraint model can lead to higher CER gains than the unconstraint model among most predictors. In particular, when investors use our constraint approach, the higher CER gains can be yielded compared to the unconstraint model among all predictors, which means that eliminating outliers in stock return forecasts can help investors to allocate their assets and create more substantial economic gains.
Following Dai and Zhu [25], we learned that the mixed model can also significantly strengthen the predictive performance. erefore, based on the model of interest in this paper, we construct three mixed models, namely, CDA&CT, CT&IQR, and CDA& IQR. e first model is to impose CT constraint on the prediction information of CDA strategy, the second model is to impose IQR constraint on the prediction information of CT constraint, and the last model is to impose IQR constraint on the prediction information of CDA strategy. e results show that the mixed model can yield more accurate prediction information and better economic gains than the unconstrained model and the competition model in this paper. Particularly, the mixed model connected with our constraint model shows the best prediction ability.
In addition, in terms of the business cycle, different risk aversion coefficients, transaction cost equal to 50, alternative weight ranges, and alternative out-of-sample periods, our results with the IQR constraint are all robust. In a word, our constraint model shows a stronger out-of-sample forecasting ability and creates considerable economic gains. e arrangement of this paper is as follows: Section 2 describes the empirical data with summary statistics. Section 3 presents the constraint model and forecast strategy. Section 4 reports the empirical results for out-of-sample forecasting performance based on unconstraint model and constraint models. Section 5 shows some extensions and robustness test. Section 6 makes the conclusion.

Empirical Data
To predict out-of-sample stock returns, we employ 14 popular macroeconomic variables (we thank Amit Goyal for providing the latest data, which can be obtained from http://www.hec.unil.ch/agoyal/), which are mentioned in Welch and Goyal [6] who will update the data for the latest year. ese variables are often used in many recent pieces of literature to predict stock returns [26][27][28][29]. We also use Amit Goyal's dataset to obtain the excess stock returns (we use the monthly return of S&P 500 index minus the short T-bill rate). Next, we will give the description of 14 predictors.
(i) Log dividend-price ratio (DP): log of a 12-month moving sum of dividends paid on the S&P 500 index minus the log of stock prices (S&P 500 index) (ii) Log dividend yield (DY): log of a 12-month moving sum of dividends minus the log of lagged stock prices (iii) Log earnings-price ratio (EP): log of a 12-month moving sum of earnings on the S&P 500 index minus the log of stock prices (iv) Log dividend-payout ratio (DE): log of a 12-month moving sum of dividends minus the log of a 12month moving sum of earnings (v) Stock return variance (SVAR): sum of squared daily returns on the S&P 500 index (vi) Book-to-market ratio (BM): book-to-market value ratio for the Dow Jones Industrial Average (vii) Net equity expansion (NTIS): ratio of a 12-month moving sum of net equity issues by NYSE-listed stocks to the total end-of- e summary statistics for 14 predictors and excess stock return from January 1927 through December 2018 are given in Table 1. We compute that the monthly average excess stock return can achieve 0.498%, which is a considerable return. is is also an important reason why we are committed to studying the stock return forecasts. All data stem from Amit Goyal's dataset.

Forecasting Regression Model.
According to the relevant literature, we utilize the following individual variable forecasting regression model to predict excess stock returns: where r t+1 represents the excess stock return at month t + 1 over risk-free bill at month t, x i,t denotes the available variable for the ith predictor at month t, and ε i,t+1 is a disturbance term which is assumed to be independent and identically distributed, and the mean equals zero. Based on the literature of Devpura et al. [30], we obtain out-of-sample excess stock return forecasts by applying a recursive estimation window. Particularly, the whole sample consisting of T observations is divided into the first M observations and T-M observations as the inner part of the sample and the outer part of the sample, respectively. Next, the following procedure is utilized to obtain out-of-sample forecast: where α i,t and β i,t are the ordinary least squares (OLS) estimates of regression coefficients αand β and each estimated result is generated by regressing r t m t�2 on a constant and x i,t m−1 t�1 . Continuing this step to make recursive prediction, we can obtain the out-of-sample forecast sequence, r i,t T t�m+1 , a total of T-M. Following the paper of Goyal and Welch [31] and Hammerschmid and Lohre [32], after implementing the predictive regression model, we further apply the historical average benchmark model to evaluate the prediction accuracy of the original model and constraint model. e benchmark model actually assumes that β i in formula (1) is zero and is given as follows:

Constraint Models.
To compare our constraint model, we also consider two prevailing constraint models as competition model. e first one is proposed by Campbell and ompson [17] who believe that it is difficult for investors to hold stocks when the forecast signal is negative (that is, the stock excess return forecast is negative), thus constraining the negative return forecasts to zero. e specific model is given as follows: e second one is put forward by Yi et al. [20]. ey decompose the predictor into the long-term mean part and the short-term deviation part and then carry out the least square regression for the decomposed content. Finally, compute stock return forecast based on the corresponding return forecast of the two parts by utilizing the simple average combination strategy. e specific operation process is given as follows: where C is the lag order of the long-cycle term and equals 12.
x i,τ denotes the available variable for the ith predictor.
x C i,t,mean denotes the long-term mean of the ith at month t. x C i,t,dev denotes the short-term deviation of the ith at month t.
where α i,t,mean , β i,t,mean , α i,t,dev , and β i,t,dev are the ordinary least squares (OLS) estimates of regression coefficients αand β and generated by regressing r t m t�2 on a constant and is the ultimately out-of-sample predictive content we need.
We advance a new constraint model which applies the concept of statistical quartile and uses the quartile difference (that is, the difference between the third quartile and the first quartile, also known as Interquartile Range) to deal with abnormal forecasts. erefore, we name this approach the Interquartile Range constraint (IQR, hereafter). More specifically, we establish a judgment mechanism for the new constraint model, which sets an interval for the original prediction value without outliers. We give its mathematical form as follows: where r IQR i,t+1 is the excess stock return forecast of our new constraint model for month t + 1 and ith predictor. IQR t denotes the Interquartile Range of actual stock return at month t. r i,t+1 represents the excess stock return forecast at month t + 1 based on regression in (2). Notably, we choose three times the Interquartile Range as the criterion for truncating outliers to acquire economic constraint more

Empirical Results for Equity Premium Forecasts
In this part, we evaluate the forecasting performance of the predictive model from an economic and statistic point of view. e out-of-sample period is 1948 : 01-2018 : 12.

Assessing the Out-of-Sample Statistics Performance.
We employ a widely accepted approach to evaluate the outof-sample forecasting performance of constraint and unconstraint models. is approach generally compares the forecasting models with the historical average benchmark model to weigh the accuracy of return forecasts (see, e.g., [32][33][34][35][36][37]). e method is R 2 (R 2 OS ) statistic which was proposed by Campbell and ompson [17]. R 2 denotes as the percentage reduction for the mean squared prediction error of the forecasting model relative to the mean square prediction error of the benchmark model. Next, we calculate the out-of-sample R 2 (R 2 OS ) statistic.
As can be seen from the above formula, R 2 OS is a positive value, which means that the model of interest will produce more accurate forecasting information than the benchmark model. Following previous literature of Buncic and Tischhauser [38], we further apply the Clark and West [24] test (CW hereafter) to examine whether the forecasting regression model has a lower MSPE relative to the benchmark model. e null hypothesis is that MSPE b ≤ MSPE i and the alternative hypothesis is that MSPE b > MSPE i . e CW statistic is a modification of Diebold and Mariano [39] statistic, which is convenient to evaluate statistical significance. To calculate the CW statistic, we first define the MSPE-adjusted statistic: Next, we regress f i,t T t�m to a constant so that we can get the t-statistics of the constant. en, the standard normal distribution is used to calculate the p value of one side (upper tail) test. e p value represents the statistical significance that the model of interest has a smaller mean square prediction error compared with the benchmark model.
is table shows the out-of-sample statistical performance of the original (no constraint) model and constraint models based on 14 macroeconomic predictors. We evaluate the accuracy of out-of-sample stock return forecast by employing R 2 OS . R 2 OS computes the reduction in MSPE (mean squared prediction error) for the predictive regression model relative to the benchmark model. e statistical significance of R 2 OS is examined by employing Clark and West [24] test. * * * , * * , and * mean significance at the 1%, 5%, and 10% levels, respectively. Table 2 reports the out-of-sample statistical performance with regard to the original model and constraint models. We apply R 2 OS and CW statistics to evaluate statistical performance and statistical significance, respectively. Our empirical results reveal that only SVAR, TBL, and TMS can yield positive R 2 OS based on univariate predictive regression, which once again validates the research of Welch and Goyal [6] who believed that it is difficult for individual predictor to obtain better predictive performance relative to the historical average.
Comparing the CT constraint and the CDA strategy to the original model, the results display that there are 10 and 12 variables that produce larger R 2 OS , respectively. On the other hand, 6 and 7 out of them are positive, respectively. At the same time, we find that the predictive performance of all individual predictors can effectively strengthen after implementing IQR constraint to deal with the outliers in the original model. What is more, our new constrain model can yield 11 positive R 2 OS relative to the benchmark model. When we apply CW test, DP, DY, EP, NTIS, TBL, LTY, and TMS significantly outperform the historical average model. More specifically, based on the whole sample, the new constraint model we submitted can promote R 2 OS for DP from −0.074% to 0.653%, promote it for DY from −0.389% to 0.252%, promote it for SVAR from 0.154% to 0.836%, promote it for NTIS from −0.633% to 0.313%, promote it for TBL from 0.067% to 0.844%, promote it for TMS from 0.103% to 0.994%, and promote it for INFL from −0.060% to 0.727%. We also calculate the average R 2 OS to evaluate the comprehensive predictive ability of the corresponding predictive model and find that the CT constraint and CDA strategy can promote the average R 2 OS from −0.502% to −0.272% and 0.094%, respectively, while our IQR constraint model generates the largest average of R 2 OS at 0.261%. Consequently, we believe that the IQR constraint model can effectively improve the out-of-sample predictive accuracy compared with the original model, CT constraint, and CDA strategy.

Assessing the Out-of-Sample Economic Performance.
We have testified that the constraint approach significantly outperforms the benchmark model in the prediction of stock returns from the point of view of statistics. To describe stock return forecasts with utility-related indicator from an economic perspective (Campbell and ompson [17]; Dai and Kang [40]), we further evaluate the forecasting performance of stock return by calculating the certainty equivalent return (CER). We assume that a mean-variance investor allocates their assets among the S&P 500 index and the risk-free property. To maximize the CER, we set the optimal weight of the stock in the portfolio during t + 1 as follows:

Mathematical Problems in Engineering
where c represents the risk aversion coefficient of investor, r i,t+1 is the stock return forecast based on a different predictive model, and σ 2 t+1 is the variance prediction of stock return time series during t + 1. For risk aversion coefficient, c � 3, as in Rapach et al. [26]. We followed Campbell and ompson [17] with a five-year rolling window to compute variance forecasts. Meanwhile, we assign 1 − ω * i,t+1 to riskfree bills. In addition, following the information in the relevant literature (e.g., Huang et al. [34]; Pan et al. [23]), in order to prevent more than 50% of short selling and financial leverage, ω * i,t is limited to −0.5 to 1.5. e portfolio returns of investor in stocks and risk-free bills are defined as follows: where R f,t+1 represents risk-free rate. Finally, we employ the following formula to solve CER: where μ i and σ 2 i represent the mean and variance of the portfolio of stocks and risk-free bills, respectively. CER gains are the metric of economic performance that equals portfolio management fees that mean-variance investors are willing to pay to obtain predictive formation relative to the historical average benchmark. Accordingly, the result of CER gains is the difference between the CER of the model of interest and the CER of the benchmark model.
is table shows the out-of-sample portfolio performance of the original (no constraint) model and constraint models based on 14 macroeconomic predictors. We employ the certainty equivalent return (CER) gains to remark the performance of the portfolio. CER is a utility-related indicator which is derived from investors investing their assets between stock and T-bill rate. e CER gains indicate the difference between the forecasting model we applied and the benchmark forecast model. We use a risk aversion coefficient equal to 3 and limit the weight to −0.5 to 1.5 to calculate the CER. Furthermore, the CER gains should be multiplied by 1200 as an annualized percentage value. Table 3 reports the economic performance of the stock return forecasts assessed in terms of CER gains of an investor who applies the given constraint model and original model relative to the benchmark forecasting model. In this table, the CER gains are multiplied by 1200 to represent the annualized percent values. Comparing the original model with constraint models based on economic performance, the result reveals that the CT constraint model fails to significantly strengthen CER gains. In addition, the CDA strategy can effectively promote CER gains. Delightfully, the CER gains are further promoted compared to the original model and CT constraint model when investors use our constraint model. It is noteworthy that all predictors improve their economic performance after implementing the IQR constraint model. What is more, 10 of 14 the CER gains are greater than zero, which means that our constraint model can strongly beat the benchmark model in terms of economic performance. More specially, by the application of IQR constraint model, the CER gains promote it from −0.207% to 0.216% for SVAR, from 0.530% to 0.919% for NTIS, from 1.497% to 1.755% for TBL, from 0.622% to 0.864% for LTY, from 0.591% to 0878% for LTR, and from 1.547% to 1.967% for TMS. In a word, implementing the IQR constraint approach can help investors gain more economic gains.

Mixing Existing Model.
ere are many studies that have shown that combination forecasting can effectively promote the prediction accuracy of stock returns. e combination approach includes two modes: one is to combine the prediction variables (see, e.g., Kelly and Pruitt [41]; Bahrami et al. [42]; Pettenuzzo et al. [21]), and the other is to combine the popular models (see, e.g., [43][44][45][46][47][48]). Following the literature of Dai and Zhu [25], we employ a mixed model by mixing existing forecasting models to further study stock return forecasts. In this paper, existing forecasting models refer to the CT constraint model, CDA strategy, and IQR constraint model, respectively. With regard to the first mixed model, we consider the reorganization of CT constraints and CDA strategy, namely, the CDA&CT constraint model. e specific implementation steps are as follows: e second mixed model is to impose the IQR constraint on the stock return forecasts of the CT constraint, namely, the CT&IQR constraint model. e details are as follows: For the last mixed forecasting model, we consider applying the IQR constraint to eliminate outliers obtained by the CDA strategy, namely, the CDA&IQR constraint model.
is table shows the out-of-sample statistical performance of the mixed existing constraint models based on 14 macroeconomic predictors. We evaluate the accuracy of outof-sample stock return forecast by employing R 2 OS . R 2 OS computes the reduction in MSPE (mean squared prediction error) for the predictive regression model relative to the benchmark model. e statistical significance of R 2 OS is examined by employing Clark and West [24] test. * * * , * * , and * mean significance at the 1%, 5%, and 10% levels, respectively. Table 4 presents the statistical performance of the mixed models based on the existing model; similarly, R 2 OS is used for evaluation. e results reveal that utilizing the CDA&CT model can enhance 10 of 14 R 2 OS compared to the original model. After implementing the CT&IQR model, we find that all of the 14 macroeconomic predictors are able to significantly intensify concerning out-of-sample prediction ability. What is more, the CT&IQR constraint model can yield 12 positive R 2 OS relative to benchmark model. When we apply the CDA&IQR constraint model, all of R 2 OS are positive, which means that the CDA&IQR constraint model can significantly outperform the historical average model. Notably, the prediction performance of some predictors can be greatly improved by utilizing CDA&IQR constraint, such that R 2 OS promotes it from −1.463% to 0.839% for EP, from −0.389% to 1.004% for DY, from −0.074% to 1.011% for DP, from 0.067% to 1.066% for TBL, from −0.761% to 1.411% for LTR, and most notably, from −0.677% to 2.384% for LTY. On the other hand, the average R 2 OS for the three mixed models arrives at 0.140, 0.505, and 0.850, respectively, which is significantly higher than the original model. In addition, we find that the average R 2 OS of each mixed model is effectively higher than the corresponding existing model. e above evidence shows that the mixed model owns more accurate prediction contents than the existing forecasting model in this paper.
is table shows the out-of-sample portfolio performance of the original (no constraint) model and constraint models based on 14 macroeconomic predictors. We employ the certainty equivalent return (CER) gains to remark the performance of the portfolio. CER is a utility-related indicator which is derived from investors investing their assets between stock and T-bill rate. e CER gains indicate the difference between the forecasting model we applied and the benchmark forecast model. We use risk aversion coefficient equal to 3 and limit the weight to −0.5 to 1.5 to calculate the CER. Furthermore, the CER gains should be multiplied by 1200 as an annualized percentage value. Table 5 reports the economic performance of the mixed models, assessed in terms of CER gains. e results of the mixed model show that most of the CER gains are greater than zero, which implies that the combination of existing models can effectively exceed the benchmark forecasting model. We further compare the CER gains of the mixed model with the original model and the existing models and find that the mixed model is obviously superior to not only the original model but also the corresponding existing model. In terms of the three mixed models in this paper, CDA&IQR model has the best portfolio performance. Note that CDA&IQR model has generated 12 positive CER gains and achieved more than 1% CER gains with regard to LTR, TBL, and LTY. In particular, the CDA&IQR model promotes the CER gains for LTR from 0.591% to 1.675%, promotes it for TBL from 1.497% to 2.159%, and promotes it for LTY from 0.622% to 3.379%. Overall, the mixed model has a strong predictive utility in terms of portfolios.

Extensions and Robustness Analysis
In this section, in order to further evaluate the results of Section 4, we propose several extensions and robustness tests which include the out-of-sample performance over the business cycle, CER for alternative risk aversion coefficient,   [49] and Fama and French [50] have shown that the predictability of the stock market is affected by the business cycle. One explanation is that investors take on high risk during the recession, which makes investors desire for higher risk premium, thus creating the predictability of stock returns. On the other hand, compared with the expansion period, the recession period can obtain better prediction performance (Jiang et al. [12]; Rapach et al. [26]; Wang et al. [28]). Consequently, we next examine the statistical performance and economic gains of the IQR constraint model and the original model over the business cycle. R 2 OS based on recession and expansion periods are given as follows: where I c t (c � EXP, REC) is an indicator variable equal to one when month t is within the scope of the expansion (recession) period and zero otherwise. e expansion and recession period of the NBER-dated business cycle are derived from the FRED database.
is table shows the out-of-sample statistical performance and economic performance of the original (no constraint) model and quartile difference constraint model based on 14 macroeconomic predictors over business cycles. We evaluate the accuracy of out-of-sample stock return forecast by employing R 2 OS . * * * , * * , and * mean significance at the 1%, 5%, and 10% levels, respectively. We employ the certainty equivalent return (CER) gains to remark the performance of the portfolio. CER is a utilityrelated indicator which is derived from investors investing their assets between stock and T-bill rate. We use risk aversion coefficient equal to 3 and limit the weight to −0.5 to 1.5 to calculate the CER. Furthermore, the CER gains should be multiplied by 1200 as an annualized percentage value. Table 6 reports R 2 OS and CER gains for the original model as well as the IQR constraint model during the recession and expansion period. As can be seen from the second to fifth columns of this table, the forecasting performance in the period of recession is much better than that in the period of expansion, which is consistent with the conclusion of the previous article (Neely et al. [10]; Pan et al. [23]). Furthermore, in either the recession or expansion period, the IQR constraint model can yield larger R 2 OS than the original model in every univariate forecasting regression. Particularly, in the period of recession, we can obtain 11 positive R 2 OS after implementing the IQR constraint model. For the period of expansion, our constraint strategy leads the number of positive R 2 OS to increase from 1 to 9. In terms of economic performance, the difference between expansion and recession is similar to that of statistical performance. During the recession, almost all individual predictor forecasting models realize positive CER gains based on the IQR constraint model, and nearly all CER gains of the IQR constraint model are much larger than the original model. When considering economic performance from the perspective of the expansion period, although the economic performance of all the predictors seems to be not as good as the recession, 13 of the 14 CER gains have made progress after the implementation of the IQR constraint model. In a word, the prediction results of the expansion period and recession period show that rational investors will achieve more accurate prediction information and better economic value compared with the original model and historical average benchmark model with the application of the IQR constraint model.

Alternative Parameters for CER.
For a portfolio exercise, the optimal portfolio weight is often affected by the risk aversion coefficient, the weight range allocated to the stock investment, and the transaction cost. Consequently, we employ some alternative parameters to create portfolio performance for stock return prediction. Specifically, we apply c � 2, 4, and 6 instead of c � 3 to verify the stability of portfolio return. In addition, we consider limiting the weight range to 0 to 1.5 instead of −0.5 to 1.5 as previously used. Furthermore, we followed Neely et al. [10] who consider a proportional transaction cost; that is, the value of each transaction is 50 basis points rather than 0 bases.
is table shows the out-of-sample portfolio performance of the original (no constraint) model and constraint models based on 14 macroeconomic predictors. We employ the certainty equivalent return (CER) gains to remark the performance of the portfolio. CER is a utility-related indicator which is derived from investors investing their assets between stock and T-bill rate. In this table, we use risk aversion coefficient (i.e., c � 2, 4, 6), the weight range (i.e., ω � [0, 1.5]) of substitution, and the alternative transaction cost (i.e., TC � 50) to calculate the CER. e specific  Table 7, we report the main results for portfolio performance linked to the alternative parameters. It is worth noting that we need to keep the remaining parameters consistent with the previously mentioned parameters when one of the parameters changes. From the first column to the sixth column of this table, we can see that although the CER gains of the original model and the IQR constraint model have obvious differences under different risk aversion coefficients, the economic performance of all of the predictors based on the IQR constraint model significantly outperforms the original model. In addition, when the risk aversion coefficient equals 2, 4, and 6, we can obtain 10, 9, and 9 positive CER gains after the application of the IQR constraint model, respectively, which is consistent with the results of Table 2. Table 7 also suggests that the economic performance of all predictors related to our constraint model has improved after limiting the weight range from 0 to 1.5, and 10 positive CER gains have been generated. Finally, when we consider the transaction cost (i.e., TC � 50), all CER gains fail to beat the ones compared with nontransaction costs. is is because of the stability of the historical average; however, economic performance is inversely proportional to the historical average. Inevitably, our forecasts increase the frequency of transactions, which in turn increases transaction costs. Nevertheless, we can still improve the CER gains of all predictors with the implementation of the IQR constraint model relative to the original model.
According to the above analyses, the economic gains are robust based on the IQR constraint model after applying the alternative parameters to CER.

Multivariate Information Analysis.
Univariate forecasting model usually has model uncertainty, even if the economic constraint model may also have model uncertainty. To solve this issue, the existing pieces of literature have developed strand multivariate combination approaches to predict stock returns (see, e.g., Neely et al. [10], Dai and Zhu [25], and Rapach et al. [26]). e combination approach that we mainly use is proposed by Rapach et al. [26], which is given as follows: where N(N � 14) denotes the number of predictors; r i,t+1 represents the out-of-sample prediction result of OLS estimation of the ith predictor; ω i,t is the combined weight corresponding to the predictor; and r com t+1 denotes the combination prediction for month t + 1. Next, following Zhang et al. [51,52], we adopt five combination methods corresponding to different weight strategy.
is table shows the out-of-sample statistics performance of combination forecast. We evaluate the accuracy of out-of-sample stock return forecast by employing R 2 OS . R 2 OS computes the reduction in MSPE (mean squared prediction error) for the predictive regression model relative to the benchmark model. e statistical significance of R 2 OS is examined by employing Clark and West [24] test. * * * , * * , and * mean significance at the 1%, 5%, and 10% levels, respectively. e out-of-sample period of this table is 1948 : 01-2018 : 12. Table 8 shows that the prediction performance derived from combination forecast approaches is significantly superior to the historical average, which is consistent with antecedent literature related to stock return prediction (see, e.g., Rapach et al. [26]; Wang et al. [28]). In addition to CT constraint that fails to beat the original model and CDA strategy, our IQR constraint model and the corresponding combination model can all surpass the original model. Notably, our IQR constraint model as well as the mixed models related to IQR constraint model can generate R 2 OS of more than 1%, while other models or combination models can hardly yield R 2 OS exceeding 1%, which once again  [24] test. * * * , * * ,   and * mean significance at the 1%, 5%, and 10% levels, respectively. e out-of-sample period of this table is 1982 : 01-2018 : 12. Table 9 reports the statistical performance of the existing models and the mixed models in the alternative out-ofsample period. e second column in this table is the out-ofsample R 2 OS of the original model, which is used to compare with the model applied earlier in this paper. We find that the CT constraint model and the CDA strategy can improve the forecasting performance of 7 and 8 predictors, respectively. Notably, with the application of the IQR constraint model, all predictors can generate larger R 2 OS , and most of them are positive numbers. What is more, based on the mixed model, almost all variables achieve larger R 2 OS , especially the mixed model that combines the IQR constraint model. More specially, the CT&IQR model and the CDA&IQR model lead R 2 OS of all variables to become positive, where R 2 OS of LTY changes the most (R 2 OS increases from −0.787% to 2.758%). Overall, there is more accurate prediction information in predicting stock returns by applying the IQR constraint model and its mixed model.

Conclusion
We construct a new approach based on the principle of boxplot to constrain stock return forecasts. is approach imposes upper and lower bound to achieve reasonable return forecasts. As abnormal prediction information is usually accompanied by high risk, it is difficult for investors to believe abnormal prediction information. Fortunately, our approach not only identifies outliers but also eliminates outliers. In order to better demonstrate the predictability of our new constraint approach, we compare two constraint approaches previously proposed by Campbell and ompson [17] and Yi et al. [20]. e out-of-sample empirical results exhibit that from the perspective of statistics or economics, our new constraint approach can archive greater economic gains compared with the competition model as well as the unconstraint approach. In addition, we also analyze the predictability and economic gains of constraint and unconstraint approaches from combination model, multivariate strategy, alternative parameters (i.e., c � 2, 4, 6, ω � [0, 1.5], and TC � 0), business cycles, and alternative out-of-sample period. All empirical results demonstrate that our constraint method is robust and more stable than the unconstraint model and other competition constraint models.

Data Availability
e data or materials used to support the findings of this study will be available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.