Abstract
We suggest a new method dealing with the problem of endogeneity of the threshold variable in structural threshold regression models based on copula theory. This method enables us to relax the assumption that the threshold variable is normally distributed and to capture the dependence structure between the threshold regression error term and the threshold variable independently of the marginal distribution of the threshold variable. For Gaussian and Student’s t copulas, this dependent structure can be captured by copula-type transformations of the distribution of the threshold variable, for each regime of the model. Augmenting the threshold model under these transformations can control for the endogeneity problem of threshold variable. The single-factor correlation structure of the threshold regression error term with these transformations allows us to consistently estimate the threshold and the slope parameters of the model based on a least squares method. Based on a Monte Carlo study, we show that our method is robust to nonlinear dependence structures between the regression error term and the threshold variable implied by the Archimedean family of copulas. We illustrate the method by estimating a model of the foreign-trade multiplier for seven OECD economies.
ACKNOWLEDGMENTS
We thank the editor (Todd Clark), an associate editor, two anonymous referees, Stelios Arvanitis, Yiannis Bilias, George Kapetanios, Dimitris Karlis, Andros Kourtelos, Andrew Patton, Peter Phillips, and Thanasis Stengos for comments on an earlier version of the article.
Notes
1 Threshold models have recently been gained growing interest (see Teräsvirta, Tjøstheim, and Granger Citation2010 for a survey). Applications include studies on business cycles (e.g., Potter Citation1995), the debt-growth relationship (e.g., Reinhart and Rogoff Citation2010; Kourtellos, Stengos, and Tan Citation2013), financial markets and volatility (see Tsay Citation2010 for a survey), monetary policy models (e.g., Davig and Leeper Citation2007; Kazanas and Tzavalis Citation2015).
2 To mitigate the effects of this assumption on the threshold parameter estimates, Kourtellos et al. (Citation2017) suggested a semiparametric estimator.
3 For a concise overview of applications of copula theory to economic and financial data, see, for example, Patton (Citation2006, Citation2012) and Fan and Patton (Citation2014). Recently, this theory has been applied to capture contemporaneous dependence between the error terms of two time series modeled separately, by the self-exciting threshold autoregressive model using as threshold variables lagged values of these series (see, e.g., Wong et al. Citation2017).
4 Note that, if δ = 1, then becomes
, which gives the joint distribution of
for
, which is the set of real numbers
. Analogously, we can obtain
, when δ = 0.
5 Note that, as in Heckman’s approach (1979), the selectivity bias correction terms in Lee’s (Citation1983) model are given by the inverse Mills ratios of the transformed distribution of the selection equation error term, truncated at the condition function values of the selection equation.
6 For instance, note that, apart from the Gaussian and Student’s t copulas, a closed form solution of
, i = 1, 2, can also be derived for the case that Ci is the Farlie–Gumbel–Morgenstern (FGM) copula frequently used in practice, which assumes nonlinear dependence between
and zt. As shown in Crane and Van Der Hoek (Citation2008), this can happen under the assumptions that
and
. Then, we have,
7 A similar problem appears in binary choice selectivity models (see Vella Citation1998 for a survey). Note that, for nonnormal distributions of xkt, the nonlinear pattern of the transformation
, or
, reduces the correlation between
and xkt. For instance, if
, where U is the uniform distribution, we can easily see that
where
, is the inversion function of the standard normal distribution. Based on a Maclaurin series expansion of
around
(i.e.,
),
can be written as
The last relationship shows more clearly that the degree of correlation between and xkt is also determined by the higher order terms in the above expansion of
which mitigates the problem of multicollinearity.
8 The discontinuity structure of distribution of zt, at
, adds extra randomness which, in a nonparametric framework of estimating model (1), also helps to identify the slope coefficients of the model across its two regimes (see Yu and Phillips Citation2018).
9 Note that this approach is different from that of Hansen (Citation2000)—see also Gonzalo and Pitarakis (Citation2002) and Kourtelos et al. (2016)—assuming that the threshold effects vanish with T, that is, as
10 Note that this approach is different from that of Hansen (Citation2000)—see also Gonzalo and Pitarakis (Citation2002) and Kourtelos et al. (2016)—assuming that the threshold effects vanish with T, that is, as
11 Note that similar results are obtained if we set the value of at the 1st quantile of the distribution of
12 Similar performance to the nonparametric method of estimating and
has also found for the ECDF method.
13 Kraay (2012) used a similar relationship, without threshold effects, to obtain estimates of the government spending multiplier.
14 In particular, REER is calculated from the nominal effective exchange rate (NEER) and a measure of the relative price of country i and its trading partners j, for all (e.g., Darvas Citation2012):
, where
is the geometric mean of the nominal bilateral exchange rates of a country i with its trading partners j, wij is the weight of trading partner j,
is the consumer price index of country i, and
is the geometric mean of the CPI’s of trading partners j.
15 Our earlier working paper version of this article (European Central Bank Working Paper No. 2136) shows graphs of the variables and lists their summary statistics. These are suppressed here for reasons of space.
16 Note that Canada can reject the linear relationship at a 15% level.
17 As shown by Roth (Citation2013) and Ding (2016) for the multivariate case, is also a Student’s t distribution. See also Kotz and Nadarajah (Citation2004).
18 Note that a column of this matrix also contains the values of the quantile
, or
, based transformed variable of regressor
dealing with its possible endogeneity.