Nonlinear Co-Integration Between Unemployment and Economic Growth in South Africa Andrew Phiri North-West University, South Africa phiricandrew@yahoo.com In this paper, a momentum threshold autoregressive (mtar) model is used to evaluate nonlinear equilibrium reversion between unemployment and economic growth for South African data between the periods 2000-2013. To attain this objective we estimate the first-difference and the gap model variations of Okun's specification. For the latter model variation, we employ three de-trending methods to obtain the relevant 'gap' data; namely, the Hodrick-Prescott (hp) filter, the Baxter-King (bk) filter and the Butterworth (bw) digital filter. A common finding from our empirical analysis is that Okun's law holds concretely for South African data regardless of the model specification or the de-trending technique that is used. Moreover, our analysis proves that unemployment granger causes economic growth in the long-run, a result which may account for the jobless-growth phenomenon experienced by South Africa over the last decade or so. Key Words: unemployment, economic growth, Okun's law, South Africa, mtar model, nonlinear unit root tests, nonlinear co-integration, nonlinear Granger tests, Hodrick-Prescott filter, Baxter-King filter, Butterworth high-pass filter jel Classification: C22, C51, E23, E24 Introduction High economic growth in conjunction with low unemployment under a low inflation environment can be deemed as the ultimate objective of macroeconomic policy in South Africa. Over the last decade or so, two prominent macroeconomic policy frameworks have embodied these objectives, those being, monetary policy's 'inflation-targeting' regime and fiscal policy's Accelerated and Shared Growth Initiative of South Africa (asgisa). Implemented in February 2002 and still in use to date, the inflation-target policy rule specifies that the South African Reserve Bank (sarb) should contain inflation at levels of between 3 and 6 percent, whereas the asgisa initiative seeks to halve unemployment and attain a Managing Global Transitions 12 (4): 303-324 6% economic growth rate by the year 2014. The assumed compatibility of the aforementioned policy objectives is inevitable demonstrated as monetary policy in South Africa is designated towards manipulating nominal variables like interest rates and inflation as a means of influencing real variables such as output growth and employment. Ultimately, the success of disinflation policy is reflected in its effect on unemployment and output growth. However, up-to-date South Africa has not only to managed to achieve arguably the highest economic growth rates in Africa since 1994, but the economy simultaneously boasts one of the highest youth unemployment rates in the world. So even though the South African Reserve Bank (sarb) can be credited for containing inflation within its set target which has been accompanied with steadily improved economic growth, such acquired growth has been characterized by what is popularly referred to as a 'jobless growth' syndrome (Hodge 2009). A mystery is warranted since the 'jobless growth' phenomenon contradicts the epic rise of unemployment caused by the sharp decline of real output experienced worldwide during the great depression. Therefore, a classical challenge for academics and policymakers alike is to provide an adequate account of unemployment-growth correlations in the South African economy. The question regarding the linkage between economic growth and unemployment gained prominence after Okun (1962) depicted the extent to which the unemployment rate is negatively correlated with output growth. By analyzing data over the period of 1947 to 1960, Okun (1962) documented that unemployment in the United States tends to fall by a one percentage point for every 3 percentage point rise in output growth. Thereafter, the United States was dubbed as having an estimated 'Okun coefficient' of 3 and a plethora of subsequent authors sought to estimate Okun's coefficient by either adopting a single-country approach (see Caraiani 2010; Ahmed, Khali, and Saeed 2011), panel-data approach (see Dixon and Shepard 2002,1997; Lal et al. 2010) or a multi-regional approach (see Freeman 2000; Adanu 2005; Villaverde and Maza 2009). The appeal of Okun's relationship is attributed to its simplicity and its extensive empirical support qualifies it to belong at the core of modern macroeconomics (Jardin and Gaetan 2011). As noted by Silvapulle, Moosa, and Silvapulle (2004), estimating the Okun coefficient has important implications for the business cycle since it relates the level of activity in the labour market to the level of activity in the product market. Whilst Okun's law implies that more labour is typically required for increased productivity levels, Okun's coefficient serves as an indication of the cost of unemployment in terms of output growth (Noor, Nor, and Ghani 2007). And in consolidation with the Phillips curve; Okun's relationship assists macroe-conomic policy in determining the optimal or desirable growth rate as a prescription for reducing unemployment (Moosa 1997). Overall, Okun's law is recommended as 'a rule of thumb' which provides policymakers with an understanding of how different markets adjust, and thus allowing for correct policies to be selected when facing shocks (Pereira, Bento in Silva 2009). In reality, Okun's law is more of a statistical relationship rather than a structural feature of the macroeconomy (Knotek 2007). The development of a pure theoretical foundation for Okun's relationship has been largely neglected in the academic literature, such that empirically, no functional form has been dominantly preferred to any other on the basis of theory (Weber and West 1996). As a consequence, the empirical examination of Okun's law is typically subject to revisions with the co-movement between output growth and unemployment frequently being analyzed under different settings. So while there is no contention on the importance of Okun's law, debates have evolved on the econometric techniques used to establish this relationship; how the cyclical components are extracted; and whether a dynamic or static specification is adopted (Turturean 2007). Recently, the possibility of asymmetric behaviour between economic growth and the unemployment rate has added a new dimension in the development of the academic literature. Take for instance Jardin and Gaetan (2011) who consider asymmetries in Okun's relationship as being important because asymmetric behaviour can adequately account for the varying effectiveness of structural and stabilization policies. Other commentators, such as Geldenhuys and Marnikov (2007), consider the impact of asymmetric behaviour on policy forecasting practices. In particular, these authors argue that if Okun's relationship is indeed found to be asymmetric, forecasts based on linear estimates of Okun's coefficient can lead to biased error terms. And yet another cluster of authors can also be identified, who advocate on the necessity of incorporating asymmetries in Okun's relationship as a means of reinforcing asymmetric behaviour in the Phillips curve. The rationale behind this line of thought is that if Okun's coefficient changes between regimes, then the sacrifice ratios should also change between regimes. In other words, different degrees of gradualism in the disinflation process may imply different im- pacts on unemployment for the same reduction in inflation (Beccarini and Gros 2008). Our study contributes to the literature by addressing the economic significance of asymmetric behaviour in Okun's relationship for South African data. To this end, our study makes use of the momentum threshold (mtar) autoregressive framework of Enders and Granger (1998). The logic behind the choice of our adopted approach can be described as follows. Engle and Granger (1987) argue that evidence of unit roots between a pair of time series variables necessitates the use of co-integration analysis prior to the estimation of any regression formed by the variables. According to the authors, the presence of co-integration would then imply that the variables follow a common long-run trend and the ols estimation of the time series will not yield spurious results. This is an important implication for our case study since previous empirical works have cautioned of unit root I(i) behaviour in output growth and unemployment variables for South African data (see Hodge 2006; Burger and Marnikov 2006; Gupta and Uliwingiye 2010). And yet it should also be noted that these conclusions are based on studies which assume a linear data generating process (dgp) among the series. Of recent, it has become widely accepted that standard unit root tests, suffer from low power when a linear approximation of an otherwise nonlinear time series is used to evaluate the integration properties of a time series (Enders and Granger 1998). A similar contention has risen for co-integration analysis, in which researchers like Enders and Dibooglu (2001) prove that the implicit assumption of symmetric adjustment is problematic if the adjustment towards long-run equilibrium is not linear. In particular, the authors argue that the presence of nonlinearities between a pair of time series signifies a high probability of nonlinear adjustment processes towards the long-run equilibrium for the data. With this in mind, our paper probes into the possibility of asymmetric behaviour between the unemployment rate and output growth using the mtar model. We choose this model because it represents a simple yet flexible framework that can simultaneously facilitate for (1) nonlinear unit root tests, (2) nonlinear co-integration analysis; and (3) nonlinear causality analysis. Therefore, against this backdrop, we present the remainder of the paper as follows. The following section of the paper presents the empirical framework of the study whereas section four presents the empirical results of the study. The paper is concluded in section five by providing policy recommendations and suggesting avenues for future research. Empirical Framework Our paper uses two classes of Okun's law specifications; namely, the first differences model and the gap model. To ensure that we obtain a balanced, robust view on the estimation results, we specify the Okun's specifications on both the direct and the reverse regressions of unemployment on output growth. For instance, in specifying the 'first differences' version of Okun's law, the link between the unemployment rate (ur) and economic growth (gdp) is represented as: ' Agdpt^ / \ ßi 0 / \ Aurt + / \ 61 K Aurt / , 0 ß2 , , AgdPt, where A is the first difference operator such that Agdpt = gdpt - gdpt-1 and Aurt = urt - urt-1. On the other hand, the 'gap model' measures these variables in terms of their deviations from long-run trends and is specified as: 'gdPts / \ ß1 0 ( C \ urt + / \ 61 ,urt, , 0 ß2 , , gdPt , where urct = urt - ur*t and gdpct = gdpt - gdp* are representative of the cyclical components of the unemployment rate and real output, respectively; with gdp* denoting a measure of potential output gap and ur* the unemployment gap variable. Having specified our baseline theoretical models, we can proceed to introduce co-integration analysis amongst the variables. We, therefore, take heed of Enders and Granger (1998) and model asymmetric adjustment between the unemployment and real output growth variables by allowing the residual deviations (i. e. gti) from the long-run equilibrium of regressions (1) and (2) to behave as a tar process. Formally, these residuals are modelled as follows: p 8%ti = kp&-i + (1 - It Wt-i + J]PiAgt-i + et. (3) i=1 In our paper, we identify four types of co-integration relations which govern the asymmetric dynamics within Okun's law, namely; tar with a zero threshold; consistent tar with a nonzero threshold; mtar with a zero threshold; and consistent mtar with a nonzero threshold. In the tar model with a zero threshold, the indicator function, It, is set according to: Under the tar model with a nonzero threshold, we set It, as: It 1, if \Çt-i > T 0, if Çt-i < T (5) where t is the value of the threshold variable. Enders and Granger (1998) suggest the use of a grid search procedure, as demonstrated in Hansen (1997), to derive a consistent estimate of the threshold i. e. the threshold estimate yielding the lowest rss is considered the true threshold estimate. The tar models are designed to capture potential asymmetric deep movements in the residuals if, for example, positive deviations are more prolonged than negative deviations (Enders and Dibooglu 2001). Enders and Granger (1998) and Caner and Hansen (2001) suggest that by permitting the Heaviside indicator function, It, to rely on the first differences of the residuals, A£t-1, a mtar version of equation (11) can be developed. The implication of the mtar model is that correction mechanism dynamic since by using A£t-1, it is possible to access if the momentum of the series is larger in a given direction relative to the direction in the alternative direction. In other words, the mtar model can effectively capture large and smooth changes in a series whereas the tar model shows the 'depth' of the swings in equilibrium relationship. In modelling mtar threshold co-integration with a zero threshold, the indicator function Mt, is set as: While in the mtar model with a nonzero threshold, Mt, is set as: For both tar and mtar specifications, Enders and Silkos (1998) demonstrate that a sufficient condition for stationary of gt-1 is that p1, p2 < 0. If gt-1 is found to be stationary, the least squares estimates of p1 and p2 have an asymptotic multivariate normal distribution for any given value of a consistently estimated threshold. Moreover, the null hypothesis of no (6) (7) co-integration (i. e. h01: p1 = p2 = 0) can be formally tested using a standard F-statistic for both tar and mtar models. If the null hypothesis of no co-integration is rejected, it is possible to test for the null hypothesis of symmetric adjustment (i.e. h02: p1 = p2) against the alternative of asymmetric adjustment (i. e. h12: p1 * p2) using a similar F-test. The empirical F-distribution for the null hypothesis; p1 = p2 = ois tabulated in Enders and Dibooglu (2001) whereas Enders and Siklos (2001) report critical values for testing the null hypothesis of p1 * p2. If both null hypotheses of no co-integration and no asymmetric co-integration can be simultaneously rejected, the granger representation theorem is satisfied and thus an associated error correction model can be estimated for the pair of time series variables. Thus in validating the presence of threshold co-integration, the error correction model can be modified to take into account asymmetries as in Blake and Fomby (1997). In our study we augment each of our threshold co-integration regressions with thresholds error correction specifications. In particular, the tar-tec model can be expressed as: Agdpt Aurt = AnIt fr-i + All (l - It )&-! P P + ^ aiiAgdpt-i + ßiiAurt-i. (8) i=i Whereas the mtar-tec model is specified as: Agdpt Aurt = AliMt£-i + An(i - Mt)£-1 PP + ^ aiiAgdpt-i + ^ßiiAurt-i, (9) i=i where the indicator functions for the tar and mtar model specifications are represented by It and Mt respectively. Through the above described systems of error correction models, two types of joint hypotheses can be tested. Firstly, the presence of asymmetries between the variables could initially be examined by examining the signs on the coefficients of the error correction terms. This involves testing the null hypothesis of H03: A^t-i = ^i2^t-i against the alternative H13: Ai&-i * ¿¿2&-i. The second type of hypothesis tested is that of granger causality effects which relatively examines whether all Agdpt-k and Aurt-k are statistically different from zero. Granger tests are used to examine whether the lagged values of one variable do not improve on the explanation or 'granger-cause' another variable. In particular, the null hypothesis that urt does not lead to gdpt can be denoted as: h04: a = o, i = 1,..., k; whereas the null hypothesis that gdpt does not lead to urt is: h05: fii = o, i = 1,..., k. All aforementioned hypotheses are based on a standard F-test. Furthermore, three types of joint hypotheses can be formed from the tec model. Firstly, granger causality tests can be implemented by testing whether all Agdpt-k and Aur— are statistically different from zero based on a standard F-test and if the A coefficients of the error correction are also significant. Empirical Analysis empirical data The data used in the empirical analysis consists of the annual percentage change in the real gross domestic product which is gathered from the South African Reserve Bank (sarb) online database whereas the unemployment rate for all persons aged above 15 years of age is collected from various issues of the quarterly labour force surveys (qlfs) as complied by Statistics South Africa (statssa). Our empirical analysis uses quarterly adjusted data obtained for the periods extending from 2000 to 2014. The choice of our sample period and periodicity reflects the limitations in the availability of the time-series data on unemployment and economic growth for South Africa. Although it would be desirable to employ a longer span of data, the available data provides the advantage of avoiding the issue of potential structural breaks related to South Africa's political and structural reforms such as those experienced in 1994. Moreover, we take note that while our data is relatively short, it is, however, up-to-date and further eliminates the problem of data unreliability associated with the South African unemployment series before 2000. Further given that gross domestic product is available on a quarterly basis and the unemployment rate is limited to half-yearly data, we use cubic spline interpolation to convert the half-yearly unemployment data into quarterly data over the same time period. We favour the use of cubic spline interpolation over other time series data conversion techniques due to its computational accuracy and stability of computation. Moreover, cubic spline interpolations satisfy the further condition at the end point. As a part of our data construction, we introduce the de-trending meth- ods used to extract the 'potential output' and 'unemployment gap' variables necessary to estimate the gap version of Okun's specification. The construction of these 'gap variables' is necessary since there exists no observable data on the trend components of the unemployment and output growth variables. Also taking into consideration that a majority of these de-trending techniques are not without scepticism, it is standard practice to apply a variety/different de-trending techniques to ensure robustness in the regressions analysis. Therefore in following along this course of reasoning, our study considers three alternative de-trending techniques, namely the Hodrick-Prescott (hp) filter; the Baxter-King (bk) filter and the Butterworth (bw) digital filter as respectively introduced by Hodrick and Prescott (1997), Baxter and King (1999); and Pollock (2000). The purpose of using these three de-trending techniques is to enable a robust analysis concerning the sensitivity of the estimated Okun's coefficient to the different choices of our gap variable estimates. UNIT ROOT TESTS In testing for unit roots, we begin on the simple premise of subjecting a univariate time series, yt, to the following generalized autoregression: Yt = 0) + ut, (11) ayt = st (st-1 < t) + st (st-1 > t) + ut, (12) Whereas the mtar version of the unit root test regression with a zero threshold and a consistent threshold estimate threshold are, respectively, specified as: ayt = st(ast-1 < 0) + st(ast-1 > 0) + ut, (13) table 1 Nonlinear Unit Root Tests Variable Model Lag Asymmetry test Unit root test Decision (i. e. pl = p2) (i. e. p1=p1 = 0) gdp TAR 2 0.94 12.63*** Linear I(o) (3-32)* (16.46)*** Nonlinear I(o) C-TAR 2 3.94* 15.59*** Nonlinear I(o) (7.87)* (21.28)*** Nonlinear I(0) MTAR 2 0-95 12.13*** Linear I(0) (9.46)** (22.89)*** Nonlinear I(0) C-MTAR 2 4.90* 16.03*** Nonlinear I(0) (6.67)* (19.96)*** Nonlinear I(0) ur TAR 0 2.45 2.86* Linear I(0) (4.96)* (7.22)** Nonlinear I(0) C-TAR 0 2.37 2.81* Linear I(0) (5.21)* (7.40)** Nonlinear I(0) MTAR 0 2.59 2.94* Linear I(0) (3.44)* (6.17)** Nonlinear I(0) C-MTAR 0 2.70 3.00* Linear I(0) (3.37)* (6.12)" Nonlinear I(0) notes Significance level codes: ***, ** and * denote the 1%, 5% and 10% significance levels respectively. Tests statistics for the first differences of the variables, i. e. Agdpt and Aurt are given in parenthesis. Ay = St (Ast- < t) + St (Ast- > t) + u. (14) Thereafter, two hypotheses can be formed from regressions (ii)-(i4). The first hypothesis tests for asymmetries within the time series. To this end, we test the null hypothesis of no asymmetric effects as h00: p1 = p2 against the alternative hypothesis of an asymmetric data generating process (i.e. h01: p1 ± p2). Subsequent to testing for asymmetric effects, we then proceed to test for unit root behaviour within the time series. Pragmatically, the null hypothesis of a unit root is tested as h10: p1 = p2 = 0 against the alternative hypothesis of an otherwise stationary asymmetric process (i.e. h10:p1 + p2 ± 0). The aforementioned tests of asymmetry and unit root behaviour are performed on time series variables of economic growth and the unemployment rate. The lag length of the threshold models which facilitate these tests are determined by the aic information criterion. As is evident from table 1, the empirical test results obtained for the time series in their levels are quite mixed. For instance, in scanning through the model tests conducted on the unemployment variable, we find that we cannot reject the null hypothesis of a symmetric process and yet we are able to reject the null hypothesis of a unit root process for same time series. Thus for the unemployment variable in its levels, we conclude a linear, stationary data generating process for the series. However, for the output growth variable in its levels, we conversely find that the c-tar and c-mtar versions of the employed tests simultaneously reject both null hypotheses of symmetry and unit root behaviour. This particular result implies a nonlinear, nonstationary data generating process for the output growth variable in its levels. And yet, in turning to the empirical results obtained for the time series in their first differences, our analysis reveals a common finding of a nonlinear yet stationary process for all variables under all model specifications. All in all, we can conclude that all utilized time series appear to be both nonlinear yet stationary processes in their first differences. Therefore, the results obtained from our preliminary unit root analysis paves the way for the threshold co-integration analysis which we conduct next. co-integration analysis Having investigated the integration properties of the unemployment and economic growth variables, we proceed to investigate threshold co-integration and error correction effects amongst the times series. However, prior to estimating any threshold models, we must first test a number of hypotheses to select which models best capture asymmetric behaviour in Okun's specification. To this end, we employ three threshold tests which have been previously discussed previously discussed. To recall, (1) we test for co-integration effects; (2) we test for threshold co-integration effects and (3) we test for threshold error correction effects. The results of these tests are reported in table 2. In referring to these results, we find that at least one type of threshold model manages to reject all three hypotheses at least a 10 percent significance level for all variations of Okun's law. This is quite an encouraging result since it implies that the data displays at least one form of nonlinearity for each version of Okun's specification. Another interesting result is that the mtar specification is most suitable for modelling nonlinear behaviour between unemployment and economic growth for South African data. The only exception holds for the cf filter estimates which favour a tar model specification. Furthermore, all estimated versions of Okun's law unveil significant asymmetric co-integration behaviour only when output growth is placed as the dependent variable in the regression. In summing up the test results reported in table 2, we can draw two broad conclusions thus far. Firstly, our analysis infers significant asymmetric behaviour between unemployment and economic growth for South African data. In this respect, our results adhere with those obtained in Geldenhuys and Marnikov (2007). However, in slightly differing from Geldenhuys and Marnikov (2007), we find smooth nonlinear adjustment behaviour in the data as opposed to an abrupt one. This result is expected since the otherwise abrupt nonlinearity is most suited for data containing structural break periods. Seeing that our data does not cover such periods, it therefore becomes reasonable that we detect smooth nonlinear behaviour among the data. Our second conclusion is that we establish economic growth as being the driving variable in the asymmetric relationship detected between the time series. This is worth observing since it serves as a guideline on how to estimate each of the selected threshold regressions. In our instance, we specify the mtar models under the assumption that economic growth is regressed on the unemployment rate. This is of course with the exception of the cf filter regression in which we model tar nonlinearity and yet retain economic growth as the dependent variable in the regression. Our estimation results of the first difference model specifications are reported below in table 3 whereas the results obtained for the gap model versions are reported in table 4. Starting with the results reported in table 3 for first differences model, we take note of a long-run coefficient estimate of -0.09. Technically speaking, the magnitude of this coefficient estimate as obtained under both first difference models implies that a 1 percent decrease in the unemployment rate is associated with a -0.09 percent increase in productivity output. This result is seemingly plausible as it does not violate traditional theory of a negative unemployment-growth co-relationship as initially postulated by Okun (1962). Furthermore, the magnitude of this relationship is consistent with some of the Okun coefficient estimates obtained in previous studies. Among these previous studies are the works of Adanu (2005) who obtain a similar estimate of -0.09 for Alberta province in Canada; Villaverde and Maza (2009) who find a -0.08 estimate for a regional group of Spanish data and also Geldenhuys and Marnikov (2007) who obtain an estimate of -0.11 for South African data. In moving on to examining the regime switching behaviour among the co-integration error terms, we firstly note that all threshold estimates table 2 Threshold Cointergation and Error Correction Tests Model (1) (2) tar-tec mtar-tec H(l) Ho H0 „(3) Ho „(1) Ho „(2) Ho „(3) Ho First differences Agdpt Aurt 25-36 (0.00)*** 4.10 (0.05)* 0.47 (0.50) 32.71 (0.00)*** 9.16 (0.01)** 2.47 (0.13) Aurt Agdpt 41.82 (0.00)*** 0.68 (0.42) 0.01 (0.91) 50.82 (0.00)*** 1.66 (0.2l) 0.01 (0.95) hp filter gdPC urt 6.84 (0.0l)** 1.07 (0.31) 0.66 (0.43) 6.15 (0.01)** 0.16 (0.69) 0.25 (0.62) urC gdPC 4.36 (0.02)* 0.22 (0.64) 2.78 (0.11)* 4.46 (0.02)* 1.19 (0.28) 0.49 (0.49) bk filter gdpc urt 28.51 (0.00)*** 3.56 (0.07)* 2.94 (0.11)* 33.43 (0.00)*** 6.70 (0.02)* 1.59 (0.23) urc gdPC 27.28 (0.00)*** 0.01 (0.91) 0.23 (0.64) 32.79 (0.00)*** 0.09 (0.76) 1.10 (0.32) bw filter gdpc urt 26.51 (0.00)*** 4.34 (0.05)* 0.65 (0.43) 34.03 (0.00)*** 9.29 (0.01)** 3.51 (0.08) urt gdPC 54.27 (0.00)*** 1.06 (0.3l) 0.01 (0.94) 55.93 (0.00)*** 0.96 (0.34) 0.66 (0.43) c-tar-tec c-mtar-tec „(1) h(2) "o „(3) „(1) „(2) „(3) First differences Agdpt Aurt 29.08 (0.00)** 6.84 (0.02)* 0.79 (0.39) 32.75 (0.00)*** 9.19 (0.01)** 2.78 (0.1l)' Aurt Agdpt 42.23 (0.00)*** 0.86 (0.36) 0.96 (0.34) 67.86 (0.00)*** 8.18 (0.01)** 1.85 (0.19) hp filter gdPC urt 6.84 (0.01)** 1.06 (0.31) 0.01 (0.98) 10.04 (0.00)*** 5.27 (0.03)* 3.75 (0.07) urt gdPC 5.20 (0.01)* 1.47 (0.24) 3.64 (0.07)* 6.85 (0.01)* 4.81 (0.04)* 5.26 (0.03) bk filter gdPC urt 28.74 (0.00)*** 3.71 (0.07)* 1.08 (0.32) 33.91 (0.00)*** 7.01 (0.01)* 1.82 (0.20) urt gdPC 27.71 (0.00)*** 0.27 (0.61) 0.23 (0.64) 32.79 (0.00)*** 0.09 (0.76) 1.10 (0.32) bw filter gdPC urt 32.08 (0.00)*** 8.35 (0.01)** 1.27 (0.28) 33.28 (0.00)*** 8.77 (0.01)** 2.22 (0.15) urt gdPC 56.83 (0.00)*** 1.99 (0.17) 0.24 (0.63) 60.65 (0.00)*** 2.58 (0.12) 0.44 (0.52) notes Column headings are as follows: (1) dependent variable, (2) independent variable. Significance level codes: ***, ** and * denote the 1%, 5% and 10% significance levels respectively. table 3 Threshold Co-Integration and Error Correction Estimates for First Difference Model Specification MTAR-TEC C-MTAR-TEC y X y X Agdp Aur Agdp Aur ßi -0.09 -0.09 (0.00)*** (0.00)*** P&-1 -O.72 -0.72 (0.0l)** (0.01)** P&-1 -1.76 -1.76 (0.00)*** (0.00)*** T 0 0.11 AAgdp+-k -0.39 -1.18 -0.38 -1.26 (0.47) (0.31) (0.47) (0.27) AAgdp--k -0.30 -0.50 -0.29 -0.47 (0.36) (0.47) (0.36) (0.50) AAm+-k -0.04 -0.80 -0.04 -0.80 (0.64) (0.00)*** (0.66) (0.00) AAw--k -0.09 -0.99 -0.09 -0.99 (0.28) (0.00)*** (0.29) (0.00) 1 0.21 2.39 0.19 2.53 (0.83) (0.27) (0.85) (0.23) 1 -1.82 -1.05 -1.81 -1.06 (0.00)*** (0.14)* (0.00)*** (0.13) R2 0.80 0.86 0.80 0.85 DW 1.61 2.42 1.61 2.39 p-value 0.37 0.31 0.35 0.31 LB 0.31 0.55 0.27 0.59 JB 3.59 3.82 3.65 3.98 notes Significance level codes: ***, ** and * denote the 1%, 5% and 10% significance levels respectively. dw and lb respective denote the Durbin Watson and Ljung-Box test statistics for autocorrelation whereas jb denotes the Jarque-Bera normality test of the residuals. are encouragingly close to zero in value. Moreover, the threshold error term estimates satisfy the convergence condition of error term stationar-ity i. e. p1,p2 < 0 and (1 - p1)(i - p2) < 1. In further diagnosing these co-integration threshold error terms, we observe that negative deviations are eliminated quicker than positive ones. We can make such inference since the estimate of 1 is of a lower absolute value in comparison to its 2 counterpart. Notably, Harris and Silverstone (2001) make similar inferences in their study for both us and uk data. In addition, our estimates of the threshold error correction terms also bear a slight resemblance to those obtained in Harris and Silverstone (2001), in the sense of producing correct negative estimates in the lower regimes of the estimated models. However in differing from these authors, we are able to obtain significant values for the estimates of the threshold error correction terms and thus we can draw meaningful interpretations of the error correction coefficients. In this respect, we not only discover that the long-run error correction terms for both m tar-tec and c-mtar-tec models are almost identical in magnitude, but we more importantly note that the speed of adjustment in both models is quicker when there is a shock to economic growth as opposed to a shock to the unemployment rate. Meanwhile, we are only able to identify significant short-run effects for the lagged coefficients of the economic growth variable when shock has been induced on the unemployment rates, whilst we are find no short-run effects for shocks to economic growth variable. In diverting our attention to the empirical results of the estimated gap versions of Okun's law as reported in table 5, we generally observe that the regression estimates, more or less, bear close resemblance to those attained for the first difference models. For instance, the long-run regression coefficient obtained from the gap version models produce similarly negative estimates, albeit the magnitude of these estimates vary between 0.09 and 0.98 for the different de-trending methods employed. In further considering the absolute coefficient values of the threshold error terms formed by the long-run regressions, we note that the gap model estimates also bear similarities to those obtained for the first difference models. Specifically, we observe that the absolute values ofp1 are significantly higher when the unemployment rate is the driving variable, whilst the values ofp2 are higher when the unemployment rate is the dependent variable in the co-integration system. As previously explained, this result infers that negative shocks are eliminated quicker when economic growth is the driving variable, whereas positive shocks are eliminated quicker when the unemployment rate is the dependent variable. However, after scrutinizing through the threshold error correction model estimates, we find the estimates from the gap models to be less encouraging. This especially becomes apparent when mainly considering the long-run error correction terms, from which we observe that only two table 4 Threshold Co-Integration and Threshold Error Correction Estimates for First De-Trended Model Specification hp filter bk filter bw filter c-MTAR- -TEC c-MTAR- -TEC TAR-TEC MTAR- TEC y X y X Y X Y X Agdp Aur Agdp Aur Agdp Aur Agdp Aur P -0.2 (0.02)** -0.15 (0.01) -0.09 (0.03)* -0.10 (0.01)** P16-1 -0.13 (0.66) -0.88 (0.01)** -0.97 (0.01)** -0.73 (0.01)** P16-1 -0.98 (0.00)*** -0.16 (0.48) -1.68 (0.00)*** -1.77 (0.00)*** T -0.286 -1.747 0 0.254 AAgdp+-k 0.32 (0.27) 0.37 (0.55) 0.44 (0.45) 0.48 (0.11) -0.30 (0.40) -0.56 (0.34) -0.32 (0.55) -1.55 (0.19) AAgdp--k 0.32 (0.61) -0.26 (0.84) -1.09 (0.10)* -1.29 (0.00)*** -1.22 (0.05)* 0.32 (0.74) -0.30 (0.35) -0.30 (0.66) AArnlk 0.08 (0.63) -1.07 (0.00)*** -0.48 (0.23) 0.13 (0.54) 0.06 (0.73) -1.14 (0.00)*** -0.02 (0.82) -0.78 (0.00)*** AAur--k -0.09 (0.42) -0.39 (0.14) -0.48 (0.07)* -0.19 (0.16) 0.12 (0.50) -0.36 (0.22) -0.10 (0.23) -0.99 (0.00)*** ¿+6-1 0.09 (0.83) -1.63 (0.07)* -0.64 (0.03) -0.05 (0.73) -0.09 (0.91) 1.82 (0.19) 0.08 (0.94) 3.03 (.017) -0.88 (0.07)* -0.53 (0.61) 0.12 (0.59) 0.18 (0.12)* -0.44 (0.57) -1.54 (0.24) -1.83 -1.16 (0.00)*** (0.12) R2 0.54 0.57 0.61 0.50 0.46 0.80 0.80 0.85 dw 2.10 1.56 1.43 1.68 2.42 1.85 1.60 2.63 p-value 0.89 0.23 0.09 0.28 0.39 0.62 0.31 0.13 lb 0.54 0.62 0.25 0.18 0.50 0.58 0.23 0.44 jb 3.56 4.10 3.89 4.86 3.79 4.26 3.98 4.58 notes Significance level codes: ***, ** and * denote the 1%, 5% and 10% significance levels respectively. dw and lb respective denote the Durbin Watson and Ljung-Box test statistics for autocorrelation whereas jb denotes the Jarque-Bera normality test of the residuals. models manage to produce negative and significant estimates i. e. the hp and bw filter specifications with economic growth placed as the driving variables in both models. Therefore, we are restricted to interpreting the error correction coefficient estimates solely for these two model specifi- cations. In drawing inference from these estimates, we conclude equilibrium reverting behaviour over the business cycle for the hp filter model when a shock has been induced on either the economic growth or the unemployment variables. Similarly, for the bw filter estimates, long-run equilibrium reversions occur only in the event of a shock to economic growth. It is also interesting to find that for both cases of the first difference models, we obtain significant short-run coefficient estimates of the lagged unemployment variable when a shock has been induced on the unemployment rate. Thus we collectively observe a distinct pattern over the business cycle, in which the unemployment rate is a driving factor of equilibrium conversions over the short-run whilst economic growth is responsible for equilibrium adjustment over the long-run. Having established various forms of threshold co-integration within Okun's law for the data implies that there must exist some form of causality between the variables in the granger sense. However, the direction of causality cannot be assumed a prior and thus should be investigated through a formal analysis. We are permitted to examine causality effects amongst the variables via a standard F-test. The construction of these tests has been adequately discussed in the previous section of the paper. Table 7 reports the results of the causal analysis. The most striking feature of our obtained results is that, in all cases save one, we are able to reject the null hypothesis of unemployment not causing output growth at conventional levels of significance. Conversely, we fail to reject the null hypothesis of economic growth not leading the unemployment rate. We have noted an exceptional case for the hp filter model with economic growth as the driving variable, in which we detect no causal effects within the data. In summing up these results, we can safely assume that our results depict unidirectional causality running from the unemployment rate to economic growth for the data as a whole. This result is plausible seeing that we have already established that economic growth is regressed as being the dependent on the unemployment rate but not vice versa. Concluding Remarks The goal of this paper was to examine nonlinear co-integration and causality effects in Okun's law for South African data dating between the periods of 2000 and 2013. This objective was facilitated through the use of mtar modelling techniques. We favour this approach on the premise of allowing for unit root testing, co-integration analysis and causality analysis under a single, comprehensive framework. Moreover, our study dif- table 5 Granger Causality Tests Model Y X Decision First differ- mtar-tec gdP ur 1.11 34.71 ur ^ gdP ences (0.35) (0.00)*** gdP + ur c-mtar-tec gdP ur 1.19 35.24 ur ^ gdP (0.33) (0.00)*** gdP + ur hp filter c-mtar-tec gdP ur 0.65 0.36 ur ^ gdP cf filter (0.54) (0.70) gdP + ur Bw fiter c-mtar-tec ur gdP 3.97 1.50 ur ^ gdP (0.04)** (0.25) gdP + ur c-tar-tec gdP ur 0.49 12.61 ur ^ gdP (0.62) (0.00)*** gdP + ur c-mtar-tec gdP ur 1.27 32.37 ur ^ gdP (0.31) (0.00)*** gdP + ur notes Significance level codes: ***, ** and * denote the 1%, 5% and 10% significance levels respectively. Definitions of notations: « and + represent unidirectional causality, bi-directional causality and no causality, respectively. fers from previous South African case studies as we are able to introduce nonlinearity in a strict co-integration sense. Having applied the mtar framework to South African unemployment and economic growth data has produced a number of interesting policy considerations. First of all, in quantifying the long run correlation coefficient, we find negative Okun coefficients ranging from -0.09 to -0.20 for all estimated threshold models. Clearly, these observations have far reaching ramifications as they give rise to the intriguing possibility of a long run trade-off between unemployment and economic growth. However, the aforementioned observations are of limited policy value in absence of knowing the causal relations amongst the variables. In examining the empirical results obtained from the causal analysis, we discover that during abrupt shocks to the economy there are no causal effects between the variables. This essentially means that in the event of sharp or anticipated shocks to the economy there is very little that policy intervention can do for long-run equilibrium restoration between unemployment and economic growth. However, during smooth shocks, unemployment granger causes economic growth thus allowing for direct labour policies to have an impact on output productivity. We substantiate these smooth shocks as carefully implemented and monitored policies directives which are aimed at narrowing the existing gap between the demand and supply within South African labour markets. Inclusive of such shocks are policy programmes aimed at improving the higher education system through intensifying further education and training (f et) programmes and the recently proposed 'target wage subsidy' programme which is intended to facilitate for the school-to-work transition within the youth population. We also note that under no circumstance does economic growth granger cause unemployment thus insinuating that policies aimed directly at improving economic growth such as foreign exchange policies would exert little or no influence on eradicating unemployment over the long run. This is particularly worth noting since it has been previously assumed that the stability of the exchange rate would lead to a direct improvement of employment growth in import-competitive and export-oriented sectors, especially the manufacturing sectors. Our study implies that, whilst these macro-policies may create a sustainable environment for improved economic growth, they are of little use with regards to directly eradicating unemployment. Therefore, the overall finding of uni-directional causality from unemployment to economic growth provides an adequate explanation for the 'job-less' growth pandemic experienced in South Africa over the last two decades or so. In recent South African recession periods, unemployment has continued to rise despite economic growth seemingly returning to its previous long-run trend. Deriving from our study, there exist two rational explanations to this pandemic. Firstly, negative shocks to economic growth are eradicated quicker than negative shocks to unemployment. This implies that in the event of smooth shocks to output productivity, it should be expected that economic growth should return back to its long-run steady state at a quicker rate than its unemployment counterpart. Secondly, our general finding of causality running from unemployment to economic growth highlights the ineffectiveness of macroeconomic policies aimed at reducing unemployment through improved productivity growth. Specifically, our empirical estimates suggest that smooth unemployment shocks, in the form of structural labour policies, would help stabilize the structural and cyclical components of unemployment over both the short and the long run. Therefore, for the specific case of South Africa, the job-less growth pandemic can be attributed to structural factors underlying the overwhelming unemployment rates facing the economy. Specifically, structural unemployment in South Africa is a result of a mismatch between jobs offered by employers and potential workers. Thus factors which can minimize the extent of structural unemployment within the economy need to be addressed in policies. These factors include the rigidity of labour markets, real wage rigidity, high minimum wages compared to relatively low productivity and other factors which would lead to job creation and job security. Overall, we conclude that labour policies aimed at stabilizing and eradicating structural unemployment within the economy may be a panacea towards simultaneously reducing overall unemployment and boosting economic growth over the long run. References Adanu, K. 2005. 'A Cross-Province Comparison of Okun's Coefficient for Canada.' Applied Economics 37 (5), 561-70. Ahmed, K., S. Khali, and A. Saeed. 2011. 'Does There Exist Okun's Law in Pakistan?' International Journal of Humanities and Social Science 1 (12): 293-99. Baxter, M., and R. King. 1999. 'Approximate Band-Pass Filter for Economic Time Series.' The Review of Economics and Statistics 81 (4): 575-93. Beccarini, A., and D. Gros. 2008. At What Cost Price Stability? New Evidence about the Phillips Curve in Europe and the United States.' ceps Working Document 302, Centre for European Policy Studies, Brussels. Blake, N. S., and T. B. Fomby. 1997. 'Threshold Co-integration.' International Economic Review 38 (3): 627-45. Burger, P., andM. Marinkov. 2006. 'TheSouthAfrican PhillipsCurve: How applicable is the Gordon model?' South African Journal of Economics 74 (2): 172-89. Caner, M., and B. Hansen. 2001, 'Threhsold Autoregression with a Unit Root.' Econometrica 69 (6): 1555-96. Caraiani, P. 2010. 'Bayesian Linear Estimation of Okun Coefficient for Romania: Sensitivity to Prior Distributions.' The Romanian Economic Journal 13 (38): 53-65. Dixon, R., and D. Shepard. 2002. 'Economic Growth and Unemployment in Europe.' European Business Journal 12 (4): 185-98. Enders, W., and S. Dibooglu. 2001. 'Long-Run Purchasing Power Parity with Asymmetric Adjustment.' Southern Economic Journal 68 (2): 43345. Enders, W., and C. Granger. 1998. 'Unit-Root Tests and Asymmetric Adjustment with an Example Using the Term Structure of Interest Rates.' Journal ofBusiness and Economic Statistics 16 (3): 304-11. Enders, W., and P. Silkos. 2001. 'Co-Integration and Threshold Adjustment.' Journal ofBusiness and Economic Studies 19 (2): 166-76. Engle, R., and C. Granger. 1987. 'Co-Integration and Error Correction Representation, Estimation and Testing.' Econometrica 55 (2): 251-76. Freeman, D. 2000. 'A Regional Test of Okun's Law.' International Advances in Economic Research 6 (3): 557-70. Geldenhuys, J., and M. Marnikov. 2007. 'Cyclical Unemployment and Cyclical Output: An Estimation of Okun's Coefficient for South Africa.' South African Journal of Economics 75 (3): 373-90. Gupta, R., and J. Uliwingiye. 2010, 'Dynamic Time Inconsistency and the South African Reserve Bank.' South African Journal of Economics 78(1): 76-88. Hansen, B. 1997. 'Inference in tar Models.' Studies in Nonlinear Dynamics and Econometrics 2 (1): 1-14. Harris, R., and B. Silverstone. 2011, 'Testing for Asymmetry in Okun's Law: A Cross Country Comparison.' Economics Bulletin 5 (2): 1-13. Hodge, D. 2006. 'Inflation and growth in South Africa.' Cambridge Journal ofEconomics 30(2): 163-80. -. 2009. 'Inflation, Unemployment and Economic Growth in South Africa.' South African Journal of Economics 77 (4): 488-504. Hodrick, R., and E. Prescott. 1997. 'Postwar us Business Cycle: An Empirical Investigation.' Journal of Money, Banking and Credit 29 (1): 1-16. Jardin, M., and S. Gaetan. 2011. 'How Okun's Law is Non-Linear in Europe: A Semi-Parametric Approach.' Working paper, University of Rennes, Rennes. Knotek, E. 2007. 'How Useful is Okun's Law?' Economic Review Q473-103. Lal, I., S. Muhammad, A. Jalil, and A. Hussian. 2010. 'Test of Okun's Law in Some Asian Countries: Co-Integration Approach.' European Journal of Scientific Research 40 (1): 73-80. Moosa, 1.1997. 'A Cross Country Comparison of Okun's Coefficient.' Journal of Comparative Economics 24 (3): 335-56. Noor, Z., N. Nor, and J. Ghani. 2007. 'The Relationship between Output and Unemployment in Malaysia: Does Okun's Law Exist.' International Journal ofEconomics and Management 1 (3): 337-44. Okun, A. 1962. 'Potential gnp: Its Measurement and its Significance.' Cowles Foundation Paper 190. http://cowles.econ.yale.edu/P/cp/p01b/ p0190.pdf Pereira, E., C. Bento, and R. Silva. 2009. 'Co-Integration and Asymmetric Adjustment between Output and Unemployment.' Management and Industrial Engineering Working Paper 52/2009, University of Aveiro, Department of Economics. Pollock, D. 2000. 'Trend Estimation and De-Trending via Rational Square-Waves Filters.' Journal of Econometrics 99 (2): 317-34. Silvapulle, P., I. Moosa, and M. Silvapulle. 2004. 'Asymmetry in Okun's Law.' The Canadian Journal of Economics 37 (2): 353-74. Turturean, I. 2007. 'Okun's law in Romania during 1992-2004.' mpra Working Paper 10638,, University Library of Munich, Munich. Villaverde, J., and A. Maza. 2009. 'The Robustness of Okun's Law in Spain, 1980-2004: Regional Evidence.' Journal of Policy Modeling 31 (2): 28997. Weber, C., and J. West. 1996. 'Functional form in Regression Models of Okun's Law.' Applied Economics Letters 3 (9): 607-9. This paper is published under the terms of the Attribution-NonCommercial-NoDerivatives 4.0 International (cc by-nc-nd 4.0) License (http://creativecommons.org/licenses/by-nc-nd/4-0/).