Global COE Hi-Stat Discussion Paper Series 026

Research Unit for Statisticaland Empirical Analysis in Social Sciences (Hi-Stat)

Job Search, Bargaining, and Wage Dynamics

Hi-Stat Discussion Paper

Shintaro Yamaguchi

January 2009

Hi-Stat
Institute of Economic Research
Hitotsubashi University
2-1 Naka, Kunitatchi Tokyo, 186-8601 Japan
http://gcoe.ier.hit-u.ac.jp

JOB SEARCH, BARGAINING, AND WAGE DYNAMICSSHINTARO YAMAGUCHIA BSTRACT. This paper constructs and estimates a model of strategic wage bargaining with on-the-job search
to explore three different components of wages: general human capital, match-specific capital, and outside
option. As the workers find better job opportunities, the current employer has to compete with outside firms toretain them. This between-firm competition results in wage growth even when productivity remains the same.
The model is estimated by a simulated minimum distance estimator and data from the NLSY79. The results
indicate that the improved value of outside option raises wages of ten-year-experienced workers by 16%.

1. I NTRODUCTION
Wages of young workers grow rapidly. For example, wages of a typical white male high-school graduate
in the U.S. grow by 50% during the first ten years of his permanent transition to the labor market. This
early wage growth accounts for about two thirds of the entire wage growth over the lifetime. Given thatnew entrants to the labor market can look forward to about 40 years of work, this fact implies that a majorpart of wage growth occurs in a relatively short period of time. The objective of the paper is to empiricallyunderstand the sources of wage growth of young workers.
One standard interpretation is based on the accumulation of human capital. Workers learn skills more
quickly while they are young. The rapid wage growth can be interpreted as a direct consequence of the rapid
accumulation of human capital. Another explanation is job turnover. Young workers change jobs frequently.
Workers start their careers without detailed knowledge about the nature of their jobs or where they can findbetter jobs. They gradually find a better place to work. Wages grow when a worker finds a better job andswitches from the current job to a new job. Topel and Ward (1992) find that about a third of early career
wage growth can be explained by job switching. This paper considers another source of wage growth: thevalue of outside option of workers. This effect can increase wages even when a worker stays with the same
employer and productivity remains the same.
In the model employed workers are occasionally contacted by firms other than the current employer.
When the value of the offer is better than that of the current employment contract, the worker wants to quit
Date: revised in July 2008 (first version: November 2005).
Address: Department of Economics, McMaster University, 1280 Main St. West, Hamilton, ON., Canada L8S 4M4, URL:
socserv.mcmaster.ca/yamtaro, Email: yamtaro@mcmaster.ca.
1

the current job to accept the offer. However, the current employer may be willing to pay the worker a higher
wage to retain him. In this way, offers from outside firms improve the bargaining position of a worker when
he renegotiates wages with the current employer. Wages can grow due to an increase in outside option value,even when productivity remains the same. A worker is in a position of renegotiation, because rent from amatch exists due to search frictions. Workers do not have alternative job opportunities in hand, while firms
are not able to immediately find new workers that replace the current employees. Finding an alternative
partner who is as good as the current partner is more difficult when agents are heterogeneous, because theagents similar to the current partner can be scarce in the market. This labor market friction is a source of therent from a worker-firm match.
Understanding the sources of wage growth is important for active labor market policies. If a large fractionof wage growth is due to the labor market frictions, policies that attempt to accelerate matching process will
have significant effect on individual wages. Public Employment Services that assist job search not only
reduce unemployment, but also enhance the wage growth of employed workers.
To consider the effects of outside option value, match quality, and human capital accumulation on wages,
I develop and estimate a model of on-the-job search with strategic bargaining. The model is based on acontinuous-time search model in which unemployed and employed workers search for jobs that are characterized by idiosyncratic match quality. On the one hand, employed workers accumulate human capitalthrough learning-by-doing. On the other hand, they occasionally receive job offers from new potential employers. When the match quality with a new potential employer is high enough, a bargaining game begins
among three-players: the worker, the current employer, and a new potential employer who attempts to recruit him. In this three-player bargaining, a worker first bargains with a firm that is strategically chosen
by the worker. He then bargains with the other firm if the bargaining with the first firm is not successful.
The worker exploits these multiple job opportunities and uses a firm as a threat in bargaining with anotherfirm. Wages grow over the duration of a job, not only because workers accumulate human capital, but also
because they receive job offers, which improves their value of outside option.1
My work is related to the paper by Postel-Vinay and Robin (2002) that analyzes wage dispersion in France
using a model of on-the-job search with counteroffers. Although their model has some similar features to mymodel, they focus on wage distribution in a stationary equilibrium, instead of wage dynamics. To address

1

Barron, Berger, and Black (2005) provide empirical analysis for counteroffers using the data from the 2001 Small Business
Administration Survey. They find that factors that are likely to raise match rents indeed increase the likelihood that a firm will
extend a counteroffer to an employee with a competing offer. This empirical evidence supports the theoretical argument in thispaper.
2

the question of sources of wage growth, my model departs from their previous contributions in two ways.
First, the model allows for productivity growth of workers through learning-by-doing. The productivity
growth contributes to wage growth, because it increases not only the value of a match, but also the valueof outside option such as unemployment. The intuition is that experienced workers are more likely to get
employed and that they would have higher value of a match if they are employed. My empirical results
suggest that the growth of worker productivity actually explains the largest part of wage growth. Second,
wage renegotiation following idiosyncratic productivity shocks is considered in this paper. The model by
Postel-Vinay and Robin (2002) does not predict any wage decreases on the job and it assumes exogenous
destruction of a match. In contrast, in my model productivity shocks may cause a wage decrease on the job
and a separation endogenously occurs if a match becomes inefficient due to negative productivity shocks.
The price of these extensions is that an explicit solution to the equilibrium wage and its distribution are toocomplicated to be derived analytically because of nonstationarity. I use a numerical solution for simulationand estimation of the model. These extensions are important to match the observed dynamics of wages andworker flows.
The main contribution of the paper is to empirically explore the effect of labor market frictions on wagesof young workers. The structural parameters of the model are estimated by a simulated minimum distanceestimator using a sample of white male high school graduates from the National Longitudinal Survey ofYouth 1979. The parameter estimates are used to simulate three counterfactuals. The first exercise evaluates
the effect of the improved outside option value on wages. I find that wages of ten-year-experienced workers
would be 16% lower, if their outside option is unemployment, not the employer with the second best match.
This result indicates that the effect of the improved value of outside option is substantial and it accounts forabout 30% of the observed wage growth during the first ten years. In the second simulation, labor market
friction is shut off so that workers are able to find the best matches immediately. In the absence of frictions,wage growth is solely driven by the accumulation of human capital. In this simulation, wages grow by
25% during the first ten years of their careers, which implies that only about a half of the observed wage
growth reflects the accumulation of human capital. The third simulation imposes that human capital does not
grow. I find that human capital accumulation affects wage profile not only because it directly changes laborproductivity, but also because it alters job search behavior. A lack of human capital accumulation reduces
the value of a match, and thus, it increases the reservation match quality. As a result, the initial wage inthe career is raised by 15%, but wages grow only by 15% in the first ten years. The mean unemployment

3

duration is also increased by 20%. The results suggest that human capital accumulation has important
implication on job search and that they should be considered together.
The plan of the paper is as follows. The following section reviews the related literature. The theoreticalmodel is presented in section 2. Section 3 describes the data. The estimation strategy is detailed and theresults are presented in section 4. To evaluate the effect of human capital and labor market frictions on
wages, the parameter estimates are used to simulate counterfactuals in section 5. Section 6 concludes.
Proofs and technical details are gathered in appendices.

2. T HE M ODEL
2.1. Setup. I consider the labor market for infinitely lived ex ante heterogeneous workers and ex post heterogeneous firms. Both the workers and the firms are risk neutral, have the same discount rate r, and havecomplete access to the capital market. Workers are heterogeneous with respect to innate ability γ0 and general human capital h. Innate ability is fixed for each individual’s lifetime, while general human capital grows
with work experience. Both are precisely observed by all agents in the economy.
Unemployed workers receive utility value b, including both pecuniary and non-pecuniary unemployment
benefits, and meet a firm at a rate of λ0 . Employed workers receive a flow wage payment w and produce
a flow output y. They are also contacted by firms other than the current employer at a rate of λ . When aworker is contacted by a firm, match quality θ and an idiosyncratic shock to the match εθ are drawn from
distributions with CDFs F and G, respectively. Both θ and εθ are bounded above θMAX , εθ ,MAX and below
θMIN , εθ ,MIN . All workers and firms precisely observe the values of the match quality θ and the idiosyncratic
shock εθ . The accumulation of match-specific human capital is captured by a stochastic growth of match
quality. This modeling saves the amount of computation because no additional state variable such as tenureis used, although an alternative (and more common) way is that match-specific capital grows with tenure.
The match quality grows from θ to θ + (θ ) at the Poisson rate of δ1 (θ ). The idiosyncratic shock changes
occasionally at the Poisson rate of δ .
Firms are ex post heterogeneous in the sense that match quality varies across firms. However, they are ex
ante homogeneous because neither workers nor firms know their match quality until they meet. A match isexogenously destroyed at a rate δ0 due to family problem, health problem, and etc.
General human capital is accumulated through learning-by-doing while workers are employed and it
does not depreciate. This implies that human capital does not change while workers are unemployed. Hence
I have h(t + a) ≥ h(t) ∀a > 0 where t is the index for time and equality holds only when workers are
4unemployed from period t to period t + a. I also assume that human capital does not grow once it reaches acertain level h̄ < ∞.
The output y produced by a worker-firm match is strictly increasing in innate ability γ0 , general humancapital h, match quality θ , and idiosyncratic shock to the match εθ . For example, the production function
can take a multiplicative form y = exp(γ0 + γ1 h + γ2 h2 ) · θ · εθ where γ1 > 0 and γ2 < 0.
The assumption that all workers meet firms at the same rate and draw match quality from the common
distribution may seem restrictive, because high-skilled workers should receive more job offers and face abetter distribution of job opportunities. This observation is plausible and it is also true in the model. The
intuition is that I distinguish between meeting rates (λ and λ0 ) and offer rates. In the model, the highskilled workers receive more job offers than the low-skilled workers when they meet the same number offirms. Firms are less likely to employ unskilled workers because they are unproductive. Moreover, highskilled workers face a better distribution of job opportunities, even though all workers have the common
match quality distribution F. To see this, consider the distribution of productivity of a match, rather thanthe distribution of match quality. The distribution that high-skilled workers face is better than that for lowskilled workers in the sense of first order stochastic dominance because the productivity is increasing in
innate ability and human capital. For example, if the production technology is additive in workers’ skill andmatch quality, workers’ skill shifts the productivity distribution to the right.2
2.2. Wage Bargaining.
2.2.1. Bilateral bargaining between an unemployed worker and a firm. Unemployed workers receive unemployment benefit b and meet a firm at the Poisson rate of λ0 . When an unemployed worker meets afirm, the match quality θ and idiosyncratic shock εθ to the match are randomly determined. The worker
and the firm bargain over an employment contract given all relevant information such as skill and match
quality. Once the contract is signed, the firm pays a constant flow of wage w and the worker supplies flow oflabor service until the contract is renegotiated or separation occurs. The contract is renegotiable with mutual
consent. Also, workers can quit a job and firms can fire workers at any time.
Let Θ = {θ , εθ } be the vector of match quality and idiosyncratic shock to the match. The continuation
values of employed workers and firms are denoted by U(h, Θ, w) and V (h, Θ, w), respectively. Let U0 (h)be the continuation value for unemployed workers. I omit the innate ability γ0 in the notation, because it isconstant over time. The firms’ value of a vacancy is zero because of the free entry condition. For notational
2Strictly speaking, the value of a job is not measured by the productivity at the time when a worker and a firm meet. But, as I willsee below, the value of a job is strictly increasing in the initial productivity and the argument here still works.
5

convenience, I define the joint value J(h, Θ),J(h, Θ) = U(h, Θ, w) +V (h, Θ, w).
Observe that w does not change the joint value J, although it certainly changes the allocation. This will beclear when I see the condition for job turnover and the equilibrium payoffs in the case of job turnover, which
are both independent of wages.
The equilibrium payoffs and wage are determined by the Nash bargaining solution. The equilibrium ischaracterized as follows(1) J(h, Θ) < U0 (h): If the match is inefficient, they immediately separate to look for other potential
partners. The resulting payoffs are {U0 (h), 0},(2) J(h, Θ) ≥ U0 (h): If the match is efficient, they sign a contract immediately and the resulting payoffs
are{U(h, Θ, w),V (h, Θ, w)} = {β J(h, Θ) + (1 − β )U0 (h), (1 − β )(J(h, Θ) −U0 (h))}
where β is a measure of bargaining power of a worker.
The worker’s threat point is U0 and the firm’s threat point is 0. They split the surplus J −U0 and the worker
takes his fraction of β as well as U0 . The contract is signed only when the match is efficient.
2.2.2. Multilateral bargaining between a worker and two firms. Employed workers are contacted by firms
other than the current employer at a rate λ . When a worker meets a firm, the worker draws a match quality
with the new firm and an idiosyncratic shock to the match. The outcome is precisely observed by the worker,
the new firm, and the incumbent firm.
Dey and Flinn (2005) and Cahuc, Postel-Vinay, and Robin (2006) consider the equilibrium of this threeplayer bargaining. Dey and Flinn (2005) impose an axiomatic Nash bargaining solution assuming thatthe worker’s threat point is the match value with the losing firm. Cahuc, Postel-Vinay, and Robin (2006)
analyzes this bargaining problem using a strategic model and find a particular bargaining procedure thatcan support the assumption used by Dey and Flinn (2005). Given these two previous papers, I assume thatthe threat point of the worker is the match value with the firm that has a lower match quality. Thus, thedifference between the two match values is the surplus.
Let ΘC and ΘN be the vectors of match quality and idiosyncratic shocks to the match with the current
employer and the new potential employer, respectively. In the following, I assume that the match with the6current employer is efficient. This is actually true in equilibrium because no inefficient matches exist. Also,
the worker prefers the current employer to the new potential employer and he prefers not bargaining when
he is indifferent. The equilibrium outcome is summarized as follows(1) J(h, ΘC ) ≥ J(h, ΘN ) and U(h, ΘC , w) ≥ β J(h, ΘC ) + (1 − β )J(h, ΘN ): the worker does not initiate
bargaining. The worker stays with the current employer and the current contract remains effective.
(2) J(h, ΘC ) ≥ J(h, ΘN ) and β J(h, ΘC ) + (1 − β )J(h, ΘN ) > U(h, ΘC , w): the worker bargains with thecurrent employer to induce a counteroffer. The outcome is that the worker stays with the current
employer and he is promoted. The equilibrium payoffs of the worker and the current employer are
β J(h, ΘC ) + (1 − β )J(h, ΘN ) and (1 − β )[J(h, ΘC ) − J(h, ΘN )], respectively. The new firm fails to
recruit the worker and receives the value of vacancy, which is zero.
(3) J(h, ΘN ) > J(h, ΘC ): the worker bargains with the new potential employer. The outcome is that theworker quits the current employer to take an offer from the new potential employer. The equilibrium
payoffs of the worker and the new potential employer are β J(h, ΘN ) + (1 − β )J(h, ΘC ) and (1 −β )[J(h, ΘN ) − J(h, ΘC )], respectively. The incumbent firm loses the worker and receives the valueof vacancy, which is zero.
2.2.3. Renegotiation in response to a productivity shock. An idiosyncratic productivity shock to a match εθoccasionally changes at a rate δ . The new shock is drawn independently from the previous shocks. Firms
and workers may want to renegotiate the contract in response to the shock. I assume the following wage
renegotiation rule.3(1) J(h, Θ) ≥ U0 (h) and V (h, Θ, w) ≥ 0: When the match remains efficient and the participation constraint of the firm is still satisfied, they keep the current contract.
(2) J(h, Θ) ≥ U0 (h) and V (h, Θ, w) < 0: When the match remains efficient but the firm’s participation
constraint is not satisfied, wage is cut so that the firm’s participation constraint is satisfied withequality. The updated wage w∗ gives the worker the continuation value J(h, Θ) and it gives the firmzero.
(3) J(h, Θ) < U0 (h): When the match becomes inefficient, separation occurs. The worker’s value is
U0 (h) and the firm opens vacancy and receives zero.
Notice that a positive productivity shock does not cause a renegotiation because it does not make a firm ora worker worse off and thus, the participation constraints do not get violated due to the shock. Also, the3Postel-Vinay and Turon (2006) applies the same wage rule to demonstrate that i.i.d. productivity shocks generate persistent wage

dynamics. Assuming serially uncorrelated shocks, they are able to characterize wage dynamics analytically.
7

participation constraint of the worker does not get violated due to the shock unless it is already violated
before the bargaining is initiated. Of course, this does not happen in equilibrium.
Harris and Holmstrom (1982) find that this wage contract is optimal for risk-averse workers who cannot
commit to remaining in the job. The intuition is that these wage dynamics minimize the fluctuation of wages
given a lack of commitment, and thus they are preferred by risk-averse agents. Another interpretation is thatthe proposed wage dynamics minimize the deviation from the ex ante optimal wage dynamics in the case offull commitment. Although I assume that agents are risk-neutral for analytical simplicity, this argument isstill applicable if workers are even slightly risk-averse.
Endogenous separation is an important property from empirical viewpoint. In the model, matches withhigh-skilled or experienced workers are more likely to survive negative productivity shocks than those withlow-skilled or inexperienced workers. Matches are destroyed when they become inefficient. Workers who
have high labor productivity are unlikely to become unemployed. This prediction of the model is consistent
with many empirical findings. These wage dynamics and worker mobility are useful for a quantitative
analysis.4

2.3. Value Function. I have not discussed how value functions are defined in terms of the primitives of themodel. The next proposition provides the existence and the uniqueness of the value functions.

Proposition 2.1. The value functions U, V, U0 , and J are uniquely determined by the following contraction
mappings. The value of unemployment U0 is given by(r + λ0 )U0 (h) = b + λ0 EX max[β J(h, X) + (1 − β )U0 (h), U0 (h)]where X is a vector of the new match quality and idiosyncratic shock. The joint value of the match J is givenby

4Wages are assumed to be constant until renegotiation. Alternative assumption is that wage is a product of efficiency unit of laborand its price and that the price of efficiency unit of labor is constant until renegotiation. However, the model does not fit the observed
patterns in the data. More specifically, the model does not match the employment hazard rate and the wage variance simultaneously.
In this specification, all variations in productivity are translated into variations in wages. If I allow for endogenous separation andthe model matches the observed employment hazard rate, the model generates too large wage variance compared with the observed
wage variance. In contrast, my original specification in which agents negotiate wages fits both the employment hazard rate and thewage variance simultaneously, because productivity is not directly translated into wages. The detailed discussion and the simulation
results are provided in the appendix.
8

(r + λ + δ0 + δ1 (θ ) + δ )J(h, Θ)= y(h, Θ) + λ Pr(J(h, X) > J(h, Θ)) ·EX [β J(h, X) + (1 − β )J(h, Θ)|J(h, X) > J(h, Θ)]˙ Θ).
+δ0U0 (h) + δ1 (θ )J(h, Θ+ (θ )) + δ EX max[J(h, X), U0 (h)] + J(h,The continuation value for the worker U and the continuation value for the firm V are given by(r + λ + δ0 + δ1 (θ ) + δ )U(h, Θ, w)= w + λ EX max[β J(h, X) + (1 − β )J(h, Θ),β J(h, Θ) + (1 − β )J(h, X), U(h, Θ, w)] + δ0U0 (h)+δ1 (θ )U(h, Θ+ (θ ), w) + δ EX max[min{U(h, X, w), J(h, X)}, U0 (h)] + U̇(h, Θ, w)(r + λ + δ0 + δ1 (θ ) + δ )V (h, Θ, w)= y(h, Θ) − w + λ Pr(J(h, X) ≤ J(h, Θ)) ·EX min[(1 − β )(J(h, Θ) − J(h, X)), V (h, Θ, w)|J(h, X) ≤ J(h, Θ)]+δ1 (θ )V (h, Θ+ (θ ), w) + δ EX max[V (h, X, w), 0] + V̇ (h, Θ, w).
Proof. The value of unemployment U0 is determined by the following asset pricing equation,(r + λ0 )U0 (h) = b + λ0 EX max[β J(h, X) + (1 − β )U0 (h), U0 (h)]where the first term on the right hand side is the flow of unemployment benefit and the second term captures the value when the worker meets a firm. Notice that there is no human capital accumulation while
unemployed. Given J, the equation defines a contraction mapping for U0 .
The value of employment U is determined by the following asset pricing equation,(r + λ + δ0 + δ1 (θ ) + δ )U(h, Θ, w)= w + λ EX max[β J(h, X) + (1 − β )J(h, Θ),β J(h, Θ) + (1 − β )J(h, X), U(h, Θ, w)] + δ0U0 (h)+δ1 (θ )U(h, Θ+ (θ ), w) + δ EX max[min{U(h, X, w), J(h, X)}, U0 (h)] + U̇(h, Θ, w)

9where U̇ is a time-derivative term, which is given by U̇ = dU/dh · dh/dt. The second term on the righthand side captures the change of the continuation value when the worker is contacted by a firm other thanthe current employer. Specifically, the first element in the bracket is the value from job turnover, the second
element is the value from being promoted, and the third element is the value of remaining in the current
employment contract. The third term captures the value of exogenous separation. The value of the growthof match quality is captured by the fourth term. The fifth term captures the continuation value when an
exogenous shock hits the match. The first element is the value of keeping the current contract, the second
element is the value from renegotiation (wage cut), and the third element is the value from being fired(separation.) Given J and U0 , the recursive equation above defines a contraction mapping for U.
Similarly, the firm’s value of filling a position is determined by the asset pricing equation,(r + λ + δ0 + δ1 (θ ) + δ )V (h, Θ, w)= y(h, Θ) − w + λ Pr(J(h, X) ≤ J(h, Θ)) ·EX min[(1 − β )(J(h, Θ) − J(h, X)), V (h, Θ, w)|J(h, X) ≤ J(h, Θ)]+δ1 (θ )V (h, Θ+ (θ ), w) + δ EX max[V (h, X, w), 0] + V̇ (h, Θ, w).
where V̇ is a time-derivative term. The first term on the right hand side is the income from output, thesecond term is the labor cost, the third term captures the continuation value when the employee receives
outside offers. The value of match quality growth is captured by the fourth term. The fifth term captures thecontinuation value when a productivity shock hits the match. In the third term, the first element is the valueof promoting the employee to retain him and the second element is the value from remaining in the current
contract. In the fourth term, the first element is the value when the contract remains the same and the second
element is the value from firing (separation) or the value of employment when wage cut occurs. Given J,
the Bellman equation defines a contraction mapping for V .
Finally I solve the value function J. Combining the worker’s value function U and the firm’s valuefunction V , I find the joint value function J,(r + λ + δ0 + δ1 (θ ) + δ )J(h, Θ)= y(h, Θ) + λ Pr(J(h, X) > J(h, Θ)) ·EX [β J(h, X) + (1 − β )J(h, Θ)|J(h, X) > J(h, Θ)]˙ Θ).
+δ0U0 (h) + δ1 (θ )J(h, Θ+ (θ )) + δ EX max[J(h, X), U0 (h)] + J(h,10Notice that the joint value J is independent of wages as I mentioned at the beginning of the section, because
the wage only affects the allocation, not the size of the rent. This Bellman equation defines a contraction
mapping for J given U0 .

The equilibrium wage is set so that the worker’s continuation value U(h, Θ, w) equals his payoff fromthe bargaining. Suppose that the wage is determined through bargaining among three players, and that the
type-Θ firm wins and the type-Θ0 firm loses. The equilibrium wage w∗ satisfies(2.1)

U(h, Θ, w∗ ) = β J(h, Θ) + (1 − β )J(h, Θ0 )where the left-hand side is the payoff in a three-player bargaining which is independent of wage. The
continuation value for a worker U is strictly and monotonically increasing in wage w, which is clear fromthe Bellman equation shown above. The equilibrium wage is uniquely determined so that the continuation
value for the worker equals the payoff obtained from the bargaining. I find it intractable to derive an explicit
solution to the equilibrium wage w∗ . The wage is numerically derived in the following empirical sections.
3. DATA
3.1. Sampling Criteria. The data are taken from the National Longitudinal Survey of Youth (NLSY) 1979
which includes information on the weekly work history of individuals from 1978 through 2002. The survey
consists of individuals who were 14-21 years old as of January 1, 1979. There are some advantages fromusing the NLSY to analyze career. First, the NLSY has a rich work history. In particular, more than one
job is recorded in a survey year. Young workers frequently switch jobs and having more than one job ina year is not uncommon. To understand how young workers build their careers, it is important to record
job switching as accurately as possible. Second, the survey period covers the initial transition to the labormarket. I am particularly interested in the wages and labor mobility of young workers. The NLSY focuseson this part of the life cycle.
I use a sample of white male high school graduates for the following reasons. First, I consider the labormarket to be heterogeneous across race, sex, and education. Second, this is the largest demographic groupif I take a subset of individuals in terms of the conditions given in the first point. Third, this group is themost commonly used for wage and worker mobility analysis in the literature and thus, my results can becompared to the previous contributions to some extent.
Following Farber and Gibbons (1996), I limit my sample to jobs that were held after an individual firstmade a primary commitment to become a full-time worker. I say that an individual has made a long-term
11transition to the labor market if he works in full-time jobs for six or more quarters during three years after
graduating. A full-time job is defined as a job that consists of 30 or more hours per week worked. I also
omit self-employed jobs and jobs without pay. Details of the construction of the data set are contained inappendix B.3. The final sample consists of 9,462 observations for 662 individuals and 3,882 full-time jobs.
3.2. Variable Definitions.
3.2.1. Wages. Hourly wages are computed by dividing weekly earnings by hours worked per week. If
earnings or hours per week worked are recorded on a monthly or an annual basis, I convert them to weekly
data. Wages are deflated by the monthly consumer price index of the year of 2004.
3.2.2. Job-to-Job Transition Indicator. I say that an individual makes a job-to-job transition if (1) he stays
unemployed or has a part-time job for three weeks or less before he starts working for a new job and (2) heleaves the job voluntarily. Otherwise, I say that an individual changes jobs via non full-time employment.
For example, if a worker is unemployed for four weeks between jobs, I say that he changes jobs via nonfull-time employment. Bowlus, Kiefer, and Neumann (2001) define job-to-job transition as the condition(1) only. But condition (2) is also relevant for my analysis because between-firm competition must occur inthe model whenever job-to-job transition occurs.

4. E STIMATION
4.1. Solution and Simulation Method.
4.1.1. Empirical Specification. I allow for worker heterogeneity in terms of productivity. The productivityof a type i worker who has h years of experience, match quality θ , and an idiosyncratic productivity shockεθ is given by(4.1)

yi (h, θ , εθ ) = exp(γ0,i + γ1 h + γ2 h2 + γ3 h3 + γ4 h4 ) · θ · εθ .

The distribution of individual types is assumed to have two support points: low-skilled and high-skilled.
I normalize the low-skilled type as γ0,LOW = 0. The parameters to be estimated are the support point for
high-skilled type γ0,HIGH and its weight pHIGH . The sampling distribution of match quality is assumed to belognormal, ln θ ∼ N(µ, σθ2 ). I discretize match quality into M support points so that each point is equally
distant. Let θm be the match quality at the m-th level (1 ≤ m ≤ M). When the match quality with the current
12employer is θm , match quality grows θm+1 at the Poisson rate δ1 (θm ),δ1 (θm ) = δ1a + δ1b m i f 1 ≤ m < M
0 otherwise.
The distribution of productivity shocks is assumed to have two support points: bad state and good state.
Because it is imposed that Eεθ = 1 for normalization, I only have to estimate the support point of bad stateεθ ,BAD and its weight pBAD . I also assume wages are measured with error but this error is unbiased relative
to the true log wage ln w∗ ,ln w = ln w∗ + εε ∼ N(0, σε2 )where ln w is the observed wage. This measurement error is necessary for estimation because I will match
not only the conditional mean of wages, but also the conditional variance of wages. The interest rate r is
set to 0.075. The parameter for bargaining power β is fixed5 at 0.5 for identification purpose, which will be
discussed in the section 4.4.

4.1.2. Model Implementation. The model is numerically solved using a contraction mapping in discretizedstate space. Match quality is discretized to 20 levels. They are equally distant and their weights are givenby the method proposed by Kennan (2006). The lowest match quality is at 2.5% quantile and the highest
is at 97.5% quantile. Wages are discretized to 35 levels. The lowest hourly wages on the grid is $1.0 andthe highest is $45.0. This grid turns out not to be restrictive because any simulated wage is on either ofthese boundaries in the estimation. For a numerical solution of a continuous time model, the model must be
converted into a discrete time model. The length of a decision period is assumed to be a quarter (13 weeks.)This relatively fine grid is especially relevant while agents are young, because they change jobs frequently
and have multiple jobs during a year. If I use a year as a decision period, I would lose more information
about the early careers of young workers. All events described by Poisson process with the rate a in thecontinuous time model are approximated by a Bernoulli process with the rate a∆h where ∆h (= 0.25) is thelength of the period in the discretized model. Each individual is simulated for K (= 10) times from his firsttransition to the labor market until his first exits from the labor market or the survey ends.
5As an additional result, I also estimated this parameter. The results can be found in the appendix.
13

It is important to match the length of survey and the length of simulation so that I can re-generate the actual sampling process. Suppose I am interested in the conditional mean wage given ten years of experience.
If the agents are simulated for ten years, more productive workers are likely to be included in the simulated
sample of workers with ten years of experience, because high-skilled workers tend to be more experienced.
On the contrary, if the simulation length is as long as one hundred years, probably all types of workers will
be included in the simulated sample of workers with ten years of experience. Due to the sampling mechanism, the conditional mean wage will be higher in the former simulated sample. The simulation length canaffect the auxiliary model parameter estimates and thus, matching it with the survey length is important in
the estimation.
4.2. Model Estimator. The parameters of the model are estimated by a simulated minimum distance estimator, which is closely related to the simulated method of moments, the efficient method of moments, andindirect inference.6 The model is simulated to generate an artificial data set of employment history and wage
path. This simulated data is summarized by descriptive statistics such as hazard rate of job and unemployment, wage growth rate, and wage variance. These summary statistics of simulated data are compared withthe corresponding statistics of actual data. The objective of estimation is to look for structural parametersthat give the best fit. Summary statistics used in a simulated minimum distance estimator are called auxiliary
parameters. Auxiliary parameters are often moments of data, but they can be parameters of a reduced form
model. For example, coefficients of a wage regression are used as auxiliary parameters.
There are a few reasons to use simulated minimum distance estimator to estimate the model. First,
implementation is straightforward. For example, maximum likelihood is computationally too demanding
for my model because it requires me to integrate out serially correlated unobserved state variables. In
particular, the match quality in one job is correlated with those in other jobs. The match quality withan alternative job opportunity is also correlated over time, even when the worker does not change jobs. In
addition, idiosyncratic shocks also need to be integrated out, because they change only occasionally. Second,
I directly match the important features of the data. By looking at the gap between auxiliary parameter of thedata (i.e., summary statistics of the data) and their simulated counterparts, I can tell which features of thedata the model fails to fit.
The simulated minimum distance estimator minimizes the criterion function that is defined as the weighted
average of the squared distance between the sample auxiliary parameters and the corresponding simulated
auxiliary parameters. Let φ be the vector of structural parameters. The vector of sample auxiliary parameter
6See Carrasco and Florens (2002) for a survey.
14

estimated from the NLSY is denoted by ρ̂. The vector of simulated auxiliary parameter computed from asimulated data set with K simulations is denoted by ρK . The simulated minimum distance estimator φ̂ isdefined asφ̂= arg min[ρK (φ ) − ρ̂]0W [ρK (φ ) − ρ̂]φwhere W is a positive definite matrix. The variance of the estimated structural parameters is minimized
when the weighting matrix W is the inverse of the covariance matrix of the auxiliary parameters. Noticethat an auxiliary parameter is a function of structural parameters. The choice of auxiliary parameters affects
efficiency and they should be good descriptions of the data. One interpretation is that I see the data throughthe “lens” of a descriptive statistical model. When the optimal weighting matrix W ∗ is used, the covariancematrix for the structural parameters is given byVar(φ̂ ) = (1 + 1/K)


∂ g0 (φ ) ∗ ∂ g(φ ) −1W∂φ∂φ

where W ∗ = Var(ρ̂)−1 , g(φ ) = ρ(φ ) − ρ̂, and K is the number of simulation.

4.3. Auxiliary Models. The auxiliary parameters that describe employment history data are hazard rates of
40 , full-time employment spells {p }40 , and non full-time employment spells
full-time job spells {pJ,t }t=1
F,t t=1
5 . Non full-time employment includes part-time employment and unemployment. A spell of non{pN,t }t=1

full-time employment ends if a worker finds a new full-time job. The hazard rate of non full-time employment captures the transition from unemployment or part-time employment to full-time employment. A spellof full-time employment ends if a worker loses a full-time job and gets unemployed or switches to a new
part-time job. Thus, the hazard rate of full-time employment captures the transition from full-time employment to unemployment or part-time employment. Lastly, a spell of full-time job ends if a worker loses
full-time employment status or he switches to another full-time employment. The hazard rate of full-time
job captures the job-to-job transition as well as a transition to non full-time employment. Quarterly hazard
rates of full-time job spells (up to 40 quarters), full-time employment spells (up to 40 quarters), and non
full-time employment spells (up to 5 quarters) are nonparametrically estimated by a Kaplan-Meier estimatorusing both censored and non-censored spells. The length of the periods is chosen so that roughly 90% of
spells die during that time interval. For example, the hazard rate of a full-time job at quarter t is given by
pJ,t

=

Nt
C )
∑τ≥t (Nτ + Nτ+1
15

where Nt is the number of jobs that end at quarter t and NtC is the number job spells that are censored atquarter t.
The auxiliary parameters that describe the wage data are the average log wage conditional on general work
experience, within-job wage growth rate, and within-individual wage growth rate. Let wi jt , GXi jt and T Ni jt
be log wage, experience, and tenure of individual i when he works for employer j in period t. Following
Murphy and Welch (1990), the average log wage conditional on general work experience is estimated by
OLS with quartic term,
wi jt= β0,OLS + β1,OLS GXi jt + β2,OLS GX 2 + β3,OLS GX 3 + β4,OLS GX 4 + εi jt

2where εi jt is a statistical error. The coefficients βOLS and the residual variance σOLS
are used as the auxiliary

parameters. This auxiliary model provides a parametrically estimated conditional mean and variance oflog wage given experience. Notice that the estimated coefficients should not be interpreted as returns toexperience because of endogeneity. In particular, experience and unobserved innate ability are positively
correlated. This motivates me to estimate within-individual wage growth rate using a fixed effect model.
Denote the average log wage and experience of individual i by wi and GX i . Let w̃i jt = wi jt − w̄i be thedeviation of log wage from the average over the entire job history of individual i. I also define the averages
and the deviations for GX and ε in a similar way. The fixed effect model for within-individual wage growthis given byw̃i jtg4 i jt + ε̃i jt .
g3 i jt + β4,FE GX
g2 i jt + β3,FE GX
f i jt + β2,FE GX= β1,FE GX

Notice that I use wage variation within individual job histories to identify returns to experience. The coeffi2 are used as the auxiliary parameters. The residual variance of thiscients βFE and the residual variance σFE
2 is expected to be smaller than the residual variance of the OLS σ 2 , because thefixed effect estimator σFE
OLS

individual fixed effect is eliminated in this model. Thus, the difference between two is useful to identify thedispersion of individual innate ability.
The estimated returns to experience include both within-job wage growth and the wage growth throughjob switching. To understand the wage growth structure more deeply, I estimate within-job wage growthusing a fixed effect estimator. Denote the average log wage and tenure of individual i at job j by wi j and
T N i j . Let w̃˜ i jt = wi jt − wi j be the deviation of log wage from the average over the duration of job j. I also
define the averages and the deviations for other variables such as T N and ε in a similar way. The fixed effect16

model for within-job wage growth is given byw̃˜ i jt

g4
g3
g2
f
N i jt + ε̃˜i jt .
N i jt + β4,OJFE Tg
N i jt + β3,OJFE Tg= β1,OJFE Tf
N i jt + β2,OJFE Tg

4.3.1. Estimates of Auxiliary Parameters. The estimated hazard rates are plotted in figure A.4 through A.6.
Some selected numbers are also presented in table 1 for detailed information. The quarterly hazard rate ofnon full-time employment at the first quarter is 0.51. It decreases to 0.46 in the fifth quarter. The estimates
imply that the average time to find a full-time job is about 2 quarters. The quarterly hazard rate of full-time
employment at the first quarter is 0.16. It quickly decreases in the first few years and drops to 0.05 at thetenth quarter. The rates seem to be stable around 0.03 after the fifteenth quarter. The hazard rate of full-time
job at the first quarter is 0.20, which is higher than the employment hazard rate because some jobs end with
job-to-job switching. Similar to the full-time employment hazard rate, the job hazard rate decreases quickly
in a first few years and it drops to 0.07 at the tenth quarter. The pattern found here is very similar to the one
reported by Farber (1999).
The estimated sample auxiliary parameters for wages are presented in table 2. The first two columns are
the OLS estimates. The results imply that the average initial log wage is 2.1 (= $8.1) and the average log
wage with ten years experience is 2.6 (= $13.6.) The next two columns present the regression estimates forwithin-individual wage growth. The results imply that the average annual return to experience in the firstten years is 0.054. Also notice that the residual variance is smaller than that of OLS, because individual
effect is eliminated in this model. The last two columns show the regression estimates for within-job wage
growth. The results imply that the average annual wage growth rate on the job is 0.033. Again, the residual
variance is smaller than those from OLS and within-individual wage regressions, because both individualand employer fixed effects are eliminated.
The distance between the sample auxiliary parameters and the simulated counterparts is optimally weighted
by the estimated variance-covariance matrix of the sample auxiliary parameters. I estimate it by a nonparametric bootstrap with 10,000 replications. To preserve the dependent structure of the data, I sample individuals with blocking over the time dimension. If a particular individual i is selected, his entire wage andemployment history are included in the sample. There are some advantages using a bootstrap to estimate
the optimal weighting matrix. First, I can avoid small sample bias for second order moments as is pointed
out by Altonji and Segal (1996). Second, analytical solutions to covariance between auxiliary parameters
can be complex. But I can estimate it through simulation. The estimated standard deviations are presented
in table 1 and 2. Third, no update of the weighting matrix is necessary. This simplifies computation.
17

The number of auxiliary parameters is 101 (40 from the full-time job hazard rate, 40 from the full-time
employment hazard rate, 5 from the non full-time employment hazard rate, 6 from the OLS wage regression,
5 from the fixed effect model for on-the-job wage growth, 5 from the fixed effect model for overall wage
growth.) The simulated minimum distance estimator φ̂ is formally given byφ̂= arg min[ρk (φ ) − ρ̂]0W ∗ [ρK (φ ) − ρ̂]
φ

where
φ

= {γ1 , γ2 , b, λ0 , λ , δ0 , δ , γ0,HIGH , pHIGH , εθ ,BAD , pBAD , α, θMIN , σε }

ρ =


40
40
5
2
2
2
{pJ,t }t=1
, {pF,t }t=1
, {pN,t }t=1
, βOLS , σOLS
, βFE , σFE, βOJFE , σOJFEW ∗ = Var(ρ̂)−1 .
Remember that ρ̂ is the sample auxiliary parameter estimated by the observed sample.

4.4. Identification. Flinn and Heckman (1982) show that search models are under-identified if only unemployment duration and initial accepted wages are used. Because employment duration and the subsequent
wages after the initial accepted wages are also used in the estimation, I might be able to fully identify themodel. However, identification through such information is not fully studied.7 Thus, following Flinn andHeckman (1982), I impose three parametric assumptions to satisfy the recoverability condition to assure
the identification of the structural parameters. The additional information such as employment duration andwages other than the initial ones improves the efficiency of the estimator. First, I fix the annual discount rate
r at 0.075. They show that it is not separately identified from the value of unemployment benefit b. Second,
the match quality distribution is assumed to be lognormal. This shape restriction is necessary to identify
the fraction of unaccepted match quality, and thus the rate at which an unemployed worker meets a firm. In
other words, the reason for unemployment cannot be decomposed into the reservation match quality and job
availability without this shape restriction. Third, the bargaining power of workers β is set to 0.5, because
it is weakly identified without additional structure.8 Intuitively, it is difficult to distinguish the bargaining
power from the size of the ‘pie.’ These parametric assumptions are commonly used in the literature such as
paper by Dey and Flinn (2005).
7Barlevy (2005) finds nonparametric identification conditions for a search model without bargaining.
8The main result is reported by assuming β =0.5. I also estimated the model parameters including β and find that the results arelittle different from the case with fixed bargaining power. See the appendix B.2 for the details.
18

Although the model is identified using the parametric assumptions, it is probably worthwhile to intuitively
discuss the identification of the bargaining effect on wages. For simplicity, let us suppose that the matchspecific capital (e.g. tenure effect) does not grow on the job. More general case will be discussed below.
Consider two groups of workers with same general work experience. Workers in the first group have worked
for years in the same employer, and thus their tenure is long. Workers in the second group were unemployed,
but they have just found a new job. The bargaining position and the match quality of workers in the formergroup are expected to be high, because they stay with the same employer for years for good match quality
and because they should have received many job offers from outside on the job. Thus, the wage growthof these workers is solely driven by the accumulation of general human capital. On the other hand, onthe-job wage growth of workers with zero tenure is also driven by the change in bargaining position. The
difference in the average on-the-job wage growth between two groups identifies the wage growth throughbargaining. Notice that workers in both groups are highly experienced, because long tenure implies highlabor market experience. Nevertheless, the estimated effect of bargaining can be applied to inexperienced
workers, because the match quality distribution and the rate at which a worker meets a firm are same acrossworkers with different experience, which is also a common identification assumption in the search literature.
When match-specific capital grows over time, its effect on wages must be taken into account in order to
identify the bargaining effect on wages. This effect may not be negligible for workers with zero tenure, while
it probably has little effect for workers with long tenure. The growth of match-specific capital is included in
the model so that match quality with the current employer stochastically grows. This productivity component
is identified using job-to-job transition rate. When the growth of match quality is fast, job turnover rate isstrongly decreasing in tenure (i.e., the slope of job hazard rate will be steep.) The contribution of thiscomponent to wage growth is derived from the model, because it explicitly relates productivity to wages.
Thus, the effect of bargaining on wage growth is still identified by comparing the two groups.
Although I do not have a theorem about the uniqueness of the estimated parameters, I address this potential problem numerically by using the simulated annealing algorithm9 in the estimation that is one of global
optimization algorithms. Unlike local optimization algorithms such as Nelder-Mead (simplex) method andquasi-Newton algorithm, a global optimization algorithm does not stuck in local optima. The probability thatthe simulated annealing algorithm terminates with the global optimal solution approaches 1 as the annealing

9

See Goffe, Ferrier, and Rogers (1994) for a detailed discussion. They also provide a FORTRAN software.
19schedule is extended. The parameters are estimated by the simulated annealing with eight different starting parameters. Each run includes about 80,000 function evaluations and converges to the same parametervalues, which strongly supports that the reported parameters uniquely minimize the objective function.

4.5. Estimation Results. The estimation results are presented in table 3. A worker’s productivity growsby 25% during the first ten years of his career through accumulation of general human capital. Noticethat this is not the entire productivity growth of a typical worker with ten years work experience, because
the productivity also grow through match quality change. The value of unemployment benefit amounts to$5.1 per hour in 2004 dollar. This value is lower than the minimum wage during early 80s, which was
between about $7.00 and $5.50. Those who do not have a full-time job meet 3.8 firms during a year. Fulltime job holders are contacted by 1.2 firm during a year. Notice that these numbers are not same as thenumber of offers actually made by firms. In many cases, they immediately separate because they realize
they do not have a good match quality as soon as they meet. The match quality distribution is assumed to
be lognormal. The mean and the standard deviation of log of match quality are 0.52 and 0.32. The rate at
which the match quality grows is about 0.2 for workers whose match quality is at the median of the sampling
distribution. When match quality changes, it improves by 11% (calculated by taking log difference.) Given
these estimates, the average annual match quality growth rate is about 2%. A negative productivity shock
occurs once in every 1.5 years. This negative shock causes a significant productivity drop. When a workerfirm match is in a bad state, their output is almost zero. Also a match is destroyed due to an exogenous shock
once in 8.4 years. The proportion of high-skilled workers is 67% and they are 56% more productive thanlow-skilled workers. Lastly, the lowest wage and the highest wage in the simulation are $4.59 and $39.6,
respectively. This implies that my wage grid is not restrictive, because the grid covers wages between $1.00
and $45.00.

4.6. Model Fit.

4.6.1. Hazard Rates. The sample hazard rates and simulated hazard rates are plotted in figures A.4, A.5,
and A.6. The results indicate that the model captures the basic features of the hazard rates. The hazard
rate of full-time employment is declining because of both duration dependence and heterogeneity. Workers
accumulate human capital through learning-by-doing. This process makes a worker-firm match more andmore productive. As a consequence, the transition rate from full-time employment to unemployment (orpart-time employment) decreases over time. Matches are heterogeneous in terms of the innate ability and20

the experience of the worker and the match quality. Because less efficient matches are destroyed by negative
productivity shocks, the selection mechanism generates declining hazard rate over time.
The full-time job hazard rate decreases over time for a similar reason, but there is one important difference
because this includes job-to-job transition. Observe that human capital and innate ability do not affect the
job-to-job transition rate, because these are transferable across jobs. However, a worker is more likely toseparate from the match and switch to a new job if the match quality is low. This is another source thatgenerates a declining hazard rate.
4.6.2. Wages. The sample and simulated auxiliary parameters of wage regression are presented in table
4. The prediction of average log wage given experience is reasonably close to the data, as is shown infigure A.1. The difference between the simulation and the observation is 0.02 when experience is ten years.
Notice that the slope quickly decreases around the point where experience is five years. This feature cannot
be captured by wage regression if only linear and quadratic terms are used, which is also emphasized by
Murphy and Welch (1990). The prediction of average within-individual wage growth rate also captures thebasic features of the data (see figure A.2), although the simulated wage growth path is less concave thanthe observation. The model generates a similar wage growth path to the simulation, partly because humancapital accumulation function includes a quartic term. The difference between the simulated path and theobservation is largest at five years of experience and it is 0.03. Figure A.3 presents within-job wage growthrate. The model captures the basic patterns of on-the-job wage growth. However, the simulated on-thejob wage growth path is less concave than the observation: the simulated wage growth is slower than theobserved wage growth when tenure is short, while it is too fast when tenure is long. A possible explanation
for the gap between the model and the observation is that match quality distribution and the rate at which
a worker meets a firm are assumed to be the same across workers with different experience and tenure. If
inexperienced workers (or workers with short tenure) meet firms more frequently than experienced workers(or workers with long tenure), the wages of inexperienced workers will grow more rapidly, which helps themodel fitting the observed on-the-job wage growth path.
5. C OUNTERFACTUALS
5.1. Is Improving Outside Option Important? The most interesting feature of the model is that not only
labor productivity, but also the value of outside option changes over time and contributes to wage growth.
The effect of the outside option value on wages can be evaluated by comparing wages with a counterfactual
wage when the outside option is unemployment. The actual (predicted) wage profile and the counterfactual
21

wage profile are presented in figure A.7. At the initial wage (zero experience) the actual and the counterfactual wages are same because the actual outside option at the initial wage is unemployment. However,workers quickly improve their outside options, particularly for the first five years. The actual wage is 15%
higher than the counterfactual wage when workers have five years of experience. After five years of experience, the bargaining position little contributes to wage growth. The cumulative wage growth due tothe outside option value is 16% when experience is ten years. Because wages grow by about 50% during
the first ten years, the improved value of outside option can account for about 30% of this wage growth.
This simulation exercise indicates that improving outside option substantially contributes to wage growth ofyoung workers, particularly in the very early stage of their careers.
5.2. Reducing labor market frictions. Another empirical question of the paper is what the effect of thelabor market frictions on wages is. Due to the labor market frictions, workers have to search for better jobs
and outside options, and bargain wages. In the absence of frictions, workers are able to immediately findthe best match and an alternative job opportunity which is as good as the best one. Thus, the wage growthis solely driven by the growth of human capital in this economy. By comparing the actual wage profile
with the counterfactual wage profile in the absence of frictions, I can quantify the effect of the labor market
frictions on wages. Since the simulation of the model is carried out in a discrete time model, it is impossible
to simulate an economy with absolutely zero friction. But it is still worth simulating an economy with thefeasibly lowest possible frictions.
To implement the idea, the model is parameterized in the following way. First, workers always havethe best match quality with all firms. This “best” match quality is at the 97.5% quantile in the estimatedunderlying match quality distribution and the idiosyncratic productivity shock is in a good state. In thesimulation of the estimated model, 18.6% of workers have this match quality when they have ten years ofwork experience. Second, workers meet one firm in every quarter regardless of their employment status.
So the parameters are set so that λ0 = λ = 4. Remember that Poisson processes are approximated by thecorresponding Bernoulli processes and that the decision period is a quarter. This parameterization implies
that the probability that a worker and a firm meet in a quarter is 1. Third, workers do not lose a job. In alabor market without frictions, workers are able to find a job immediately even if they lose one. To bringabout this situation in a discrete time model, the parameters are set so that δ0 = δ = 0. All other parameters
remain the same as the estimated model.
The simulated wage profiles are presented in figure A.8. There are two interesting features. First, the wage
growth in the frictional labor market is more rapid than the wage growth without frictions. In the absence of22

frictions, wage growth is solely driven by the growth of human capital. To the contrary, in a frictional labormarket, workers are improving their match quality and the value of outside option, which also increases
their wages. The wage gap between two wage profiles is quickly decreasing during the first five years by
about 25% (measured by log difference) and the remaining gap is little reduced after that. Because theobserved wage growth rate during the first ten years is about 50%, labor market frictions account for about
a half of the observed wage growth during this period. This estimate is greater than the number claimed
by Topel and Ward (1992), which is about a third, because they do not take into account the effect of the
improved value of outside option. The second interesting feature of the comparison of two wage profiles
in figure A.8 is that the wage gap remains unfilled even after ten years of job search. This is explained by
the average productivity of worker-firm matches. In a labor market without frictions, workers are always
assigned to the best match. But, in a frictional labor market, workers are occasionally taken away from a job
due to productivity shocks even if they once find the best match quality. Thus, the average match quality ina frictional labor market is strictly lower than the best match quality even after several years of job search.
5.3. No Human Capital Growth. Yet another counterfactual is simulated to examine the effect of humancapital accumulation on wages and job search behavior. In this simulation, workers do not accumulate their
human capital, which has two implications. First, wages, as well as labor productivity, grow slowly because
only job search and bargaining improve the continuation value for workers. Second, the reservation match
quality is increased because a lack of human capital growth reduces the value of a worker-firm match whilethe value of unemployment is little reduced. As a result, unemployment spells become longer than thosein the standard case. This counterfactual experiment is implemented by setting γi = 0 for i = 1, · · · , 4. Allother parameters remain the same as the estimated model.
There are two interesting features in the simulated log wage profiles in figure A.9. First, the initial wages
are increased by about 15% if human capital does not grow. This high initial wages are consequences of ahigh reservation match quality. Remember that the continuation value for the worker is proportional to thesum of the current wage and the expected future income gains as shown in proposition 2.1. When humancapital does not grow, the value of a match is reduced for a given level of match quality. Because the valueof unemployment changes little, workers have to wait for a higher match quality to start a job. Since theinitial match quality is high, there is not much room for the match quality to be improved through on-the-job
search. This high reservation match quality and a lack of human capital accumulation decrease the expected
future income gains. When the changes of continuation value for the worker are small and the expected
future income gains decrease a lot, the current wage has to increase. In other words, when workers cannot
23expect high income in the future, they must be compensated by a high wage today. This is why the initial
wages are higher when human capital does not accumulate.
Second, wages grow slowly in the absence of human capital growth. Specifically, wages grow only by
15% on average during the first ten years of careers (see table 5.) This is partly due to my specification ofmultiplicative production technology, but a high reservation match quality is also the reason. When workers
start with a high match quality, there is not much room to improve match quality through job search. These
explanations also account for the difference in wage growth rate between low-skilled workers and highskilled workers. High-skilled workers are not only more productive, but also have a lower reservation match
quality than low-skilled workers, as shown in table 6.
The increased reservation match quality also results in a longer unemployment duration. In the standard
simulation with human capital growth, the predicted unemployment hazard rate in the first quarter is 0.48.
But it decreases to 0.39 if human capital does not accumulate. In my simulations, the average unemployment
duration increases from 2.1 quarters to 2.5 quarters.
Through this simulation exercise I find that human capital accumulation affects wage growth not only byits direct effect on the current labor productivity, but also by its indirect effect on job search behavior due tolow future productivity.

6. C ONCLUSION
The main contribution of this paper is to investigate wage dynamics using a model of wage bargaining
with on-the-job search. This is a departure from the conventional wage growth literature in the sense that
I consider a new component of wage growth: the improved value of outside option. The model improves
on previous work in two respects: it allows for productivity growth, and it deals with wage renegotiation
following productivity shocks. These changes make the model analytically complex and the equilibrium
wage is numerically solved rather than using the explicit wage solution used in previous work. However,the first change allows me to examine the effect of general human capital growth, which makes a substantial
contribution to wage growth. I find that human capital accumulation affects wage profile by changing job
search behavior as well as labor productivity itself. The model also generates a wage decrease on the job
and endogenous separation due to productivity shocks. Both of these are necessary to match the observed
wage dynamics and worker mobility.
The model captures important features of wage dynamics and worker mobility in the NLSY 79. The
empirical results show that the contributions of job search and bargaining to wage growth is significant. The
24

improved value of outside option accounts for about 30% of the observed wage growth during the first ten
years. This paper also quantifies the effect of human capital on wage growth and search behavior. Without
human capital accumulation, the reservation match quality is increased, which results in 15% higher initial
wage and only 15% of wage growth in the first ten years. The results suggest a tight relationship between
human capital accumulation and job search.
This paper provides a detailed econometric model of bargaining with on-the-job search that allows forrich forms of heterogeneity and time-varying productivity components. Nevertheless, there are important
extensions that will be addressed in the future work. First the role of fellow workers in wage bargaining
is ignored. Large firms can easily replace their employees from their internal labor market. Thus, thebargaining position of the worker tends to be low in larger firms. Second, parametric assumption for match
quality distribution is strong. The assumption of lognormality is common in the literature, but there are anumber of possible alternative assumptions. Nonparametric identification might be possible using the ideaof record statistics, which is proposed by Barlevy (2005) for a job search model without bargaining.

25

A PPENDIX A. TABLES AND F IGURESTABLE 1. Sample Hazard Rates

Note: The hazard rates are estimated by the Kaplan-Meier estimator. Non full-time employment includes part-time employment
and unemployment. The standard deviation of the estimated hazard rates is obtained by a nonparametric bootstrap with blocking
for time dimension at the individual level. The number of replications in bootstrap is 10,000. The survival rate of full-time jobs
and full-time employment at 40th quarter is 0.07 and 0.14, respectively. The survival rate of non full-time employment at 5th
quarter is 0.04.

TABLE 2. Sample Auxiliary Parameters

Note: Log of hourly wages deflated by 2004 price are regressed on experience (Exp) or tenure (Ten.) Eighteen observations ofhourly wages are higher than $100 and they are excluded from regressions. The fixed effect estimator is used to estimate
within-individual wage growth and within-job wage growth. The standard deviation is obtained by a nonparametric bootstrap withblocking for time dimension at the individual level. The number of replications in bootstrap is 10,000.
26

TABLE 3. Parameter Estimates
Parameter
Notation Estimates Std. Dev.
Human Capital
γ1
0.448
0.061
γ2
-0.346
0.094
γ3
0.150
0.046
γ4
-0.0020.001
Unemp. Benefit
b
5.140
0.657
Poisson Rate
Unemp.
λ0
3.869
0.322
Emp.
λ
1.261
0.088
Separation
δ0
0.119
0.004
Match Q. Change
δ1a
0.300
0.047
δ1b
-0.0090.002
Prod. Shock
δ
1.811
0.107
Type Dist.
High-skilled
γ0,HIGH
0.566
0.079
Weight
pHIGH
0.674
0.175
Shock Dist.
Badεθ ,BAD-2.7390.218
Weight
pBAD
0.385
0.020
Match Q. Dist.
Mean
µ
0.523
0.034
Std. Dev.
σθ
0.320
0.096
Measurement Error
σε
0.250
0.014
Interest Rate: r = 0.075
Bargaining Power: β = 0.5
No. of Moment Conditions: 101
Function Value at the Minimum: 211.648
Note: The coefficients of production function are denoted by γ. The hourly value of unemployment compensation measured in
2004 dollar is b. The annual arrival rate of an offer for those who do not have a full-time job is λ0 . The annual offer arrival rate on
the full-time job is λ . The annual arrival rate of an exogenous separation is δ0 . The parameters for the annual rate of match qualitychange are δ1a and δ1b . The annual arrival rate of productivity shock is δ . The support point for high-skilled individuals and its
weight are denoted by γ0,HIGH and pHIGH . The support point for low-skilled workers are normalized to zero. The support point
for bad (negative) shock and its weight are εθ ,BAD and PBAD . Those for good (positive) shock are εθ ,GOOD and PGOOD . I impose a
restriction that the mean of the shock be one for normalization. The location and the scale parameter of match quality distribution,
which is lognormal, are µ and σθ , respectively. The standard deviation of measurement error is σε .

TABLE 4. Sample and Simulated Auxiliary Parameters (Wage Regression)

Note: Hourly log-wages are regressed on experience (EXP) or tenure (TEN). For the detailed information on the results for sample
auxiliary parameter estimation, see table 2.
27

TABLE 5. Log-Wage Profile with No Human Capital Growth.

TABLE 6. Reservation Match Quality.

Note: The reservation match quality depends on the state of idiosyncratic shock. The ranks in thediscretized match quality space are in parenthesis (20 is the highest.)
28

F IGURE A.1. Average Log-wage by Experience

F IGURE A.2. Within-Individual Wage Growth Rate

F IGURE A.3. Within-Job Wage Growth Rate

29

F IGURE A.4. Non Full-time Employment Hazard Rate

F IGURE A.5. Full-time Employment Hazard Rate

F IGURE A.6. Full-time Job Hazard Rate

30

F IGURE A.7. Wage Profile When Outside Option is Unemployment

F IGURE A.8. Log-Wage Profile in a Labor Market with Few Frictions.
31

F IGURE A.9. Log-Wage Profile with No Human Capital Growth.

F IGURE A.10. Unemployment Hazard Rate

32

A PPENDIX B.
B.1. Alternative Bargaining Assumption. In this paper wages are assumed to be constant until renegotiation. However, this is not the only possible assumption. Alternative assumption is that wage is a productof efficiency unit of labor and its price and that the price of efficiency unit of labor is constant until renegotiation. In this alternative specification, wages grow over time without renegotiation as workers accumulate
their human capital. Let q be the price of efficiency unit of labor and y be the labor productivity. The flow
wage payment is then w = qy. In bargaining, the price is set so that the continuation value of the worker
coincides with the equilibrium payoff in the bargaining. Suppose that the wage is determined through bargaining among three players, and that the type-Θ firm wins and the type-Θ0 firm loses. The equilibrium
price q∗ satisfies
U(h, Θ, q∗ ) = β J(Θ, h) + (1 − β )J(Θ0 , h).
Similar to the equation (2.1), this equilibrium price q∗ is also determined uniquely.
One advantage of this specification is that wage decomposition is straightforward. Assuming the production technology is given by equation (4.1), wages can be decomposed asln w = γ0,i + γ1 h + γ2 h2 + γ3 h3 + γ4 h4 + θ + εθ + ln q(h, Θ, Θ0 )where q only captures the effect of change in outside option value. For normalization, the lowest price q is
set 0.01 and the highest price q̄ is set 1.00, because price and labor productivity are not separately identified.
I have estimated this model using the same set of auxiliary parameters. The estimation results (available
upon request) indicate that the model does not fit the observed patterns in the data. More specifically, themodel does not match the employment hazard rate and the wage variance simultaneously. In this specification, all variations in productivity are translated into variations in wages. If I allow for endogenous separationand the model matches the observed employment hazard rate, the model generates too large wage variancecompared with the observed wage variance. In contrast, my original specification in which agents negotiate
wages fits both the employment hazard rate and the wage variance simultaneously, because productivity is
not directly translated into wages.

B.2. Additional Estimation Results. For the main results of the paper, the bargaining power parameter isfixed at 0.5. To be able to estimate the bargaining power in a reliable way, employer side information should
be used. But, it is not included in the NLSY. However, I also attempt to estimate the parameters including
33

bargaining power to see how other parameter estimates but bargaining power are affected. The results are
summarized in table 7. The estimated bargaining power is 0.49, which is close to the parameterized valueof 0.5. Hence, other parameter estimates are also close to those in the main results.
TABLE 7. Parameter Estimates
Parameter
Human Capital

Notation Estimates Std. Dev.
γ1
0.437
0.057
γ2
-0.404
0.060
γ3
0.184
0.026
γ4
-0.001
0.001
Unemp. Benefit
b
4.620
1.174
Poisson Rate
Unemp.
λ0
3.673
0.270
Emp.
λ
1.143
0.084
Separation
δ0
0.117
0.005
Match Q. Change
δ1a
0.312
0.037
δ1b
-0.0100.001
Prod. Shock
δ
1.969
0.114
Type Dist.
High-skilled
γ0,HIGH
0.653
0.105
Weight
pHIGH
0.772
0.101
Shock Dist.
Badεθ ,BAD-2.7550.338
Weight
pBAD
0.392
0.028
Match Q. Dist.
Mean
µ
0.486
0.024
Std. Dev.
σθ
0.402
0.074
Measurement Error
σε
0.252
0.018
Bargaining Power
β
0.494
0.060
Interest Rate: r = 0.075
No. of Moment Conditions: 101
Function Value at the Minimum: 205.847
Note: See table 3.

B.3. Details of Data Set Construction.
B.3.1. Weekly Job Information. First I collect information about each job. In particular, I am interestedin (1) the start week of the job, (2) the end week of the job, (3) wages, and (4) hours worked per week.
To determine the start and the end week of a job, I collect information from the most recent survey. I use
white male high school graduates in this analysis. High school graduates are defined as individuals who
have high school diploma but did not pursue further formal education until the most recent survey year.
There are 950 white male high school graduates in my sample and they had 10,128 jobs in total. Because anindividual may have a job for more than one survey year, I have multiple observations for each job. I have
24,565 observations in my sample. To reduce the initial condition problems, I need information about thetransition from school to work. I omit individuals who graduated high school before January 1978, which34

is earlier than the first survey year. This condition leaves me 18,114 observations for 685 individuals and7,780 jobs. Because I am interested in jobs after the transition to labor market is made, I omit jobs that
ended before graduating. I have 16,364 observations for 684 individuals and 6,555 jobs after imposing thiscondition. Next I omit self-employment jobs and jobs without pay. Then I have 15,695 observations for 684
individuals and 6,417 jobs.
I determine whether a job is full-time or not by taking the average hours over all observations for each job.
A full-time job is defined to be a job that individuals work at for 30 hours per week or more. The numberof hours per week worked is topcoded so that it does not exceed 96 hours before the average is calculated.
My model precludes a transition between part-time and full-time without changing employer. I expect that
this does not cause a serious bias, because more than 90% of jobs (5,914 jobs) did not have such switching.
Omitting all part-time jobs, I have 13,686 observations for 678 individuals and 5,388 full-time jobs in thesample.
Dual Jobs. Several individuals have more than one job at the same time. In particular, there are some shortterm full-time jobs that are started and ended while an individual holds a primary full-time job. I omit theseshort term jobs for two reasons. First, these jobs may be misclassified as full-time because individuals areunlikely to have two or more full-time jobs. Second, even if they are correctly classified, those secondaryjobs are temporary and less important for an individual’s career. There are 601 secondary jobs and the
median duration of a job is 17 weeks (the mean duration is 34 weeks.) On the other hand, there are 4,787
primary jobs and the median duration of a job is 51 weeks (the mean duration is 126 weeks.)
Some jobs are started before old jobs are ended. To reconcile the gap between the data and the model,
I move the start week of a job, while the end week of a job is not edited. If there are two jobs that haveoverlapping periods, I assume that the second job started right after the first job is ended. There are 250
jobs (about 5% of all jobs) of which start weeks have to move. The median change is 15 weeks (mean is 34
weeks.)
After all the operations described above, I have 12,846 observations for 678 individuals and 4,787 jobs.
The information is used to construct a quarterly job array in the followings.
Job-to-Job Transition Indicator. As is mentioned in the main text, I say that an individual makes a job-to-job
transition if (1) he spends three weeks or less before he starts working for a new job and (2) he leaves thejob voluntarily. Otherwise, I say that an individual changes jobs via non full-time employment. Under thisdefinition, I find that 1,612 jobs are started as job-to-job transitions, while 2,480 jobs are not. If I change the35TABLE 8. The number of observations.

search duration from three weeks to two weeks, 1,461 jobs are started as job-to-job transition. If I change
the search duration from three weeks to four weeks, 1,694 jobs are started as job-to-job transition.
B.3.2. Quarterly Work History Array. Using the information constructed in the previous subsection, I construct a quarterly work history array for each individual. A quarter consists of 13 weeks. The quarterly labor
force status of an individual is determined by the labor force status of the first week in a given quarter. Thisimplicitly assumes that the labor force status does not change during a quarter, although it may for someshort term jobs. It is true that I could avoid this discretization bias by making the decision period shorter,
but that may cause another problem. Weekly work history data would contain more noise. As Neal (1999)
points out, NLSY work history file includes measurement errors. As I make the decision period shorter, I
may pick up every noise contained in the data. Aggregating the weekly information to quarter may moderate the bias. In addition, using shorter decision period can be computationally expensive and make myestimation intractable.
Entry and Exit from the Full-time Labor Market. I trace each individual from his initial long-term transition
to the full-time labor market until he exits from the full-time labor market or the survey ends. The definitionof an initial long-term transition to the full-time labor market is similar to that of Farber and Gibbons (1996).
I consider an individual has made this transition if he works for six or more quarters during three years.
Exit from the full-time labor market occurs when an individual does not work at a full-time job for three
years or more. If this is observed in the data, I record his separation from a job and drop the information
after that. Re-entry to the full-time labor market is not included in my sample, because I do not have any
information about when individuals re-entered the labor market, i.e., the time job search is re-started.
Wages. Hourly wages are computed by dividing weekly earnings by hours worked per week. If those are
not recorded on a weekly basis, such as monthly or annual, I convert them to weekly data. Wages are
deflated by monthly CPI (the average in 2004 is set to 100.) If an individual had a job that ended before theinterview week, I use the wage is observed in the last week of the job. Otherwise, I use the wage observed36

in the interview week. In both cases, I assume that the wage is constant during the quarter that includes
the week when wage is observed. When I need to move the starting date of the job as I explained above,
the wage observation is also moved. I keep the link between employment record and corresponding wage
observations.
Variables for Auxiliary Models. My auxiliary parameters are hazard rates and conditional mean and varianceof wage given experience and tenure. To compute the auxiliary parameters, I need full-time job spells, fulltime employment spells, non full-time employment spells, wages, and the corresponding experience andtenure.
Job spells, Full-time Employment spells, and Non Full-time Employment spells. Using the quarterly work
history array made above, I construct full-time job spells, full-time employment spells, and non full-time
employment spells. A full-time job spell ends when an individual switches to a new job or loses full-time
employment status (his new status is non full-time employment.) On the other hand, a full-time employment
spell ends only when an individual does not have a full-time job any more. Thus, the difference between job
spells and employment spells is whether they include job-to-job transitions. A non full-time employment
spell ends if an individual who does not have a full-time job finds a new full-time job.
Work Experience. All individuals start their career with zero experience and zero tenure, regardless of theiractual experience recorded in original NLSY work history file. I construct these variables using only myquarterly work history array. It is important to make my variables consistent with my model. Otherwise
I would be unable to interpret the data in terms of my model. As it is assumed in my theoretical model,
experience does not depreciate. So, experience is the cumulative number of quarters that an individual had
a full-time job in the past.

B.3.3. Sensitivity Checks. One major concern is a potential bias introduced by ignoring all part-time job
experience. I constructed another quarterly work history array that includes part-time jobs. The differencesfrom the original array are the following. First, if individuals have both a full-time job and a part-time job ina given week, only the full-time job is recorded in weekly work history array. Then I convert it to a quarterly
array. Second, I assume part-time experience accounts for a half of full-time experience.
Using this quarterly work history array I construct two different measures of work experience: one includes both full-time and part-time job experience, and another includes full-time job experience only. The
paths of the two different measures are reported in figure B.1. If I include part-time jobs, the average work
experience is 7.4 years ten years after graduating, while it is 6.9 years if I do not include part-time jobs.
37

F IGURE B.1. The difference in work experience when part-time jobs are not counted.

F IGURE B.2. Ratio of Part-time Jobs

Even 20 years after graduating, the difference is small and it is 0.8 years. This result suggests that the bias
from ignoring part-time jobs is negligibly small.
I also present the proportion of all jobs that are part-time in figure B.2. Part-time jobs are quite common
among the young. About a half of job holders worked in part-time jobs right after graduating. After a year,
the proportion of all jobs that are part-time quickly drops to 23% and it continues to decrease. This patternis consistent with the result found by Klerman and Karoly (1994). This result implies that part-time jobs arecommon only within a couple of years after graduating.
B.4. Numerical Issues.
38

B.4.1. Optimal Weighting Matrix. The optimal weighting matrix used in a simulated minimum distanceestimator is given by the covariance matrix of the sample auxiliary parameters. I first introduce someTinotations. Let Wi = {wi,t }t=1
be the wage history of individual i from the initial period (t = 1) to the lastperiod included in my sample (t = Ti .) Although subscripts for a job are omitted, this wage history includes
iwages in different jobs. Let diJ = {di,J j }Jj=1
be the vector of job spells that individual i have worked at.

F }Fi be the vector
The number of jobs that individual i had in my sample is denoted by Ji . Let diF = {di,kk=1

of full-time employment spells for individual i where Fi is the number of full-time employment spells ofN }Ni be the vector of non full-time employment spells for individual i where N
individual i. Let diN = {di,lil=1

is the number of non full-time employment spells of individual i. Finally, let Di = {Wi , diJ , diF , diN } be thework history of individual i and let D = {Di }Ni=1 be the whole sample.
The covariance matrix of the sample auxiliary parameters is estimated using nonparametric bootstrap
with blocking at the individual level. Specifically, I draw Di for N times with replacement. This is importantto preserve the dependent structure of the data. For example, I assume that wages are serially correlated dueto the individual fixed effect and the employer fixed effect. Another example is that full-time job spells andfull-time employment spells are positively correlated, because a long job spell necessarily implies a long
employment spell. Also, the model predicts that wages and job spells are positively correlated, because
workers with high match quality earn high wages and they are likely to remain in the job longer.
The B-th bootstrapped data set DB is used to compute the auxiliary parameters ρ̂B . I replicate the bootstrap
for 1,000 times and compute the covariance matrix of the auxiliary parameters. Although this computation
takes about an hour for a 2-GHz CPU machine, it is worth doing because the estimated covariance matrix
can be used for many different specifications of the structural model. This is very useful because I haveto estimate models with different prefixed parameters such as bargaining power parameter β for sensitivity
checks.
There are a few reasons that a nonparametric bootstrap should be used to estimate the optimal weighting
matrix. First, this matrix can be computed even when the analytical form is very complex. One such example
is the covariance between the coefficients of wage regression (e.g., βOLS ) and the hazard rates (e.g., pJ,t .)However I can still compute it through simulation in a straightforward way. Second, the weighting matrix
does not have to be updated over during the iteration in the estimation of structural parameters. This reduces
the computational burden.
B.4.2. Minimizing the Objective Function Value. The objective function is not smooth in the structuralparameters, because actions of simulated agents are discrete such as whether or not to change jobs. The
39

equilibrium wages of agents are also chosen from the discretized wage grid. Suppose a simulated agent
does not change his action if I change the structural parameter values from φ to φ + ε, but he changes theaction if I further change the parameter to φ + 2ε. This discrete response makes the objective function
non-smooth through the auxiliary parameters. It is true that this non-smoothness vanishes as the numberof simulations is increased up to infinity, but there are limitations for computational time and for computer
memory to store actions of simulated agents. I could not make the objective function smooth enough to use
a quasi-Newton method.
For the minimization of the objective function, I use a simulated annealing.10 The method is robust toa non-smooth objective function and it works well for my model. Alternative algorithm that can be used
for a non-smooth function is the Nelder-Mead simplex method. But, it converges to one of local minima,presumably because the surface of the objective function is too noisy. This problem still happens when thenumber of simulation is 100, but I could not increase it given my computational resources. The additional
advantage of simulated annealing is that it searches for the global minimum. The global optimum can befound as the number of iteration goes to infinity. The main disadvantage is that simulated annealing often
requires several ten thousands of iterations to achieve convergence. In particular, about 80,000 iterations areneeded for my model, which takes me three days.
B.4.3. Standard Deviation of the Structural Parameters. Non-smoothness of the objective function is also
problematic when I estimate asymptotic standard deviations of the structural parameters, because Jacobianof the objective function has to be computed. Remember that the objective function is not differentiable due
to discontinuous responses of simulated agents. The most formal method for this problem is to use a nonparametric bootstrap. However, it would take me a year for about 100 replications, which is not tractable. A
practical solution here is to compute the Jacobian using a finite difference method with assuming smoothness. This should work well when the number of simulations is large enough, because the discontinuity
vanishes as the number of simulations is increased to infinity. If this is not the case, the estimated standard
deviations can be extremely small particularly when the step size for derivative is small. The reason is thatthe finite difference method captures a noise on the surface of the objective function and that the absolutevalue of the computed derivative is very large.
Different step sizes are tested to see how the estimated standard deviations change. As a reference point,
the final step size of simulated annealing is used because it reflects the surface of the objective function
around the minimum. I will also check the validity of the method in a little informal way. I will estimate
10Goffe, Ferrier, and Rogers (1994) implemented a FORTRAN software that is downloadable from http://www.netlib.org.
40the structural parameters using nonparametric bootstrap with a small number of replications, say, ten times.
This should not be used to make a formal inference for the estimates. But it is worth doing to compare theresults obtained by a finite difference method and those obtained by a bootstrap. If these two sets of results
are very different, a different step size should be used.

R EFERENCES
A LTONJI , J., AND L. S EGAL (1996): “Small Sample Bias in GMM Estimation of Covariance Structures,”
Journal of Business and Economic Statistics, 14(3), 353–66.
BARLEVY, G. (2005): “Identification of Search Models using Record Statistics,” Federal Reserve Bank ofChicago.
BARRON , J. M., M. C. B ERGER , AND D. A. B LACK (2005): “Selective Counteroffers,” Syracuse University.
B OWLUS , A., N. K IEFER , AND G. N EUMANN (2001): “Equilibrium Search Models and The Transition
from School to Work,” International Economic Review, 42(2), 317–43.
C AHUC , P., F. P OSTEL -V INAY, AND J.-M. ROBIN (2006): “Wage Bargaining with On-the-Job Search:
Theory and Evidence,” Econometrica, 74(2), 323–64.
C ARRASCO , M., AND J.-P. F LORENS (2002): “Simulation Based Method of Moments and Efficiency,”
Journal of Business and Economic Statistics, 20(4), 482–92.
D EY, M., AND C. F LINN (2005): “An Equilibrium Model of Health Insurance Provision and Wage Determination,” Econometrica, 73(2), 571–627.
FARBER , H. (1999): “Mobility and Stability: The Dynamics of Job Change in Labor Markets,” in Handbookof Labor Economics, ed. by O. Ashenfelter, and D. Card. Elsevier Science Publisher (North-Holland),vol. 3.
FARBER , H. S., AND R. G IBBONS (1996): “Learning and Wage Dynamics,” Quarterly Journal of Economics, 111.
F LINN , C., AND J. H ECKMAN (1982): “New Methods for Analyzing Structural Models of Labor Force
Dynamics,” Journal of Econometrics, 18, 115–168.
G OFFE , W., G. F ERRIER , AND J. ROGERS (1994): “Global Optimization of Statistical Functions with
Simulated Annealing,” Journal of Econometrics, 60(1-2), 65–99.
H ARRIS , M., AND B. H OLMSTROM (1982): “A Theory of Wage Dynamics,” Review of Economic Studies,
49(3), 315–33.
41

K ENNAN , J. (2006): “A Note on Discrete Approximations of Continuous Distributions,” University of
Wisconsin-Madison.
K LERMAN , J. A., AND L. A. K AROLY (1994): “Young Men and the Transition to Stable Employment,”
Monthly Labor Review, 117(8), 31–48.
M URPHY, K. M., AND F. W ELCH (1990): “Empirical Age-Earnings Profiles,” Journal of Labor Economics,
8(2), 202–29.
N EAL , D. (1999): “The Complexity of Job Mobility among Young Men,” Journal of Labor Economics, 17,
237–61.
P OSTEL -V INAY, F., AND J.-M. ROBIN (2002): “Equilibrium wage dispersion with worker and employer
heterogeneity„” Econometrica, 70(6), 2295–2350.
P OSTEL -V INAY, F., AND H. T URON (2006): “On-the-job Search, Productivity Shocks, and the Individual
Earnings Process,” University of Bristol.
T OPEL , R., AND M. WARD (1992): “Job mobility and the careers of young men,” Quarterly Journal of
Economics, 107(2), 439–79.

42