Detecting Cyclical Turning Ponits Anas Ferrara 2004
Detecting Cyclical Turning Ponits Anas Ferrara 2004
Detecting Cyclical Turning Ponits Anas Ferrara 2004
2003 EDITION
Luxembourg: Office for Official Publications of the
European Communities, 2003
ISBN 92-894-3414-7
ISSN 1725-4825
Cat. No. KS-AN-03-033-EN-N
Detecting cyclical
turning points:
the ABCD approach and
two probabilistic indicators -
26th CIRET Conference,
Taipei, October 2002
E U R O P E A N
COMMISSION
1 THEME 1
General
statistics
Europe Direct is a service to help you find answers to your questions about the European Union
New freephone number:
00 800 6 7 8 9 10 11
A great deal of additional information on the European Union is available on the Internet.
It can be accessed through the Europa server (http://europa.eu.int).
ISBN 92-894-3414-7
ISSN 1725-4825
Introduction ................................................................................................................ 2
1 Detecting cyclical turning points: the ABCD approach ........................................... 2
1.1 The concept of detection................................................................................. 2
1.1.1 Detecting past turning points: dating the cycles
(turning points chronology) ................................................................... 3
1.1.2 Detecting present turning points: real time detection ............................ 3
1.1.3 Detecting future turning points: predicting turning points ...................... 4
1.2 The ABC approach ......................................................................................... 5
2 A couple of probabilistic indicators......................................................................... 9
2.1 A leading probabilistic indicator of the growth cycle........................................ 9
2.1.1 The Neftçi turning point detection method ............................................ 10
2.1.2 Construction of the leading probabilistic indicator................................. 13
2.2 A coincident probabilistic indicator of the business cycle ............................... 15
2.2.1 The Markov-Switching model of Hamilton (1989) ................................. 15
2.2.2 Construction of the coincident probabilistic indicator ............................ 18
3 Detection of the US and Eurozone cyclical turning points ..................................... 21
3.1 Decision rules and signals interpretation ........................................................ 21
3.2 Performance of the probabilistic indicators over the US cycle ........................ 22
3.2.1 Historical review of past cycles ............................................................. 22
3.2.2 Validation stage..................................................................................... 23
3.2.3 Real-time results on the last US cycle .................................................. 25
3.3 Performance of the probabilistic indicators over the Eurozone cycle.............. 29
3.3.1 Historical review of past cycles ............................................................. 29
3.3.2 Detecting the Eurozone growth cycle.................................................... 31
3.3.3 Detecting the Eurozone classical cycle................................................. 35
References ................................................................................................................. 40
26th CIRET Conference, Taipei, October 2002
Session:
Composite and leading Indicators / New methods
Abstract
The intricate issue of detecting and forecasting turning points of macroeconomic cycles
has been one more time well illustrated recently with the global downturn experienced by
most countries around the world in 2000-2001. Governments and Central Banks are very
sensitive to economic indicators showing signs of deterioration in order to be able to
adjust their policies sufficiently in advance to avoid more deterioration or a recession.
Those indicators require at least two qualities: they must be reliable and they must
provide a readable signal as soon as possible.
In this paper, we first discuss the concept of detection and propose the ABCD strategy of
the COE to identify the relevant cyclical turning points. Second, we introduce a couple of
indicators able to nowcast and to forecast those turning points. Both indicators are
probabilistic and are based on two different approaches. The first one is computed by
using the turning point detection algorithm of Neftçi (1984) and aims to forecast the
fluctuations of the growth cycle. The second one is grounded on the Markov-Switching
model proposed by Hamilton (1989) and is used to detect in real time peaks and troughs
of the classical cycle. The paper will review the performance of those indicators which
have been disseminated into the public by the COE since 1996. The analysis of those
leading and coincident indicators will particularly focus on the United States and the
Eurozone cyclical turning points.
JEL Classification:
∗
Centre d’Observation Economique, 27 avenue de Friedland, 75382 Paris Cedex 08,
France ([email protected] and [email protected])
26th CIRET Conference, Taipei 2
Introduction
in the short term. The present document will focus on the two last aspects: TPs
detection and prediction.
But let us review some specific issues when identifying successively past, present
and future TPs:
1.1.1 Detecting past turning points : dating the cycles (turning points
chronology)
In the United States, the NBER’s Business Cycle Dating Bureau’s Committee is
widely recognized as the authority for determining the peaks and troughs of the
classical business cycle (points B and C). However, there is a substantial delay before
the announcement of those dates. For example, the July 1990 peak was announced in
April 1991 and the March 1991 trough only in December 1992. More recently, the
March 2001 peak was announced in November 2001.
In other countries, there is no official dating of the classical business cycle. The
main issue is the definition of criteria used to recognize an economic fluctuation as a
cycle. The Conference Board refers to the 3D’s rule (diffusion, deepness, duration). If
dating the classical business cycle is not so easy, then dating the growth cycle is even
more difficult since the series must first be de-trended. Moreover, the way the series is
seasonally adjusted (directly or indirectly for geographic aggregates like Eurozone
indicators) and previously adjusted for calendar effects may impact on the datation (see
for example Astolfi et al., 2001 and Lommatzsh and Stephan, 2001). It may therefore
happen that different estimates are available on the market, we refer, for instance, to
Anas (2000) or Krolzig (2001) for Eurozone business cycle datations.
There have been many attempts to create an algorithm which would establish the
TPs dates. The most famous one is the Bry and Boschan (1971) procedure still in use
in many countries in order to estimate a series’ TPs (see, for example, Kim, Buckle and
Hall, 1995). Apart from those non parametric approaches, a great number of parametric
models has been developed lately, which could be useful to date the TPs of the
classical business cycle, based mainly on the Hamilton’s (1989) Markov-Switching
model. But a decision rule is still needed to identify the TPs as discussed by Harding
and Pagan (2001). In the case of switching regime models, the identification is
undertaken with a “natural” decision rule made on estimated smoothed probabilities of
the “hidden” regimes (probability higher than 50%). In this respect, a controversy has
recently emerged with Hamilton regarding the usefulness of using those sophisticated
models versus more transparent and simple methods for dating cycles.
If the use of quarterly GDP series may be sufficient to provide a dating of the past
TPs, it is clearly not operational in real time. GDP is only available on a quarterly basis
26th CIRET Conference, Taipei 4
with a delay of one to three months, sometimes with significant revisions. Thus, GDP is
not a good candidate to assess TPs in real time and the use of other series is
unavoidable. A solution is to use a GDP proxy commonly called a coincident index
(estimated by use of diverse linear methods). Stock and Watson (1989) were the first to
revive consideration on comovement of variables along the cycle by introducing a
dynamic factor model in order to extract a common factor. In this case, methods have
to be determined to estimate the probability of a TP of this common factor. In this
respect, Diebold and Rudebush (1996) recently proposed to mix together dynamic
factor models and regime switching (see also Kim and Nelson, 1998). If no coincident
indices are used, other generally non linear methods may directly produce the
probability of a TP in real time. For example, the multivariate Markov-Switching Vector
Autoregressive (MS-VAR) model proposed by Krolzig (1997) or the univariate Markov-
Switching model combined with a probability aggregation method proposed by the COE
and developed in this document (see section 2).
The timing of the prediction is very important. It is quite difficult to predict TPs in
the medium or long term (over 9 months). Even if economic imbalances sometimes
make an adjustment plausible or necessary in the future, it is difficult or even
impossible to predict when this adjustment will occur. In the short-term, however, the
TP prediction is, or should be, easier, except for important and sudden external shocks
(like the September 11th terrorist attack in New-York). Indeed, foreseeable changes in
economic policies should not reverse the course of economic development due to the
impact delay of these measures and the inertia of economic evolution. This is why it
may be useful to complement macroeconomic modeling with short-term leading or
coincident indicators.
We may distinguish three different ways to predict a TP:
1. A coincident index is elaborated and then projected through different techniques
(for example, VAR models in the Stock and Watson approach, comovement in
dynamic factor models). A method and a rule to detect the TPs of the projected
coincident indicator have then to be used: probit, Hamilton, Neftçi or ad-hoc rules
are possible candidates.
2. A leading index is elaborated and a method and a rule are applied to detect the
TP of that leading indicator. In this case, however, the average lead has to be
known.
3. A third approach consists in using multivariate models: detecting the TPs of
various leading indicators and aggregating the corresponding probabilities in order
to provide a signal for a future TP (Anas, 1997, and Anas and Nguiffo-Boyom,
2001).
26th CIRET Conference, Taipei 5
Figure 2 - Evolution of the US cycle over the 1986-1992 period and the ABCD
approach
26th CIRET Conference, Taipei 8
Figure 3 - Evolution of the Eurozone cycle over the 1987-1994 period and the
ABCD approach
26th CIRET Conference, Taipei 9
forecast the fluctuations of the growth cycle. First, we present the Neftçi’s sequential
algorithm.
The aim of this algorithm is to detect cyclical turning points in real time, which
mark the beginning or end of a cyclical downturn. For this purpose, Neftçi (1982)
developed a stochastic model for macroeconomic time series, based on the ingenious
work of Shiryayev (1978), to detect probability changes over processes. It is based on
the assumption that the series behaves differently depending on the downward or
upward regime in which it evolves.
Let us consider the stochastic process (Xt)t, where for all t, Xt represents the
observation on increments of the macroeconomic time series considered. According to
the finite sample (xt)t=1,…T, we will infer the occurrence or non-occurrence of a change in
the economic regime. Let Ζ (respectively Ζ’) be an integer-valued random variable
denoting the date following a peak (respectively trough)1. Let us suppose Z = i (or
Z’ = i), for i = 2,…,t, with T ≥ t ≥ 2, which means that a turning point has appeared
between dates i-1 and i. With the two following assumptions, we will be able to
characterize the cumulative distribution function.
• Assumption 1. The probability distribution of (Xi+j)j=0,1,2,…, is different and
independent of the distribution of (Xi-j)j=1,2,….
• Assumption 2. The realizations of the stochastic process (Xt)t between and within
regimes are independent.
If we consider that a peak appeared between dates i-1 and i, i.e. Z = i, with
T ≥ t > i ≥ 2, then we get:
P ( X 1 ≤ x 1 , K , X i ≤ x i , K , X t ≤ x t ) = F 1 ( x 1 , K , x i − 1 ) F 0 ( x i , K , x t ) , (1)
where F1(.) and F0(.) are the two cumulative distribution functions for the upward
and downward regime respectively. Generally, F1(.) and F0(.) are chosen to be
Gaussian cumulative distribution functions. The variable Z is not directly observable.
Based on historical values of (Xt)t, we intend to determine, at any date t, whether a
turning point has already occurred ( Z ≤ t ) or not (Z > t ) .
1
We suppose that Z refers to the date of peaks, but the results are diametrically symmetric for
troughs.
26th CIRET Conference, Taipei 11
Suppose that the practitioner has gathered some experience from the study of
past turning points and has subjectively defined a priori probabilities. Let T t be the a
priori transition probability of the change from upward to downward regime, i.e.
T t = P ( Z = t Z > t − 1) , (2)
and T t ' the a priori probability of the change from downward to upward regime, i.e.
Let us note x t = ( x 1 , K , x t ) the historical values of (Xt) since the last trough.
Given x t , let us evaluate at any date t the probability of occurrence of a turning point in
the recent past. Let Pt ( Pt ' ) denote the a posteriori probability of occurrence of a peak
(trough) at or before date t based on observations x t , that is:
Pt = P ( Z ≤ t x t ) . (4)
P (x Z ≤ t)P (Z ≤ t)
Pt = , (5)
P(xt )
and by extension:
P (xt Z ≤ t)P (Z ≤ t)
Pt = , (6)
P ( xt Z ≤ t)P (Z ≤ t) + P ( xt Z > t)P (Z > t)
Thus, by using this previous equation, the following Neftçi's formula is recursively
derived (for peaks (t ≥ 1 ) ):
Pt =
[Pt − 1+ (1 − Pt − 1 ) T t ] f 0 ( x t )
, (7)
[Pt − 1 + (1 − Pt − 1 ) T t ] f 0 ( x t ) + [(1 − Pt − 1 )( 1 − T t ) ] f 1 ( x t )
where f 0 (.) is the density function of (Xt)t during a downward regime and f 1 (.)
during an upward regime, and where P1=0.
We can see from equation (7) that Neftçi's formula allows to compute the a posteriori
probability of occurrence of a turning point, incorporating current information into the posterior
probabilities estimated over previous periods. As described in Niemira (1991): ''(...) the Neftçi
method accumulates probabilities from the start of the previous turning point. This particular
dynamic characteristic of the Neftçi method is a major improvement over its predecessors.'' This
26th CIRET Conference, Taipei 12
is an advantage over, for instance, a Probit approach which has poor dynamic contents and
may therefore be less powerful if the lead times are unstable.
The transition probabilities Tt and T’t in the Nefçi’s formula indicate the degree of
persistence of the process. For instance, Hamilton (1989) assumed that these
probabilities are constant overtime. However, recent works (see, for instance, Filardo,
1994 or Diebold, Lee and Weinbach, 1994) propose time-varying transition probabilities
as a function of the phase age or based on a leading indicator. Neftçi (1982)
considered that the transition probabilities were non constant and estimated them from
past experience.
Parameter estimation
The parameters of the probability distribution function of (Xt)t are estimated over
samples made of upward and downward regimes. The a priori transition probabilities
denoted T t ' and T t must also be estimated.
The first step consists in an a priori dating of the cycle peaks and troughs
(classical business cycle or growth cycle) of (Xt)t. The data are split into upward and
downward regimes in order to obtain two separate samples made respectively of
observations belonging to upward and downward regimes. This ex-ante determination
of peaks and troughs is done visually or by using automatic techniques based on ad-
hoc rules. For example, in Artis et al. (1995a), a method called ALT is designed and
applied to the reference series to verify its performance. In Artis et al. (1995b), a variant
of the maximum distance criterion used in discriminant analysis is developed. Those
techniques of « dating » a cycle are numerous (see, for example, Harding and Pagan
(2001) for an extensive review).
In the second step, the parameters of the probability density functions f0(Xt) and
f1(Xt) and the a priori probabilities (Tt’ and Tt) must be estimated. A few assumptions
are needed to estimate these parameters.
The probability density functions are estimated using an empirical distribution of
(Xt)t or by fitting a tabulate density function to observations of (Xt)t in each regime. In
the paper of Neftçi (1982), the density functions f0(Xt) and f1(Xt) are estimated by using
the empirical frequency distribution of (Xt) during upward and downward periods. On
the contrary, in the papers of Diebold and Rudebush (1989, 1991), Artis et al.(1995a)
and Anas (1997) a Normal distribution is fitted. In Diebold and Rudebush, the densities
are calculated first in a static way (1989) then dynamically (1991) to have an ex ante or
real-time evaluation of the performance of the famous CLI index. In Artis et al. (1995b),
a similar rolling-estimation technique is used .
Concerning a priori transition probabilities Tt’ and Tt, we may suppose that the
probability of a turning point is an increasing function of the age of the regime. In this
case, a priori transition probabilities are duration-dependent, as supposed in the paper
of Neftçi (1982). However, in Diebold and Rudebush (1989, 1991) evidence was
26th CIRET Conference, Taipei 13
Aggregation procedure
Suppose we selected N leading time series (Xkt)t for k=1,…,N (see the next
subsection for the choice of the series). For k=1,…,N, we associate a latent variable
(Skt)t such that, for all t, Skt = 1 if a turning point of the series (Xkt) has occurred before
date t and Skt = 0 otherwise. Moreover, consider a forecast horizon h, we note (Rt)t the
variable such that Rt = 1 if a cyclical turning point of the global economy occurs
between t and t+h (a peak for example) and Rt = 0 otherwise. We want to estimate the
value P(Rt = 1), for all t.
For each leading time series (Xkt)t the probability of an upcoming cyclical turning
point can be developed by using the bayesian formula as follows:
P ( R t = 1 ) = P ( R t = 1 S tk = 1 ) P ( S tk = 1 ) + P ( R t = 1 S tk = 0 ) P ( S tk = 0 ) . (8)
The two risks2 associated with this approach are, first α tk the risk of a false signal
(or type I error), defined as:
α tk = P ( R t = 0 S tk = 1 ) , (9)
2
Both risks are widely discussed in the seminal article of Okun (1960).
26th CIRET Conference, Taipei 14
and second, β tk the risk of missing the cyclical turning point3 (or second type
error), defined as:
β tk = P ( R t = 1 S tk = 0 ) . (10)
We assume that both risks are constant overtime, i.e. for all t; α tk = α k
and
β tk = β k . An estimate of P(Rt = 1) is Pk(Rt = 1) defined by:
P k ( R t = 1) = (1 − α k ) Pt k + β k
(1 − Pt k ) (11)
= β k
+ (1 − α k
− β k ) Pt k , (12)
where αk and βk are empirical estimates and where Ptk is the a posteriori
probability of an upcoming cyclical turning point given by the Nefçi’s formula (equation
(7)) applied to the variable (Xkt)t.
In the same way the diffusion indices are computed, we may consider a kind of
diffusion index of these probability estimates Pk(Rt = 1) given by equation (12) through
an aggregation procedure over the k leading series:
1 N
1 N
N
∑ [Pk ( R t
k =1
= 1) ] =
N
∑
k =1
β k
+ (1 −α k
− β k
) Pt k , (13)
N
(1 − α k
− β k
)P
= β + ∑
k =1 N
t
k
. (14)
Lastly, we decide to normalize the formula (14) so that it would equal 1 as soon as
all a posteriori probabilities equal 1. The final index we use, that we call IARC (in
French : “Indicateur Avancé de Retournement Conjoncturel”), which is an estimate of
P(Rt = 1) for all t, is such that:
β N
(
1−α k − β k
)
IARC = + ∑ N Pt k , (15)
1−α k =1
∑ 1−α
k
( )
k =1
where α and β are the averages of the type I and type II risks.
3
Either because the leading indicator missed the general economic TP or because the signal
was too late.
26th CIRET Conference, Taipei 15
Note that the minimum value of the IARC indicator is not 0 but β / (1 - α ),
because there always exist a risk of missing the turning point. For communication
purposes, the IARC indicator is put negative when in search for a trough (see section
3).
We use several criteria for selecting the leading components of the IARC
indicator. As outlined by the OECD, an economic rationale is needed to avoid taking
only the statistical performance over a period of time into account. Also, the series
need to be available for a long period of time to allow for estimation. But the main
selection criteria relates to the degree and stability of the operational lead as well as to
the low degree of first and second type risks. There is generally a trade-off: when the
lead increases, the risks increase at the same time.
where, for k=0,…,p, ak,St = ak,1 when St = 1, and ak,St = ak,2 when St = 2 and
where (εt)t is a white noise process with finite variance σ2.
Moreover, the whole specification of the Markov-Switching model needs the
specification of (St)t, as a 2-state first order Markov chain. That is, the value of the time
series St, for all t, depends only on the last value St-1, i.e.,for i,j=1,2,:
Lastly, it can be shown that the average length Lj of both regimes, for j=1,2, is
given by:
Lj = 1 / ( 1-pjj ). (19)
Parameter estimation
somewhat classical in the statistical literature, but in this case the main difficulty stems
from the fact that the latent process St cannot be observed and has therefore to be
estimated, for all date t. The MLE method aims to find the parameter θ such that the
conditional log-likelihood L(θ) is maximum, with L(θ) expressed as :
T
L(θ ) = ∑ log f ( x t / Ft −1 , θ ) , (21)
t =1
where, for all t, Ft denotes the vector of observations obtained through date t and
where f(xt / Ft-1, θ) is the conditional density of the MS(2)-AR(p) model, which can be
written as :
2
f ( x t / Ft −1 , θ ) = ∑ f ( x t / S t = i, Ft −1 , θ ) P ( S t = i / Ft −1 , θ ) , (22)
i =1
where f(xt / St=i, Ft-1, θ) is the conditional density of xt, assuming the current state
is known for each date t, given by, under the Gaussian assumption, :
p
1
f ( x t / S t , Ft −1 , θ ) = exp[( x t − a 0, St − ∑ a k , St x t − k ) 2 / 2σ 2 ] , (23)
2π σ k =1
Thus, by using equations (21) to (23) the log-likelihood L(θ) can be evaluated for a
given parameter θ. However, according to equation (22), the evaluation of the
conditional log-likelihood L(θ) asks for the knowledge of P(St = i / Ft-1,θ), for i=1,2. This
estimation is computed by using properties inherent to Markov chains: this is the
forecast of being in the state i given the information through date t-1. This estimated
probability P(St = i / Ft-1, θ), for i=1,2, is referred to as the filtered probability of being in
state i. This filtered probability will be saved in output to build our coincident business
cycle indicator. Note also that another conditional probability of being in the state i can
be computed, given all the available information through date T. This probability P(St = i
/ FT, θ), for i=1,2, is referred to as the smoothed probability of being in state i, often
used for recession dating procedures.
Model extensions
We describe now the way used to build our coincident probabilistic indicator,
starting from the filtered probabilities given by the Markov-Switching model applied to
the increments of diverse time series carefully chosen.
Aggregation procedure
The aggregation procedure is basically the same as the one developed to build
the leading probabilistic indicator IARC in the section 2.1.2. Assume we selected N
coincident time series (Xkt)t, for k=1,…,N (see the next subsection for the choice of the
series). For k=1,…,N, we associate a latent variable (Skt)t such that, for all t, Skt = 1 if
the series Xkt belongs to a low regime corresponding to a recession regime and Skt = 0
otherwise. Moreover, we define the variable (Rt)t such that Rt = 1 if the economy is in
recession and Rt = 0 otherwise. We want to estimate P(Rt = 1), for all t, which will
constitute our recession indicator.
26th CIRET Conference, Taipei 19
Similarly as in the subsection 2.1.2, for each coincident time series (Xkt)t the
probability of a recession can be developed by using the bayesian formula given by
equation (8). The two risks αk and βk associated with this approach are respectively the
risk of a false signal (type I error) and the risk of missing the business cycle turning
point (type II error), given respectively by equations (9) and (10). However, in this case,
we generally get P(Rt = 1 / Skt = 0) = 0, that is recessions are never missed. This is
understandable insofar as a recession is a main macroeconomic event widely diffused
all over the series. Thus, we estimate P(Rt = 1) by Pk(Rt = 1), for all t, defined by :
Pk ( R t = 1 ) = (1 − α k ) Pt k , (24)
(1 − α k )
N
SERI = ∑ N Pt k . (25)
k =1
∑
(1 − α k )
k =1
One of the main issue of this kind of indicator is the way to choose the different
components to include. Indeed, we search for series with a pretty strong persistence,
because volatility can lead unreliable signal, as well an ability of reaction in case of
recession, to provide a signal as soon as possible. For instance, regarding the United-
States, the series considered by the Business Cycle Dating Committee of the NBER
(industrial production, employment, real income and wholesale-retail sales) seem to be
potential candidates. These latter series are also integrated in some other economic
composite indicators (see for instance Stock and Watson (1993)). Moreover, as our aim
is to develop a monthly indicator to detect recession, the considered series have, of
course, to be sampled on a monthly basis. Therefore, series such as GDP or Eurozone
employment cannot be included in the indicator.
To discriminate among the huge set of economic monthly time series available in
data bases, we consider a criterion able to measure the goodness of recession real-
time detection of the series. The chosen criterion is the quadratic probability score
26th CIRET Conference, Taipei 20
(QPS) of Brier (1950), suggested for example by Diebold and Rudebusch (1989),
defined as follows:
1 T
QPS =
T
∑ (R
t =1
t − Pt ) 2 , (26)
4
See www.nber.org/cycles
5
see www.conference-board.org
26th CIRET Conference, Taipei 21
In this section, we apply our both indicators, IARC and SERI, for real-time turning
points detection to the US economy, then to the Eurozone economy. For each zone,
we validate in a first step the indicators through an historical analysis over the past,
then, in a second step, we present the results on the last cycles. First, we precise the
decision rules associated with both indicators and how the given signals have to be
interpreted by practitioners.
A decision rule is important to decide whether a turning point has occurred or not.
In our case, the decision rule is the determination of a threshold over which the
probability of a turning point is understood as a signal. For example, in the case of
Markov-Switching models with 2 regimes, the “natural” 50% threshold is recommended
by Hamilton (1989). However, it does not generally allow to avoid false signals and
must be therefore under-evaluated. Obviously, the signal is thus more reliable but also
more lagged. In the case of the US GDP series over the period 1952-84 studied by
Hamilton (1989), the low GDP volatility and the quasi absence of growth cycles in the
USA avoid this problem. Therefore, the question of the definition of a statistically robust
threshold is still open. In the Neftçi approach, a threshold of 90% or 95% is
recommended.
These decisions imply first and second type risks and they should be constructed
in view of the costs generated by these errors. Unfortunately, these errors are difficult
to assess so that, in practice, no decision rule is determined. It may be possible to
directly evaluate the probability estimates without using any decision rule. For example,
Diebold and Rudebush (1989) evaluate turning point forecasts on a number of
attributes like accuracy, calibration, resolution and sharpness, based on the works of
Winkler (1969). The accuracy statistics is the Brier’s (1950) QPS, analog to the mean
square error, given by equation (26). As another example, the diverse recession
indexes computed by Stock and Watson, and released monthly on their web site6, are
raw probabilities without decision rules associated with. Therefore, when the recession
probability reached 73% (as it was the case in June 2001), what is the conclusion to
draw ? In the COE approach, an empirical threshold is determined on the aggregated
probability based on past performances, comparatively to reference datations, in order
to make the understanding of the signal easier for decision-makers.
6
see Jame’s Stock web site : http://ksghome.harvard.edu/~.Jstock.Academic.Ksg/
26th CIRET Conference, Taipei 22
7
These crossings are interpreted as a signal of overheating or excessive slowdown.
26th CIRET Conference, Taipei 23
9600
8400
7200
6000
4800
3600
1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000
200
150
100
50
-50
-100
-150
-200
Baxter-King
-250
1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000
Let’s first focus on the start-end recession index (SERI). Over the period January
1972 – December 2000, the automatic selection procedure based on the QPS criterion
and an economic consideration lead us to keep the four following series:
• Unemployment rate of civilian workers
• Manufacturing industrial production index
26th CIRET Conference, Taipei 24
We consider that a series gives a signal of recession when the filtered probability
crosses the threshold value of 50%. None of the signals provided by these four series
miss a recession, but the signals are slightly lagged (around 1 month), excepted the
one stemming from construction spending which is advanced. Indeed, as a recession is
a main economic event, resulting in a significant decline widely diffused over the whole
economy, all the series considered always react when a recession occurs. In fact,
these four series present a high degree of persistence. However, they send sometimes
false signals. Especially, the IPI series send three false signals of business cycles
troughs. Therefore, this latter series possesses the weakest weight in the composite
indicator, according to the aggregation procedure. The most reliable component is the
unemployment rate, which never gives a false signal. However, this series is the most
lagged towards the reference business cycle. This phenomenon illustrates well the
trade-off existing between the lead and the reliability. According to our experience, it
appears extremely difficult to develop a reliable indicator able to predict with a
reasonable horizon when recessions will begin. The prediction of the end of a
recession should be easier, insofar as the economy presents the property of duration-
dependence during this phase (see Diebold et al., 1993). This means that the
probability of a recovery increases as the recession phase ages. In that case, a TVTP
Markov-Switching model could be fruitfully used.
The results provided by the SERI indicator are strongly coherent with the
reference recession dates given by the NBER. The average lag over a recession start
is 1.4 months, lags varying between 0 and 2 months, and the average lag over a
recession end is 1.6 months, lags varying between 1 and 2 months. We note that lags
are stable overtime. The SERI never emits a false signal of a recession start, but a
false signal of a recession end has been given during the 1974 recession, where all the
series switched almost simultaneously to a non-recession regime.
Let’s now consider the IARC index. The components have been chosen by
examining a set of series commonly considered as leading series, according to their
advance and their correlation with the reference growth cycle. We selected the six
following series:
• S&P’s 500 stock index
26th CIRET Conference, Taipei 25
• Yeld curve
• ISM survey
• Consumers expectations
• Stocks in the manufacturing industry
• New construction permits
Over the considered period, these series present cyclical evolutions in advance
over the GDP growth cycle, with an average lead of ten months, but with a detection
delay varying between four and seven months. All the six growth cycles have been
experienced by the series, although the movement is sometimes lagged. The financial
components and consumers expectations exhibit a large advance along with a high
volatility, thus generating some false signals. On the opposite, the three other
components are less advanced, but have a higher coefficient of correlation with the
reference growth cycle. In comparison with the series used in the SERI indicator, these
series present a higher degree of volatility.
The IARC indicator possesses an average lead of 2.5 months over the reference
growth cycle, but the signal is lagged at three times. Especially, the 1994-1995 cycle is
badly anticipated by the components. At the 80% level, the IARC do not provide any
false signal.
We are now interested in the last US cycle and how the IARC and SERI indicators
have helped to detect in real time and to predict TPs of the growth and classical cycles.
We focus first on the US growth cycle, that is the detection of points A and D of
the ABCD approach. In April 2000, the IARC crossed the 80% value which means that
a turning point of the growth cycle, a peak A in this case, would probably occur in the
next three months (see Table 1). Indeed, by using the most up-to-date GDP data, the
US growth cycle, estimated by applying a classical Baxter-King filter, shows a peak in
May 2000. When the peak was announced, it meant that the US growth rate would
decrease below its trend growth rate estimated at more than 3% at this time.
Knowing that a point A has been detected, let’s focus now on the peak of the
classical cycle, that is the detection of point B of the ABCD approach. In November
2001, the NBER Datation Commitee determined officially that a peak in business
activity occurred in March 2001, ending a period of 10 years of expansion beginning in
March 1991. This latter period is the longest expansion period identified by the NBER
and has been characterized by strong growth rates (around 3.5 % on the average over
this period).
26th CIRET Conference, Taipei 26
Table 1 - Evolution of the IARC indicator in the search for a peak of the
growth cycle from January 1999 to May 2000
janv-99 0,14
févr-99 0,16
mars-99 0,22
avr-99 0,26
mai-99 0,24
juin-99 0,22
juil-99 0,23
août-99 0,38
sept-99 0,50
oct-99 0,59
nov-99 0,50
déc-99 0,42
janv-00 0,44
févr-00 0,52
mars-00 0,68
avr-00 0,82
mai-00 0,85
The signal provided by the SERI indicator was very clear (see Table 2 and Figure
5). In March 2001, the SERI reached the value of 60% and thus crossed the threshold
value of 50%. This signal meant that the estimated peak in business activity (point B)
occurred during January 2001, taking into account the average delay of 1.4 months.
This signal has been confirmed the month after, in April 2001. This result is very close
to the NBER datation, however the SERI index of March 2001 has been released in
April 2001, thus 7 months before the NBER announcement. Of course, we do not
pretend to point out the lack of swiftness of the NBER, which aims to date and not to
detect in real time and which is under strong political and public pressures, but this
result underlines the reliability and the timeliness of our indicator.
The various series which compose the indicator provided signals at different
times. First, the industrial production index crossed the threshold value of 50% in
January 2001, which is not surprising insofar as the IPI is well known as a leading
indicator of the classical cycle. Then, the two series related to employment have
indicated the start of the recession in March and April 2001, which points out that
employment is coincident with the classical cycle. Lastly, the construction spending
series attained 55% in July 2001 and did not stay a long time in a recession regime.
This phenomenon reflects partly the resilience of the American households to decrease
their spending during this recession, which avoided a deeper crisis.
26th CIRET Conference, Taipei 27
Figure 5 - Evolution of the SERI indicator from January 2001 to May 2002
1 00 %
P ro b ab ility o f recessio n
5 0%
0%
a vr-0 1
ju il-0 1
o ct-0 1
a vr-0 2
ju in -0 1
n ov-0 1
d éc-0 1
ja nv-0 1
fé vr-0 1
m a i-01
ja nv-0 2
fé vr-0 2
m a i-02
a oû t-0 1
sep t-0 1
m a rs-0 1
m a rs-0 2
Once the point B was detected, we were looking for the point C, that is the end of
the recession period. The SERI crossed the 50% value in January 2002, which meant
that the trough of the classical cycle can be estimated to be in November 2001.
Assuming this date will be confirmed later by the NBER, this implies that this recession
lasted 9 months, a length close to the average length over the eight last recessions (11
months). The US recessions length is lying between six (in 1980) and sixteen (in 1974)
months.
In parallel, the IARC was used to detect the point D, i.e. the trough of the growth
cycle (see Table 3 and Figure 6). A signal was emitted in July 2001 (the IARC crossed
the threshold value of 80%), implying a trough in the next three months. This signal
was cancelled out by an unpredictable shock : the September 11th terrorists attack. As
a result, the values of the IARC from September to November 2001 declined under
80%. Finally, a signal of a trough was emitted by the IARC with the data of December
2001, implying thus that the point D would occur during the first quarter of 2002.
Indeed, if we look at the most up-to-date data concerning the US growth rate, it seems
that the trough of the US growth cycle can be dated to the first quarter of 2002.
26th CIRET Conference, Taipei 29
Table 3 - Evolution of the IARC indicator in the search for a trough of the
growth cycle from January 2001 to December 2001
janv-01 -0,44
févr-01 -0,50
mars-01 -0,53
avr-01 -0,68
mai-01 -0,73
juin-01 -0,80
juil-01 -0,84
août-01 -0,84
sept-01 -0,76
oct-01 -0,72
nov-01 -0,77
déc-01 -0,94
Source: COE, May 2002
Figure 6 - Evolution of the IARC indicator in the search for a trough (point D)
from January 2001 to December 2001
In the case of the Eurozone area, the datation of the business cycle (points B and
C) is not easy for different reasons.
• For some countries of the Eurozone, homogeneous series of quarterly GDP are
not available for a long period of time. The series of german GDP was starting in
1992 because of the reunification but a retropolation should be soon available.
26th CIRET Conference, Taipei 30
• In some series, the methodology has been changing overtime; for example GDP
for France is corrected for trading days since 1995. With new base years,
estimates are not strictly comparable along time.
• The method for seasonal adjustment or correction of trading days is not
homogeneous across countries. The adjustment may be done before or after the
aggregation (indirect versus direct method)
For dating the growth cycle (points A and D), there is a need to select the de-
trending method and to perform it before or after the aggregation over the countries
(indirect versus direct approach). We do not want here to discuss those issues. Let
mention the recent work of the Conference Board (see Zarnowitz and Ozyildirim, 2001)
showing that the TPs identification show great similarity when using the PAT approach
or the Hodrick-Prescott or band-pass (like Baxter-King) filtering methods, at least for
the United States.
In this paper, we will use a datation on a direct GDP aggregate for the Eurozone
derived from Eurostat data for the period starting in 1995 and using COE’s calculations
before. Therefore, this series is quite provisional. For estimating the growth cycle, we
use both above mentioned filters. Different turning points datation techniques are
available. An easy one is the Bry and Boschan (1971) approach, applicable directly on
the classical cycle or on the estimated growth cycle. The dating of the last growth cycle
is never easy because of edge effects. If we locate graphically the peaks and troughs
of the cycles, criteria of length, intensity and “trend crossing” should be used. The mid-
80’s and the 1998-99 cycle are mild ones because they do not cross the trend.
Zarnowitz and Ozyildirim (2001) indicated in their study that the 1998-99 cycle was
identified by the Bry and Boschan algorithm but not accepted by them because they
were uncertain that it could be qualified as a growth cycle contraction. Therefore, we
also could consider that the ascending phase lasted from the end of 1997 to mid-2001
without a real slowdown due to the Asian countries. In fact, the growth cycle is clearly
detected in Italy and in Germany while it is not in France, where it is known as the “trou
d’air”.
Over the period January 1972 – March 2002, the Eurozone economy experienced
eight growth cycles (see Figure 7). Three of them have been followed by a business
cycle. The growth cycle peaks of 1977, 1986, 1995, 1998 and 2000 (points A in the
ABCD approach) have not been followed by business cycle peaks (points B in the
ABCD approach). The average delay between points A and B was about 3 quarters in
the first recession (first oil shock) and around one quarter in the two following
recessions starting in 1980 and 1992. Regarding the last cycle, the peak may be
located in the third quarter of 2000 while the trough is not yet clearly established but,
according to our projections, should have occurred in the second quarter of 2002.
Despite a decrease of GDP in the last quarter of 2001, there was no recession in the
Eurozone during that period of time.
26th CIRET Conference, Taipei 31
120
112
104
96
88
80
72
64
56
1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002
-1
-2
Hodrick-Prescott
Baxter-King
-3
1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002
We first need to select the leading economic time series. Each leading series
used for the detection of the Eurozone turning points growth cycle is a weighted
average of member countries leading series. We have selected five leading series:
• a Eurozone stock market index
• the spread between long and short term interest rates
26th CIRET Conference, Taipei 32
• the first factor of a principal component analysis of business survey results in the
intermediate goods industry in the Eurozone
• manufacturer prices expectations concerning wholesale trade in the Eurozone
• the COE leading indicator of the American growth cycle.
For each series, a calculation of the operational mean lead has been performed
over the 1990-1999 period. The 1980's do not provide a good period for estimating the
performance of indicators because of the relative heterogeneity of economic growth
cycles in the zone.
The performance can be evaluated in Figure 8: a strong signal is always emitted
in the three months before the turning point.
26th CIRET Conference, Taipei 33
The Eurozone growth cycle upward phase ended at the end of the last quarter
2000. Actually, the peak can be found in the third quarter of 2000 according to the last
estimate of the cycle done in May 2002 (because of edge effect it takes at least one
year to locate a peak). A strong signal was given by the COE leading indicator in
October 2000 (see Figure 9), anticipating a peak in the three next months. This was
quite a good performance since at that time most economists were rather optimistic
about economic growth in 2001 and did not anticipate any peak in the growth cycle.
Figure 9 - COE leading indicator for the Eurozone: search for the next peak
in December 2000
More recently, in February 2002, the COE leading indicator for the Eurozone
entered the range indicating a strong probability of recovery within the next three
months, meaning that the Eurozone growth rate should be above its trend -which is
now assessed to be 2 per cent- in the third quarter of 2002 (see Figure 10). According
to the methodology used to build this indicator, it means that the centred three quarters
moving-average of quarterly changes should overpass 2% in annual terms in the third
quarter of 2002.
26th CIRET Conference, Taipei 35
Figure 10 - COE leading indicator for the Eurozone: search for the next
trough in April 2002
The direct use of GDP to detect the turning points of the classical cycle is not
operational because of the availability delay and the degree of revisions of GDP data.
Therefore, we must use a GDP proxy. We want to use a set of coincident series to get
a coincident signal of recession. The idea is to develop a Eurozone coincident
probabilistic indicator similar to the coincident indicator which was developed by the
COE for the American business cycle (see section 3.2). The series to be included in the
indicator must have an economic meaning. But other statistical restrictions also have to
be mentioned, such as fast availability, weak revisions, long enough history, etc...
As mentioned before, the employment series is a widely recognised coincident
indicator. The Eurozone employment series does not provide any false signal.
However, the signal is given with a certain delay. An inconvenient is the quarterly
frequency of the series and its delay of publication. Monthly series are only available for
some European countries.
Another series of interest, close to the employment series, is the unemployment
rate series. This series appears to give more or less the same results than the
employment series. The advantage is its monthly availability and the reasonable delay
of publication. Two periods of low regimes for the unemployment series are identified
by the MS model since 1982 (see Figure 11).
26th CIRET Conference, Taipei 36
11.0
10.5
10.0
9.5
9.0
8.5
8.0
7.5
1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002
Filtered probability
1.0
0.8
0.6
0.4
0.2
0.0
1982 1984 1986 1988 1990 1992 1994 1996 1998 2000 2002
We also have tried other series derived from some European surveys. Below are
presented the filtered probabilities of being in a low regime, stemming from a MS
model, for the industrial survey, restricted to the intermediate goods sector (see Figure
12), and for the construction sector (see Figure 13). Regarding the first survey, a
common factor has been derived from a principal component analysis on major
questions (tendency, orders and stocks). As concerns the construction sector, the
synthetic index calculated by the European commission has been used.
We observe that those surveys are very reactive to the economic climate.
However, regarding the classical business cycle, these series provide too much false
signals of recession. Actually, they are informative of the growth cycle.
26th CIRET Conference, Taipei 37
100
50
-50
-100
-150
-200
1985 1987 1989 1991 1993 1995 1997 1999 2001
Filtered probability
1.0
0.8
0.6
0.4
0.2
0.0
1985 1987 1989 1991 1993 1995 1997 1999 2001
26th CIRET Conference, Taipei 38
Construction Survey
8
-8
-16
-24
-32
-40
-48
-56
1977 1979 1981 1983 1985 1987 1989 1991 1993 1995 1997 1999 2001
Filtered probability
1.0
0.8
0.6
0.4
0.2
0.0
1977 1979 1981 1983 1985 1987 1989 1991 1993 1995 1997 1999 2001
would be the low-growth regime (positive growth but below the trend growth rate) and
an expansion phase where the growth rate is above the trend growth rate.
Finally, more series should be added to those series in order to produce an
aggregated signal. The aggregation rule described in section 2.2.2 will be applied. The
Eurozone coincident probabilistic indicator, under development in the COE, should be
soon available.
IPI Eurozone
120
110
100
90
80
70
60
1971 1973 1975 1977 1979 1981 1983 1985 1987 1989 1991 1993 1995 1997 1999 2001
Filtered probability
1.0
0.8
0.6
0.4
0.2
0.0
1971 1973 1975 1977 1979 1981 1983 1985 1987 1989 1991 1993 1995 1997 1999 2001
26th CIRET Conference, Taipei 40
References
Anas J. (1997), “Forecasting turning points of the French growth cycle”, paper
presented at 22nd Conference Ciret, Helsinki, Finland.
Anas J. (2000), “Le cycle économique de la zone euro: datation et détection”, in Les
Mutations de l’Economie Mondiale, ed. C. de Boissieu, Economica, Paris.
Artis M. J., Bladen-Hovell R., Osborn D., Smith G., and Zhang W. (1995), « Turning
point prediction for the UK using CSO leading indicators », Oxford Economic
Paper, vol. 47, pp. 397-417.
Astolfi R., Ladiray D. and Mazzi G.L.(2001), “Business cycle extraction of eurone GDP:
direct versus indirect approach”, Quarterly Journal of Economic Research, vol 3,
pp 377-398.
Bry G. and Boschan C. (1971), « Cyclical analysis of time series: Selected procedures
and computer programs », NBER, Technical Paper, n° 20.
Diebold F. X., Rudebusch G. D. and Sichel, D.E. (1993), ”Further evidence on business
cycle duration dependence”, Business Cycles, Indicators and Forecasting, ed. J.
Stock and M. Watson, University of Chicago Press for the NBER, pp. 255-280.
reprinted in : Business Cycles: Durations, Dynamics and Forecasting, Princeton
University Press, Princeton, 1999.
Diebold F. X., Lee J.-H. and Weinbach G. C. (1994), « Regime switching with time-
varying transition probabilities », in: Nonstationary Time Series Analysis and
Cointegration, C. Heargreaves (ed.), pp. 283-302, Oxford University Press,
reprinted in : Business Cycles: Durations, Dynamics and Forecasting, Princeton
University Press, Princeton, 1999.
Filardo A. J. (1994), « Business cycle phases and their transition », Journal of Business
and Economic Statistics,
vol. 12, pp. 299-308.
Hall, R. (2002), “The NBER’s Business Cycle Dating Procedure”, Business Cycle
Dating Commitee, NBER, January 10th.
Kim C. J. and Nelson C. R. (1998), « Business cycles turning points, a new coincident
index and tests of duration dependence based on a dynamic factor model with
regime switching », Review of Economics and Statistics, vol. 80, pp. 188-201.
Krolzig H.M., (2001), ‘Markov-Switching procedures for dating the Eurozone business
cycle’, Quarterly Journal of Economic Research, vol 3, pp 339-351.