Next Article in Journal
Financial Development and Countries’ Production Efficiency: A Nonparametric Analysis
Next Article in Special Issue
Bond Risk Premia and Restrictions on Risk Prices
Previous Article in Journal
Nonparametric Estimation of a Conditional Quantile Function in a Fixed Effects Panel Data Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stationary Threshold Vector Autoregressive Models

1
Department of Economics, University of Western Ontario, Social Science Centre, London, ON N6A 5C2, Canada
2
Department of Statistical and Actuarial Sciences, University of Western Ontario, Western Science Centre, London, ON N6A 5B7, Canada
*
Author to whom correspondence should be addressed.
J. Risk Financial Manag. 2018, 11(3), 45; https://doi.org/10.3390/jrfm11030045
Submission received: 15 June 2018 / Revised: 20 July 2018 / Accepted: 3 August 2018 / Published: 5 August 2018
(This article belongs to the Special Issue Financial Econometrics)

Abstract

:
This paper examines the steady state properties of the Threshold Vector Autoregressive model. Assuming that the trigger variable is exogenous and the regime process follows a Bernoulli distribution, necessary and sufficient conditions for the existence of stationary distribution are derived. A situation related to so-called “locally explosive models”, where the stationary distribution exists though the model is explosive in one regime, is analysed. Simulations show that locally explosive models can generate some of the key properties of financial and economic data. They also show that assessing the stationarity of threshold models based on simulations might well lead to wrong conclusions.

1. Introduction

Correct theoretical and empirical modelling of financial time series remains challenging. First of all, the usual linear framework often falls short of properly describing the data, which, instead, exhibit important nonlinear features. Secondly, economic theory regularly results in models with multiple equilibria and asymmetries which the time series model should be able to accommodate. Finally, data is often interconnected and hence simple univariate models generally fall short of appropriately describing the complex nature of the data. The Global Financial Crisis in 2007–2008 demonstrated this very clearly and reinforced the need to use a multivariate nonlinear framework in economic models, in general, and in empirical finance, in particular.
Among the many possible candidate nonlinear models, threshold models are particularly interesting and they have been extensively used in the existing empirical literature. These models are straightforward generalizations of linear models. For example, the simple two regime Threshold Autoregressive (TAR) model specifies a different autoregressive structure for each of the regimes and a threshold variable that determines which regime is active. These models are therefore relatively simple to estimate, and, since at time t the regime state is known, they are more suitable for forecasting than other nonlinear models, in particular hidden Markov models. Finally, TAR models allow for reasonably simple tests of the linear structure against nonlinear alternatives and to test the number of regimes. The multivariate generalization of the TAR model instead uses vector autoregressive structures in the regimes and is therefore naturally referred to as the Threshold Vector AutoRegressive (TVAR) model (Hubrich and Teräsvirta 2013; Tsay 1998).
Empirical studies have used threshold models to explore the asymmetry of shocks and nonlinear relationship between variables in financial markets and data from the real and monetary economy. For instance, TVAR models are widely used to study the asymmetric effect of fiscal and monetary policies in different credit, interest rate and inflationary regimes (Balke 2000; Fazzari et al. 2015; Shen and Chiang 1999). Balke (2000), for example, studies the propagation of shocks to output growth, the Federal Funds rate, inflation and measures of credit conditions during “tight” and “normal” credit market conditions using a TVAR framework with two regimes. The results suggest that shocks have a larger effect on output in “tight” credit regimes and that contractionary monetary shocks are more effective than expansionary ones. A similar approach is followed by Calza and Souza (2006) to study the transmission of monetary shocks across two credit regimes in the EU area and by Li and St-Amant (2010) to evaluate the effect of financial stress conditions on monetary policy effectiveness in Canada.
Another important application of threshold models has been to study the business cycle. For example, Altissimo and Violante (2001) study the joint dynamics of US output and unemployment using a bivariate TVAR model for recessions and expansions. Here, the lagged feedback variable, which measures the depth of the recession, defines the regime. The resulting model is a VAR with a fixed number of lags when the economy is in expansion and a time varying lag order when the economy is in recession. The authors find that nonlinearities are statistically significant only for unemployment, but it transmits to output through cross-correlation. Further evidence on the usefulness of threshold models for analysing the business cycle can be found in Koop and Potter (1999), Peel and Speight (1998), Koop et al. (1996), and Potter (1995), amongst others.
Threshold models are also popular in financial markets studies to explore the asymmetric relation between variables. In particular, a common application of TAR models includes determining the threshold effect in price movements related to transaction cost (Yadav et al. 1994). The threshold autoregressive conditional heteroskedastic class of models has been applied to study the nonlinear effect in volatility processes (Rabemananjara and Zakoian 1993). Finally, multivariate threshold models have been extensively used in studying the dynamics in stock prices, returns, volatilities, inflation and economic activity (Barnes 1999; Griffin et al. 2007; Huang et al. 2005; Li et al. 2015). Griffin et al. (2007), for example, study the joint dynamics of stock market turnover, returns and volatility in 46 countries using a TVAR model with two regimes that are separated by the sign of the past return. The authors conclude that small negative return shocks, rather than large ones, are the drivers for the decrease in turnover after a decrease in returns. Li et al. (2015) study the interaction between Shanghai and Shenzhen stock markets in a bivariate three regime TVAR model where the threshold variable is the average difference of the log returns between the two markets. Their results suggest that the strength of interaction between markets is regime dependent. In particular, the Shanghai market leads most of the time, except for the third regime, where both markets interact simultaneously. A detailed review of the application of threshold models in empirical economics can be found in Hansen (2011).
One challenge with nonlinear time series models, in general, and by extension therefore also with threshold models, in particular, is to assess model stationarity. Establishing stationarity is important as it is a fundamental assumption in most theoretical research. Indeed, the asymptotic properties of estimators in threshold models are generally established under a set of standard regularity conditions, which include the existence of finite higher order moments and the strict stationarity of the data generating process (Tsay 1998). Moreover, existing inference approaches assume stationarity of the data generating process (Hansen 1996, 2000; Tsay 1998) and violation of this assumption might lead to spurious nonlinearity (Calza and Souza 2006) and could invalidate the use of, e.g., Hansen (1996) simulated p-values for inference.
While significant progress has been made to establish conditions which ensure stationarity for the univariate threshold case under realistic assumptions (Brockwell et al. 1992; Chan and Tong 1985; Chen et al. 2011; Knight and Satchell 2011; Petruccelli and Woolford 1984), to the best of our knowledge, very little is known about the multivariate extension. See Chen et al. (2011) for an extensive review about recent findings regarding the stationarity of TAR models. If one was to use the general approach from this literature to establish the stationarity of TVAR models, it would require proving the convergence of an infinite sum of products of random matrices. This is clearly difficult and likely explains the absence of theoretical results for the TVAR model.
In this paper, we fill this gap in the existing literature and analyse the properties of the TVAR model in detail. To achieve this, we assume that the trigger variable is exogenous and that the regime process follows a Bernoulli distribution. We first derive necessary and sufficient conditions for second order stationarity, which are not present in the previous literature, when the variance–covariance matrices of the random vector and the error process are assumed to have full rank. Next, we characterize the joint conditional distribution of the data generating process when the error vector follows a multivariate normal distribution. Finally, we derive the unconditional distribution for a special case of the TVAR model, and we demonstrate that, in this case, the distribution of the threshold model is an infinite mixture of normals. This shows that TVAR models are very general and can accommodate many of the stylized features of financial data.
As a first interesting application of our results, we consider the special case where the elements of the random vector are positively correlated and we describe a model that is explosive in one regime, but still allows for the existence of steady state distribution. A similar idea was introduced in Knight et al. (2014) in the univariate case as a so-called “locally explosive model”. In particular, they study the univariate threshold autoregressive model with exogenous trigger and its application to bubble formation. We expand the notion of locally explosive models to the bivariate TVAR model. The derived conditions for the existence of the stationary distribution have simple economic intuition and are easy to interpret. In particular, our results show that, in the stationary model, there is a trade-off between autoregressive dependence in the regime and the probability of the regime.
Next, we conduct an empirical analysis of the locally explosive models. In the absence of explicit theoretical conditions that guarantee stationarity of the model, the previous literature suggested to establish stationarity indirectly by demonstrating, using a simulation study, that the estimated model does not appear to contradict the stability assumption. To assess this procedure, we simulate the bivariate locally explosive TVAR model for different distributions of the regimes. Our results show that a simulation study aimed at verifying stability of a particular model might give inconclusive or even wrong results. Specifically, we show that the simulation exercise may very well fail to reject stability of non-stationary TVAR models when the probability of the explosive regime is low.
Finally, we empirically document that the locally explosive TVAR model can be associated with bubble formation processes. In fact, our simulated locally explosive models appear to possess explosive and unit root behaviour while overall remaining stationary. These properties are implied by the definition of bubbles prevailing in the current literature and formally described by Evans (1991) and Phillips and Yu (2011). Our results, therefore, should encourage further research into multivariate threshold models and their use to study the formation of and existence of bubbles in financial data.
The structure of the paper is as follows: In Section 2, we derive the necessary and sufficient conditions for second order stationarity and for the existence of a stationary distribution for the TVAR model. This section also derives closed form solutions for the stationary distribution. In Section 3, we consider the so-called locally explosive models, in which the TVAR model is explosive in one regime, while overall remaining stationary. This section also presents some interesting special cases and reports the results from a simulation study. Finally, Section 4 concludes the paper. All proofs can be found in Appendix A and Appendix B contains additional figures.

2. The Threshold Vector Autoregressive Model

Throughout this paper, we consider the threshold vector autoregressive model given by
Y t = Φ 1 I ( X t 1 R 1 ) Y t 1 + Φ 2 I ( X t 1 R 2 ) Y t 1 + ϵ t ,
where Y t is a ( n × 1 ) random vector, Φ 1 and Φ 2 are ( n × n ) parameter matrices, where n is the number of time series, I ( ) is the indicator function, X t is an i i d random variable, which determines the regime, and ϵ t is a sequence of independent multivariate random vectors, such that E ( ϵ t ) = 0 and V a r ( ϵ t ) = Σ , t , where Σ is positive definite with full rank. We assume that E ( ϵ t | X s ) = 0 for all s t and that the sequence ( ϵ t , X t ) , t 1 , is i i d .
The regime process is defined as S t = I ( X t R 2 ) , t , where P r o b ( X t R 2 ) = π and P r o b ( X t R 1 ) = 1 π , with R 1 R 2 = R and R 1 R 2 = . From this, it follows that S t is an i i d Bernoulli variable with S t = 0 with probability 1 π and S t = 1 with probability π . Using S t , Equation (1) can be rewritten as Y t = ( Φ 1 + S t 1 Φ 0 ) Y t 1 + ϵ t , where Φ 0 = Φ 2 Φ 1 . If we further denote by B t = S t Φ 0 π Φ 0 , where E ( B t ) = 0 , t , the model in Equation (1) can be rewritten as a Random Coefficient Model (RCM) (see Nicholls and Quinn (1982)) given by
Y t = ( Φ + B t 1 ) Y t 1 + ϵ t ,
where Φ = Φ 1 + π Φ 0 = ( 1 π ) Φ 1 + π Φ 2 .
In the following sections, we examine the TVAR model specified above in detail. First, we provide the necessary and sufficient conditions under which the TVAR model is second order stationary. We also derive expressions for the moments and the stationary solution to the model given in Equation (1). Secondly, we derive the distribution associated with this data generating process. For simplicity, we assume only two regimes in Equation (1). However, the theoretical results obtained here can easily be generalized to multiple regimes.

2.1. Stationarity of the TVAR Model

Theorem 1 provides conditions under which the TVAR model above is second order stationary, i.e., that E ( Y t ) is constant and C o v ( Y t , Y t + h ) depends only on the lag h.
Theorem 1.
The process Y t , t = 0 , 1 , 2 , defined in Equation (1) is second order stationary with positive definite covariance matrix V = V a r ( Y 0 ) if and only if:
1.
μ = 0 , where μ is a mean of the initial vector, μ = E ( Y 0 ) ,
2.
the covariance matrix V solves V Φ V Φ E ( B t 1 V B t 1 ) = Σ , and
3.
λ < 1 , where λ is the maximum eigenvalue of the matrix ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 .
Proof. 
See Appendix A. ☐
Condition 3 of Theorem 1 provides a simple and intuitive eigenvalue condition for establishing second order stationarity of the TVAR model in Equation (1). We note that conditions similar to what we present here have been derived previously, though this has been done using alternative and, we believe, less realistic and empirically interesting assumptions (see, e.g., Nicholls and Quinn (1981), Feigin and Tweedie (1985) and Saikkonen (2007)). In particular, Nicholls and Quinn (1981), Feigin and Tweedie (1985) as well as Saikkonen (2007) assume independence of the error term, ϵ t , and the threshold variable, X t . We do not have this assumption and instead make a, we believe, far less restrictive assumption of the exogeneity of the threshold variable, X t , i.e., E ( ϵ t | X t ) , which does not rule out, for example, the possibility that the threshold variable X t could be a part of the ϵ t process.
Condition 2 of Theorem 1 provides an expression for calculating the covariance matrix of the second order stationary process Y t . Notice that, after vectorization of this expression, we can obtain a closed form formula for this. Remark 1 provides this formula.
Remark 1.
From vectorization of the expression V Φ V Φ E ( B t 1 V B t 1 ) = Σ , the equation for the variance of Y t can be obtained from
v e c V = ( I Φ Φ π ( 1 π ) Φ 0 Φ 0 ) 1 v e c Σ .
Theorem 2 provides an expression for the stationary solution to the model in Equation (1) and the corresponding conditions for the existence of this solution. Theorem 2 also shows that this solution is unique and strictly stationary.
Theorem 2.
Assume that V is positive definite with full rank. Then, the TVAR model in Equation (1) has a unique stationary solution given by
Y t = ϵ t + n = 1 k = 1 n Φ + B t k ϵ t n ,
if and only if λ < 1 , where λ is the maximum eigenvalue of the matrix ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 .
Proof. 
See Appendix A. ☐
In Remark 2, we provide the restriction on the eigenvalues of the matrix Φ , which is necessary for the stationary model in Equation (1) and follows from Theorems 1 and 2. This condition is more tractable, and it is used in Section 3 to simplify the analysis of the stationary TVAR model with one explosive regime.
Remark 2.
Let the process Y t , t = 0 , 1 , 2 , defined in Equation (1) be stationary with positive definite covariance matrix V. Then, the maximum eigenvalue of the matrix Φ is less than 1.
Proof. 
See Appendix A. ☐
The results of Theorems 1 and 2 can be extended to TVAR models with more than one lag. Corollary 1 presents the conditions for the stationarity of the TVAR model, which contains more than one lag.
Corollary 1.
Consider the following two-regime TVAR model with p lags in each regime
Y t = I ( X t 1 R 1 ) j = 1 p Φ 1 j Y t j + I ( X t 1 R 2 ) j = 1 p Φ 2 j Y t j + ϵ t ,
where the properties of X t and ϵ t are those following Equation (1). This model has a unique stationary solution given by
Z t = η t + n = 1 k = 1 n A + D t k η t n ,
where Z t and η t are n p × 1 vectors given by Z t = [ Y t , Y t 1 , Y t 2 , Y t ( p 1 ) ] and η t = [ ϵ t , 0 , 0 , , 0 ] , respectively, and D t = ( S t π ) A 2 + ( π S t ) A 1 , with A i , i = 1 , 2 , defined as
A i = Φ i 1 Φ i 2 Φ i 3 Φ i ( p 1 ) Φ i p I n 0 0 0 0 0 I n 0 0 0 0 0 I n 0 0 0 0 0 I n 0 ,
if λ < 1 , where λ is the maximum eigenvalue of the matrix ( 1 π ) A 1 A 1 + π A 2 A 2 , and only if λ 1 < 1 , where λ 1 is the maximum eigenvalue of the matrix A = ( 1 π ) A 1 + π A 2 .
Proof. 
See Appendix A. ☐
The distinctive feature of the TVAR model is that it is a linear Vector Autoregresive Model (VAR) in each of the regimes and an interesting question therefore is how the stability of each regime contributes to the stationarity of the whole TVAR model. Knight and Satchell (2011) investigate this question in detail for the univariate TAR model and Niglio et al. (2012) provide evidence that, when the univariate TAR model is stationary in both regimes, the whole TAR model cannot explode. The more interesting situation, however, occurs when the model in Equation (1) is explosive in one of the regimes.
The results of Theorems 1 and 2 can be used to analyse this particular situation, one in which the TVAR model in Equation (1) is explosive in one of the regimes. For example, the following example shows that the TVAR model can still be stationary in that case provided the probability to be in the explosive regime is not too large. See also Section 3 for further analysis.
Example. Consider the model in Equation (1), where Φ 1 = 0.70 0.21 0.31 0.80 , Φ 2 = 0.20 0.32 0.10 0.25 and π = P r o b ( X t R 2 ) = 0.3 . Since one of the eigenvalues of Φ 1 is equal to 1.01 , the model is not stationary in regime one. On the other hand, ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 = 0.36 0.12 0.12 0.0 0.16 0.41 0.06 0.14 0.16 0.06 0.41 0.14 0.07 0.18 0.18 0.47 , and its maximum eigenvalue λ = 0.78 . Thus, overall the model is stationary.

2.2. The Stationary Distribution

In this section, we describe the stationary distribution associated with the model in Equation (1). Throughout this section, we assume that ϵ t N ( 0 , Σ ) are independent random vectors. Let Y t be defined by Equation (4) and let S n ( t ) = k = 1 n ( Φ + B t k ) , n 1 , with S 0 ( t ) = 1 . It then follows that
Y t = ϵ t + n = 1 S n ( t ) ϵ t n .
From this, we have that
Y t | S n ( t ) N ( 0 , Σ + n = 1 S n ( t ) Σ S n ( t ) ) ,
and from the definition of S n ( t ) we notice that the stationary distribution of Y t is a complicated mixture of Normal distribution. Since it is difficult to establish the distribution of Y t in general, we will derive it under the assumption that Φ 1 = 0 . In this case, if Y t denotes returns, then prices follow a random walk in regime 1.
From Equation (8), we see that the characteristic function of Y t conditioned on S n ( t ) is given by
ϕ ( t , Y t | S n ( t ) ) = e x p 1 2 t Σ t 1 2 t n = 1 S n ( t ) Σ S n ( t ) t .
Notice that, when Φ 1 = 0 and Φ 2 = Ψ , then B t = ( S t π ) Ψ and Φ = π Ψ , and hence S n ( t ) = k = 1 n S t k Ψ . Note also that k = 1 n S t k Ψ Σ k = 1 n S t k Ψ = k = 1 n S t k Ψ Σ k = 1 n Ψ . The conditional characteristic function in Equation (9) therefore becomes
ϕ ( t , Y t | S n ( t ) ) = e x p 1 2 t Σ t 1 2 t n = 1 k = 1 n S t k Φ 2 Σ k = 1 n Φ 2 t .
Given the conditional characteristic function and the distribution of S n ( t ) , we can obtain the unconditional characteristic function and the marginal stationary distribution of Y t . The results are presented in Theorem 3.
Theorem 3.
The stationary distribution of the TVAR process with Φ 1 = 0 and Φ 2 = Ψ has the following characteristic function
ϕ ( t , Y t ) = ( 1 π ) K = 0 π K e x p 1 2 t n = 0 K Ψ n Σ Ψ n t .
Moreover, the probability distribution function is given by
f ( Y t ) = ( 1 π ) K = 0 π K N 0 , n = 0 K Ψ n Σ Ψ n ,
where N ( A , B ) is the multivariate normal distribution function with mean A and covariance matrix B.
Proof. 
See Appendix A. ☐
Theorem 3 is a generalization of a result for the univariate threshold autoregressive process developed in Knight and Satchell (2011) and shows that when Φ 1 = 0 the distribution function of Y t given in Equation (12) is an infinite mixture of multivariate Normals. This type of distribution can generate excess kurtosis. Such distributional characteristics are interesting when it comes to analysing financial markets and economic problems, since it can accommodate the special features of this type of data. For instance, the distributions of equity returns and typical measures of the realized volatility are characterized by large kurtosis. Thus, the theorem shows that TVAR models can be used to study these processes.

3. Locally Explosive TVAR Models

Threshold autoregressive models where one regime is non-stationary are related to the so-called locally explosive models. Knight et al. (2014) defines the locally explosive model as a model in which some regimes may be explosive, but the whole model has a stationary distribution. They study univariate threshold models and apply the idea of locally explosive models to investigate the formation of bubbles. In this section, we generalize the notion of locally explosive models to the bivariate setting. In order to do so, we need to link the stationarity of the whole model in Equation (1) provided in Theorems 1 and 2 to the stability of the model in each particular regime. Notice that the locally explosive models considered in this paper are models, which are state explosive. When X t = t instead the TVAR model is related to the models derived in Phillips and Yu (2009) and Phillips et al. (2011) where the explosive behaviour is defined in the time series context.
The derived conditions for the existence of a stationary solution are conditions on the matrix ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 , and it is not in general possible to relate the eigenvalues of this matrix to the eigenvalues of the parameter matrices Φ 1 and Φ 2 without adding extra structure. In the following section, we therefore consider a bivariate TVAR model, where the parameter matrices Φ 1 and Φ 2 have either positive entries only or are upper triangular. We first obtain the conditions on the parameter matrices under which the locally explosive TVAR model remains stationary. We next provide a simulation study to examine the characteristics of this model and show that graphically it is very difficult to assess model stationarity using simulated data.

3.1. Special Cases of the TVAR Model

We consider the special case where Y t in Equation (1) is bivariate and the parameter matrices have positive entries. We introduce the following additional notation for Φ 1 = ϕ 11 1 ϕ 12 1 ϕ 21 1 ϕ 22 1 and Φ 2 = ϕ 11 2 ϕ 12 2 ϕ 21 2 ϕ 22 2 . The following corollary to Theorems 1 and 2 provides conditions in terms of the individual ϕ ’s, under which the TVAR model is second order stationary. These conditions do not rule out the possibility of an explosive regime, and, if we assume that one regime is explosive, we derive the conditions on the coefficient matrix of the stationary regime.
Corollary 2.
Let the matrices Φ 1 and Φ 2 have positive entries. If ( 1 π ) ( ϕ j 1 1 + ϕ j 2 1 ) 2 + π ( ϕ i 1 2 + ϕ i 2 2 ) 2 < 1 , i , j = 1 , 2 , then the model in Equation (1) is stationary. Moreover, if the model in Equation (1) is explosive in one of the regimes x { 1 , 2 } , then ( ϕ i 1 x + ϕ i 2 x ) < 1 , i = 1 , 2 , where x { 1 , 2 } \ { x } .
Proof. 
See Appendix A. ☐
Corollary 2 shows that if the model in Equation (1) is explosive in one regime, the persistence of the variables in this regime is restricted by the probability of the regime and the persistence, defined as the column sum of the coefficients of Φ 1 or Φ 2 respectively, of the variables in the other regime. In other words, Corollary 2 states that there is a trade-off between how persistent a given regime can be and the probability of this particular regime. In addition, when the conditions of Corollary 2 hold and one of the regimes is explosive, the sum of the coefficients of the other regime’s matrix is naturally bounded by one.
Corollary 2 provides sufficient conditions for stationarity of the model, even when the underlying relationship is explosive in one of the regimes. We believe that the above finding might be useful for a number of financial and macroeconomic models. In fact, the assumption of positive entries only in Φ 1 and Φ 2 is not very restrictive for economic research and there are a variety of well documented cases with positive relationships between variables and their lags. For example, it is documented to be the case for asset returns and asset market illiquidity, consumption and GDP, volatility and trading volume and inflation and stock volatility, among many other pairs (Amihud and Mendelson 1986; Engle and Rangel 2008; Jagannathan et al. 2000; Wang and Yau 2000).
When we add slightly more structure and assume that Φ 1 and Φ 2 are triangular matrices with nonnegative diagonal entries, we can derive the necessary conditions directly in terms of the eigenvalues of Φ 1 and Φ 2 . Corollary 3 summarizes these findings.
Corollary 3.
Let the process Y t , t = 0 , 1 , 2 , defined in Equation (1) be stationary. Then, the following conditions hold:
1.
λ 1 2 λ 2 2 1 π ,
2.
λ 1 1 λ 2 1 1 ( 1 π ) ,
3.
λ 1 1 λ 2 2 1 ( 1 π ) π , and
4.
λ 1 2 λ 2 1 1 ( 1 π ) π ,
where λ 1 i and λ 2 i are the eigenvalues of the matrix Φ i , i = 1 , 2 .
Proof. 
See Appendix A. ☐
Since the eigenvalues of a triangular matrix are its diagonal entries, Corollary 3 could equivalently be stated as follows.
Corollary 4.
Let the process Y t , t = 0 , 1 , 2 , defined in Equation (1) be stationary. Then, the following conditions hold:
1.
ϕ 11 2 ϕ 22 2 1 π ,
2.
ϕ 11 1 ϕ 22 1 1 ( 1 π ) ,
3.
ϕ 11 1 ϕ 22 2 1 ( 1 π ) π , and
4.
ϕ 11 2 ϕ 22 1 1 ( 1 π ) π .
Corollaries 3 and 4 illustrate explicitly that there is a trade-off between how persistent a regime in the TVAR model can be and the probability of that regime while ensuring the overall stationarity of the process. Again, it is noteworthy that the stationarity of the TVAR model does not rule out the possibility of an explosive regime, but it restricts the value of the own autoregressive coefficients.

3.2. Simulation

Second order stationarity implies that means, variances and covariances are time-invariant and finite. If stationarity is not satisfied, however, it could be that shocks to the data generating process could lead to a time series that have unbounded moments. Previously, and in the absence of explicit stationarity conditions such as the ones derived in this paper, the literature instead suggested to verify that the estimated model does not contradict stability assumptions by use of simulation studies (Hubrich and Teräsvirta 2013). Specifically, the literature proposed to switch off the noise and simulate the estimated model for different histories. If the generated series converge to the same point, the natural conclusion would be that the simulated model is stationary. In contrast, finding at least one starting point that leads to an explosive time series would be sufficient to invalidate the stationarity assumption.
In this section, we perform a graphical analysis to “test” the stationarity of the TVAR model as suggested in the existing literature for different locally explosive TVAR models. Our results show that this “test” does not always allow us to draw the correct conclusion and the outcome of it is affected by the distribution of the explosive regime and the persistence of this regime. To be specific, we simulate the bivariate TVAR model in Equation (1) with different parameter values. We generate time series from the model of length equal to n = 250 , which is equivalent to one year of daily observations. The number of simulations is equal to m = 200 . The initial values of the time series, Y 0 , are equally distributed over the interval given by [ 0.15 , 0.24 ] for the first series and equally distributed over the interval given by [ 0.17 , 0.23 ] for the second series. In Appendix B, we report additional results when n = 2000 to check the robustness of our result.
The parameter values used in the simulation study are shown in Table 1. As the table shows, regime 2 is by construction always (locally) explosive and we vary the value of π , the probability of regime 2, such that the overall TVAR model can be stationary or non-stationary. This is indicated by the maximum eigenvalue of the matrix ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 , which is reported in column six labelled λ m a x . In particular, we define three groups of models, such that models within each group have the same coefficient matrices, but the probability to be in the explosive regime 2, π , varies.
Models 1–6 are more strongly related to lags in the explosive regime 2 than in regime 1. We contrast our models such that the persistence of the models in the second regime is stronger in group 2 than group 1. When the second regime is mildly explosive, like the models from group 1, this regime has to occur very frequently, in order to make the whole TVAR model non-stationary. In contrast, model 6 is unstable even when the probability to be in the explosive regime is as low as 30%. This confirms numerically that, when one regime is not stable, the distribution of the regimes is crucial for the stationarity of the whole TVAR model.
Figure 1 shows the simulated paths from models 1–3. When π is fairly low (Panel a), the time series appear stationary. When π gets larger and λ is closer to 1, the simulated model looks like a unit root (Panel b). While models 1 and 2 generate spikes, the simulated series return to the initial level all the time, a characteristic of stationary processes. When π = 0.7 (Panel c), the series is no longer stable and this is also evident from the figure. This conclusion is also valid when n = 2000 (see Figure A1 in Appendix B).
Figure 2 shows the simulated paths from models 4–6. These models are very persistent in regime 2 and they can generate huge spikes even when the probability of this regime is low (Panels a and b). Both models 4 and 5 look like unit root models, which explode, though they return to the initial level afterwards. Thus, the simulation exercise cannot reject stability of model 5, even though it is non-stationary by construction. The simulation study though does reject stability of model 6, when the probability to be in the explosive regime increases to 50% (Panel c). Thus, the results of the simulation might be misleading about non-stationary TVAR model with low probability of the explosive regime. The result of the simulation of models 4–6 prevails when n = 2000 (see Figure A2 in Appendix B). It is of course impossible to check all starting points in a simulation study and we might simply not be lucky enough to have a starting point that allows rejecting stability of the model.
Models 7–9 describe a type of relationship, where a particular time series is more strongly related with its own lag in regime 1 and with the other time series in regime 2. These models are quite persistent in regime 1, but still remain stationary in this regime. Figure 3 shows the simulated paths from these models. The explosive performance of model 9 is evident from Panel c. The conclusion, however, is not clear about model 8. This model is quite persistent in both regimes, thus it can generate growing series, like those shown in Panel b, and simulation of model 8 may in fact lead to rejecting the stability of a stationary model. However, we cannot reject stability of the model when n = 2000 (Panel b of Figure A3 in Appendix B). In fact, when the length of the simulated time series is increased to n = 2000 , the series from model 8 grows first but then returns to the initial level later on. Thus, the results from simulating model 8 show that the conclusion from this type of simulation study may also be sensitive to the sample size used in the simulation.
Figure 4 shows the simulated Y 1 t from TVAR models 4 and 5 specified in Table 1. We end this section by noting that the simulated series of Y t shown could be associated with data generating processes of financial or economic bubbles. Evans (1991) defines periodically collapsing explosive processes of bubbles such that the explosive behaviour of this process prevails through the whole sample, with non zero probability to collapse when it faces some threshold level. Phillips and Yu (2011) suggest a locally explosive process of bubbles, where asset prices transit from a unit root regime to an explosive regime and claim that this approach is consistent with other propagation mechanisms in financial markets like rational bubbles, exuberant responses to economic fundamentals and herd behaviour. Our simulation exercise shows that a simple bivariate locally explosive yet globally stationary TVAR model can generate unit root or explosive behaviour, which is consistent with these existing definitions of bubbles. An open question in the literature relates to how one can test for bubbles. In a recent paper, Ahmed and Satchell (2018) examine the performance of the Generalized Sup Augmented Dickey Fuller test proposed by Phillips et al. (2015) for the detection of explosive roots in univariate TAR models. They show that the power of the test drops considerably even though locally explosive regimes continue to be present when the process has a stationary distribution. We conjecture that this conclusion generalizes to the multivariate setting used in our paper.

4. Conclusions

This paper derives the necessary and sufficient conditions for the existence of a stationary distribution of the TVAR model with two regimes, when the regime process follows a Bernoulli distribution. These results are, to the best of our knowledge, unavailable in the existing literature. We further derive a closed form solution for the stationary distribution in the special case when there is no autoregressive structure in one of the regimes.
When the variables of interest are positively related, we describe a bivariate TVAR model, which is explosive in one regime, but allows for a stationary distribution along with finite moments. These results are related to so-called locally explosive models and our results extend the notion of locally explosive univariate processes to the bivariate case. We show that such models may remain stationary and, to ensure this, there is a trade-off between the persistence in a given regime and the probability of this regime.
In an empirical application, we simulate from various bivariate TVAR models, which are explosive in one of the regimes. We show how these models can capture the unit root and explosive behaviour, usually implied by the literature on bubble formation. We also demonstrate that a simulation study may fail to reject the stability of non-stationary TVAR models, when the probability of the explosive regime is low.

Author Contributions

The authors contributed equally to all aspects of this paper.

Funding

This research was funded by the Danish National Research Foundation grant number [DNRF78].

Acknowledgments

The authors would like to thank John Knight (deceased), Sergii Pypko, Charles Saunders, Tim Conley, Adlai Fisher, participants at the 2016 NFA Meeting, the 4th Annual Doctoral Workshop in Applied Econometrics at the University of Toronto and the University of Western Ontario Econometric Workshop and three anonymous referees for valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proofs

Proof of Theorem 1.
We first prove necessity. Let Y t be second order stationary such that E ( Y t ) = μ . Taking expectation of Equation (2), we have that μ = ( 1 Φ ) 1 0 = 0 . From Equation (2), we have that
( Y t μ ) ( Y t μ ) = ( B t 1 μ + ( Φ + B t 1 ) ( Y t 1 μ ) + ϵ t ) ( B t 1 μ + ( Φ + B t 1 ) ( Y t 1 μ ) + ϵ t ) .
Taking expectations on both sides of Equation (A1) and noticing that the expectation of cross products are zero, we have
V a r ( Y t ) = V a r ( B t 1 μ ) + V a r ( ( Φ + B t 1 ) ( Y t 1 μ ) ) + V a r ( ϵ t ) ,
which can be rewritten as
V a r ( Y t ) = V = V a r ( B t 1 μ ) + Φ V a r ( Y t 1 ) Φ + E ( B t 1 V a r ( Y t 1 ) B t 1 ) + V a r ( ϵ t ) ,
or
V = Φ V Φ + E ( B t 1 V B t 1 ) + Σ .
It then follows that,
V Φ V Φ E ( B t 1 V B t 1 ) = Σ .
By definition, Σ is a positive definite matrix of full rank. Conlisk (1974) and Conlisk (1976) show there is a unique positive definite V if and only if the maximum of the moduli of Φ Φ + E ( B t 1 B t 1 ) is less than 1. Notice that E ( B t 1 B t 1 ) = π ( 1 π ) Φ 0 Φ 0 and Φ Φ + π ( 1 π ) Φ 0 Φ 0 = π Φ 2 Φ 2 + ( 1 π ) Φ 1 Φ 1 . Thus, the conditions used in Conlisk (1974) and Conlisk (1976) transform to λ < 1 , where λ is the maximum eigenvalues of the matrix π Φ 2 Φ 2 + ( 1 π ) Φ 1 Φ 1 .
We now show sufficiency. Let conditions 1–3 hold. Taking expectation of Equation (2) at t = 1 shows that E ( Y 1 ) = E ( Y 0 ) = μ . Iterating further, it is possible to show that E ( Y t ) = μ , t . Similarly, calculating the variance of Equation (2) at t = 1 shows that V a r ( Y 1 ) = Φ V Φ + E ( B t 1 V B t 1 ) + Σ = V = V a r ( Y 0 ) . Iterating further, it is possible to show that V a r ( Y t ) = V , t . Since λ < 1 , it follows from Conlisk (1974) and Conlisk (1976) that V is positive definite. Premultiplying Equation (2) by Y t + h and taking expectations, we have
c o v ( Y t , Y t + h ) = ( π Φ 2 Φ 2 + ( 1 π ) Φ 1 Φ 1 ) c o v ( Y t , Y t + h 1 ) .
Iterating further, we have
c o v ( Y t , Y t + h ) = ( π Φ 2 Φ 2 + ( 1 π ) Φ 1 Φ 1 ) h c o v ( Y t , Y t ) = ( π Φ 2 Φ 2 + ( 1 π ) Φ 1 Φ 1 ) h V .
Thus, the process Y t , t = 0 , 1 , 2 , is second-order stationary. ☐
Proof of Theorem 2.
Let Y t be stationary and defined by Equation (4), i.e., the moments of Y t exist and they are finite. Then, from Equation (7), it follows that
E ( Y t Y t ) = E n = 0 S n ( t ) ϵ t n ϵ t n S n ( t ) .
We may rewrite this in v e c form as
v e c E ( Y t Y t ) = v e c E n = 0 S n ( t ) ϵ t n ϵ t n S n ( t ) = E n = 0 S n ( t ) S n v e c ( ϵ t n ϵ t n ) = E v e c ( ϵ t ϵ t ) + n = 1 k = 1 n ( Φ + B t n ) ( Φ + B t n ) v e c ( ϵ t n ϵ t n ) .
Since k = 0 n A k k = 0 n B k = k = 0 n A k B k for any matrices A k and B k whenever the matrix product exists, the later can be rewritten as
v e c E ( Y t Y t ) = n = 0 ( Φ Φ + E ( B t n B t n ) ) n v e c Σ .
Then,
v e c V = n = 0 ( π Φ 2 Φ 2 + ( 1 π ) Φ 1 Φ 1 ) n v e c Σ .
Furthermore, since ( π Φ 2 Φ 2 + ( 1 π ) Φ 1 Φ 1 ) v e c V = n = 1 ( π Φ 2 Φ 2 + ( 1 π ) Φ 1 Φ 1 ) n v e c Σ = v e c V v e c Σ , we have that
v e c V ( π Φ 2 Φ 2 + ( 1 π ) Φ 1 Φ 1 ) v e c V = v e c Σ ,
or
v e c V ( Φ Φ + E ( B t n B t n ) ) v e c V = v e c Σ ,
which is equivalent to
V Φ V Φ E ( B t n V B t n ) = Σ .
Since V and Σ are both positive definite, the maximum eigenvalue of Φ Φ + E ( B t n B t n ) is less than 1 (Conlisk 1974, 1976). Thus, the the maximum eigenvalue of ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 , λ , is less than 1.
We now prove sufficiency. Let all the eigenvalues of the matrix ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 be less than 1. Following Nicholls and Quinn (1982), we consider
W r ( t ) = ϵ t + n = 1 r k = 1 n ( Φ + B t k ) ϵ t n = n = 0 r S n 1 ( t ) ϵ t n .
Given that the eigenvalues of the matrix ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 are less than 1, the limit W ( t ) of W r ( t ) exists in mean square and thus in probability. Moreover, W ( t ) = ϵ t + n = 1 ( k = 1 n Φ + B t k ) ϵ t n satisfies Equation (2) and W ( t ) is stationary. Now, suppose U ( t ) is another stationary solution of Equation (2) and define
X ( t ) = W ( t ) U ( t ) .
By definition, X ( t ) = ( Φ + B t 1 ) X ( t 1 ) , E ( X ( t ) ) = 0 and X ( t ) is stationary. Then E ( X ( t ) X ( t ) ) = Φ E ( X ( t 1 ) X ( t 1 ) Φ ) + E ( B t 1 E ( X ( t 1 ) X ( t 1 ) ) B t 1 ) . Since X ( t ) is stationary, and Φ Φ + E ( B t 1 B t 1 ) = ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 , we have ( I ( ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 ) ) v e c E ( X ( t ) X ( t ) ) = 0 . However, since the eigenvalues of ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 are less then 1, E ( X ( t ) X ( t ) ) =0. Thus, W ( t ) = U ( t ) , and W ( t ) is the unique solution of Equation (2). Since W ( t ) is the same for all t, it is also the strictly stationary solution of Equation (2). ☐
Proof of Remark 2.
Following Theorems 1 and 2, the maximum eigenvalue of Φ Φ + E ( B t n B t n ) is less than 1. Consider S = Σ + E ( B t n V B t n ) and K = V . S and K are positive definite, since Σ and V are positive definite. From Equation (20), we have that
K Φ K Φ = S .
From Barnett and Storey (1970), it follows that the maximum eigenvalue of Φ is less than 1. ☐
Proof of Corollary 1.
Given the definitions of Z t , η t , and A i , i = 1 , 2 , we can rewrite the model in Equation (5) in its companion form
Z t = A 1 I ( X t 1 R 1 ) Z t 1 + A 2 I ( X t 1 R 2 ) Z t 1 + η t .
Now, define A = ( 1 π ) A 1 + π A 2 . Then, the model in Equation (A12) can be rewritten as a Random Coefficient Model (Nicholls and Quinn 1982) given by
Z t = ( A + D t ) Z t 1 + η t ,
where D t = ( S t π ) A 2 + ( π S t ) A 1 , such that E D t = 0 .
The proof of sufficient conditions are similar to the proof of Theorem 2 and it suffices to show necessity. Define Ω = v a r Z t and assume it exists and that it is finite. Notice that η ( t ) = ( 1 , 0 , 0 , . . , 0 ) ϵ ( t ) = l ϵ ( t ) . Define H = l l . Then, following the first part of the proof of Theorem 2, we have
v e c Ω = n = 0 ( A A + E D t n D t n ) n v e c ( H Σ ) .
Following the same strategy as in the proof of Theorem 2, it is straightforward to show that
Ω = A Ω A + E D t n Ω D t n + H Σ .
Let z = [ z 1 , . . z p ] be the left eigenvector of the matrix A with corresponding eigenvalue λ and z i are n × 1 vectors. Then,
( 1 λ 2 ) z Ω z = z 1 Σ z 1 + z E D t n Ω D t n z .
Since Ω is positive semidefinite, E D t n Ω D t n is positive semidefinite. Since Σ is positive definite, λ < 1 when z 1 0 . Now let z 1 = 0 . Since z is the left eigenvector of A with eigenvalue λ , we have the following system of equations:
z 1 Φ i + z i + 1 = λ z i , i = 1 , . . , p 1 ,
and
z 1 Φ p = λ z p .
Since λ 0 , z p = 0 . Thus, z i = 0 , i = 1 , . . , p 1 . However, since z 0 , this contradicts that z 1 = 0 . ☐
Proof of Theorem 3.
The characteristic function of Y t is
ϕ ( t , Y t ) = E S n ( t ) e x p 1 2 t Σ t 1 2 t n = 1 S n ( t ) Σ S n ( t ) t .
Thus, it is defined by the distribution of S n ( t ) . Since n = 1 S n ( t ) = n = 1 k = 1 n S t k Ψ , the probability space of n = 1 S n ( t ) is { 0 , n = 1 K Ψ n , K 1 } , and n = 1 S n ( t ) has a Geometric distribution with P ( n = 1 S n ( t ) = 0 ) = ( 1 π ) and P ( n = 1 S n ( t ) = n = 1 K Ψ n ) = ( 1 π ) π K . It follows that
ϕ ( t , Y t ) = E S n ( t ) e x p 1 2 t Σ t 1 2 t n = 1 k = 1 n S t k Ψ Σ ( k = 1 n S t k Ψ ) t = ( 1 π ) e x p 1 2 t Σ t + ( 1 π ) K = 1 π K e x p 1 2 t n = 1 K Ψ n Σ Ψ n t = ( 1 π ) K = 0 π K e x p 1 2 t n = 0 K Ψ n Σ Ψ n t .
Integrating the above expression, we have that the probability distribution function of Y t is a weighted average of normal distributions
f ( t , Y t ) = ( 1 π ) K = 0 π K N 0 , n = 0 K Ψ n Σ Ψ n ,
where N ( A , B ) is the multivariate normal distribution with mean A and covariance matrix B. ☐
Proof of Corollary 2.
The Perron–Frobenius theorem states that, for a matrix X with positive entries, there is a unique maximum eigenvalue λ such that min i j x i j λ m a x i j x i j . Let λ 2 be the maximum eigenvalue of matrix Φ 2 . Then, 1 λ 2 max i { 1 , 2 } ( ϕ i 1 2 + ϕ i 2 2 ) . Let λ be the maximum eigenvalue of ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 . Then,
λ max i { 1 , 2 } F i ( ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 ) ,
where F i ( ) denotes the column sum for each row i. From the last equation, we have that λ max i { 1 , 2 } F i ( ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 ) ( 1 π ) max i { 1 , 2 } ( ϕ i 1 1 + ϕ i 2 1 ) 2 + π max i { 1 , 2 } ( ϕ i 1 2 + ϕ i 2 2 ) . Since the condition of the corollary holds for any i , j = 1 , 2 , we have that
( 1 π ) max i { 1 , 2 } ( ϕ i 1 1 + ϕ i 2 1 ) 2 + π max i { 1 , 2 } ( ϕ i 1 2 + ϕ i 2 2 ) < 1 .
Thus, λ < 1 and, from Theorems 1 and 2, the model in Equation (1) is stationary. Now, suppose the model is explosive in regime 2 and let λ 2 be the maximum eigenvalue of matrix Φ 2 . Then, from the Perron–Frobenius theorem,
1 λ 2 max i { 1 , 2 } ( ϕ i 1 2 + ϕ i 2 2 ) .
Equations (A18) and (A19) show that 1 ( max i { 1 , 2 } ( ϕ i 1 2 + ϕ i 2 2 ) ) 2 < 1 π ( 1 π ) π max i { 1 , 2 } ( ϕ i 1 1 + ϕ i 2 1 ) 2 and thus ( ϕ i 1 1 + ϕ i 2 1 ) < 1 for any i = 1 , 2 . ☐
Proof of Corollary 3.
From Remark 2, we know that, if the model is stationary, then the eigenvalues of the matrix Φ are less than 1. Let λ 1 and λ 2 be the eigenvalues of matrix Φ = ( 1 π ) Φ 1 + π Φ 2 where Φ 1 = ϕ 11 1 ϕ 12 1 0 ϕ 22 1 , Φ 2 = ϕ 11 2 ϕ 12 2 0 ϕ 22 2 . Then it follows that,
λ 1 λ 2 = det Φ = ( 1 π ) 2 λ 1 1 λ 2 1 + π 2 λ 1 2 λ 2 2 + π ( 1 π ) ϕ 11 1 ϕ 22 2 + π ( 1 π ) ϕ 22 1 ϕ 11 2 .
Since ϕ 11 i and ϕ 22 i , i = 1 , 2 are nonnegative,
λ 1 λ 2 ( 1 π ) 2 λ 1 1 λ 2 1 + π 2 λ 1 2 λ 2 2 .
Now suppose λ 1 2 λ 2 2 1 and π 1 λ 1 2 λ 1 2 . Since Φ 1 and Φ 2 are upper triangular matrices with nonnegative diagonal entries, we have that
λ 1 λ 2 π 2 λ 1 2 λ 2 2 1 .
Thus, there is an eigenvalue of Φ , which is greater than 1. From Theorem 2, it follows that the process Y t , t = 0 , 1 , 2 , is not stationary. Similarly, it can be shown that the process Y t is not stationary when condition 2 of Corollary 3 does not hold. ☐

Appendix B. Additional Figures

Figure A1. Simulated paths represented by different colors from models 1–3 for a different set of histories over an n = 2000 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Figure A1. Simulated paths represented by different colors from models 1–3 for a different set of histories over an n = 2000 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Jrfm 11 00045 g0a1
Figure A2. Simulated paths represented by different colors from models 4–6 for a different set of histories over a n = 2000 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Figure A2. Simulated paths represented by different colors from models 4–6 for a different set of histories over a n = 2000 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Jrfm 11 00045 g0a2
Figure A3. Simulated paths represented by different colors from models 7–9 for a different set of histories over a n = 2000 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Figure A3. Simulated paths represented by different colors from models 7–9 for a different set of histories over a n = 2000 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Jrfm 11 00045 g0a3

References

  1. Ahmed, Muhammad Farid, and Stephen Satchell. 2018. What Proportion of Time Is a Particular Market Inefficient?... A Method for Analysing the Frequency of Market. Efficiency when Equity Prices Follow Threshold Autoregressions. Journal of Time Series Econometrics 10: 1–22. [Google Scholar] [CrossRef]
  2. Altissimo, Filippo, and Giovanni L. Violante. 2001. The nonlinear dynamics of output and unemployment in the U.S. Journal of Applied Econometrics 16: 461–86. [Google Scholar] [CrossRef]
  3. Amihud, Yakov, and Haim Mendelson. 1986. Asset pricing and the bid-ask spread. Journal of Financial Economics 17: 223–49. [Google Scholar] [CrossRef]
  4. Balke, Nathan. 2000. Credit and economic activity: Credit regimes and nonlinear propagation of shocks. Review of Economics and Statistics 82: 344–49. [Google Scholar] [CrossRef]
  5. Barnes, Michelle. 1999. Inflation and returns revisted: A TAR approach. Journal of Multinational Financial Management 9: 233–45. [Google Scholar] [CrossRef]
  6. Barnett, Stephen, and Colin Storey. 1970. Matrix Methods in Stability Theory. New York: Barnes and Noble. [Google Scholar]
  7. Brockwell, Peter J., Jian Liu, and Richard L. Tweedie. 1992. On the existence of stationary threshold autoregressive moving-average processes. Journal of Time Series Analysis 13: 95–107. [Google Scholar] [CrossRef]
  8. Calza, Alessandro, and Joao Sousa. 2006. Output and Inflation Responses to Credit Shocks: Are There Threshold Effects in the Euro Area? Studies in Nonlinear Dynamics & Econometrics 10: 1–19. [Google Scholar]
  9. Chan, Kung S., and Howell Tong. 1985. On the use of the deterministic Lyapunov function for the ergodicity of stochastic difference equations. Advances in Applied Probability 17: 666–78. [Google Scholar] [CrossRef]
  10. Chen, Cathy W. S., Mike K. P. So, and Feng-Chi Liu. 2011. A review of threshold time series models in finance. Statistics and Its Interface 4: 167–81. [Google Scholar] [CrossRef]
  11. Conlisk, John. 1974. Stability in a random coefficient model. International Economic Review 15: 529–33. [Google Scholar] [CrossRef]
  12. Conlisk, John. 1976. A further note on stability in a random coefficient model. International Economic Review 17: 759–64. [Google Scholar] [CrossRef]
  13. Engle, Robert F., and Jose Gonzalo Rangel. 2008. The spline-garch model for low-frequency volatility and its global macroeconomic causes. Review of Financial Studies 21: 1187–222. [Google Scholar] [CrossRef]
  14. Evans, George. 1991. Pitfalls in testing for explosive bubbles in asset prices. American Economic Review 81: 922–30. [Google Scholar]
  15. Fazzari, Steven M., James Morley, and Irina Panovska. 2015. State dependent effects of fiscal policy. Studies in Nonlinear Dynamics and Econometrics 19: 285–315. [Google Scholar] [CrossRef]
  16. Feigin, Paul D., and Richard L. Tweedie. 1985. Random coefficient autoregressive processes: A Markov chain analysis of stationarity and finiteness of moments. Journal of Time Series Analysis 6: 1–14. [Google Scholar] [CrossRef]
  17. Griffin, John M., Federico Nardari, and René M. Stulz. 2007. Do investors trade more when stocks have performed well? Review of Financial Studies 20: 905–51. [Google Scholar] [CrossRef]
  18. Hansen, Bruce. 1996. Inference when a nuisance parameter is not identified under the null hypothesis. Econometrica 64: 413–30. [Google Scholar] [CrossRef]
  19. Hansen, Bruce. 2000. Sample splitting and threshold estimation. Econometrica 68: 575–603. [Google Scholar] [CrossRef]
  20. Hansen, Bruce. 2011. Threshold autoregression in economics. Statistics and Its Interface 4: 123–27. [Google Scholar] [CrossRef]
  21. Huang, Bwo-Nung, Ming-Jeng Hwang, and Hsiao-Ping Peng. 2005. The asymmetry of the impact of oil price shocks on economic activities: An application of the multivariate threshold model. Energy Economics 27: 455–76. [Google Scholar] [CrossRef]
  22. Hubrich, Kirstin, and Timo Teräsvirta. 2013. Thresholds and smooth transitions in vector autoregressive models. Advances in Eonometrics 32: 273–326. [Google Scholar]
  23. Jagannathan, Ravi, Ellen R. McGrattan, and Anna Scherbina. 2000. The declining U.S. equity premium. Quarterly Review, Federal Reserve Bank of Minneapolis 24: 3–19. [Google Scholar]
  24. Knight, John, and Stephen Satchell. 2011. Some new results for threshold AR(1) models. Journal of Time Series Econometrics 3. [Google Scholar] [CrossRef]
  25. Knight, John, Stephen Satchell, and Nandini Srivastava. 2014. Steady state distributions for models of locally explosive regimes: Existence and econometric implications. Economic Modelling 41: 281–88. [Google Scholar] [CrossRef]
  26. Koop, Gary, and Simon M. Potter. 1999. Dynamic asymmetries in U.S. unemployment. Journal of Business and Economic Statistics 17: 298–312. [Google Scholar]
  27. Koop, Gary, M. Hashem Pesaran, and Simon M. Potter. 1996. Impulse response analysis in nonlinear multivariate models. Journal of Econometrics 74: 119–47. [Google Scholar] [CrossRef]
  28. Li, Fuchun, and Pierre St-Amant. 2010. Financial Stress, Monetary Policy, and Economic Activity, Working Paper. Ottawa: Bank of Canada.
  29. Li, Johnny Siu-Hang, Andrew C. Y. Ng, and Wai-Sum Chan. 2015. Managing financial risk in chinese stockmarkets: Option pricing and modeling under a multivariate threshold autoregression. International Review of Economics and Finance 40: 217–30. [Google Scholar] [CrossRef]
  30. Nicholls, Desmond F., and Barry Quinn. 1981. The stability of random coefficient autoregressive models. International Economic Review 22: 741–44. [Google Scholar]
  31. Nicholls, Des F., and Barry G. Quinn. 1982. Random Coefficient Autoregressive Models: An Introduction. Lecture Notes in Statistics. New York: Berlin/Heidelberg: Springer, vol. 11, pp. 1–154. [Google Scholar]
  32. Niglio, Marcella, Francesco Giordano, and Cosimo Damiano Vitale. 2012. On the Stationarity of the Threshold Autoregressive Process: The Two Regimes Case, Working Paper. Rome: Italian Statistical Society.
  33. Peel, David A., and Alan E. H. Speight. 1998. Modelling business cycle nonlinearity in conditional mean and conditional variance: Some international and sectoral evidence. Economica 65: 211–29. [Google Scholar] [CrossRef]
  34. Petruccelli, Joseph D., and Samuel W. Woolford. 1984. A threshold AR(1) model. Journal of Applied Probability 21: 270–86. [Google Scholar] [CrossRef]
  35. Phillips, Peter, and Jun Yu. 2009. Limit Theory for Dating the Origination and Collapse of Mildly Explosive Periods in Time Series Data, Working Paper. Singapore: Singapore Management University.
  36. Phillips, Peter, and Jun Yu. 2011. Dating the timeline of financial bubbles during the subprime crisis. Quantitative Economics 2: 455–91. [Google Scholar] [CrossRef]
  37. Phillips, Peter, Yangru Wu, and Jun Yu. 2011. Explosive behavior in the 1990’s NASDAQ: When did exuberance escalate asset values? International Economic Review 52: 201–26. [Google Scholar] [CrossRef]
  38. Phillips, Peter, Shuping Shi, and Jun Yu. 2015. Testing for multiple bubbles: Historical episodes of exuberance and collapse in the S&P 500. International Economic Review 56: 1043–78. [Google Scholar]
  39. Potter, Simon. 1995. A nonlinear approach to US GNP. Journal of Applied Econometrics 10: 109–25. [Google Scholar] [CrossRef]
  40. Rabemananjara, Roger, and Jean-Michel Zakoian. 1993. Threshold ARCH models and asymmetries in volatility. Journal of Applied Econometrics 8: 31–49. [Google Scholar] [CrossRef]
  41. Saikkonen, Pentti. 2007. Stability of mixtures of vector autoregressions with autoregressive conditional heteroskedasticity. Statistica Sinica 17: 221–39. [Google Scholar]
  42. Shen, Chung-Hua, and Thomas Chi-Nan Chiang. 1999. Retrieving the vanishing liquidity effect—A threshold vector autoregressive model. Journal of Economics and Business 51: 259–277. [Google Scholar] [CrossRef]
  43. Tsay, Ruey. 1998. Testing and modeling multivariate threshold models. Journal of the American Statistical Association 93: 1188–202. [Google Scholar] [CrossRef]
  44. Wang, George, and Jot Yau. 2000. Trading volume, bid-ask spread, and price volatility in futures markets. Journal of Futures Markets 20: 943–70. [Google Scholar] [CrossRef]
  45. Yadav, Pradeep K., Peter F. Pope, and Krishna Paudyal. 1994. Threshold autoregressive models in finance: The price differences of equivalent assets. Mathematical Finance 4: 205–21. [Google Scholar] [CrossRef]
Figure 1. Simulated paths represented by different colors from models 1–3 for different set of histories over a n = 250 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Figure 1. Simulated paths represented by different colors from models 1–3 for different set of histories over a n = 250 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Jrfm 11 00045 g001aJrfm 11 00045 g001b
Figure 2. Simulated paths represented by different colors from models 4–6 for different set of histories over a n = 250 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Figure 2. Simulated paths represented by different colors from models 4–6 for different set of histories over a n = 250 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Jrfm 11 00045 g002
Figure 3. Simulated paths represented by different colors from models 7–9 for different set of histories over a n = 250 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Figure 3. Simulated paths represented by different colors from models 7–9 for different set of histories over a n = 250 period using m = 200 simulated paths. The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Jrfm 11 00045 g003
Figure 4. Simulated paths of Y 1 t of model in Equation (1) represented by different colors with parameters from models 4 and 5 for over n = 250 period, where ϵ t N { 0 , 1 } The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Figure 4. Simulated paths of Y 1 t of model in Equation (1) represented by different colors with parameters from models 4 and 5 for over n = 250 period, where ϵ t N { 0 , 1 } The parameters are those from Table 1 and the probability to be in the explosive regime 2 is equal to π .
Jrfm 11 00045 g004
Table 1. Parameter values used in simulating the bivariate TVAR model.
Table 1. Parameter values used in simulating the bivariate TVAR model.
GroupModelRegime 1, Φ 1 Regime 2, Φ 2 Probability of Regime 2, π λ max
110.20.30.30.80.30.72
0.30.40.50.8
20.20.30.30.80.50.95
0.30.40.50.8
30.20.30.30.80.71.18
0.30.40.50.8
240.20.31.11.20.10.84
0.30.41.21.05
50.20.31.11.20.31.8
0.30.41.21.05
60.20.31.11.20.52.8
0.30.41.21.05
370.90.050.30.80.10.95
0.70.30.20.8
80.90.050.30.80.30.99
0.70.30.20.8
90.90.050.30.80.51.03
0.70.30.20.8
Notes: This table shows the parameter values used in the simulated TVAR models. The distribution of the regimes is Bernoulli with probability to be in regime 2 equal to π . Notice that regime 2 is not stable in any of the models. In the right-hand column, we report the maximum eigenvalue of the matrix ( 1 π ) Φ 1 Φ 1 + π Φ 2 Φ 2 , λ m a x .

Share and Cite

MDPI and ACS Style

Grynkiv, G.; Stentoft, L. Stationary Threshold Vector Autoregressive Models. J. Risk Financial Manag. 2018, 11, 45. https://doi.org/10.3390/jrfm11030045

AMA Style

Grynkiv G, Stentoft L. Stationary Threshold Vector Autoregressive Models. Journal of Risk and Financial Management. 2018; 11(3):45. https://doi.org/10.3390/jrfm11030045

Chicago/Turabian Style

Grynkiv, Galyna, and Lars Stentoft. 2018. "Stationary Threshold Vector Autoregressive Models" Journal of Risk and Financial Management 11, no. 3: 45. https://doi.org/10.3390/jrfm11030045

APA Style

Grynkiv, G., & Stentoft, L. (2018). Stationary Threshold Vector Autoregressive Models. Journal of Risk and Financial Management, 11(3), 45. https://doi.org/10.3390/jrfm11030045

Article Metrics

Back to TopTop