D I S A G G R E G A T E D Y N A M I C S A N D E C O N O M I C G R O W T H IN CANADA by ELIZABETH CLARE WAKERLY B.A., University of Sheffield, 1987 M.A., Carleton University, 1988 A THESIS SUBMITTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in THE FACULTY OF GRADUATE STUDIES Department of Economics We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA January 1997 ©Elizabeth Clare Wakerly, 1997 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of the thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of the thesis for financial gain shall not be allowed without my written permission. Department of Economics The University of British Columbia # 997-1873 East M a l l Vancouver, B . C . , Canada V 6 T 1Y2 Date: Abstract This thesis takes the form of three essays in which I use disaggregate and aggregate information to examine Canadian economic growth. In the first essay, I present evidence that the process of economic growth differs for low income per capita provinces and industries. This contrasts with results from traditional studies of economic convergence. In those papers, estimates of a rate of convergence suggest that poor provinces eventually "catch up" to rich provinces by growing faster. Unfortunately, this approach ignores the pattern of economic growth within the cross-section distribution. Explicitly modelling the evolving distribution, I find little mobility in the cross-sectional ordering and some evidence of divergence. In the long run, the poor stay (relatively) poor and the rich remain (relatively) rich. In the second essay, I examine the dynamic effects of aggregate and disaggregate disturbances on both economic growth and the interaction between disaggregates. The approach is motivated by the class of models which predict two-way interaction between aggregate and disaggregate behaviour, such as Durlauf [28]. The disaggregate disturbance is identified as having no long-run impact on aggregate economic growth. I find that the aggregate shock has a large impact on aggregate income in both the short and long run; and accounts for most of its variation. The disaggregate shock contains some information for aggregate activity at business cycle horizons. Most interaction is explained by the disaggregate disturbance; the aggregate shock contributes little. In the third essay, I present results from a variety of unit root tests on provincial and manufacturing industry panel income data. Standard DickeyFuller unit root tests applied to panels require averaging of data across the cross-section. More powerful tests allow pooling of cross-section and time-series information. Using these methods, I find that the null hypothesis of a unit root is rejected—strongly contrasting with results obtained using the standard Dickey-Fuller methodology. 11 Contents Abstract ii Contents m List of Tables v List of Figures vii Acknowledgement ix Chapter 1 Introduction 1 Chapter 2 The Evolving Distribution of Canadian Economic Growth . 6 2.1 Introduction 6 2.2 Related Literature 8 2.3 Econometric Model 14 2.4 Cross-Section Dynamics 15 2.4.1 Provincial 16 Personal income 16 Wages, salaries and supplementary labour income . . . . 19 2.4.2 Industrial 21 2.5 Robustness 24 2.5.1 Provincial data adjustment 25 Relative to Ontario (1) 25 Relative to Ontario (2) 26 Aggregating some provinces 26 2.5.2 Industry data 27 Ontario 27 Western provinces 28 Eastern provinces 28 2.6 Conclusions 29 Chapter 3 The Dynamic Effects of Aggregate and Disaggregate Disturbances 37 3.1 Introduction 37 3.2 Identification 39 3.3 Interpretation 40 3.4 Results 41 3.4.1 Provincial data 41 3.4.2 Industrial data 45 3.5 Conclusions 46 Chapter 4 Testing For Canadian Unit Roots: A Panel Data Approach 57 4.1 Introduction 57 iii 4.2 Unit Roots and Panel Data 4.3 Methods 4.3.1 The panel Dickey-Fuller test 4.3.2 The <5-bar test 4.3.3 The i-bar test 4.4 Results 4.4.1 Provincial data Panel Dickey-Fuller test 6-hax test: base case £-bar test: index adjustment t-bar test 4.4.2 Industrial data Panel Dickey-Fuller test 6-bar test: base case S-bax test: index adjustment t-bar test 4.4.3 Dickey-Fuller test 4.5 Conclusions Chapter 5 Conclusions Bibliography Appendices A The Data A . l Data Sources A . 2 Manufacturing Industries B Fractile Transition Probability Matrices C Further Estimations C . l Provincial Economic Growth per Capita C.2 Provincial W S S L I per Capita C.3 Provincial Employment per Capita C.4 Industrial G D P per Employee D Tests For Robustness E Results Using Coulombe and Lee Data E . l Data E.2 Mobility Dyanmics E.2.1 Output and income per capita E.2.2 Productivity measures E.2.3 Longer sample period E . 3 Conclusions F Quah and Sargent's Dynamic Index Model F . l Theoretical Model F.2 Application iv 58 62 62 63 67 69 69 69 70 71 71 73 73 73 73 74 76 77 79 83 89 90 90 93 94 95 95 96 97 98 107 117 117 118 118 119 120 120 122 122 124 List of Tables 2.1 Transition Probability Matrix—Personal Income per Capita by Province 16 2.2 Fractile Transition Probability Matrix—Personal Income per Capita by Province 18 2.3 Transition Probability Matrix—Wages, Salaries and Supplementary Labour Income per Employee, by Province 20 2.4 Transition Probability M a t r i x — G D P per Employee, 21 Manufacturing Industries 23 3.1 Variance Decomposition of Aggregate Income per Capita and Provincial Interaction 44 3.2 Variance Decomposition of Aggregate G D P per Employee and Manufacturing Industry Interaction 47 4.1 £-Bar Test Statistics for Provincial Data 70 4.2 £-Bar Test Statistics (with Index Adjustment) for Provincial Data 71 4.3 t-Statistics (t-Bar Test) for Provincial Data 72 4.4 6-Bar Test Statistics for Manufacturing Industry Data 74 4.5 6-bar Test Statistics (with Index Adjustment) for Manufacturing Industry Data 74 4.6 t-Statistics (tf-Bar Test) for Manufacturing Industry Data . . . . 75 4.7 D F Test Statistics for Provincial and Manufacturing Industry Data 76 B l W S S L I Per Employee B y Province 94 B2 G D P Per Employee, 21 Manufacturing Industries 94 C I Transition Probability Matrix—Income Growth by Province . . . . 95 C2 Transition Probability M a t r i x — W S S L I Per Capita by Province . . 97 C3 Transition Probability Matrix—Employment Per Capita by Province 98 C4 Transition Probability M a t r i x — G D P Per Person at Work by Industry Sector 99 C5 Transition Probability M a t r i x — G D P Per Person Hour at Work by Industry Sector 99 D I Transition Probability Matrix—Income by Province, Relative to Ontario 107 D2 Transition Probability Matrix—Income by Province, Subtracting Ontario 108 D3 Transition Probability Matrix—Income by Province, Aggregating some Provinces 108 D4 ing D5 ing D6 ing El E2 E3 Transition Industries, Transition Industries, Transition Industries, Data Used Data Used Data Used Probability M a t r i x — G D P Per Employee, ManufacturOntario Probability M a t r i x — G D P Per Employee, ManufacturWestern Provinces Probability M a t r i x — G D P Per Employee, ManufacturEastern Provinces in Coulombe and Lee [21] in Coulombe and Lee [47] in Coulombe and Lee [22] vi 109 109 110 117 118 118 List of Figures Figure 1 Distribution Dynamics 31 Figure 2 Log of Provincial Personal Income per Capita (1926-1994) . . 32 Figure 3 Quantiles for Log of Provincial Personal Income per Capita (1928-1992) 33 Figure 4 Log of Provincial W S S L I per Employee (1966-1994) 34 Figure 5 Quantiles for Log of Provincial W S S L I per Employee (19681993) 35 Figure 6 Log of Manufacturing Industry G D P per Employee (1961:11994:11) 36 Figure 7 Impulse Response: Aggregate Disturbance on Aggregate Income Per Capita 49 Figure 8 Impulse Response: Disaggregate Disturbance on Aggregate Income per Capita 50 Figure 9 Impulse Response: Aggregate Disturbance on Interaction . . 51 Figure 10 Impulse Response: Disaggregate Disturbance on Interaction 52 Figure 11 Impulse Response: Aggregate Disturbance on Aggregate G D P per Employee 53 Figure 12 Impulse Response: Disaggregate Disturbance on Aggregate G D P per Employee 54 Figure 13 Impulse Response: Aggregate Disturbance on Interaction Measure 55 Figure 14 Impulse Response: Disaggregate Disturbance on Interaction Measure 56 Figure C l Growth of Provincial Personal Income per Capita (1927-1993) 100 Figure C2 Quantiles for Growth of Provincial Personal Income per Capita (1929-1992) 101 Figure C 3 Log of Provincial W S S L I per Capita (1961-1992) 102 Figure C 4 Quantiles for Log of Provincial W S S L I per Capita (1963-1991) 103 Figure C5 Provincial Employment per Capita (1966-1993) 104 Figure C6 Index of G D P per Employee 1986=100 (1961-1994) 105 Figure C7 Index of G D P per Person Hour Worked 1986=100 (1961-1994) 106 Figure D I Quantiles for Log of Provincial Personal Income per Capita, Divided by Ontario Data (1928-1993) Ill Figure D2 Quantiles for Log of Provincial Personal Income per Capita, Subtracting Ontario Data (1928-1993) 112 vii Figure D3 Quantiles for Log of Aggregated Provincial Personal Income per Capita (1928-1991) Figure D4 Quantiles for Log of Manufacturing Industry G D P per E m ployee: Ontario (1973-1990) Figure D5 Quantiles for Log of Manufacturing Industry G D P per E m ployee: Western Provinces (1973-1990) Figure D6 Quantiles for Log of Manufacturing Industry G D P per E m ployee: Eastern Provinces (1973-1990) Figure F l Sample Standard Deviations: Including Employment (19681988) Figure F2 Sample Standard Deviations: Including G N P (1968-1988) . Figure F3 Sample Standard Deviations: Including Employment vs including G N P (1968-1988) Figure F4 Sample Standard Deviations: Two-Index Model vs Total Employment (1968-1988) Figure F5 Sample Standard Deviations: Two-Index Model vs Total G N P (1968-1988) viii 113 114 115 116 125 126 127 128 129 Acknowledgement I would like to thank my thesis committee, John Cragg, John Helliwell and James Nason for helpful comments and guidance. I am indebted to my fellow graduate students and my parents for their support. In particular, I would like to thank Shaun Vahey, without whose support, encouragement and patience this research would not have been possible. ix Chapter 1 Introduction This thesis takes the form of three essays in which I use aggregate and disaggregate information to exarriine Canadian economic growth. The importance of disaggregate information in explaining aggregate activity is emphasised i n the papers by Durlauf [28], Galor and Zeira [35], Long and Plosser [53] and Quah [68]. I find that, in general, averaging out disaggregate effects causes a loss of information and generates misleading results. In the first essay, Chapter 2, I examine the contribution of disaggregates to economic convergence. Traditionally, empirical work focuses on estimating a single rate of convergence across disaggregates. Using Canadian provincial data and techniques Barro [3] and Barro and Sala-i-Martin [4] develop, various authors, that include Coulombe and Lee [21] and [47], Helliwell [40] and Sala-iM a r t i n [77], regress growth rates on initial income per capita levels. They find convergence of around 1-3% per annum for selected post-WWI samples. A n interpretation of this result takes the form of evidence that, i n the long run, poor regions "catch up" rich regions (or the differential between the rich and poor decreases) by growing faster. Unfortunately, the approach used in these studies reveals little about the evolving pattern of economic growth within the cross-section because averaging annihalates the disaggregate information. A convergence rate of say 2% per annum across Canadian provinces, sheds no light on whether rich and poor economies possess different short-run growth dynamics. Nor does it provide any information about the probability of a relatively poor province improving its rank i n the cross-section. (And arguably, both policy makers and individuals are more concerned with rankings than with absolute economic growth.) A n average convergence rate cannot tell us whether it is the already rich which are converging towards each other or whether the poor are becoming relatively poorer. 1 Existing studies also neglect the relationship between provincial and industrial patterns of growth. If, for example, the pattern of economic growth for the major manufacturing industries i n Newfoundland and Alberta—food and chemicals respectively—is very different, then the finding that income per capita i n the two provinces is converging should be treated with some caution. Since the growth patterns of the leading industries differ, the path of provincial income per capita most likely differs as well. As an alternative to the traditional approach to measuring convergence, I explicitly study the relationships among Canadian disaggregates. To do this, I model the cross-section dynamics of Canadian income at a disaggregate level. The approach I follow builds on work by Quah [69]. I use disaggregate provincial and industrial data, and examine the sensitivity of the results to different measures of economic growth, to the inclusion of Ontario data, and to asymmetries in province size. I show that although the long-run cross-section distribution of provincial income per capita is generally unimodal, there exists virtually no mobility i n the cross-sectional ordering of income levels. I also report evidence to support the view that the process of economic growth is distinct for poorer regions. Put differently, poor provinces remain relatively poor and rich provinces remain relatively rich. Adjusting the data to account for the size differences of the provinces, and for the disproportionate role played by Ontario, does not affect the main results. It seems that it is those factors which determine a region's comparative advantage (geographic location, climate, endowment of natural resources, etc.) which determine the steady-state distribution of income across provinces. Disaggregation by manufacturing industry across Canada, using monthly data on gross domestic product ( G D P ) per employee, reveals a bi-modal or divergent pattern. However, the manufacturing sectors in the top part of the distribution are evenly spread across provinces; and, i n general, those i n the bottom contribute only a small proportion of provincial G D P . The bi-modal pattern is also evident in Ontario manufacturing industry data, but not for aggregations of western and eastern provinces. This suggests that the pattern in Ontario is driving the national picture of divergence. The different results for provincial income per capita and manufacturing 2 industry G D P per employee data are striking. The income per capita data— which include transfer payments—show a unimodal long-run distribution; the G D P per employee data—which exclude transfer payments—exhibit divergence or a bi-modal distribution. This highlights the importance of transfer payments: in their absence, the long-run tendency of Canadian economic growth is towards divergence. B y ignoring the disaggregate information, existing studies that average across disaggregates present an incomplete picture of convergence. These studies fail to capture the shape and mobility dynamics within the cross-section distribution. In the second essay, Chapter 3, I examine the impact of two types of disturbance on Canadian economic growth. Building on the finding i n the previous Chapter that disaggregate dynamics contain important information for explaining the pattern of convergence, I directly identify the dynamic effects of aggregate and disaggregate disturbances. I use the long-run, just-identified structural vector autoregression ( S V A R ) method of Blanchard and Quah [14]. I estimate an unrestricted vector autoregression ( V A R ) system that includes an aggregate measure of economic growth and a measure of disaggregate interaction. To transform the unrestricted V A R to the S V A R , I assume that there exists an aggregate shock and a disaggregate or regional shock. The aggregate disturbance has permanent effects on aggregate income. This could be, for example, a technology shock. The disaggregate disturbance is identified as having no long-run impact on the level of aggregate income. For example, a provincial government fiscal stimulus might have short-run aggregate and disaggregate effects, but in the long-run leave aggregate economic activity unchanged. The two types of disturbance are assumed to be orthogonal. Interaction, or mobility, is assumed to be unaffected by either disturbance in the long run. If disturbances had long-term effects, it would imply that the level of interaction could drift over time. However, as shown in Chapter 2, this is not the case. The short-run dynamics of the two shocks on interaction are unrestricted. I examine annual provincial and monthly manufacturing industry data. I find that the aggregate shock has a large positive impact on provincial income per capita in both the short and the long run; and accounts for most of its variation. The disaggregate disturbance initially increases income per capita; but the 3 economy quickly adjusts. The effect of this disturbance becomes insignificant within 5 years. The short-run impact of the two types of disturbance on the provincial interaction measure is imprecisely estimated. The disaggregate shock provides the strongest impulse; and explains most of the variation. But its effect decays to zero within 3 years. The aggregate shock has a larger medium-term impact, which disappears within 7 years. A similar pattern emerges with the industry data. The aggregate disturbance has a positive long-term effect on the aggregate measure, and explains most of its variation. The effect of the disaggregate disturbance is close to zero. Neither disturbance has much of an effect on the interaction measure. W i t h i n only a few months, the impacts are indistinguishable from zero. The results i n this Chapter confirm that disaggregate dynamics contain i m portant information for explaining the pattern of economic growth, at least at business cycle horizons. B y ignoring disaggregate disturbances, researchers ascribe too much importance to aggregate disturbances in explaining this pattern. The role of the aggregate disturbance in explaining most of aggregate economic activity is consistent with the neoclassical growth model, where only technological shocks are assumed to bring about long-run changes i n output. In the third essay, Chapter 4, I present the results from a variety of unit root tests on Canadian income panel data. I compare the results of the standard Dickey-Fuller method with those from three tests designed specifically for panels. Applying standard unit root tests to panel data loses information, since data must be aggregated across cross-sections. More powerful tests, which allow cross-section, as well as time-series information to be included, have recently been developed. Unlike the Dickey-Fuller test, these have limiting normal distributions. The first panel data test is similar to the standard approach, but the initial observation for each disaggregate time-series is subtracted from each subsequent observation i n that series to allow for fixed effects. The unit root hypothesis is then tested by applying a conventional t-test to the coefficient on the lagged dependent variable, using data which is averaged across the cross-section. I refer to this as the "panel Dickey-Fuller" test; see Breitung and Meyer [16]. In the second test, the data are adjusted by subtracting the cross-section mean from all observations. The test statistic is an average of individual co- 4 efficients on the lagged dependent variable. I refer to this as the "£-bar" test (Levin and L i n [51]). I extend this methodology to allow for potentially more complex common effects across individual series. Using a dynamic index model developed by Quah and Sargent [76], I identify multiple common factors and subtract the resulting indices (in place of the cross-section average) from the data. In the third approach, the "f-bar" test (Im, Pesaran and Shin [42]), the initial observation is subtracted—in a similar manner to the panel Dickey-Fuller case. The test statistic is a (small sample adjusted) average of the Dickey-Fuller tests for the individual series. I examine Canadian annual data on gross provincial product and monthly data on gross domestic product per employee for 19 manufacturing industries. Using the standard Dickey-Fuller test, applied to data averaged across the crosssection, I cannot reject the null hypothesis of a unit root i n either the provincial or industrial data. In contrast, none of the three panel unit root tests suggest that the data are characterised by a unit root. This suggests that discarding cross-section information can mislead the researcher on the dynamic properties of the model. This result is i n the spirit of Granger's [36] hypothesis, i n finding that aggregating a dynamic multivariate process can lead to a univariate process with fundamentally different properties. In the final chapter, Chapter 5, I draw some conclusions. 5 Chapter 2 The Evolving Distribution of Canadian Economic Growth 2.1 Introduction Canadian economic growth has long concerned both economists and politicians. Recent empirical work has focused on estimating the "rate of convergence" across provinces. Coulombe and Lee [21] and [47], Helliwell [40] and Salai-Martin [77] regress growth rates on initial income per capita levels. This technique was developed by Barro [3] and Barro and Sala-i-Martin [4] for looking at convergence across U.S. states. Most studies find convergence of around 1-3% per annum for various post-WWI samples. A n interpretation of this result takes the form of evidence that, i n the very long run, poor economies "catch up" rich economies by growing faster. Unfortunately, this approach reveals little about the evolving pattern of economic growth within the cross-section. A convergence rate of, say, 2% per annum across Canadian provinces neither sheds light on whether rich and poor economies possess different short-run growth dynamics nor does it provide any information about the probability of a poor province improving its rank in the cross-section. A n average convergence rate cannot tell us whether it is the already rich which are converging towards each other or whether the poor get (relatively) poorer. A n d arguably, policy makers and individuals are more likely to be concerned with the position of a province relative to others, rather than with absolute changes over time. Another issue existing studies neglect is the relationship between provincial and industrial patterns of growth. Clearly, the two are closely related. Suppose that the pattern of economic growth for the major manufacturing industries in Newfoundland and Alberta—food and refined petroleum products respectively—is very different. Then, the finding that income per capita i n the two provinces is converging should be treated with some caution. Since the 6 growth patterns of the leading industries differ, the path of income per capita growth is also likely to differ. In this Chapter, I study the relationships among Canadian provincial and industrial disaggregates, explicitly modelling the dynamics. I examine both the shape of the cross-section distribution over time and the degree of mobility within it. The approach I follow builds on work by Quah [69]. I use disaggregate data on provincial income per capita and G D P per employee for manufacturing industries, to give a more complete picture of the mechanisms driving economic growth. I also examine the sensitivity of the results to different measures of economic growth (in Appendix C ) , to the inclusion of Ontario data, and to asymmetries i n province size. I show that although the cross-section distribution of provincial income per capita is generally unimodal, there is virtually no mobility i n the cross-sectional ordering of income levels. Put differently, poor provinces remain poor and rich provinces remain rich. Newfoundland remains at the bottom of the distribution and Ontario maintains its place at the top. Adjusting the data to account for the size differences of the provinces, and for the disproportionate role played by Ontario does not affect the main results. These findings suggest that the initial distribution of resource endowments is an important determinant of the steady-state distribution of income across Canadian provinces. In contrast to the unimodal provincial distribution, disaggregation by manufacturing industry across Canada, using monthly data on G D P per employee, reveals a bi-modal or divergent pattern. However, the manufacturing sectors i n the top part of the distribution are evenly spread across provinces; and those in the bottom contribute a small proportion of provincial G D P . The bi-modal pattern is also evident i n Ontario manufacturing industry data, but not for aggregations of western and eastern provinces. This suggests that the pattern i n Ontario could be driving the national picture of divergence. The different results for the provincial and manufacturing industry data are striking. The provincial income per capita data—which include transfer payments—show a unimodal long run distribution; the manufacturing industry G D P per employee data—which exclude transfer payments—exhibit divergence or a bi-modal distribution. This suggests that transfer payments play an important role: in their absence, the long-run tendency of Canadian economic growth 7 is likely to be towards divergence. For policy makers concerned with provincial inequalities, this study provides a mixture of good and bad news. On the one hand, there is evidence of provincial unimodality—consistent with long-run convergence. O n the other hand, there is no evidence to suggest that poorer provinces can overtake richer regions in the short or medium run. •' The remainder of this Chapter is organized as follows. In Section 2.2,1 review the related literature on economic convergence. In Section 2.3, I present the econometric model. I examine cross-section mobility dynamics i n Section 2.4, using provincial measures of personal income per capita, and wages and salaries per employee; and manufacturing industry data on G D P per employee. In Section 2.5, I investigate the robustness of the results. I draw some conclusions in Section 2.6. 2.2 Related Literature Following Barro [3], Barro and Sala-i-Martin [4] and [5] and Mankiw et al [55], empirical studies of economic growth and convergence have tended to concentrate on estimating rates of convergence across economies. The neoclassical growth model (Cass [19], Koopmans [44] and Solow [80]) predicts that income per capita converges to a steady-state, where per capita output, capital stock and consumption grow at a given rate of technological progress. Diminishing returns to capital imply that additions to the capital stock generate more output when the initial stock is relatively small. Researchers begin by assuming a Cobb-Douglas production function, (2.1) Y(t) = K(ty(A(t)L(t))i-° where 0 < a < 1, Y is output, K is capital, L is labour and A is labouraugmenting technological progress. L and A are assumed to grow at rates n and g respectively, with laws of motion given by, L(t) = L(0)e nt A(t) = A(0)e<*. Denning s to be the rate of saving, the dynamic equation for the stock of capital per unit of effective labour, fc, is given by, (2.2) fc(<) = sy(t) -(n + g + 6)ic(t) = sk(i) a 8 -(n + g + S)k(t), where y is output per unit of effective labour and 6 is the rate of depreciation. k converges to its steady state value given by, Substitution gives the following expression for steady state per capita income, (2.3) log{y) = logA(0) + gt + -^—log(s) 1—a where y(t) = Y(t)/L(t) - -^—log(n 1 —a is income per capita. + g + 6) Given differing resource en- dowments, climate and institutions across countries, researchers assume that logA(0) = a + e where a is a constant and e is a country-specific shock term. Approximating around the steady state level of income per effective worker, y*, gives, (2.4) y(t) = /3[log(y*) - log(y(t))]. where f3 = (n + g + 6)(1 — a). Solving this equation and rearranging gives the average growth rate of per capita income y between dates 0 and T , - (2 5) 1, \y(T)] l-e-P , T r y* " T [m= -^ [w). io9 9+ lo9 where f3 is the rate of convergence. If f3 > 0, poor economies grow faster than rich economies—as predicted by neoclassical growth models. This is known as ^-convergence. Denning cr to be the variance of y, the data show cr-convergence 2 if ax < o~t-i for all t. That is, there is a decline i n dispersion, over time, i n the (entire) cross-section. Researchers estimate equations similar to equation (2.5). The coefficient on log declines i n magnitude with the length of the sample, for a given /?. As T increases, the effect of the initial position on the average growth rate gets smaller; as T —* oo the coefficient tends to zero. f3 is estimated non-linearly to take account of the associated value of T in equation (2.5). In this way, similar estimates of /? should be obtained, regardless of the length of the sample. Estimated on a wide range of data—across and within countries—values of f3 derived from this model suggest a uniform convergence rate of 2% per annum. For example, Barro and Sala-i-Martin [5] estimate /? at 0.0175 (standard error 0.0046) for per capita personal income i n the US states over the period 18801988. In an earlier paper ([4]), the authors show that the dispersion of US state 9 personal income has declined from about 0.50 i n 1880 to around 0.20 in the 1980s. In [4], Barro and Sala-i-Martin also examine convergence across U . S . nonagricultural sectors over the period 1963-1986. 1 Estimates of /? range from 0.0093 in Wholesale and Retail Trade, to 0.0460 in Manufacturing (standard errors are 0.0064 and 0.0082, respectively). Barro and Sala-i-Martin [5] obtain similar results from cross-country studies. For a sample of 98 countries, they find evidence of "conditional" convergence for the period 1960 to 1985. The neoclassical growth model implies only that economies with the same preferences and technology will reach a common steady state. When countries are diverse, it is necessary to allow for heterogeneity i n certain institutional variables—for example, initial school enrolment rates, the ratio of government consumption expenditure to G D P , and proxies for political stability—in order to apply the framework in cross-country studies. This is the sense in which the convergence is "conditional". Several studies recognise that the transfer of technology is an important part of the story. Using annual data (1963-1989) on technical progress i n O E C D countries, Helliwell [39] shows that there has been significant international convergence in the rate of technical progress: initially poorer countries have faster technical progress. In Helliwell and Chung [41], analysis of the Solow residuals for O E C D countries for 1960-1985 shows convergence in rates of technical progress. The authors conclude that convergence i n per capita G D P is not just a function of differences in investment rates. Recent Canadian studies include Coulombe and Lee [21] and [47], Helliwell [40], Lee [46] and Sala-i-Martin [77]. Coulombe and Lee ([21] and [47]) examine Canadian provincial income and output per capita data during 19611991. They split the sample into three sub-periods, 1961-1971, 1971-1981 and 1981-1991, and regress each province's growth rate, relative to the Canadian average, on the initial level of income, relative to the Canadian average, for each time period. This approach helps to mitigate the effect of sample-specific time trends. Coulombe and Lee report estimates of /3 that range from 0.0105 for gross provincial product per capita to 0.0289 for personal disposable income The eight sectors are: mining; construction; manufacturing; transportation; wholesale and retail trade; finance, insurance and real estate; services; and government. l 10 per capita. Lee [46] uses an extended version of the neoclassical model to account for differences i n various investments and industrial structures across Canadian provinces. Annual growth i n output per worker is regressed on the initial level of output per worker, the sum of technical progress, labour growth and the capital depreciation rate, and the different investment rates (human, private physical and public physical). Over the period 1966-1992, Lee [46] finds that the dispersion of productivity declined steadily. Helliwell [40] shows that provinces which were poorer in 1961 have faster average growth rates over the period 1961-1989 and that there is a significant downward trend in the interprovincial variation of income levels during 19261990. Sala-i-Martin [77] estimates f3 for Canadian income per capita between 1961 and 1991 at 0.024 (standard error of 0.008). Traditional methods of measuring convergence are also related to growth patterns among trading partners (Ben-David [8]) and in explaining educational and financial development (Berthelemy and Varoudakis [11]). Ben-David [8] finds that grouping countries according to their primary trade affiliations tends to produce significant income convergence within the groups. Such significant convergence is uncommon among these countries when they are grouped randomly. Berthelemy and Varoudakis [11] find that two groups of countries, separated by a financial development threshold, form distinct convergence clubs to different steady-state growth paths. Several recent papers question the methodology used in the traditional approach. Evans and Karras [32], for example, show that the traditional approach to estimating convergence is only valid if economies have identical first-order autoregressive dynamic structures and all permanent cross-country differences are completely controlled for. They develop an alternative approach, based on Levin and Lin's [51] unit root test. 2 However, the traditional approach and Evans and Karras' [32] methodology lead to similar conclusions: the 48 contiguous US states and a group of 54 countries show conditional convergence. Lee, Pesaran and Smith [48] consider the predictions of the Solow theory using panel data on output per capita for 102 countries for the period 1960-1989. They examine three notions of convergence: /J-convergence; cr-convergence; and See Chapter 4 for more details of this test. 2 11 whether, i n the time-series dimension, each country is converging to its own steady-state equilibrium, i.e. whether the data contain a unit root. They strongly reject a common steady-state growth rate: across the sample, the variance of the log of per capita output increases with time. Time-series estimates of the speed of convergence are about 20%, but it is argued that these are b i ased upwards. This is confirmed by a range of unit root tests which suggest that the series are non-stationary. The authors cannot reject the hypothesis of no convergence, even to country-specific steady states. Unfortunately, the estimation methods used by all these studies ignore the evolving pattern of the cross-section distribution. No account is taken of the nature of interaction across economies. As Quah ([67] and [75]) argues, /?- convergence, or looking at the coefficients of a cross-section regression, gives no information about the changing nature of that distribution over time. These coefficients represent only average behaviour. A convergence rate neither provides information about transitory growth dynamics across economies nor does it indicate anything about the relative rankings of different provinces. Conditional convergence only shows whether each economy converges to its own steady state, which is different from that of other countries. Similarly, a measure of 3 cross-section variance (<r-convergence) cannot distinguish between a situation in which economies are criss-crossing within the cross-section and a situation i n which the rich remain rich and the poor remain poor. Suppose the income distribution across economies can be represented i n time t as shown by the density distribution in Figure 1. This shows that the majority of economies are i n the middle of the distribution. For some time t + s, s > 0, another density distribution is shown. This is bi-modal: rich economies group together; poor economies group together; and the middle ground is disappearing. The figure illustrates both the external "shape" of the cross-section distribution, and the internal "mobility" or churning within the distribution. Traditional convergence studies do not address either of these data characteristics. But these dynamics provide important information on the dynamics of the poor catching up the rich and on convergence clubs. Conditioning on explanatory variables leads researchers to conclude that it is these variables that determine an economy's position. But, as pointed out by Quah [70], it is the factors determining "club" membership that are important. When different convergence clubs form, factor inputs and social characteristics endogenously align around values determined by each economy's convergence club. 3 12 Quah [69] develops an empirical model that incorporates mobility and shape dynamics. The econometric method uses information from the entire distribution, and imposes no structure, on either growth trends or the expected nature of convergence. Using this model enables investigation of short-run .growth dynamics. It can, for example, yield information on whether there are distinct, identifiable groups (clubs) of rich and poor. W i t h i n the cross-section distribution, the nature of rankings of observations can also be examined, such as whether provinces which were relatively poor yesterday remain so today. Quah [64] uses this model to investigate movements in G D P per capita income in 118 countries, relative to the world average, over the sample period 1962-63 to 1984-85. Steady-state distributions suggest cross-country incomes tend towards extremes of the cross-section. That is, middle-income countries are disappearing in favour of the very rich and very poor. In a later paper, Quah [75] studies mobility across US states—the results suggest convergence, or a unimodal steady-state. However, persistence remains high at the ends of the cross-section distribution. Using similar methods to Quah, Andres and Lamo [1] analyse G D P per capita data for 24 O E C D economies from 1960 to 1990. They find evidence of inertia in income rankings and one group of economies remains persistently at the bottom of the cross-section distribution. 4 Bianchi [12]) studies data on G D P per capita for 119 countries i n 1970, 1980 and 1989. He takes each cross-economy income distribution at time t in isolation and estimates it non-parametrically. He then applies a bootstrap test of multimodality to each. The data show unimodality i n the early part of the sample (early 1960s); but by the end of the sample (late 1980s), the data reject unimodality i n favour of bimodality. In this paper, I use the techniques developed by Quah [69] to examine Canadian provincial and manufacturing industry data. B y considering both provincial and industrial data, I gain a more complete picture of the factors driving the pattern of economic growth and convergence. In an early study of regional aspects of Canada's economic growth, Green [37] finds that, i n 1956, high income The bottom 25% in 1960 fell in a range of 0.26 to 0.54 of the OECD average. By 1990, 25% of OECD countries still had incomes between 0.3 and 0.8 of the OECD average. By contrast, countries in the second and third quantiles fluctuated about the mean, while the richest 25% started with an upper limit of 1.56 and ended at 1.25 times the average. 4 13 provinces tended to have either a large manufacturing sector (where output per worker is high) or a large share of output concentrated in a relatively high productivity resource output. He concludes that not only was the distribution of resources between agricultural and non-agricultural output important, but the composition and productivity within these sectors appeared to be important determinants of income differences. 2.3 E c o n o m e t r i c M o d e l In this Section, I outline the econometric framework developed by Quah (see, for example, [69]), which I subsequently use to examine the evolving nature of the distribution of economic growth across Canada. First I assume that the cross-sectional distribution of the n-dimensional vector of disaggregate variables, Xt, is generated by disaggregate-specific shocks, conditional on aggregate effects. aggregates at time t by Tf 5 Denote the cross-sectional distribution of dis- Then, following Quah [69], I employ a stochastic kernel equation to represent the dynamics of the cross-section distribution. The 6 dynamics of the distribution in Tt are parameterized by (Af, Q), where M is a transition probability matrix and Q is a sequence of quantiles. Each quantile-set pair (Q(t), Q(t+1)) defines a transition probability matrix M of transitions from T\ to Tt+\. The matrix M summarizes information on cross-section mobility; the sequence Q describes the shape of the cross-section distribution. Each element (j, k) of the transition probability matrix, M , indicates the probability of an observation starting i n state j in period t, ending up in state k i n period t +1. The probability of arriving at state j in period t depends only on the state the variable was in one period earlier—the history of arriving at a * These disaggregate shocks have only transitory effects on the aggregate. They could include, for example, a regional demand shock, or a drop in fish stocks or a bumper harvest in the prairie provinces. A kernel function K is usually a symmetric probability density function, such as the normal density. The kernel estimator, with kernel K is defined by, 6 i=l where X{ is a random variable, and h is the window width (or smoothing parameter). The kernel estimator is a series of "bumps" placed at the observations. The kernel function K determines the shape of the bumps while the window width h determines their width. In this estimation, a squared Epanechinikov kernel is used (see Silverman [79]). 14 certain point is ignored. 7 S-period ahead transition probabilities are obtained by multiplying the transition probability matrix M by itself s-times (i.e., M"). Varying the transition period indicates how mobility in the cross-section distribution behaves over time. These s-period transition probabilities contain information about both short and long-run disturbances. A s s changes, the model imposes identifying restrictions in different ways on the transition probability matrix, M. A s a consequence, the information contained i n the transition probability matrix is a function of the transition period. Very short periods provide details of short-run fluctuations, but are relatively uninformative for long-run shocks. A s s increases, a more accurate picture of the impact of long-run disturbances is provided. I present 8 results for different transition periods, enabling a comparison of the relative importance of short and long-run shocks to the data. I also compute the ergodic or steady-state distribution implied by the transition probability matrix M. This is given by the eigenvector associated with the (single) unit eigenvalue. 9 If it exists, the ergodic distribution provides a summary of intradistribution mobility—with aggregate effects removed, it represents movements of disaggregates relative to each other. 2.4 Cross-Section D y n a m i c s In this Section, I investigate the shape of and mobility in the distribution of disaggregate measures of Canadian economic g r o w t h . 10 I examine disaggregation by province and by manufacturing sector. Appendix A lists the sources of data. More precisely, suppose Xt is a random variable that can assume only an integer value {l,2,...,JV},then, , 7 P(X t = k\X -i = j, X _ t t 2 = i,...) = P(X t = k\X -i = j) = t p. ih This is an N-state Markov chain with transition probability pj*. Markov chains form the simplest time-series model for a discrete-valued random variable. This must be weighed against the accompanying loss in "power". In estimating s-year transition probabilities, the true sample size is approximately T divided by s. That is, the number of non-overlapping samples helps to determine the power of the test. As s increases, the number of non-overlapping samples decreases. A distribution is ergodic if one of the eigenvalues of the transition probability matrix M is unity, and all other eigenvalues are inside the unit circle. The long-run forecast for an ergodic Markov chain is independent of the current state. That is, the vector of ergodic probabilities indicates the unconditional probability of each of the different states. Estimation in this thesis is carried out using SHAZAM, RATS and Danny Quah's time series random fields (tsrf) package. 8 8 l0 15 1-year transitions (Number) (130) (129) (125) (126) (128) Steady-state distribution 5-year transitions -0.365 0.90 0.05 0.01 0.01 (Number) (128) (115) (111) (113) (121) Steady-state distribution -0.367 0.78 0.03 0.02 0.02 0.149 0.073 Upper endpoint -0.207 -0.071 0.074 0.09 0.01 0.83 0.10 0.02 0.06 0.14 0.77 0.02 0.16 0.76 0.09 0.205 0.245 0.229 Upper endpoint -0.215 -0.071 0.074 0.02 0.20 0.01 0.76 0.15 0.05 0.06 0.70 0.20 0.22 0.62 0.03 0.02 0.17 0.264 0.175 0.285 0.290 0.02 0.05 0.91 0.172 0.290 0.01 0.02 0.12 0.82 0.204 Table 2.1: Transition Probability Matrix—Personal Income per Capita by Province 2.4.1 Provincial I look at two different measures of economic growth at the provincial level of disaggregation: personal income per capita, and wages and salaries per employee. Personal income Figure 2 plots annual real personal income per capita, relative to the national average, for the Canadian provinces over the period 1926-1993. In 1926 P E I ' s real income per capita was 33% below the national average; by 1993, the difference was only 10%. In contrast, B C ' s real income per capita exceeded the national average by more than 46% i n 1926; by 1993, this figure had fallen to 15%. The shape of the cross-section distribution has also altered, with the dispersion of incomes declining markedly over the sample. Using the (log) ratio of provincial real per capita income to the national average, I estimate a transition probability matrix, M, following the methods outlined i n Section 2.3. Table 2.1 presents the results for two different assumptions about the period of transition. The endpoints of the intervals (quantiles) across 16 the cross-section are designed to ensure that approximately equal numbers of observations begin i n each interval. 11 The upper endpoints of each interval are shown on the top row of each panel. The first column gives the total number of observations with starting points in that income state. for n — b, where n is the number of quantiles. 13 12 I present the results Recall that each (j, k) entry of the transition matrix is the probability of a province in state (or quantile) j in period t transiting to state k in period t + s. The first panel of Table 2.1 contains the annual transition matrix, estimated by averaging the observed one-year transitions over every year, from 1927 to 1993. 14 This shows high persistence for the lower and upper quantiles: the corner entries are at least 90%. The ergodic distribution implied by the oneyear transition function shows an accumulation in the middle states, and a thinning at the upper and lower t a i l s . 15 The second panel of Table 2.1 shows the transition probability matrix produced from estimation with a five-year transition p e r i o d . 16 A s expected, mo- bility increases relative to the one-year transition period: diagonal elements are smaller, and off-diagonals larger. The ergodic distribution becomes more skewed—away from the lower states and towards the upper states (3 and 4)— but remains u n i m o d a l . 17 The changing shape of the steady-state distribution provides an indication of the errors in the short transition period estimates. A T h e intervals are constructed from uniformly distributing the data overits observed range. An alternative approach would be to specify the end-points of each interval, but this would impose (possibly incorrect) structure on the data. Where numbers in the first column differ from each other, it is by at most the number of years in the sample. This is due to rounding errors. I experimented with n ranging from 2 to 6, but found the results varied little. ^Entries in the lower right-hand section of each panel in Table 2.1 show the probability of an economy in a rich state staying in that state in the next period. Entries showing 0 to 2 decimal places are left blank. Estimation of a second-order transition matrix (X(t) depends on X(t — 1) and X(t — 2)) produces a similar distribution in the steady state. In this case, the corner diagonal entries of M are both 0.95. The entries in each column of each panel should sum to one. Although this holds for the one-year transition period (given minor rounding errors), it does not always apply in the longer transition periods. One possible explanation for this discrepancy lies with the sparse nature of the transition probability matrix—there are a non-trivial number of zero entries. I investigate this by experimenting with several different matrices, and find that several zero entries lead to summation errors. With no zero entries, the sum of column entries is equal to one, even for 20-year transition periods. With zero entries, errors begin to appear after only a few years. Estimation with a ten-year transition period produces similar results, with slightly increased mobility over the five-year transition period, but a similar ergodic distribution. I also restrict the sample to the period 1961-1992. Mobility and ergodic distributions vary little from the results based on the longer sample period. 11 1 2 1 3 15 18 17 17 1-year transitions (Number) 0.20 0.92 (109) (131) 0.07 (130) (129) (130) 1-year transitions (Number) (130) (174) (130) (195) 0.25 0.93 0.05 0.40 0.08 0.77 0.15 0.02 Quantile 0.60 0.80 1.00 0.02 0.12 0.84 0.02 0.02 0.98 0.15 0.73 0.12 Quantile 0.50 0.75 0.06 0.01 0.81 0.14 0.20 0.68 0.08 1.00 0.01 0.12 0.92 Table 2.2: Fractile Transition Probability Matrix—Personal Income per Capita by Province one-year transition period is too short to provide an accurate indication of the steady-state outcome. However, it yields similar qualitative results to the longer transition periods. A n alternative method of estimating transition probabilities is to examine the time-invariant fractile matrix and the sequence of quantiles. 18 A fractile matrix possesses elements in which every quantile set has equal measure. When M is fractile, the uniform distribution is always an ergodic limit. Table 2.2 . presents fractile transition probability matrices for provincial income per capita data for two different quantile specifications (which are shown in the top row of each panel). These matrices are approximately the average of a sequence of time-varying fractile transition probabilities. In comparison with the top panel of Table 2.1, the corner entries exhibit more persistence, particularly at the top of the distribution. However, mobility between states appears more pronounced at the lower end of the distribution. 1 8 A matrix M(t) is fractile, if, n (^ lm(t))<t>y,t(ll(t)) M = 4>vA<li(t)) = <t>yAli(t)) Vi. 18 Differences between the matrices in Table 2.1 and Table 2.2 can be attributed to the different methods of estimation. In Table 2.1, the transition probability matrix imposes a uniform distribution at the outset, but loses precision at the ends of the distribution. B y contrast, the fractile matrix forces the data to be uniform in the limit, pushing the mid-lower and mid-higher income groups into the tails of the distribution. This characteristic also helps to explain the relatively large entries in quantiles 2, 3 and 4 in the off-diagonal (Table 2.2): they reflect a tendency for the distribution to become more uniform over time. Figure 3, a chart of the quantiles (associated with Table 2.1), reveals that dispersion and variance across the cross-section is declining. The chart plots the log of provincial income per capita in quantile i over average income per capita against year. The bottom two quantiles display relatively steep gradients, whereas the top quantile appears much flatter: economies i n the bottom of the cross-section distribution have moved up the distribution faster than those at the top have moved down—the gap between richest and poorest is narrowing. However, the ranking of provinces within the cross-section distribution remains largely unchanged over the sample—poor provinces (Newfoundland and P E I ) remain at the bottom of the distribution, rich provinces (Ontario and Alberta) remain at the top. It is this persistence at the ends of the cross-section distribution which is ignored in the traditional studies of convergence. 19 Wages, salaries and supplementary labour income I now study the distribution of wages and salaries per employee across provinces. 20 Annual data on wages, salaries and supplementary labour income (WSSLI) are available for the period 1966-1994. 21 I examine W S S L I per em- ployee, relative to the national average. Fluctuations in the data are concentrated around the early 1980s (Figure 4). Unlike the income per capita data, Appendix C contains results for personal income per capita growth data. The transition probability matrices show mobility is higher than in the levels data, and there is no pattern of persistence. The ergo die distributions are unimodal and symmetric. I also study data on provincial annual employment per capita for 1966-1992: cross-section dispersion changes little over the sample period, and mobility within the distribution is limited. The ergo die distributions suggest convergence towards the bottom half of the cross-section distribution. Appendix C contains results for estimation of transition probability matrices using data on provincial wages and salaries per capita. The pattern of mobility is very similar to that found in the personal income per capita data: persistence at the ends of the cross-section distribution, and a unimodal ergodic pattern. WSSLI, unlike personal income, excludes interest and dividend earnings. 1 9 2 0 2 1 19 1-year transitions (Number) (54) (54) (54) (54) (54) Steady-state distribution 3-year transitions -0.112 1.00 (Number) (48) (49) (47) (47) (49) Steady-state distribution -0.112 1.00 Uppe r endpoint -0.023 0.019 0.118 0.78 0.20 0.264 0.06 0.80 0.15 0.264 0.15 0.85 0.363 Upper endpoint -0.023 0.019 0.118 0.279 0.052 0.65 0.32 0.226 0.22 0.74 0.06 0.279 0.073 0.057 0.33 0.55 0.15 0.079 0.02 0.13 0.64 0.22 0.238 0.21 0.78 0.384 Table 2.3: Transition Probability Matrix—Wages, Salaries and Supplementary Labour Income per Employee, by Province Newfoundland now fluctuates about the average rather than remaining at the bottom of the distribution, reflecting the relatively high unemployment rate i n Newfoundland. 22 The transition probability matrices (Table 2.3) show that economies starting in the bottom quantile are unlikely to move up the distribution. 23 However, there is a very high probability of leaving the middle states, particularly with a 3-year transition period. The ergodic distributions display a very different pattern from the per capita data. In this case, economies tend to diverge to the ends of the cross-section distribution. 24 Three main groupings are visible in the quantiles generated from the W S S L I per employee data (Figure 5). The bottom group displays an upward drift; the I n 1976, Newfoundland's unemployment rate averaged 13.1%; the average across Canada was 7.2%. By 1994, Newfoundland's unemployment rate had risen over 50% to 20.4%, while the national average had risen only 45% to 10.4%. Newfoundland's unemployment rate remained the highest of the ten provinces over the period 1976-1994. Fractile transition probability matrices for this and all subsequent variables studied in this Section are in Appendix B. T h e differences are not attributable to the different sample periods: re-estimation using data for 1966-1992 for WSSLI per capita and WSSLI per employee did not produce significantly different results. 22 23 24 20 middle group is fairly flat (centred about zero); and the top group exhibits a slow decline. Dispersion has fallen over time. This suggests convergence in the traditional sense: a narrowing of the cross-section range over time. However, the ergodic distributions suggest divergence within the cross-section. Although the range across provinces of W S S L I per employee is falling, provinces with relatively low W S S L I tend to remain at the bottom of the distribution; those with a relatively high value stay at the top. This points to a bimodal pattern. In the middle of the cross-section distribution, the mass is approximately zero. Differences i n the ergodic distributions between per capita and per employee data can be attributed to the gap between the number employed and the total population i n the different provinces. In general, those provinces with low population have a high unemployment rate; those with high population have a lower unemployment rate. Migration to the richer provinces may have increased the gap between rich and poor (in terms of income and employment opportunities). This is consistent with Helliwell [40], who finds that income per capita and unemployment rate gaps are strong determinants of m i g r a t i o n 25 and Lee [46], who finds that inter-provincial migration appears to slow down the convergence process. Lee and Coulombe [47] find that output per worker converges faster than output per capita, and attribute this to regional disparities i n unemployment rates. 26 Although I find that although dispersion has fallen, W S S L I per em- ployee diverges within the cross-section distribution. The discrepancy between relatively rich and poor provinces falls. However, at the same time, the probability that the relatively poor can become rich shrinks. 2.4.2 Industrial I now consider disaggregation by manufacturing group. 27 Combining the results from this Section with those for the provincial disaggregation provides a clearer Migration has primarily been directed towards the richer provinces. Appendix E contains results from estimations using data employed by Coulombe and Lee in their studies of convergence across Canadian provinces (see [21], [47] and [22]). In most cases, my results are consistent with the reduction in dispersion of different measures of output and income found by Coulombe and Lee. However, cross-section dynamics differ across variables, and for one measure—output per hour worked—there appears to be divergence within the distribution. Appendix C contains results from estimation for disaggregation by seven industry sectors, examining GDP per person at work and GDP per person hour for the period 1961-1994. Both data sets produce unimodal ergodic distributions. 25 2 6 2 7 21 picture of the evolving pattern of economic growth. It enables an analysis of whether particular industries can explain the pattern of provincial income. I study gross domestic product ( G D P ) per employee (relative to the national average) for 21 manufacturing industries, using monthly data for the period 1961:1 to 1994:11. 28 G D P per employee varies both across industries and over time, as illustrated i n Figure 6. The tobacco and tobacco products industry has the highest G D P per employee, followed by the refined petroleum and coal products, and beverages industries. A t the bottom of the distribution are the leather and allied products, clothing, and furniture industries. The refined petroleum and coal products industry possesses the highest value of capital stock per employee; the lowest value is held by the clothing industry. G D P per employee may, to a large extent, be reflecting the capital intensity of the different industries. Alternatively, it may be an indication of the low value-added in or productivity of those industries at the bottom of the distribution. Table 2.4 presents the transition probability matrices. One-month transition period probabilities show surprising amounts of movement among states, particularly in the middle ranges. However, industries do not move further than one state away from their initial position, and i n the middle ranges, are equally likely to move up or down. The ends of the cross-section distribution show more persistence. A s the transition period is increased to three years, mobility increases across the board. The steady-state distributions reveal an interesting pattern. In the one- month transition period, the steady-state is unimodal. A s the transition period increases, a bi-modal pattern becomes evident—industries tend to group at the bottom or top of the cross-section distribution. A t first glance, this result seems to contradict the provincial results which reveal that the ergodic distribution is unimodal. However, the manufacturing industries which are at the bottom of the cross-section distribution (leather, clothing and furniture and fixtures) do not, i n general, contribute a great deal to provincial G D P . Those industries at the top of the distribution (refined petroleum and coal products, tobacco, beverages and food) are evenly spread across provinces. This suggests that, at least in the long run, there is some support for the traditional theory of convergence. 2 8 Appendix A lists the manufacturing industries used in the estimation. 22 1-month transitions (Number) (1648) (1648) (1651) (1649) (1646) Steady-state distribution 6-month transitions -0.594 0.94 0.06 (Number) (1624) (1623) (1634) (1629) (1627) Steady-state distribution 12-month transition -0.594 0.92 0.09 (Number) (1593) (1594) (1619) (1606) (1599) Steady-state distribution 36-month transition -0.594 0.92 0.09 (Number) (1371) (1382) (1435) (1415) (1400) Steady-state distribution -0.594 0.86 0.14 0.202 0.205 0.214 0.214 Upper endpoint -0.202 0.068 0.325 0.06 0.85 0.09 0.09 0.80 0.10 0.10 0.81 0.09 0.197 0.199 0.201 Upper endpoint -0.202 0.068 0.325 0.08 0.79 0.12 0.11 0.71 0.17 0.17 0.67 0.01 0.16 0.194 0.202 0.201 Upper endpoint -0.202 0.068 0.325 0.08 0.81 0.10 0.09 0.79 0.11 0.11 0.79 0.10 0.202 0.201 0.188 Upper endpoint -0.202 0.068 0.325 0.12 0.69 0.17 0.15 0.66 0.19 0.18 0.70 0.13 0.178 0.199 0.210 1.140 0.09 0.90 0.201 1.140 0.01 0.16 0.83 0.198 1.140 0.09 0.90 0.195 1.140 0.12 0.87 0.199 Table 2.4: Transition Probability M a t r i x — G D P per Employee, 21 Manufacturing Industries 23 In summary, for all disaggregates studied, there is considerable persistence at the extremes of the cross-section distribution. I find that there has been a reduction in cross-section dispersion for measures of provincial income per capita over the period 1926-1992. Moreover, the relative increase in income per capita experienced by the provinces at the bottom of the rankings is larger in absolute value than the reduction at the top. However, there is also widespread persistence—little cross-section mobility—at the extremes of the cross-section distribution. In particular, Newfoundland remains at the bottom of the distribution; Ontario and Alberta remain at the top. B y contrast, there is considerable mobility among those provinces in the middle of the distribution. Data on W S S L I per employee reveal a different pattern. The ergodic distribution of W S S L I exhibits divergence within the cross-section. Those provinces with higher than average unemployment rates and lower than average wage rates diverge from those at the top of the cross-section distribution. Disaggregation by manufacturing industry, using monthly data on G D P per employee from 1961 through 1994, exhibits relatively high mobility even over a one-month transition period. In this case, the ergodic distribution is bi-modal. This could reflect different technologies and varying^ capital intensities across manufacturing industries. 2.5 Robustness In this Section, I study the sensitivity of the results to province size (measured in terms of contribution to national G D P ) . I adjust the data i n different ways before estimating the transition probability matrices. First, I adjust the data for the asymmetric size of the provinces; in particular, I test for any possible bias in the results of Section 2.4 contributed by the size of Ontario relative to the other provinces. Ontario accounts for approximately 40% of national G D P and personal income, and one third of the total population. 29 Second, I study the impact on the results of aggregating data for some of the smaller provinces. Finally, I investigate the role of Ontario i n influencing the manufacturing indusThe shares of GDP and income attributable to Ontario dip slightly in the early 1980s. This coincides with an upturn in Alberta's share of Canadian GDP and income, as the refined petroleum and chemicals industries increase in importance in that province. 2 9 24 try data by aggregating the provincial data by region and comparing the results for different regions. Appendix D contains the charts and transition probability matrices from the estimations performed in this Section. 2.5.1 Provincial data adjustment In the first adjustment to deal with any potential bias introduced by the relatively large size of Ontario, I divide all observations by the corresponding observation for Ontario. In the second adjustment, I subtract Ontario data from the national total. These two adjustments allow a comparison of absolute and relative effects. The absolute adjustment emphasises the size differential; the relative adjustment focuses on any differences in mobility. Thirdly, given the disproportionately small size of some provinces, I aggregate data from some of the provinces to form six different geographic regions: the Maritimes (Newfoundland, P E I , Nova Scotia and New Brunswick), Quebec, Ontario, the Prairies (Manitoba and Saskatchewan), Alberta and British Columbia. Relative to Ontario (1) The first adjustment—dividing the data by the corresponding observation for Ontario—displays a similar pattern of persistence i n the personal income levels data, but produces a slightly different steady-state pattern. The adjustment produces smaller actual numbers—but the fluctuations follow a similar pattern. 30 The adjusted data exhibit more mobility, particularly in the bottom half of the cross-section distribution. The adjusted data also suggest convergence (a reduction i n dispersion) i n income per capita over time, insofar as the one-year ergodic distribution shows a similar pattern to the unadjusted data—it is broadly unimodal. A s the transition period increases, the unadjusted data produces a graph which peaks i n state 4; the adjusted data peak in the final state. This is reflected in the transition probability matrices—there is a higher probability of moving down a state than of staying put with the adjusted data compared with the unadjusted data. Graphing GDP growth for Ontario and Canada between 1962 and 1994 shows that peaks and troughs are largely correlated and tend to occur in the same years. The recession in Ontario in 1981 shows clearly in the adjusted data set, but is barely visible in the unadjusted data. The variance for the Canadian data is slightly higher than that for Ontario. 3 0 25 This is because the adjusted data is divided by data from the highest value province, rather than by the average. Relative to Ontario (2) In the second adjustment, I subtract Ontario data from the observations for other provinces. 31 I find (like the first adjustment) that the steady-state distri- butions suggest a shifting up of the cross-section distribution as the transition period increases. The quantiles for the unadjusted and adjusted data (see Figures 3 and D2) display very different patterns. In the unadjusted data, there is a clear ten- dency for the bottom quantiles to display faster growth than the others. In the adjusted data, no such pattern is discernible. Instead, fluctuations occurring around 1981 (peaks i n the top two quantiles, and troughs i n the bottom three) are greatly magnified. This reflects the marked dip in personal income per capita data i n Ontario during the 1981 recession. These results confirm the importance of Ontario i n driving the national results. Absolute differences produce different quantile patterns, but ergodic distributions do not differ significantly from the unadjusted data. Aggregating some provinces In the third adjustment to the provincial income data, I aggregate observations for some of the smaller provinces to form larger regions. In this case, the results closely resemble those found in the unadjusted case. However, there is greater persistence at the bottom end of the cross-section distribution, which is reflected in the ergodic distributions: observations converge towards the middle of the cross-section distribution, rather than the top end. This is because, when aggregated, the poorer provinces have a more significant effect on the results, contributing to a downward shift i n the unimodal peak. Figure D3, which charts the different quantiles, shows that they divide into three different groups. The ranking of provinces within these quantiles is as follows: the Maritimes and the Prairies at the bottom; Quebec and Alberta in 3 1 I do not convert this data into logs. 26 the middle; and British Columbia and Ontario at the top. Both dispersion and variability diminish over time. 2.5.2 Industry data I present estimates of transition probability matrices for G D P per employee in 21 manufacturing industries, using three different provincial industry aggregates. The data set covers Ontario, the western provinces ( B . C . , Alberta, Saskatchewan and Manitoba), and the eastern provinces (Quebec, New Brunswick, Nova Scotia, P E I and Newfoundland). 32 These estimations enable me to examine more closely the role played by Ontario i n the industry data. Ontario I use annual data on G D P per employee, relative to the manufacturing industry average i n Ontario, over the period 1971 to 1991. The transportation equipment industry contributes the highest value to Ontario G D P (more than 15%), followed by industries manufacturing electrical and electronic products (10%), primary and fabricated metal (9%), and food products (9%). Employee numbers follow a similar pattern. A wide range of industries is represented i n the province. The transition probability matrices (Table D4) show considerable mobility even within a one-year transition period. Once again, most persistence is seen at the ends of the cross-section distribution: observations tend to accumulate in one of two tails. Over different transition periods, the ergodic distribution is generally U-shaped. 33 This suggests that, as for manufacturing industries across Canada, G D P per employee for manufacturing industries i n Ontario tends to diverge to the ends of the distribution. Some sectors, the winners, remain at the top of the cross-section distribution; others, the losers, stay at the bottom. Figure D4, showing the quantiles, looks similar to the quantile chart for manufacturing industries across Canada. Most of the volatility appears i n the top and bottom quantiles. These quantiles are also further apart from the remaining quantiles. See Appendix A for a list of the industries used in the estimation. A six-year transition period produced a similar pattern, with increased mobility at the bottom of the cross-section distribution. 3 2 3 3 27 Western provinces Next I aggregate annual data (1971-1991) for G D P per employee i n the same 21 manufacturing industries for the four western provinces of B . C . , Alberta, Saskatchewan and M a n i t o b a . 34 One- and three-year transition probability ma- trices (Table D5) show a high degree of mobility across quantiles: even over one year, most states can be reached from others. There is more persistence at the bottom of the distribution than at the top. The ergodic distributions display a largely uniform pattern—there is no tendency for G D P per employee to either converge or diverge over t i m e . 35 The chart of the quantiles (Figure D5) shows that again, most volatility is contained in the end quantiles. Eastern provinces Finally, I consider G D P per employee for the eastern provinces—Quebec, New Brunswick, Nova Scotia, P E I and Newfoundland—relative to the eastern average for 1971-1991. The data exhibit considerable volatility. This is particularly true in the chemical and refined petroleum industries. Estimation of one and three-year transition probability matrices (Table D6) shows that mobility at the bottom of the distribution is lower than for the western provinces, with a zero probability of moving up more than one quantile. Ergodic distributions show increasing concentration towards the top end of the cross-section distribution. The chart of the quantiles (Figure D6) shows, once 36 again, that it is the upper and lower quantiles which display most volatility. The degree of volatility has decreased since 1980, and this appears particularly true i n the bottom quantile. Given a few minor differences, adjusting the data to account for the size differences of the provinces, and for the disproportionate role played by Ontario, does not appear to affect the main results. Aggregate fluctuations i n Canada follow a similar pattern to those seen i n Ontario. In Ontario, ergodic patterns produced by estimation of transition probability matrices suggest divergence i n I weight the contribution of each province by the ratio of its population to total 'western' population. '"Estimation with a six-year transition period produces an ergodic distribution which follows a similar pattern, showing a shallow peak in the middle of the distribution. A six-year transition period produces similar results. 3 4 3 6 28 G D P per employee (to extremes of the cross-section distribution). B y contrast, the aggregation of western provinces produces roughly uniform ergodic distributions and the eastern provinces exhibit unimodality or convergence towards the top end of the cross-section distribution. These differing ergodic distributions suggest that the pattern in Ontario may be driving the national picture of divergence. 37 This should not be surprising. Ontario produces a large share of national output in a wide range of manufacturing industries. 38 2.6 Conclusions In this Chapter, I examine the dynamics of economic growth i n Canada, investigating the evolving nature of the underlying cross-section distribution. Most previous studies of economic growth and convergence in Canada focus on aggregate measures, and conclude that convergence occurs across countries (or provinces/states) and across industries. However, the methods used reveal little about fluctuations i n the pattern of economic growth within the cross-section. They tell us nothing about the mobility within the cross-section or the shape of the distribution. Using an econometric framework developed by Quah [69] that incorporates cross-section (disaggregate) dynamics, I explore relationships among cross-section disaggregates for the Canadian economy. Disaggregating data by province, I show that dispersion is declining, and that the cross-section distribution is unimodal for measures of income per capita (relative to the national average). However, persistence exists i n the cross-sectional ordering of relative income levels. Poor provinces remain relatively poor. Rich provinces remain relatively rich. This persistence i n rankings is consistent with findings by Green [37] in his study of economic growth in C a n a d a . 39 He finds that provinces with the highest average output per capita in 1890 still had the highest i n 1956, and those with the lowest continued to remain at the bottom. Differentials i n income per capita T h e Ontario results also indicate the robustness of the bi-modal pattern to the frequency of data. During the period 1961-1994, Ontario accounted for 40% of Canadian GDP. The next largest contributors were Quebec (24%), British Columbia (11%) and Alberta (10%). Ontario and Quebec both have broad industry bases. In these two provinces, no one industry accounts for more than 17% of provincial GDP. By comparison, the chemicals industry contributes more than 28% of Alberta's GDP in 1989. See also Melvin [56] and references cited therein. 37 38 3 9 29 between the top provinces and those at the bottom first widened (between 1890 and 1910) and then narrowed (from 1929). I find a strikingly different pattern for data on wages, salaries and supplementary labour income (WSSLI) per employee. Ergodic distributions are bi-modal, suggesting divergence within the cross-section distribution. Persistence at the extremes is combined with a higher mass at the ends; the middle of the distribution is disappearing. G D P per employee across manufacturing sectors also displays a bi-modal pattern. Manufacturing industries tend to remain either at the top or the bottom of the cross-section distribution. The long-run bi-modal pattern also appears i n Ontario manufacturing industry data, but not for aggregations of eastern or western provinces. This suggests that movements i n Ontario drive the national pattern. It reflects the wide range of manufacturing industries represented in Ontario, and the fact that it is the largest—in terms of G D P — province. A comparison of the results for the provincial income per capita data and for the manufacturing G D P per employee data suggests, however, that it could be government transfer payments that are behind the unimodal long-run distribution. In the absence of these payments, the long-run distribution is bi-modal and the underlying pattern of economic growth i n Canada is one of divergence. For policy makers concerned with provincial inequalities, this essay provides a mixture of good and bad news. O n the one hand, there is some evidence of long-run provincial convergence. O n the other hand, there is no evidence to suggest that poorer provinces or industries can approximate richer provinces or industries in the short or medium run. Even i n the presence of transfer payments, rankings remain unchanged within the cross-section distribution, for provincial and industrial disaggregations. 40 It appears that it is those factors which determine a region's comparative advantage (such as geographic location, climate, resource endowments, etc.) which determine the steady-state distribution of income across Canadian provinces. Although it is arguable as to whether transfer payments were designed to affect rankings. 30 Figure 1 Distribution Dynamics Time /I t+s Incc)me distritjutions \ ^s*" \><^ \ ^ \ \ 31 \ / / LU CO 1I « z oL to O CM O 1 AT 'i t1 m z 1 o » co < m , I o c ; < > O O - CJ *t <p CO O O O O (jeqA)u|/(A)u| 32 Figure 3 Quantiles for Log of Provincial Personal Income per Capita (1928-1992) 0.4 1 — -1 J Year quantO quantl quant2 — - quant3 — 33 quant4 quant5 Figure 4 Log of Provincial WSSLI per Employee (Relative to National Average) 1966-1994 34 Figure 5 Quantiles for Log of Provincial WSSLI per Employee (1968-1993) 0.3 . .5 0 1 — — J - Year quantO quantl quant2 — - — - quant3 — - - — quant4 35 quant5 Figure 6 L o g of Manufacturing Industry G D P per E m p l o y e e (Relative to Industry S e c t o r Average) 1961:1-1994:11 Chapter 3 The Dynamic Effects of Aggregate and Disaggregate Disturbances 3.1 Introduction The traditional approach to economic convergence, typified by Barro and Sala-iMartin [4], Helliwell [40], and Mankiw et al [55], involves studying whether poor economies grow faster than rich economies (after controlling for various institutional details). Evidence of this "catch-up" effect is interpreted as support for the neoclassical growth model, as set out by Solow [80]. According to this view, the growth of income per capita is driven by aggregate (perhaps, technology) disturbances, which determine the long-run steady state of the model. Another strand of the literature, exemplified by Durlauf and Johnson [29] and Quah [74] and [75], argues that by averaging out the cross-section information, important characteristics of the data disappear. These researchers argue that convergence studies should focus upon the degree of interaction between disaggregates. Using econometric methods which capture this interaction or cross-section information, they find that the conventional approach masks evidence of convergence clubs and polarisation. The pattern of cross-economy growth is consistent with the models developed by Durlauf [28], Long and Plosser [53] and Galor and Zeira [35]. In these studies, disaggregate disturbances play an important role in aggregate growth dynamics. In Durlauf [28], for example, local linkages across industries create sequential complementarities which build up over time to affect aggregate behaviour. In Chapter 2,1 find that Canadian disaggregate dynamics contain important information for explaining the pattern of convergence. These dynamics and aggregate economic growth may be the result of either aggregate or disaggregate disturbances—as predicted i n the models developed by Durlauf [28] and others. However, the methodology adopted in that Chapter does not allow identification of different types of disturbance. In this Chapter, I directly identify the dynamic effects of disaggregate and aggregate disturbances using a dynamic restrictions- 37 based technique similar to that of Blanchard and Quah [14] and Shapiro and Watson [78]. I estimate a vector autoregression ( V A R ) system, assuming that an aggregate measure of economic growth and a measure of disaggregate interaction (or mobility) are affected by two types of disturbance. I interpret the first dis- turbance as an index that represents a common transitory regional shock; and the second as an aggregate shock. The disaggregate disturbance is identified as having no long-run impact on the level of aggregate income. A n example of a disaggregate disturbance with this characteristic could be a provincial government fiscal stimulus. This might have short-run aggregate and disaggregate effects, but i n the long-run it leaves aggregate economic activity unchanged. The second (aggregate) disturbance has permanent effects on aggregate income. A n example of this type of disturbance might be a technology shock. The two types of disturbance are assumed to be orthogonal. The interaction variable, which is constructed from information contained in transition probability matrices, is assumed to go to its population mean in the long run. That is, mobility is unaffected by either disturbance i n the long run. If disturbances had long-term effects, it would imply that the level of interaction could drift over time. However, as shown in Chapter 2, this is not the case. The short-run dynamics of the two shocks on interaction are unrestricted. I examine annual provincial and monthly manufacturing industry data. I find that the aggregate shock has a large positive impact on aggregate income per capita i n both the short and the long run. The disaggregate disturbance initially increases income per capita; but the economy quickly adjusts. The effect of this disturbance becomes insignificant within 5 years; and accounts for little of the variation i n the aggregate. \ The short-run impact of the two types of disturbance on the provincial interaction measure is imprecisely estimated. The disaggregate shock provides the strongest impulse; and explains most of the variation. But its effect decays to zero within 3 years. The aggregate shock has a larger medium-term impact, which disappears within 7 years. A similar pattern emerges with the industry data. The aggregate disturbance has a positive long-term effect on the aggregate measure, and explains most of its variation. The effect of the disaggregate disturbance is close to zero. Neither 38 disturbance has much of an effect on the interaction measure. W i t h i n only a few months, the impacts are indistinguishable from zero. The results i n this Chapter, confirm that disaggregate dynamics contain important information for explaining the pattern of economic growth. B y ignoring disaggregate disturbances, researchers ascribe too much importance to aggregate disturbances. The remainder of this Chapter is organised as follows. In Section 3.2, I present the methodology used to identify the disaggregate and aggregate disturbances. In Section 3.3, I discuss the economic interpretation behind the identifying restrictions. I present the results i n Section 3.4; and conclude i n Section 3.5. 3.2 Identification I construct a just-identified structural V A R [14]. in the spirit of Blanchard and Quah The main steps are as follows. Assume that the two disturbances (dis- aggregate and aggregate) are uncorrelated at a l l leads and lags. The Wold Representation Theorem implies that, under weak regularity conditions, a stationary process can be represented as an invertible distributed lag of serially uncorrelated disturbances. In order to identify the underlying disturbances, it is assumed that they are linear combinations of the Wold innovations. Let lnY(t) and Pr(t) denote the log of the aggregate variable and the i n - teraction variable respectively. Assume that lnY(t) and Pr(t) have stochastic trends, but are not cointegrated. Let X(t) be the vector (AlnY(t), Pr(t))' and u(t) be a vector of disturbances (t*i(t), U2(t))'. The usual tests suggest that AlnY(t) is stationary. 41 These assumptions imply that X(t) follows a stationary process, given by oo (3.1) X(t) = A(0)u(t) + A(l)u(t - 1) + . . . = ^ A(j)u(t i=o -j) where E[u(t)u(t) \ = I. Equation (3.1) gives AY and Pr as distributed lags of r the two disturbances. These disturbances are assumed to be pair wise orthogonal and the variance covariance matrix diagonal, normalized to the identity. T h e T h e Dickey-Fuller unit root test allows rejection of the null hypothesis of a unit root. This holds for both data sets, with or without a time trend, and up to (and including) 5 lags of the variable. 41 39 contemporaneous effect of u on Y is given by A(0); subsequent lag effects are given by A(j), j > 1. A s X is stationary, neither disturbance has a long-run effect on A y or the interaction variable. To recover the disturbances from the data, a V A R is estimated and inverted to obtain the following Wold moving average representation, (3.2) X(t) = e(t) + C(l)e(t - 1) + . . . where e(i) - (e (t),e (t))', 1 C(0) = I and E(e(t)e(t)') 2 = Q. From the proof of the Wold Theorem, this moving average representation is known to be unique. Since the disturbances ( « i ( t ) , «2(0) innovations (ei(t),e2(t)), a r e linear combinations of the Wold ui and « 2 can be uniquely recovered. 42 From equa- tions (3.1) and (3.2), the following two conditions hold: e(t) = A(0)u(t) and A(j) = C(j)A(0) for all j. If A(0) is unique, « i and u can be recovered. Three 2 restrictions are imposed on A(0) by, (3.3) E[e(t)e(t)'] = E[A(0)u(t)(A(0)u(t))'], which implies that fi = A(0)A(0)'. The transitory nature of the disaggregate disturbance imposes a fourth: ^2fL 0 A(j) = A(0) 5^=o C ( J ) - 4 3 This gives the aggregate and interaction variables as functions of current and past disaggregate and aggregate disturbances. 3.3 I n t e r p r e t a t i o n The key restriction in the identification procedure is that the disaggregate disturbance has no long-run effects on aggregate economic activity. Examples of such disturbances might include any of the following types of province or industry-specific fiscal stimuli: the introduction of a lump-sum subsidy; the construction of infrastructure (roads, hospitals, etc.); or expenditure on education or social security. Alternatively, the source of the disaggregate disturbance could be consumer optimism i n a particular region. 44 I f this assumption does not hold, the underlying disturbances cannot be recovered. For further details of this non-fundamentalness issue see Lippi and Reichlin [52] and Blanchard and Quah [15]. Blanchard and Quah [14] present a formal argument to show that -4(0) is just identified. * Notice that in some circumstances, an initially localised (perhaps fiscal) disturbance might have spillover effects on other areas—in which case, the long-run impact on aggregate economic activity would be non-zero. Under my identification scheme, these are treated as aggregate disturbances. 42 43 4 40 It is likely, as the above examples make clear, that there are more than two sources of disturbances, perhaps with different effects on the economy. The maintained hypothesis is as i n Blanchard and Quah [14], where it is shown that if all disaggregate shocks are similar i n nature, then they lie i n a space that can be identified using the Blanchard-Quah decomposition. A simple alternative approach to the identification problem is to ignore interaction completely. For example, Bayoumi and Eichengreen [7] estimate the responses to different aggregate shocks in each disaggregate i n turn, ignoring the other disaggregates. Then the responses to shocks across disaggregates are compared. However, this procedure ignores the two-way interaction suggested in the models developed by Durlauf [28] and others—and allowed for in the methodology adopted in this Chapter. 3.4 R e s u l t s I examine both provincial and industrial disaggregate data. I use annual data on provincial income per capita from 1926 to 1993; and monthly data on G D P per employee for 21 manufacturing industries for the period 1974:2 to 1994:11. 3.4.1 45 Provincial data First, I create a summary measure of provincial interaction, or "churning", derived from transition probability matrices (see Chapter 2). This is estimated as follows. Following Quah [68], I use a stochastic kernel equation to represent the dynamics of the cross-section distribution. The dynamics are parameterized by (M, Q ) , where M is a transition probability matrix and Q is a sequence of quantiles or states. Each quantile-set pair (Q(t), Q(t + 1)) defines a transition probability matrix M of transitions from period t to period t +1. Each element (j, k) of the transition probability matrix, M, indicates the probability of an observation starting i n state j in period t, ending up in state k i n period t + 1. Associated with each transition probability matrix is an ergodic or steady-state distribution (where it exists), which gives the probability of being in a particular quantile of the distribution i n the long-run. 45 46 T h e industries are listed in Appendix A. The ergodic distribution can be thought of as representing the unconditional probability 4 6 ' i 41 I use the information contained i n the ergodic distribution to create the measure of cross-section interaction. Specifically, I calculate ergodic distributions for 4-state transition probability matrices, using 7-year overlapping sub-samples of provincial personal income per capita data. For example, the first sub-sample covers the period 1927 to 1933, the second covers 1928 to 1934, and so on. Each ergodic distribution consists of a set of 5 observations (one for each quantile or state boundary) for each sub-sample. I aggregate the transition information contained in each ergodic distribution as follows. Taking each set of 5 quantiles, I subtract the set for the previous year. 47 I take the modulus of these first differences, and calculate the average across the 5 moduli to create a measure of churning across each quantile set. I assign this interaction measure to the mid-year of that sub-sample. 48 In this way, I create a time-series measure of cross-section mobility for the period 1931-1988. I use this measure, a summary of the probability of moving between states within a 7-year period, as an indication of the interaction or churning among disaggregates. Next, I estimate a bivariate V A R system, using data on aggregate personal income per capita growth and this interaction measure. include a constant and time trend. 49 I allow for 2 lags, and 50 The dynamic effects of the aggregate and disaggregate disturbances on aggregate income per capita and the provincial interaction measure, along with the two-standard deviation bands, are shown i n Figures 7 to 1 0 . 51 The verti- cal axes show either the log of aggregate income per capita or the interaction measure; the horizontal axes show time in years. The aggregate disturbance has a cumulative effect on per capita personal of being in a particular quantile at any given time. F o r example, I subtract the set for the sub-sample 1927-1933 from the set for the subsample for 1928-1934, and so on. F o r example, for the sub-sample 1928-1934, I assign the measure to 1931, and so on. By assigning the measure to the mid-year of the sub-sample, rather than the beginning or end, the aim is to create a measure of mobility around a certain point. However, if there are big changes in mobility at the ends of the sub-sample, it may be misleading to assign the measure to the mid-year. I assume that, given the overlapping nature of the sub-samples, it is unlikely that such occurrences will distort the main results. Specifically, I use 100 x ln(income growth) and 100 X (the interaction measure). The interaction measure is the average of the modulus of the first differences of probability measures. I experimented with different lag lengths and found the results largely invariant. The standard deviation bands are obtained by 1,000 bootstrap replications as follows. First, I estimate the VAR and identify the two types of disturbance using the dynamic restrictions. Second, I draw with replacement from the distribution of the fitted errors. I use these to generate 1,000 sample replications. For each of these, I identify the disaggregate and aggregate disturbances, and compute the associated impulse responses. 47 48 49 S 0 51 42 income. After about 3 years, the peak response is almost twice the initial one. Then the effect decreases, and stabilizes to its long-run level after 7 years. The disaggregate disturbance initially increases income per capita (although the impact is much smaller than for the aggregate disturbance), but the effect falls sharply within 3 years, and disappears inside 5 years. The response of the interaction measure to an aggregate disturbance is imprecisely estimated. The peak response is around 3 years; its long-run impact is zero, reached after about 7 years. Initially, the impact of the disaggregate disturbance on provincial interaction is quite strong (and positive), but it decays very quickly. The response dies out completely within 3 years. In this case, the impact of the disaggregate disturbance is much larger than that of the aggregate disturbance. Since the interaction measure is stationary, neither disturbance has a persistent effect on provincial interaction. Even an aggregate disturbance, such as a productivity shock, has no long-run effect. Table 3.1 presents variance decompositions for the income per capita and provincial interaction measures. These assess the relative contribution of the different disturbances to fluctuations i n the two measures at various horizons. Two standard deviation bands are shown i n parentheses. The fc-period ahead forecast error i n aggregate income per capita is defined as the difference between the actual value of personal income per capita and its forecast from equation (3.2) as of fc years earlier. The forecast error is due to both aggregate and disaggregate disturbances i n the last fc years. The number under aggregate income per capita at horizon fc (fc = 0,1, 2,4,8,12,20, oo) gives the percentage of variance of the fc-year ahead forecast due to the aggregate disturbance. Impact and infinite-horizon shocks are given by the 0-year and oo (defined as the sample length) ahead forecasts respectively. The contribution of the disaggregate disturbance is 100 minus this number. Two-standard deviation bands are shown i in parentheses. The results suggest that the aggregate disturbance can explain nearly all the variation in aggregate income per capita. However, the disaggregate disturbance appears to explain most of the variance i n the interaction measure; the contribution of the disaggregate disturbance varies from 75.7% to 80.5% over different horizons. The aggregate disturbance explains little of the variance in the interaction variable. 43 Percentage of Variance Due to Aggregate Disturbances Horizon (years) Income per Capita Interaction 0 0.940 0.195 (0.77, 1.00) (0.02, 0.43) 1 0.928 0.195 (0.79, 1.00) (0.02, 0.43) 2 0.956 0.231 (0.87, 1.00) (0.06, 0.46) 4 0.979 0.243 (0.94, 1.00) (0.08, 0.50) 8 0.989 0.243 (0.96, 1.00) (0.08, 0.50) 12 0.992 0.243 (0.98, 1.00) (0.08, 0.50) 20 0.995 0.243 (0.98, 1.00) (0.08, 0.50) oo 0.998 0.243 (0.99, 1.00) (0.08, 0.50) Percentage of Variance Due to Disaggregate Disturbances Horizon (years) Income per Capita Interaction 0 0.060 0.805 (0.00, 0.23) (0.57, 0.97) 1 0.072 0.805 (0.00, 0.19) (0.56, 0.96) 0.044 2 0.769 (0.00, 0.12) (0.53, 0.93) 4 0.021 0.757 (0.00, 0.06) (0.49, 0.92) 8 0.011 0.757 (0.00, 0.03) (0.48, 0.92) 12 0.008 0.757 (0.00, 0.02) (0.48, 0.92) 20 0.005 0.757 (0.00, 0.01) (0.48, 0.92) 00 0.002 0.757 (0.00, 0.00) (0.48, 0.92) 1 Table 3.1: Variance Decomposition of Aggregate Income per Capita and Provincial Interaction 44 The response of the aggregate measure to an aggregate disturbance is consistent with the pattern suggested by a neoclassical growth model, where long-run changes i n the level of output are driven by technology shocks. The humpshaped response of aggregate income per capita to a disaggregate disturbance is similar to the response of output to a demand disturbance i n Blanchard and Quah [14]. The forecast error variance decompositions are similar to those found in the literature (see, for example, Blanchard and Quah [14] and Cogley and Nason [20].) 3.4.2 Industrial data 1 repeat the process outlined above using the monthly manufacturing industry data on G D P per employee. For the interaction measure, I estimate 4-state transition probability matrices for overlapping samples, each comprising 49 observations. I derive the summary statistic as before, as the average of the modulus of the first differences of the resulting ergodic distributions. Using the G D P per employee and interaction measures, I estimate a bivariate V A R allowing for 12 lags, and including a constant and time t r e n d . 52 The resulting impulse responses and associated two-standard deviation bands are shown i n Figures 11 to 14. The aggregate disturbance has a positive effect on G D P per employee. The initial response is followed by a decline to just under half its original level, recovering slightly before stabilizing at its long-run level after about 3 years. After an initial decline, the disaggregate disturbance increases G D P per employee to its starting level, but the effect is small and insignificant, and disappears within 2 years. The effect of the aggregate disturbance on interaction is insignificantly different from zero, even in the very short term. The disaggregate disturbance has a strong (positive) initial impact, but the impact declines quickly, disappearing within 2 years. These results are broadly similar to the results obtained using the provincial data: the aggregate disturbance has a long-run positive effect on the aggregate, but no long-run effect on the interaction measure. The disaggregate disturbance Again, the results are largely invariant to different lag lengths and the inclusion of a time trend. 8 2 45 i (by construction) has no long-run effect on either variable. For the industry data, it has an insignificant effect on the aggregate even in the short run. This insignificance could be because the underlying disturbances are regional rather than industry-specific. Variance decompositions for various horizons are shown i n Table 3.2. The aggregate disturbance explains nearly all the variation i n aggregate G D P per employee; the disaggregate disturbance explains most of the variation i n the interaction measure. 3.5 Conclusions In this Chapter, I propose a technique to examine the impact of aggregate and disaggregate disturbances on measures of aggregate economic growth and interaction across disaggregates. The interaction variable is constructed from information contained in transition probability matrices. I use data on provincial income per capita and manufacturing industry G D P per employee. The approach is motivated by those studies (such as Durlauf [28] and others) which argue that the interaction across disaggregates contains explanatory information for the aggregate, and vice versa. If the disaggregate disturbance can be shown to have some effect on aggregate economic activity, then this strengthens the case for including disaggregate information i n an explanation of aggregate growth. I find that the aggregate disturbance has positive long-run effects on provincial income per capita and G D P per employee, and explains most of their variation. The disaggregate disturbance only matters for fluctuations i n aggregate income at business cycle horizons (up to 5 years). It has little effect on G D P per employee. However, the disaggregate disturbance has important short and medium-run effects on the interaction measures. In contrast, the aggregate shock contributes little. In the long run (by construction), interaction is unaffected by either disturbance. I argue i n Chapter 2 that an explanation of economic growth and/or convergence should include disaggregate information. The results from this Chapter confirm that disaggregate disturbances contain important information for aggregate economic activity at business cycle horizons. Unfortunately, the traditional method of considering convergence concentrates on the explanatory power of 46 Percentage of Variance Due to Aggregate Disturbances Horizon (months) G D P per Employee Interaction 0 0.955 0.049 (0.59, 1.00) (0.00, 0.40) 1 0.941 0.056 (0.57, 0.99) (0.00, 0.41) 6 0.962 0.084 (0.75, 0.99) (0.05, 0.44) 12 0.978 0.090 (0.84, 0.99) (0.08, 0.45) 24 0.986 0.093 (0.90, 0.99) (0.08, 0.46) 36 0.094 0.990 (0.93, 1.00) (0.08, 0.46) 48 0.993 ! 0.094 (0.95, 1.00) (0.08, 0.46) 60 0.994 0.094 (0.96, 1.00) (0.08, 0.46) 0.094 120 0.997 (0.98, 1.00) (0.08, 0.46) oo 0.094 0.998 (0.98, 1.00) (0.08, 0.46) Percentage of Variance Due to Disaggregate Disturbances Horizon (months) G D P per Employee Interaction 0 0.045 0.951 (0.00, 0.39) (0.58, 1.00) 1 0.944 0.059 (0.01, 0.39) (0.56, 0.99) 6 0.038 0.916 (0.01, 0.25) (0.53, 0.95) 12 0.022 0.910 (0.01, 0.16) (0.53, 0.92) 24 0.014 0.907 (0.01, 0.09) (0.53, 0.92) 36 0.010 0.906 (0.00, 0.06) (0.53, 0.92) 48 0.007 0.906 (0.00, 0.05) (0.53, 0.92) 0.006 0.906 60 (0.00, 0.04) (0.53, 0.92) 120 0.003 0.906 (0.00, 0.02) (0.53, 0.92) oo 0.002 0.906 (0.00, 0.02) (0.53, 0.92) Table 3.2: Variance Decomposition of Aggregate G D P per Employee and M a n ufacturing Industry Interaction 47 aggregate shocks. Under the neoclassical growth model, long-run changes i n output are driven by technology shocks. To gain a more complete picture of convergence in the short to medium run however, the researcher should include disaggregate information. 48 Figure 7 Impulse Response: Aggregate Disturbance on Aggregate Income per Capita 10 11 12 Year 49 13 14 15 16 17 18 19 20 Figure 8 Impulse Response: Disaggregate Disturbance on Aggregate Income per Capita YDisagg yDisagg 95% yDisagg 5% Year 50 Figure 9 Impulse Response: Aggregate Disturbance on Interaction -4 + -6 + . J 8 Year 51 Figure 10 Impulse Response: Disaggregate Disturbance on interaction 10 tDisagg — tDisagg 95% | - - - tDisagg 5% 7.-8" 9 10 11 12 2+ Year 52 13 14 15 16 17 18 19 20 j Figure 11 Impulse Response: Aggregate Disturbance on Aggregate GDP per Employee 2 i — 1 0.4 4- 0.2 | 53 Figure 12 Impulse Response: Disaggregate Disturbance on Aggregate GDP per Employee 0.6 i l\ 0.4 i i . 11 * 0.2 I, 1 \ m i o < o « 3 r ~ r ^ o o o o c n o ) o -0.2 Q. E o w O a. a. a -0.4 yDisagg yDisagg 95% yDisagg 5% (3 -0.6 + -0.8 -H -1 + -1.2 Month 54 Figure 13 Impulse Response: Aggregate Disturbance on Interaction Measure 0.4 0.3 -§ 55 Figure 14 Impulse R e s p o n s e : Disaggregate D i s t u r b a n c e o n interaction M e a s u r e 0.7 0.6 0.5 0.4 0.3 tDisagg tDisagg 95% tDisagg 5% 0.2-tt 0.1 -H co o> i n , T»/-'r-»'*'co * o> ' w " i - -0.1 -0.2 Month 56 r~ co o> m i - r-~ Chapter 4 Testing For Canadian Unit Roots: A Panel Data Approach 4.1 Introduction To apply standard unit root tests to panel data, a researcher must aggregate across cross-sections, discarding the cross-sectional or disaggregate information. Recently, more powerful unit root tests have been developed for panel data, which allow cross-sectional information to be i n c l u d e d . 53 Levin and L i n ( L L ) [51] and Quah [75] show that implementing a unit root test which incorporates both cross-section and time-series information can dramatically improve power. In this Chapter, I present the results from a variety of unit root tests on panel data. I compare the results of the standard Dickey-Fuller method with those from three tests designed specifically for panels. Unlike the Dickey-Fuller test, the panel data unit root tests have limiting normal distributions. The first test takes the standard approach to unit root tests, but, to allow for fixed effects, the initial observation for each individual time-series is subtracted from each subsequent observation i n that series. The unit root hypothesis is tested by computing a standard t-ratio using data that is averaged across the cross-section. I refer to this as the "panel Dickey-Fuller" test (Breitung and Meyer ( B M ) [16]). The second test requires the data to be adjusted by subtracting the crosssection mean. The test statistic is the average of the individual coefficients on the lagged dependent variable. I refer to this as the "£-bar" test ( L L [51]). I extend the L L [51] methodology to allow for potentially more complex common effects across individual series. Using a dynamic index model developed by Quah and Sargent [76], I identify multiple common factors and subtract the resulting indeces (in place of the cross-section average) from the data. In the third approach, the "t-bar" test of Im, Pesaran and Shin (IPS) [42], the initial observation for each individual time-series is subtracted from each T h e importance of including cross-sectional information is illustrated in the context of economic convergence in Chapter 2. There, I show that ignoring cross-section information can mask underlying patterns of mobility and shape. 53 57 subsequent observation in that series—in a similar manner to the panel DickeyFuller test. The test statistic is a (small sample adjusted) average of the DickeyFuller test statistics for the individual series. I examine Canadian annual data on gross provincial product ( G P P ) for the period 1961-1990 and monthly data on gross domestic product ( G D P ) per employee for 19 manufacturing industries over the period 1974:2 to 1994:11. When I apply the standard Dickey-Fuller test, to data averaged across the crosssection, I cannot reject the null hypothesis of a unit root i n either the provincial or industrial data. In contrast, none of the three panel unit root tests suggest that the data are characterised by a unit root. The presence of a unit root has important implications for macroeconometric studies. If data contain a unit root, shocks will have permanent effects; where there is no unit root, shocks will have only temporary effects. 54 The results from this Chapter suggest that there are no permanent shocks in Canadian annual G P P and monthly manufacturing industry G D P per employee data. In a seminal paper, Granger [36] argues that the aggregation of a multivariate time series process can produce a univariate process that has fundamentally different properties from the multivariate process. The results reported here are consistent with Granger's finding: aggregating the data (in order to use an existing test) can radically alter the dynamics of the process. The rest of this Chapter is arranged as follows. In Section 4.2, I review the related literature. In Section 4.3,1 outline the theoretical methodology for each of the three panel data unit root tests used i n this Chapter. I present the results in Section 4.4, and conclude in Section 4.5. 4.2 U n i t R o o t s and P a n e l D a t a A linear process contains a unit root if one of the roots or eigenvalues of the autoregressive polynomial is unity and all other eigenvalues are inside the unit circle. Hypothesis tests concerning the coefficients of non-stationary variables cannot be conducted using tests based on the asymptotic t- and .F-distributions. The distributions of the t- and F-statistics are non-standard. 1 ! Dickey and Fuller [25] and [26] test for unit roots by deriving a non-standard S4 See, for example, Nelson and Plosser [57]. 58 ' distribution using Monte Carlo simulations. They consider three different models: (4.1) A y , = «5j/,_i + e, (4.2) A y , = a + «5y,_, + e, (4.3) A y , = a + 6y -i + yt + e . t t The first is a pure random walk model, the second adds an intercept, and the third includes both an intercept and linear time trend. If 6 = 0, the {y,} sequence contains a unit root. The test involves estimating one (or more) of the equations above using O L S to obtain the estimated value of 8 and its associated standard error. Comparing the resulting ^-statistic with the appropriate critical value (reported in the Dickey-Fuller tables) allows the researcher to determine whether to reject the null hypothesis that 6 = 0. The critical values depend on the form of regression and sample size. The same critical values are used to test for the presence of a unit root in higher order equations. A n augmented Dickey-Fuller ( A D F ) test is used when e, is serially correlated. The empirical model contains a number of lags of the dependent variable, p (4.4) A y , = a + <5y,_i + ^ 6 Ay L t L + ft where p is the order of autoregression. In Section 4.4.3 below, I estimate equation (4.4) (and a similar equation including a time trend) using Canadian annual data on G P P for 1961-1990, and monthly data on G D P per employee for 19 manufacturing industries over the period 1974:2 to 1994:11. In order to be able to use the Dickey-Fuller test, I average the data across the cross-section. I find that I cannot reject the null hypothesis of a unit root for either of the data series. These results are consistent with those of Otto and Wirjanto [58], for example, who test quarterly Canadian macro time series (1955:1-1988:4) for seasonal and non-seasonal unit roots. Their results suggest that most of the series, i n cluding G D P , contain a unit root. Banerjee, Lumsdaine and Stock [2], on the other hand, compare the null hypothesis of a unit root with stationarity about i 59 a broken trend. Using quarterly Canadian G N P data (1948:1-1989:2), they reject the null hypothesis in favour of the trend-shift alternative, with a break i n 1981:3. 55 Their interpretation is that the recession of the early 1980s led to a permanent downward shift in the trend growth rate. After the recovery, output is stationary around its original growth path. However, averaging across the cross-section loses cross-section or disaggregate information. L L [51] show, for example, that, i n the absence of individualspecific effects, the power of the standard Dickey-Fuller test is low for short time series ( T < 50). Allowing for both individual-specific intercepts and time trends, the Dickey-Fuller test has very low power even for longer time series ( T < 100). In a theoretical paper, Quah [71] examines the characteristics of the timeseries unit root coefficient estimator in panel data with simultaneously extensive cross-section and time-series dimensions. He finds that the asymptotic distribution of the estimator is neither Dickey-Fuller (as for standard time series analysis) nor normal and asymptotically unbiased (as for standard panel data analysis). It is consistent and asymptotically normal, but has a noh-vanishing bias i n its asymptotic distribution. Three recent papers propose alternative unit root tests, designed specifically for panel data, which evaluate the null hypothesis that each individual series has a unit root against the hypothesis that all series are stationary. Unlike the standard Dickey-Fuller test, each of these panel data tests has a limiting normal distribution. B M [16] apply their panel Dickey-Fuller test to German wage d a t a (19721987). The initial observation for each individual time series is subtracted from each subsequent observation in that series before a standard i-test is used to test the unit root hypothesis. To obtain the limiting distribution of the test statistic, B M [16] let the number of cross-sectional units tend to infinity. In this case, the null distribution becomes asymptotically normal rather than non-standard, as when the number of time periods approaches infinity. Allowing for both a i, linear time trend and random time effects, the unit root hypothesis cannot be rejected for firm and industry-level wages. L L ' s [51] cS-bar test statistic is an average of the individual coefficients on 5 8 The result is sensitive to the number of autoregressive lags used. 60 the lagged dependent variable. 56 Monte Carlo simulations indicate that the normal distribution provides a good approximation to the empirical distribution of the 6-bai test statistic, even for relatively small panels. The null 57 hypothesis imposes a cross-equation restriction on the first-order partial autocorrelation coefficients, enabling the S-bai test procedures to yield much higher power (against stationary alternatives) than performing a separate unit root test for each individual. Simulations indicate that this increase in power holds when pooling the time-series data for even a small number of cross-sectional units (JV > 10). W i t h individual-specific intercepts and trends, the £-bar test is powerful with relatively moderate-sized panels (JV = 10 and T = 50 or N = 25 and T = 25). The 6-bas test procedures are most appropriate for data with 10 — 250 individuals, and 25 — 250 time series observations per i n d i v i d u a l . 58 In an application of the S-bai test, MacDonald [54] examines a group of O E C D real exchange rates over the recent floating experience (1973-1992). Regardless of whether or not a time trend is included, he rejects the null hypothesis of a unit root. Using standard Dickey-Fuller tests, on the other hand, a unit root cannot be rejected. The t-b&i test is a small-sample adjustment to the average of the DickeyFuller test statistics for each individual series (calculated after first adjusting the data by subtracting the initial observation from each time-series observation). IPS [42] show that it is valid with heterogeneity across groups and i n the presence of residual serial correlation across time periods. Under the null hypothesis of a unit root, the t-bar statistic has a standard normal distribution for JV (the number of cross-sectional units) and T (the number of time periods) sufficiently large and N/T —* 0, and diverges to —oo under the alternative hypothesis of stationarity. Using Monte Carlo simulations, IPS [42] show that the i-bar statistic is superior i n small samples to the 6-bai test: i n general, the 6-hax test does not converge to the standard normal distribution, even as JV —• oo, T —• oo and N/T - » 0. Lee, Pesaran and Smith [48] present an application of the t-baii test using The data are first adjusted by subtracting the cross-section average from each observation. L L [51] show that the 5-bar statistic has a limiting normal distribution as both the crosssection and time-series dimensions of the panel grow large. The standard error of the estimated empirical power was less than 0.01 in all cases (with 10,000 replications). , Existing tests are preferable for samples with either longer time-series or 'more crosssections. 5 6 5 7 : 58 61 output per capita data for 102 countries over the period 1960-1989. They reject the existence of a unit root. Using a standard Dickey-Fuller test, they find that the unit root hypothesis is only rejected for 8 out of the 102 countries. In the next Section, I outline in more detail the methods used to calculate the unit root tests: the panel Dickey-Fuller test developed by B M [16], the 6-bar test developed by L L [51] (and my extension), and IPS's [42] t-bar test. 4.3 M e t h o d s 4.3.1 T h e panel D i c k e y - F u l l e r test Assume that the series yu, has a finite order autoregressive representation, where the autoregressive parameters do not vary over time and individuals i n the crosssection. Consider the AR(1) model given by, (4.5) y it = % , t - i + ( l - S)m + c it for t = 1,.. . , T and i = 1,... TV, with E(yi ) t — / i , - . The €,-« are uncorrelated errors distributed TV(0,cr ). Under the alternative hypothesis of no unit root, 2 \6\ < 1, the O L S estimate 6 is biased against 6 = 1 (leading to a loss of power). To overcome this problem, subtract the initial observation, y,o, from both sides of equation (4.5) to get, (4.6) x it = 6x t-i + T) it it for t = 2 , . . . , T , i = 1 , . . . , TV, where x it Tin = at + (1 — $)(yio — = yu - yto, = Vi,t-i - Vio and B M [16] show that under the null hypothesis of a unit root the asymptotic bias disappears. A conventional t-test on the O L S estimate 6 is used to test Ho : 6 = 1. To introduce a linear time trend, (4.7) yu = 6y -i 59 the model is modified to give, + y*t + /x? + e itt it where 7* = (1 — 6)y and /ij = (1 — 6)fii + 6y. 60 Ignoring individual effects leads to biased estimates of 8 and 7* under the alternative hypothesis of stationarity when using O L S . As for the basic model, subtract the initial observation from both sides of equation (4.7) to get, 59 O f the form E(yi ) = Pit + ft. 60 This is derived from (1 - 6L)(m + -ft) - (I- t 6)(m + ft) + 6-y. 62 J (4.8) x = <5x,-, _i + y*t + e - (1 - 6)(y it t it i0 - m). where xu and x , t - i are defined as before. Under the null hypothesis of a unit t root, the last term vanishes and O L S estimates of 8 and j* are consistent. The distribution of the t-statistic for 6 = 1 does not depend on individual effects. The limiting distribution for 7Y —• oo does not change if a time trend is included. (Although the limiting distribution for T —• oo is affected by including a time trend.) For wide variation i n individual effects, the performance of the panel DickeyFuller test is superior to the standard A D F approach (when cross-sectional information is important). However, it does not permit the coefficient on j/,,t_i or the lag order to vary across individuals: all individual effects are assumed to be incorporated by subtracting the initial time-series value from all subsequent observations. 4.3.2 The 6-bar test Let a stochastic process be represented by {yu} for individuals i = 1 , . . . , 7Y over time periods t = 1,.. . , T . Assume (as for the panel Dickey-Fuller test) that all individuals in the panel have identical first-order partial autocorrelation, but all other parameters of the error process are unrestricted across individuals. L L [51] present 6-bar test procedures for three different models: when the series yu (i) has zero mean for every individual in the panel (equation (4.9)); (ii) has an individual-specific mean (equation (4.10)); and (iii) has an individualspecific mean and an individual-specific time trend (equation (4.11)). (4.9) Ay = it Styu^+Cit .j (4.10) Ayi« - a,- + <5,t/j,_i + Cit (4.11) A j / a = a,- + jit + 6 u-i iy + Ot where —2 < 8i < 0 for all i = 1 , . . . , 7Y. The disturbance Cit is distributed independently across individuals and follows a stationary invertible A R M A process for each i n d i v i d u a l . 61 Under the null hypothesis, each individual time series Given by, oo 63 has a unit r o o t . Under the alternative hypothesis, the process {yu} is trend- 62 stationary for each individual in the p a n e l . 63 I outline the test procedures for Model 1. The (5-bar test requires the data to be generated independently across in1 dividuals. This assumption can be relaxed to allow for a limited amount of time-dependence via time-specific aggregate effects. As a first step, L L [51] remove the influence of aggregate effects by subtracting cross-section averages from the data. The next step calculates orthogonalized first differences and lagged levels for each individual i n the cross-section by estimating the following equations (using the adjusted data), Pi (4.12) i = Ax it - ^2 it e l i L Ax i t L L=l Pi (4.13) V it = Xit-l - ^2 82iL&Xit-L- L=l where xu = yu — yu and yu = j ^ l ^ L i 2/>'«- ^ag orders (p.) can vary across individuals. 64 Then, regress the orthogonalized innovations on the orthogonal- ized lagged level, (4.14) e it = 6iVit-i + e . it To control for heterogeneity across individuals, in and vu-i are normalized by dividing by the estimated residual standard error from equation (4.14), form eu and vu-i, (4.15) e it = 65 to where, |^ 0~ei Under Model 1, Si = 0 for all t = 1,... N; under Model 2, Si = 0 and cr,- is zero, for all i = 1, , N; under Model 3, tf; = 0 and 7; is zero, for all t = 1 , . . . , N. Under Model 1, Si < 0 for all t , . . . , N; under Model 2, Si < 0 and a; is non-zero for all N; under Model 3, Si < 0 and 7; is non-zero for all t,..., TV. T h e appropriate lag order p; is selected as follows: for a given sample length T, a maximum lag order Pmax is chosen, and then the t-statistics on the On are used to determine if a smaller lag can be used. The regression standard error is derived from, 62 63 64 6 5 T *li = T _ ._ (e.-t-g.'tfrt-l) t= +2 2 1 Pi 64 and, - —1 Wit (4.16) vtt-i - ——. Cei , Next, calculate the long run standard deviation (tr ,-) as the variation of A i j | x at zero frequency. If the data contain a time trend (Model 3), the trend is 66 removed before a i is calculated. 67 x In the next step, the ratios of long-run to short-run standard deviations for /C£Li *«) each individual (Si = ? " ) and the average ratio for the panel (5 = are calculated. The average ratio is used to adjust the mean of the unit root * test statistic. i Finally, estimate the following equation, = <5,-iS,-t_i + € . (4.17) e it it This enables several statistics to be calculated: the least squares estimate, 6; the standard error of the regression, cr ; the reported standard error of 6, e and the regression ^-statistic for testing HQ : 6 = 0, ts66 RSE(6); The 6-bar test statistic 6S T h e long run standard deviation is given by, R T *li = jTZJ ]C e=2 A X '' + ^2 2 W j R T ( T~{ L i=l X} \ A \ x i t=2+Z t A x i t ~ L I / where wj^ L/(K + 1). (Ax,t — Axit), where Axit is the average value of Aijt for individual i, is substituted for Ax,-, in the equation given in the previous footnote. These statistics are defined as, L 67 68 EN v^T l—ii=\ Zjt=2+ / N u Pi it-l T \ .=1 t=2+p; \ i = l t=2+p; / / and t = 6 / 1 2 -. RSE(S) 65 tg is then given by, (4.18) t f t s = - NTS &~ RSEtS) fi*. - i -±*L f 2 K where T = ( T — p — 1) is the average number of observations per individual i n the panel and p = 2~^LiP»- The mean and standard /i£ and cr£ respectively, are given in L L [51]. deviation adjustments, 69 L L [51] show that the 6-bar test statistic is distributed JV(0,1) in the limit, so the critical values of the standard normal distribution are used to test the null hypothesis that Si = 0 for all i = 1 , . . . , JV. L L [51] implicitly assume that there is a single aggregate common factor that has an identical impact on all individuals. I allow for more complex contemporaneous correlation i n the 6-bar test by, at the first step, subtracting index values (rather than the cross-section average) from the data. The indices are calculated using the methods of Quah and Sargent (QS) [76]. 70 The QS approach does not specify in advance any particular pattern of co-movements; it uses the orthogonality properties of permanent and common components to characterise the model. Each i = 1 , . . . , J V is assumed to be affected by J common factors Uj, j = 1 , J (common across all i), and an idiosyncratic disturbance zt, i = 1 , . . . , JV (specific to each i and orthogonal to all other i). The common factors are assumed to impact on j/,- with M a lags. The first differences of each Uj are pairwise orthogonal, and have a finite autoregressive representation (of order M ). g The idiosyncratic disturbance, z,-, is a finite order autoregression (of order JWj,). Translating the model into state space f o r m , ' . 71 QS [76] exploit the fact that the common factors have dimension O(J) (J < JV), which is independent of O(N). They use standard Kalman smoother calculations to obtain conditional T h e adjustments vary with each model. In calculating the adjustments, L L [51] set K - 3.21T / and p.- = 0 for each »' = 1 TV. Gaussian random numbers with unit variance are used to generate 250 independent random walks of length T +1. These data were used to construct the sample statistics S, 6, <T , RSE(6) and tg. Based on 25,000 replications, 69 1 3 £ the adjustment factors u* and a*, are estimated as the mean values of ——-4 f J T r- and [NTS&Z P-SE(S) J T the standard deviation of [t - NTS&r^RSEiS)] respectively. Appendix F outlines the Quah and Sargent [76] methodology in more detail. State space representation is a means of summarizing finite processes. More details of the translation process are in QS [76]. , i s 7 0 7 1 66 expectations. 4.3.3 72 Iteration on this scheme is the E M algorithm. T h e i-bar test Consider, as before, a sample of N cross-sectioned units observed over T time periods. For each cross-sectional unit and each time period, assume that the stochastic process yu is generated by the first-order autoregressive process, (4.19) y = (1 - <j>i)m + <£,y,\t-i + e,t it where (it = (i + 77,- and e,-( = Ai(L)eu a l 0 with A{(L) = a,o + anL + aiiL 2 + ..., = 1, ^2"j° o a.jj < 0 0 , i = 1,2,..., N and L is a lag operator. Assume that = the e,t are independently distributed across both groups and time periods with zero mean and finite variance. A s for the panel Dickey-Fuller test, subtract the initial observations from each side of equation (4.19) to get, (4.20) x = a + 6 Xi t-i where x = yu - Via, it i it i t + e it = Vi,t-i - Vio and <*• = (1 - 4>i)(l*i - 2/io). Taking first differences, (4.21) Ax it = at + 6iX -i itt + (a where Si = 0,- — 1. ' Combine equation (4.21) with the expression for the error term, e,( = A,(X)e,t. Since Ai(L) e it) 7 4 is invertible, and assuming a finite dimensional specification for produces the following augmented Dickey-Fuller equation, Pi (4.22) Ax it = a,- + SiX -i i)t + ^ 6 Axi^-L iL L=I + v it • T h e Kalman filter is a tool used to estimate the state vector in an optimal'way. Kalman smoothing is an inference about the value of the state vector based on the full set of data collected, e.g. estimate for 1980 using data from a full sample for the period 1961 to 1990. T h e E M algorithm is a method for maximizing a likelihood function in the presence of missing observations. It consists of two steps, an estimation and a maximization step, which are iterated to convergence. The maximization step calculates the maximum likelihood estimates of all unknown parameters conditional on a full data set. The estimation step constructs estimates of the sufficient statistics of the problem conditional on the observed data and the parameters. Missing observations are estimated based on the parameter values at one step of the iteration and then the likelihood function is maximized assuming that this is the full observable data set in the other (see, for example, Watson and Engle [81]). Approximate e, by, for eachi = 1,..., N, e,- = p.ie.^-i +P(2£i,t-2 + • • •+Pip (-i,t-p +vu fort = 1,2,...,T. 72 73 74 t t 67 i i i where the lag length p,- is allowed to vary across individuals. The null hypothesis of a unit root Ho : Si = 0 for all i, is tested against the stationary alternative fli : Si < 0 for all i. Where the stochastic process generating yu contains a linear trend, equation ( 4 . 2 2 ) is modified to, pt Ax it = a,- + jit + <5,a;,-,t_i + O Axi iL tt + v L it L=l where the vu are assumed to be distributed independently across both i and t. IPS [42] base the t-hax test on the group mean of the t-statistics for testing Si = 0 , 1 ( 4 . 2 3 ) t (p,p) N J2t (pi i) = — NT iT )P i—l where p,- is (unknown but) estimated as T / , p = (pi,j>2> • • -IPAT)', Pi = 1 3 (pn, Pi2, • • •, Pipi) and /> = (Pi, P2, • • •, /»AT). Assuming that N/T —* 0 as N —• oo 1 and T —> o o , the standardized t-bar statistic is given by, 75 N l [i - x 2 (4.24) z NT = NT a ] NT 1/2 b which is asymptotically distributed i V ( 0 , 1 ) . (^JVT) The mean (OTVT) a n d variance of the standard unit root distribution are given by, 1 N i=l and 1 N i=l respectively. IPS [42] present tabulated values of E\t-r{pi, 0)] and V[<T(P,,0)] for different values of T and p = 0 , 1 , . . . , 12, for cases when the underlying A D F regressions for each group of individuals are estimated both with and without a linear t r e n d . 76 For T large enough, these tabulated values are valid for normal and non-normal errors. These assumptions are necessary because in practice Z J V T depends on the parameters p,and is therefore not an operational statistic. The mean (E[tx{pi, 0)]) and variance (V[tx(pi, 0)]) are computed using stochastic simulations with 50,000 replications. For the model with no time trend, 4T(P>0) is the t-statistic 7S 76 68 4.4 R e s u l t s In this Section, I present the results from applying each of the three panel data unit root tests to Canadian data, and comparing these with the outcome when using a standard Dickey-Fuller test which is applied to data averaged across the cross-section. I use two different sets of data. First, I apply the tests to annual C A N S I M data on gross provincial product ( G P P ) for the ten Canadian provinces for the period 1961-1990. These data are deflated by a national price i n d e x . 77 Then I use monthly data on G D P per employee for 19 different manufacturing industries across Canada over the period 1974:2 to 1994:11. \: 78 4.4.1 Provincial data Panel Dickey-Fuller test After adjusting the data by subtracting the initial value from each observation, I estimate equation (4.6) using O L S . The resulting equation is given below, x = it 0.08273 + 0.9743z,,,_i (0.01118) (0.00985) where standard errors are shown in parentheses. The f-statistic for testing S = 1 is —2.61. The 95% critical value for 298 degrees of freedom is 1.645. Hence, using the panel Dickey-Fuller test, I can infer that the provincial data do not contain a unit root. of 6 in the ADF regression given by, p Ax = a + Sx -i t t + ^ ] Sj, Ax _£, + error. t ' J i; L=l Where a time trend is included, the following equation is used, P Axt = a + Sx —i +-yt + ^ ^flr,Ax _j, + error. t t L=l The underlying data generating process is Aa?t = £ti where et ~ N{0,1), t = 1,..., T. See Coulombe and Lee [21]. I use data for the following industries: food, beverage, tobacco, rubber, leather and allied products, primary textiles and textile products, clothing, wood, furniture and fixtures, paper, primary metals, fabricated metals, machinery, transport equipment, electrical and electronics, non-metallic mineral products, refined petroleum, chemicals and other. 77 78 69 8 a RSE(8) ts n t Model 2 -0.152120 0.990481 0.033338 -4.562972 -4.177988 Model 3 -0.226170 0.971396 0.043403 -5.210919 3.956520 Table 4.1: /5-Bar Test Statistics for Provincial D a t a Including a time trend does not change the result. I estimate equation (4.8) and find, x it = 0.07540 + 0.01185* + (0.01141) (0.00446) 0.94525z<;Li. (0.01465) \ In this case, the t-statistic for the hypothesis 6=1 equals —3.74. Once again, this value exceeds the 95% critical value of 1.645. The data do not support the null hypothesis of a unit root. <S-bar test: base case I assume initially that the data are generated by Model 2, i.e. the data may have an individual specific mean, but do not contain a time trend. I subtract the national average for each time period from each observation. Then I regress Axu and xu-i on the first two lags of the first, differences and a constant. This produces the residuals in and vu-i- Using the weights and lag truncation parameters recommended by L L [51], I derive the average standard deviation ratio (5). Finally, estimating equation (4.17) and computing t* as i n equation 6 (4.18) produces the test statistics presented in Table 4.1. Since t* ~ iV(0,1), I use the critical values of the standard normal distribus tion to test the null hypothesis that 6,- = 0 for all i = 1 , . . . , N. The absolute value of the 6-bar test statistic is 4.18. This clearly exceeds the 95% critical value of 1.96; the data reject the null hypothesis of a unit root. Under Model 3, the series {yu} may have an individual specific mean and a time trend. The results from estimation of this model are also given i n Table 4.1. The 6-bar test statistic of 3.96 clearly exceeds the 95% critical value of 1.96. The data continue to reject the null hypothesis of a unit root. 70 Model 2(1) 6 RSE{6) U n Model 2(2) -0.5435792 1.016973 0.0670925 -8.101936 -96.18277 -0.4985409 1.019640 0.0653325 -7.630826 -92.16263 Model 3(1) -0.550403 1.022312 0.058905 -9.343975 -96.76558 Model 3(2) -0.594886 1.020333 0.060517 -9.830071 -100.8222 Table 4.2: 6-Bai Test Statistics (with Index Adjustment) for Provincial Data (5-bar test: index adjustment I now repeat the 6-bar test, adjusting the data by subtracting an index at the first step, rather than the cross-section average. Using the QS [76] dynamic index model, I assume that M a = 1, M g = 1 arid M j = 2 for both a one- and a two-index model. The one-index model converges after 28 iterations; the two-index converges after 14 iterations. 79 The results for Model 2 (constant but no time trend) and Model 3 (constant plus time trend) are presented in Table 4 . 2 . 80 Adjusting the data by subtracting either a single index or the sum of two indices i n the two-index representation, the data continue to reject the null hypothesis that Si = 0 for both models (the - • - i l l 95% critical value is 1.96). The cS-bar test statistics vary little with the number j of indeces or the model used. t-bar test I adjust the data by subtracting the initial value from each observation, and estimate equation (4.22) (constant but no time trend) for each province. For most provinces, pi = 0; for Quebec, pi = 1; and for P E I , Nova Scotia and Saskatchewan, p, = 4. The t-statistics for testing Si '= 0 are shown i n Table 4.3. Values for the mean (ajvr) and variance (fr/vr) of the standard J unit root distribution are calculated as —1.48 and 0.83 respectively, giving a value for the standardized t-bax statistic of 3.14. The 95% critical value for the normal Appendix F includes charts comparing the explanatory power of two observable measures—aggregate employment and national GNP—with one- and two-index representations for the provincial data. The one-index (one common factor) model provides little explanatory power over the observable measures, but the two-index (two common factors) model appears to provide a better description of the underlying common comovements. Model i(j) refers to the estimation of Model t using a j-index adjustment. 7 9 80 71 Province Newfoundland PEI Nova Scotia New Brunswick Quebec Ontario Manitoba Saskatchewan Alberta British Columbia t-Statistic: ^-Statistic: No Time Trend Time Trend -0.39453 -2.9795 0.63051 -2.6298 1.0428 -1.8176 -0.37746 'I' -2.8174 -3.4542 -0.76778 -0.45904 -2.1115 -1.1032 -1.2389 -1.7646 -1.7158 -1.3243 -0.34315 -1.0687 -2.0090 Table 4.3: t-Statistics (t-Bar Test) for Provincial Data distribution is 1.96. Therefore, the null hypothesis of a unit root continues to be rejected. Table 4.3 also includes the ^-statistics for testing Si = 0 when a time trend is included i n the initial regression. The values for p,- are the same as for Model 2. The mean and variance statistics are calculated as a j v i = —2.11 and &JVT = 0.72 respectively, and the standardized t-bar statistic as 0.05. In this case, and for the first time using panel data tests, I cannot reject Ithe null hypothesis of a unit root. The addition of a significant time trend increases the probability of rejection of the null hypothesis. However, for the Prairie provinces, i.e. Manitoba, Saskatchewan and Alberta, addition of a time trend reduces the absolute value of the individual t-statistics (see Table 4 . 3 ) . 81 I investigate the significance of the coefficient on the time trend in the equation for each province: it is only significant (at all levels) for Quebec. When including a time trend, therefore, t: there is little change in the individual ^-statistics, but there is a considerable increase in the mean correction factor (ajvr), which varies with the inclusion of a time trend. These two effects combine to make the numerator of the standardized t-h&t test statistic close to zero. I repeat the t-h&x test, this time using the results from the equation with a time trend for Quebec, and from equation T h e [-statistics for these provinces, in the equation with ho time trend, are the largest of the group, and the closest to suggesting rejection of the null hypothesis at the individual level. 8l ; 72 (4.22) (excluding a time trend) for all other provincial regressions. The correction factors are adjusted accordingly. In this case, OATT = —1-54, 6/vr = 0.82 and the standardized t-bar statistic is 2.45. The data reject the null hypothesis of a unit root. ' r 4.4.2 Industrial data Panel Dickey-Fuller test I adjust the manufacturing industry data (as the provincial data) by subtracting the initial value from each time series observation. The resulting equations, for models with and without a time trend, are shown below (standard errors i n parentheses), x = it 0.27289 (0.01613) + 0.73456zi _i )t (0.01004) and xu = 0.23388 (0.01663) + 0.78065* + 0.70976a; _i. M (0.09025) (0.01037) ; ! The associated ^-statistics for testing 6 = 1 are calculated as —26.4 and —28.0 respectively. These clearly exceed the 95% critical value (1.645); the manufacturing industry data also reject the null hypothesisiof a unit root. 6-bar test: base case As for the provincial data, I subtract the cross-section average from each observation. I regress A i « and xu-i on 5 lagged first differences. Then I derive S, estimate equation (4.17) and calculate t* . 6 The results for both Model 2 (constant, no time trend) and Model 3 (constant plus time trend) are presented in Table 4.4. Once again, the 8-hax test statistics (16.2 for Model 2 and 32.8 for Model 3) exceed the 95% critical value (1.96). The data continue to reject the null hypothesis of a unit root. <5-bar test: index adjustment Using the QS [76] dynamic index model, I assume, as for the provincial data, that M = 1, M = 1 and Mb = 2. Table 4.5 contains the results for Model a g 73 6 RSE(6) ts n Model 2 Model 3 -0.03396090 0.9844993 0.006495736 -5.228184 16.18465 -0.126373 0.985250 0.011095 -11.38994 32.81161 Table 4.4: 6-Bax Test Statistics for Manufacturing Industry Data 6 RSE(6) ts Model 2(1) Model 2(2) -0.057845 0.978387 -0.007402 -7.814841 -17.92392 -0.048936 0.978379 0.006733 -7.267964 15.74611 Model 3(1) -0.219748 0.984247 0.014215 -15.45890 46.14207 Model 3(2) -0.176447 0.984239 0.012500 -14.11604 38.75622 Table 4.5: £-bar Test Statistics (with Index Adjustment) for Manufacturing Industry Data 2 and Model 3, assuming both one and two common factors (indeces). The critical value is 1.96. Again, the data reject the null hypothesis of a unit root for the one- and two-index cases, in models with and without a time trend. i-bar test '.;.'[;. I subtract the initial cross-section value from all observations and estimate equation (4.22) (model without a time trend). The ^-statistics for testing Si = 0 are given in Table 4.6. Assuming no time trend, and with pt = 12 for all industries, CITVT and &JVT are calculated as —1.48 and 0.77 respectively. The standardized t-bar test statistic is 2.17, which again exceeds the i95% The data reject the null hypothesis. critical value of 1.96. '' When including a time trend, p,- = 12 for all industries, ONT = —2.13 and b?jT = 0.61. The standardized t-bax statistic is given by —0.37. A s for the provincial data, when including a time trend, I can no longer reject the null hypothesis of a unit root. Again, I investigate the significance bf the time trend i n each individual industry regression: the coefficient is only significantly different from zero for the beverages industry. I compute a mixed test statistic, 74 t-Statistic: No Time Trend -2.4746 -0.54431 -1.1016 -2.0967 -0.77112 -0.89989 -1.3708 -1.6962 -2.0668 -1.0234 -0.65800 0.04411 -2.9196 -0.93501 1.3756, -2.5187 0.55835 -1.1536 0.35021 Industry Food Beverage Tobacco Rubber Leather Textiles Clothing Wood Furniture Paper and allied products Primary metals Fabricated metals Machinery Transportation equipment Electricals and electronics Non-metallic mineral products Petroleum Chemicals Other i-Statistic: Time Trend -2.7768 -0.79868 -1.7530 -3.6837 -4.0158 -2.6638 -2.7053 -1.5991 -2.6099 -2.2891 -2.5277 -2.3026 -2.9639 -2.5770 2.4703 -2.8119 -0.62480 -4.0168 -1.4189 Table 4.6: t-Statistics (t-Bar Test) for Manufacturing Industry Data i 1, 1 75 l;| i ,il ' Province 0.51989 (3.78) 1.6216 (3.78) 1.7672 (3.78) 3.4829 (5.34) 1.5653 (5.34) 1.3844 (5.34) No time trend Index (1) + no time trend Index (2) + no time trend Time trend Index (1) + time trend Index (2) + time trend Manufacturing Industry 0.84999 (3.78) 2.6638 (3.78) 2.5636 (3.78) 1.1083 (5.34) 2.9278 (5.34) : 2.8274 (5.34) Table 4.7: D F Test Statistics for Provincial and Manufacturing Industry D a t a using the t-statistics from the equation with a time trend for beverages, and from the equation without a time trend for all other industries. This gives ajvT = —1.52, 6JVT = 0.76 and a standardized i-bar statistic of 2.27. In this instance, the data reject the null hypothesis. 4.4.3 D i c k e y - F u l l e r test For comparison, I compute a standard Dickey-Fuller test for a unit root, for both data sets. In order to be able to use the test, I average the data over the cross-section. I subtract (i) the simple cross-section average (ii) the; average of a single index or (iii) the average sum of two indeces, from each cross-section observation, and estimate models with and without a time trend, given by the following equations, ' p (4.25) Ay = a + 6y -! + ]T B Ay -L L=l t t L t +e | t j .. : , -'i and p (4.26) Ay = a + 6y -i + 7* + Yl L&Vt-L d t t + U- \ • For equation (4.25) (no time trend), I calculate the F-statistic for testing a = 6 = 0. For equation (4.26) (including a time trend), I calculate the\F-statistic for testing 8 — j — 0. Table 4.7 summarises the results. Critical values are shown i n parentheses. I cannot reject the null hypothesis of a unit root for either of the data series. This result holds for both models (with and without a time trend), whether the data are adjusted by subtracting the cross-section average, a single index or a sum of two indeces. 76 4.5 Conclusions In this Chapter, I present the results from a variety of panel data unit root tests on Canadian data. Each panel unit root test statistic, unlike the Dickey-Fuller test, has a limiting normal distribution. I examine annual data on G P P (19611990) and monthly data on manufacturing industry G D P per employee (1974:2 to 1994:11). To use the standard Dickey-Fuller unit root test, data must be averaged across the cross-section. Using this approach, I cannot reject the null hypothesis of a unit root i n either data set. This result holds whether the data are adjusted (for individual effects) using the cross-section average or indeces, and whether or not a time trend is included. Using a panel Dickey-Fuller test ( B M [16]), neither the provincial nor the manufacturing industry data support the hypothesis of a unit root. For the 6-bai test ( L L [51]), I adjust the data by using a simple cross-section average, and by allowing for more complex patterns of contemporaneous correlation. In both cases, I again reject the null hypothesis of a unit root. Computation of IPS's [42] i-bar statistic, which explicitly takes account of small-sample bias, gives broadly complementary results. For this test, the results are sensitive to the inclusion of a time trend. However, I show that if, when calculating the test statistic, I only include a time trend in those individual regressions where the coefficient on the time trend is significantly different from zero, then the data once again reject the null hypothesis. The standard approach to testing for unit roots in panel data, which averages over individual effects, produces misleading results. Furthermore, re- searchers who first-difference panel data—in the belief that it contains a unit root—needlessly discard useful information. ( These results can be related to those of Granger [36], who finds that aggregation of a multivariate time series process can lead to a class of model that has fundamentally different properties. In particular, even i f disaggregate series are stationary, their aggregation may contain a unit root. This result explains the apparent inconsistency between the need to first difference the aggregate output data in Chapter 3, and the finding i n Chapter 4 that, i f cross-sectional information is included, the data appear to be stationary. If, however, the cross- 77 sectional effects are averaged out, as in the standard Dickey-Fuller approach, the hypothesis that the data contain a unit root cannot be rejected. 78 Chapter 5 Conclusions Three essays comprise this thesis. These essays study the relationship between aggregates and disaggregates in the context of Canadian economic growth. M y results support the findings in several recent papers that have emphasised the importance of including disaggregate information to explain aggregate activity (see, for example, Durlauf [28], Galor and Zeira [35] and Quah [68]). In the first essay, Chapter 2, I examine the dynamics of economic growth in Canada, investigating the evolving nature of the underlying cross-section distributions. Previous studies of economic growth and convergence in Canada focus on aggregate measures: a convergence rate of 2% per annum is interpreted as evidence that less well-off regions will, in the long run, catch up to rich regions. However, the methods used i n these studies reveal little about fluctuations i n the pattern of economic growth within the cross-section. They tell us nothing about the mobility within the cross-section rankings; or about the shape of the cross-section distribution. Using an econometric framework developed by Quah [69] that incorporates cross-section mobility and shape dynamics, I explore relationships among Canadian disaggregates. I show that the cross-section distribution is unimodal for measures of provincial income per capita (relative to the national average). However, persistence exists in the cross-sectional ordering of relative income levels: poor provinces remain relatively poor; and rich provinces remain relatively rich. The long-run distribution of gross domestic product per employee across manufacturing sectors displays a bi-modal pattern. Industries tend to remain either at the top or bottom of the cross-section distribution. The long-run bimodal pattern also appears in Ontario manufacturing industry data, but not for aggregations of eastern or western provinces. Movements in Ontario appear to drive the national pattern, reflecting the relative: size of the province and the wide range of manufacturing industries represented there. 79 ! For policy makers concerned with provincial income inequalities, this essay provides a mixture of good and bad news. O n the one hand, for provincial income per capita, the pattern of economic growth supports the concept of long-run convergence. O n the other hand, there is no evidence to suggest that poorer provinces can approximate richer provinces in the short or medium run. Rankings remain unchanged within the cross-section distribution, for provincial and industrial disaggregations. A comparison of the results for the provincial income per capita data and the manufacturing G D P per employee data suggests that it could be government transfer payments that are behind the unimodal long-run distribution. In the absence of these payments, the evidence suggests that the ergodic distribution is diverging. Hence, as noted by Green [37] i n an early study of regional aspects of Canada's economic growth, relying solely on national growth to eliminate the gap between rich and poor provinces is unlikely to be successful. This raises the issue of whether transfer payments can be sustained indefinitely, particularly i n the face of increasing constraints on government expenditure. However, this lies beyond the scope of this thesis, but would be an interesting area in which to extend the current work. It appears that it is those.factors which determine a region's comparative advantage (such as geographic llocation, climate, resource endowments, etc.) which determine the steady-state distribution of income across provinces. Despite directing substantial expenditures towards stimulating particular types of industries in poorer regions, the results suggest that there is little evidence that these industries end up at the top of the cross-section distribution. The pattern of economic growth for industries which receive subsidies does not appear to be distinct from those which do not. Those manufacturing industries at the top of the distribution (refined petroleum and coal products,: tobacco, beverages and food) are evenly spread across provinces. The results illustrate the difficulties inherent i n targeting assistance towards potential industry "winners". In the second essay, Chapter 3,1 propose a technique to examine the impact of aggregate and disaggregate disturbances on measures of aggregate economic growth and interaction across disaggregates. I use data on provincial income per capita and manufacturing industry G D P per employee. The approach is motivated by studies such as Durlauf [28], that argue that the interaction across dis- 80 aggregates contains explanatory information for the aggregate, and vice versa. If the disaggregate disturbance can be shown to have some effect on aggregate economic activity, then this strengthens the case for including disaggregate i n formation i n an explanation of aggregate growth. ..„: I find that the aggregate disturbance has positive long-run effects on provincial income per capita and G D P per employee, and explains most of their variation. The disaggregate disturbance only matters for fluctuations in aggregate income at business cycle horizons. It has little effect on G D P per employee. However, the disaggregate disturbance has important short and medium-run effects on the interaction measures. In contrast, the-aggregate shock contributes little. '! I argue i n Chapter 2 that an explanation of economic growth and/or convergence should include disaggregate information. The results from Chapter 3 confirm that disaggregate disturbances contain important information for aggregate economic activity at business cycle horizons. Unfortunately, the traditional method of considering convergence concentrates on the explanatory power of aggregate shocks: i n the standard neoclassical growth model, long-run changes i n output are driven by technology shocks. To gain a mOf e complete picture of convergence the researcher should include disaggregate information. Policy-makers mispecify the expected impact of different macroeconomic policies by ignoring the role played by disaggregate disturbances i n the short to medium run. Models which are based on the actions of a representative agent, for example, will miss important information, since they ignore the interaction across disaggregates. In the third essay, Chapter 4, I present the results from a variety of panel data unit root tests applied to Canadian provincial and manufacturing industry income data. In order to use a standard Dickey-Fuller unit root test, the researcher must average data across the cross-section. Using this approach, I cannot reject the null hypothesis of a unit root for either data set. This result holds whether the data are adjusted for individual effects using the cross-section average or a more complex index of contemporaneous correlation, and whether or not a time trend is included i n the estimated equation. Using a panel Dickey-Fuller test of Breitung and Meyer [16], that removes fixed effects, neither provincial nor manufacturing industry data support the hypothesis of a unit root. For the £-bar test (Levin and L i n [51]), I adjust 81 the data by using both a simple cross-section average and a dynamic index representation. In both cases, I again reject the null hypothesis of a unit root. Calculation of Im, Pesaran and Shin's [42] i-bar statistic, which explicitly takes account of small-sample bias, also suggests that i n both cases the null hypothesis can be rejected. The presence of a unit root has important implications for macroeconometric studies. If data contain a unit root, shocks will have permanent effects; where there is no unit root, shocks will have only temporary effects (see, for example, Nelson and Plosser [57]). The results from Chapter 4 suggest that when taking into account cross-section information, there are no permanent shocks in Canadian provincial and manufacturing industry panel'income data. Policy makers could misjudge the impact of policy changes if, when estimating the likely effect of a technological disturbance for example, they assume that it will have permanent effects. The main finding of this thesis is that aggregating data across the crosssection, or ignoring disaggregate information, can produce misleading results. This is illustrated i n different ways in the three essays. In Chapter 2, I show that the traditional approach to measuring economic convergence hides the underlying shape and mobility dynamics. In Chapter 3, I demonstrate that disaggregate shocks affect economic growth at business cycle frequencies. However, these shocks are ignored if the cross-sectional dimension is averaged out. Finally, in Chapter 4, I show that aggregating dynamic equations can lead to a model with fundamentally different time-series properties. Policy makers should be alerted to the contribution of disaggregate information and its possible impact on the aggregate. B y ignoring such information, they could seriously mispecify their model, and reach inappropriate conclusions. .• 82 Bibliography [1] Andres, Javier and A n a Lamo (1995): "Dynamics of the Income Distribution Across O E C D Countries", C E P R Discussion Paper No. 252, July. [2] Banerjee, Anindya, Robin L . Lumsdaine and James H . Stock (1992): "Recursive and Sequential Tests of the Unit-Root and Trend-Break Hypotheses: Theory and International Evidence", Journal of Business and Economic Statistics, V o l . 10, No. 3. [3] Barro, Robert J . (1991): "Economic Growth in a Cross Section of Countries" , Quarterly Journal of Economics, Vol. 106, No. 2, 407-443. [4] Barro, Robert J . and Xavier Sala-i-Martin (1991): "Convergence Across States and Regions", Brookings Papers on Economic Activity, V o l . 1, 107182. [5] Barro, Robert J . and Xavier Sala-i-Martin (1992): "Convergence", Journal of Political Economy, V o l . 100, No. 2, 223-251. [6] Bayoumi, T a m i m and Barry Eichengreen (1993): "Shocking Aspects of European Monetary Unification" in F . Torres and F . Giavazzi, eds., Adjustment and Growth in the European Monetary Union, Cambridge, U K : Cambridge University Press. [7] Bayoumi, T a m i m and Barry Eichengreen (1993): "Monetary and E x change Rate Arrangements for N A F T A " , I M F Working Paper, N o . 20, March. •,' ' [8] Ben-David, D . (1996): "Trade and Convergence among Countries", Journal of International Economics, V o l . 40, 279-298. r [9] Bernard, Andrew B . and Charles I. Jones (1996): "Productivity Across Industries and Countries: Time Series Theory and Evidence", The Review of Economics and Statistics, 135-146. [10] Bernard, Andrew B . and Steven N . Durlauf (1996): "Interpreting tests of the convergence hypothesis", Journal of Econometrics, V o l . 71, 161-173. [11] Berthelemy, J . C . and A . Varoudakis (1996): "Economic Growth, Convergence Clubs, and the Role of Financial Development", Oxford Economic Papers, Vol. 48, 300-328. • 83 [12] Bianchi, Marco (1995): "Testing for Convergence: A Bootstrap Test for Multimodality", Working Paper, Bank of England, May. [13] Blanchard, Olivier Jean and Lawrence F . Katz (1992): "Regional Evolutions", Brookings Papers on Economic Activity, V o l . 1, 1-76. [14] Blanchard, Olivier Jean, and Danny Quah (1989): "The Dynamic Effects of Aggregate Demand and Supply Disturbances", American Economic Review, V o l . 79, 655-673. [15] Blanchard, O. and D . Quah (1993): "Fundamentalness and the Interpretation of Time Series Evidence: A Reply to Lippi and Reichlin", American Economic Review, Vol. 83, 653-658. [16] Breitung, J . and W . Meyer (1994): "Testing.for unit roots i n panel data: are wages on different bargaining levels cointegrated", Applied Economics, V o l . 26, 353-361. [17] Campbell, J . and P . Perron (1991): "Pitfalls and Opportunities: What macroeconomists should know about unit roots", i n NBER Macroeconomics Annual, M I T Press. . • .• i [18] Carlino, Gerald and Leonard Mills (1994): "Convergence and the US States: A Time Series Analysis", Federal Reserve Bank of Philadelphia, Working Paper 94-13, July. [19] Cass, David (1965): "Optimum Growth i n an* Aggregative Model of Capital Accumulation", Review of Economic Studies, V o l . 32, 233-40. [20] Cogley, T . and J . M . Nason (1995): "Output" Dynamics i n Real Business Cycle Models", American Economic Review,/.Vol. 85, No. 3, 492-511. [21] Coulombe, Serge and Frank C . Lee (1994): "Convergence Across Canadian Provinces, 1961 to 1991", Canadian Journal of Economics, V o l . X X V I I I , No. 4a, 886-898. [22] Coulombe, Serge and Frank C . Lee (1995): "Long-run regional growth patterns and economic convergence in Canada", First Draft of paper to be presented at the meetings of the Canadian Economic Association, U Q A M , Montreal, 2 June. [23] de la Fuente, Angel (1996): "On the Sources of Convergence: A Close Look at the Spanish Regions", Centre for Economic Policy Research, Discussion Paper no. 1543, December. : [24] den Haan, Wouter J . (1995): "Convergence in stochastic growth models: The importance of understanding why income levels differ", Journal of Monetary Economics, V o l . 35, 65-82. [25] Dickey, D . and W . A . Fuller (1979): "Distribution of the Estimates for Autoregressive Time Series with a Unit Root", Journal of the American Statistical Association, V o l . 74, 427-431. 84 Dickey, D . and W . A . Fuller (1981): "Likelihood Ratio Statistics for A u toregressive Time Series with a Unit Root", Econometrica, V o l . 49, 10571072. Diebold, F . X . and M . Nerlove (1990): "Unit Roots i n Economic Time Series: A Selective Survey", Advances in Econometrics, No. 8, 3-69. Durlauf, Steven N . (1993): "Nonergodic Economic Growth", Review of Economic Studies, V o l . 60, No. 2, 349-366. Durlauf, Steven, N . and Paul Johnson (1994): "Nonlinearities i n Intergenerational Income Mobility", Working Paper, University of Wisconsin, Economics Department, Madison, W I 53706. Enders, W . (1995): Applied Econometric Time Series, Wiley Series i n Probability and Mathematical Statistics, Wiley: New York. Engle, R . F . and J . V . Issler (1995): "Estimating Common Sectoral C y cles", Journal of Monetary Economics, Vol. 35, 83-113. Evans, P . and G . Karras (1996): "Convergence Revisited", Journal of Monetary Economics, V o l . 37, 249-265. Evans, P . and G . Karras (1996): "Do Economies Converge? Evidence from a Panel of US States", The Review of Economics and Statistics, 384-388. Gali, J . (1992): "How Well Does the I S / L M Model F i t Postwar US Data?", Quarterly Journal of Economics, Vol. C V I I , 709-735. Galor, Oded and Joseph Zeira (1993): "Income Distribution and Macroeconomics", Review of Economic Studies, V o l . 60, No. 1, 35-52. ; Granger, C . W . J (1980): "Long Memory Relationships and the Aggregation of Dynamic Models", Journal of Econometrics, V o l . 14, 227-238. Green, A l a n G . (1971): Regional Aspects of Canada's Economic University of Toronto Press, Toronto and Netherlands. ,,' Growth, Hamilton, James D . (1989): " A New Approach to the Economic Analysis of Nonstationary Time Series and the Business Cycle" Econometrica, V o l . 57, No. 2, 357-384. Helliwell, John F . (1992): "Trade and Technical Progress", National B u reau of Economic Research Working Paper No. 4226. Helliwell, John F . (1994): "Convergence and Migration among Provinces", Working Paper No. 94-05, Department of Economics, Dalhousie University, July. ...-.». Helliwell, John F . and A l a n Chung (1992): "Convergence and Growth Linkages between North and South", National Bureau of Economic Research Working Paper No. 3948. 85 [42] Im, K . S . , M . H . Pesaran and Y . Shin (1995): "Testing for Unit Roots in Heterogeneous Panels", D A E Working Paper No. 9526, University of Cambridge. [43] Islam, Nazrul (1995): "Growth Empirics: A Panel D a t a Approach", The Quarterly Journal of Economics, November, 1127-1170. [44] Koopmans, T . C . (1965): " O n the Concept of Optimal Economic Growth", in The Economic Approach to Development Planning, Pontifical Academy of Sciences, Amsterdam: North-Holland. [45] Krugman, Paul (1993): "Lessons of Massachusetts for E M U " , i n F . Torres and F . Giavazzi, eds., Adjustment and Growth in the European Monetary Union, Cambridge, U K : Cambridge University Press. [46] Lee, Frank C.(1996): "Convergence in Canada?", Canadian Journal of Economics, V o l . X X I X , Special Issue, A p r i l , S331-S336. [47] Lee, Frank C . and Serge Coulombe (1994): "Regional Productivity Convergence in Canada", mimeo, Department of Finance, Federal Government of Canada, September. [48] Lee, Kevin, M . Hashem Pesaran and Ron Smith (1995): "Growth and Convergence: A Multi-Country Empirical Analysis of the Solow Growth Model", Department of Applied Economics Working Paper No. 9531, U n i versity of Cambridge. ;f f [49] Lefebvre, Mario. (1994): "Les Provinces Canadiennes et al Convergence: Une Evaluation Empirique", Bank of Canada Working Paper No. 94-10, November. [50] Lefebvre, Mario and Stephen S. Poloz (1996): "The Commodity-Price Cycle and Regional Economic Performance i n Canada", Bank of Canada Working Paper No. 96-12, September. 1 [51] Levin, Andrew, and Chien-Fu Lin (1993): "Unit Root Tests i n Panel Data: New Results", Discussion Paper No. 93-56, University of California, San Diego, December. y' [52] Lippi, M . and L . Reichlin (1993): " A Note on Measuring the Dynamic Effects of Aggregate Demand and Supply Disturbances", American Economic Review, Vol. 83, 644-52. [53] Long, John B . and Charles I. Plosser (1983): "Real Business Cycles", Journal of Political Economy, V o l . 91, No. 1, 39-69. [54] MacDonald, R . (1996): "Panel Unit Root Tests and Real Exchange Rates", Economics Letters, V o l . 50, 7-11. . [55] Mankiw, N . G . , D . Romer and D . Weil (1992): " A Contribution to the Empirics of Economic Growth", Quarterly Journal of Economics, V o l . 107, 407-37. 86 • ' • [56] Melvin, James R . (1987): "Regional Inequalities i n Canada: Underlying Causes and Policy Implications", Canadian Public Policy, V o l . XIII, No. 3, 304-317. [57] Nelson, C . R . and C . I. Plosser (1982): "Trends and Random Walks i n Macroeconomic Time Series", Journal of Monetary Economics, V o l . 10, 139-162. [58] Otto, Glenn and Tony Wirjanto (1990): "Seasonal Unit-Root Tests on Canadian Macroeconomic Time Series", Economics Letters, V o l . 34, 117120. [59] Pagan, Adrian (1984): "Econometric Issues in the Analysis of Regressions with Generalized Regressors", International Economic Review, V o l . 25, No. 1, 221-248. |f [60] Perron, P. (1989): "The Great Crash, the O i l Price Shock and the Unit Root Hypothesis", Econometrica, Vol. 57, 1361-1401. [61] Pesaran, M . H . and R . Smith (1995): "Estimating Long-Run Relationships from Dynamic Heterogeneous Panels", forthcoming i n Journal of Econometrics. [62] Pesaran, M . H . , R . Smith and K . Im (1995): "Dynamic Linear Models for Heterogeneous Panels", D A E Working Paper No. 9503, University of Cambridge, forthcoming in Matyas, L . and P. Sevestre, eds. Econometrics of Panel Data: Handbook of Theory and Applications, 2nd Edition, Kluwer Academic Publishers. ' . [63] Phillips, P. and P. Perron (1988): "Testing for a Unit Root i n Time Series Regression", Biometrica, Vol. 75, 335-346. [64] Quah, Danny (1992a): "Empirical Cross-Section Dynamics in Economic Growth", L S E Working Paper, October. [65] Quah, Danny (1992b): "International Patterns of Growth: I. Persistence in Cross-Country Disparities", L S E Working Paper, November. [66] Quah, Danny (1992c): "International Patterris'of Growth: II. Persistence, Path Dependence, and Sustained Take-Off i n Growth Transition", L S E Working Paper, November. / [67] Quah, Danny (1993): "Galton's Fallacy arid Tests of the Convergence Hypothesis", The Scandinavian Journal of Economics, V o l . 95, No. 4, 427-444. [68] Quah, Danny (1994): "One Business Cycle and One Trend From (Many,) Many Disaggregates", European Economic Review, V o l . 38, 605-613. [69] Quah, Danny (1994a): "One Business Cycle arid One Trend From (Many,) Many Disaggregates",European Economic Review, V o l . 38, 605-613. 87 [70] Quah, Danny (1994c): "Ideas Determining Convergence Clubs", mimeo, L S E Economics Department, F M G , and C E P R , September. [71] Quah, Danny (1994d): "Exploiting Cross Section Variation for Unit Root Inference in Dynamic Data", Economics Letters, Vol.44, 9-19. [72] Quah, Danny T . (1995): "Aggregate and Regional Disaggregate Fluctuations", mimeo, L S E Economics Department, March. [73] Quah, Danny T . (1996): "Twin Peaks: Growth and Convergence i n M o d els of Distribution Dynamics", The Economic Journal, V o l . 106, No. 437, 1045-1055. [74] Quah, Danny T . (1996a): "Twin Peaks: Growth and Convergence in Models of Distribution Dynamics", The Economic Journal, V o l . 106, No. 437, 1045-1055. \t [75] Quah, Danny (1996b): "Empirics for Economic Growth and Convergence", European Economic Review, forthcoming. [76] Quah, Danny and Thomas J . Sargent (1993): " A Dynamic Index Model for Large Cross Sections", in Business Cycles, Indicators, and Forecasting, Eds. J . H . Stock and M . W . Watson, The University of Chicago Press, Chicago, 1993. [77] Sala-i-Martin, Xavier (1994): "Regional Cohesion: Evidence and Theories of Regional Growth and Convergence", CEPR'Discussion Paper No. 1075, November. [78] Shapiro, M . and M . Watson (1988): "Sources of Business Cycle Fluctuations", NBER Macroeconomic Annual, Cambridge, Mass.: M I T Press, 3, 111-148. ;y [79] Silverman, B . W . (1986): Density Estimation for Statistics and Data Analysis, Monographs on Statistics and Applied Probability 26, Chapman and Hall, London. [80] Solow, Robert M . (1956): " A Contribution to the Theory of Economic Growth", The Quarterly Journal of Economics, V o l . L X X , 65-94. [81] Watson, M . W . and R . F . Engle (1983): "Alternative Algorithms for the Estimation of Dynamic Factor, M I M I C and Varying Coefficient Regression Models", Journal of Econometrics, Vol. 23, No. 23, 385-400. { [82] White, K . J . (1978): " A General Computer Program for Econometric Methods - S H A Z A M " , Econometrica, 239-240. 1 1 [83] Zivot, Eric and Donald W . K . Andrews (1992): "Further Evidence on the Great Crash, the Oil-Price Shock, and the Unit-Root Hypothesis", Journal of Business & Economic Statistics, V o l . 10, No.3, 251-270. 88 Appendices Mi 89 A The Data A . l D a t a Sources Variable Periods ; Source (CANSIM Matrix No.) Provincial Data Population Total employment Labour force 1926-1993 1966-1993 1966-1993 6967, 6968-6977 6967, 6968-6977 6967, 6968-6977 G D P in constant dollars G D P in current dollars GDP GDP/person Wages, salaries & S L I (GDP-based) Wages, salaries k S L I (Sources & disp. pers. inc) Average hourly earnings (1986=100) Personal income Personal income/person Disposable income Disposable income/person 1984-1993 1984-1991 1961-1993 1961-1993 1961-1993 V904-7913 7360-7369 6967, 6968-6977 6967, 6968-6977 2611-2619, 6949 1961-1993 ,5080-5097, 6965 1981-1993 6967, 6968-6977 1926-1993 1926-1993 1926-1993 1926-1993 6967, 6967, 6967, 6967, CPI Implicit price index 1978-1993 1971-1993 6967, 6968-6977 6967, 6968-6977 90 6968-6977 6968-6977 6968-6977 6968-6977 Variable Periods Source (CANSIM Matrix No.) A l l employees, all sizes Employment index Employees 1983:1-1995:1 1961:1-1983:3 1961-1992 1970-1982 1981-1993 4285 1432 7951-7964 7505-7689 5380-5586, 6849-6899 Total salaries & wages Wages, salaries k, S L I (15 sectors) Average hourly earnings 1981-1993 1926-1992 5380-5586, 6849-6899 6655 1961:1-1983:3 1983:1-1995:1 1435 4296, 4298 G D P at factor cost (1986 prices) G D P at factor cost (1986 prices) Real G D P / p e r s o n at work Real G D P / p e r s o n hour Index real G D P / p e r s o n Index real G D P / p e r s o n / h r 1961-1994 4670 1961:1-1995:1 4674 1961-1991 1961-1991 1961-1991 1961-1991 7920 7921 7926 7927 Industry Data 91 Variable Periods Source (CANSIM Matrix No.) Provincial Industry Data Total employees (21 manufacturing ind.) 1970-1982 1981-1993 A l l employees, all sizes 1983:1-1995:1 Average hourly earnings 1961:1-1983:3 1983:1-1995:1 '7482-7501 5379, 5401, 5406, 5409, 5413, 5419, 5424, 5429, 5439, 5458, 5473, 5482, 5496, 5504, 5515, 5540, 5548, 5567, 6848, 6865, 6869, 6883 4299, 4313,4327, 4341, ;,K4355, 4369,4383, 4397, ' 4411, 4425 1445, 1455, 1470, 1480, 1495 4310, 4324, 4366, 4380, 4422,4436 Salaries and wages (21 manufacturing ind.) 1970-1982 1981-1993 7482-7501 5379, 5401, 5413, 5419, 5439, 5458, • 5496, 5504, 5548, 5567, 6869, 6883 G D P at factor cost (current prices) 1971-1984 1984-1993 7871-7880 7360-7369 92 1460, 1465, 1485, 1490, 4338, 4352, 4394, 4408, 5406, 5424, 5473, 5515, 6848, • '•' 5409, 5429, 5482, 5540, 6865, A . 2 M a n u f a c t u r i n g Industries Food Beverages Tobacco products Rubber products Plastic products Leather and allied products Primary textile and textile products Clothing Wood Furniture and fixtures Paper and allied products Printing, publishing and allied Primary metal Fabricated metal products Machinery Transportation equipment Electrical and electronic products Non-metallic mineral products Refined petroleum and coal Chemical and chemical products Other manufacturing 93 B Fractile Transition Probability Matrices Table B l : 1 . - 4. W S S L I per Employee by Province „:J.: Quantile (Number) 0.20 0.40 0.60 (52) 1.00 (52) 0.75 0.25 (52) 0.25 0.73 (52) 0.02 (52) 1-year transitions Quantile (Number) 0.25 0.50 0.75 (52) 1.00 (78) 0.86 0.14 (52) 0.21 0.71 (78) 0.05 0.80 1.00 0.02 0.90 0.08 0.08 0.92 1.00 0.08 0.95 Table B2: G D P per Employee, 21 Manufacturing Industries 1-month transitions Quantile (Number) 0.20 0.40 0.60 0.80 1.00 (1620) 0.95 0.05 (1620) 0.07 0.86 0.09 , (1620) 0.09 0.81 0.09 (1620) 0.10 0.82 0.08 (1762) 0.01 0-08 0.92 1-month transitions Quantile (Number) 0.25 0.50 0.75 1.00 (2025) 0.95 0.05 (2025) 0.05 0.86 0.09 (2025) 0.09 0.82 0.08 (2167) 0.08 0.92 94 C Further Estimations In this Appendix, I present the results from estimation of transition probability matrices for data on provincial income per capita growth, on wages, salaries and supplementary labour income (WSSLI) per capita, and on employment per capita. I also examine G D P per person at work and G D P per person hour for seven different industry sectors across Canada. C.l Provincial Income Growth per Capita First, I estimate transition probability matrices for provincial personal income growth per capita. Figure C I shows the annual growth i n provincial personal income per capita, relative to the national average, for 1927-1993. Saskatchewan has the most volatile growth rate; by comparison, most other provinces vary relatively l i t t l e . 82 Table C I : Transition Probability Matrix—Income Growth by 1-year transitions Upper endpoint (Number) -0.022 -0.004 . Q 008 0.026 0.11 0.07 0.19 (128) 0.23 0.29 0:29 (126) 0.11 0.17 0.24 (123) 0.10 0.28 0.28 0.14 0.21 , 0.26 0.24 (125) 0.42 0.12 (126) 0.17 0.11 0.204 0.199 Steady-state distribution 0.198 0.203 5-year transitions Upper endpoint (Number) -0.022 -0.004 0.008 0.026 0.14 (125) 0.25 0.13 0.18 0.22 (113) 0.19 0.28 0.20 0.21 (107) 0.10 0.28 '0.27 (112) 0.26 0.24 0.31 0.10 (121) 0.17 0.31 0.17 6.10 Steady-state distribution 0.181 0.216 0.222 0.209 ; Province 0.699 0.41 0.13 0.11 0.15 0.17 0.196 f 0.699 0.30 0.11 0.14 0.09 0.26 0.172 ® The variance for Saskatchewan is more than six times the provincial average.' Ontario and Quebec have the least variable income per capita growth over the sample period. 2 95 i Transition probability matrices for these data (Table C l ) show much more mobility throughout the cross-section distribution than was evident in the levels data—there is no discernible pattern of persistence. Each panel in Table C l illustrates that it is quite common for economies to reach any other state from their initial position. The minimum entries in the off-diagonals are 7% and 10% respectively in each panel; in some cases, the largest entries are off the diagonal. In other words, provinces do not consistently exhibit high or low growth rates (relative to the national average) in income per capita. The ergodic distribution displays convergence towards the middle of the distribution; it is unimodal and symmetric. Plotting the quantiles (Figure C2) illustrates that much of the volatility occurs i n the first half of the sample period, and is contained in the end quantiles. Most of the variance in the top quantile is attributable to fluctuations i n income growth i n Saskatchewan and Manitoba. In''both provinces, a relatively large proportion of provincial G D P is attributable to agriculture and related products. 83 84 C.2 Provincial WSSLI per Capita Annual data on wages, salaries and supplementary labour income (WSSLI) per capita are available for the period 1961-1991. Figure C3 reveals that, relative to the national average, W S S L I per capita has remained fairly steady over the sample period. :] • Table C2 shows transition probability matricesjfor these data. The degree of mobility is very similar to that found i n the personal income per capita data. Note, however, the stronger persistence in states 2 and 3. This links with the patterns displayed by the ergodic distributions—both reach a peak i n state 3 . There appears to be a general tendency for W S S L I per capita to move towards the middle of the cross-section distribution and stay there. Figure C4, which graphs the quantiles, shows some reduction in dispersion (or the "shape") of W S S L I per capita, relative to the national averages In particular, the bottom quantile is increasing faster than the other quantiles, while, at the same time, the top quantile is falling over the sample period. ,f 85 ; Estimation with a ten-year transition period did not change the results significantly. *In Ontario, less than 10% of GDP is generated by agriculture and related products; in Manitoba and Saskatchewan, thisfigurerises to more than 25%. , A six-year transition period gives a similar pattern. , 83 8 1 8 5 96 ;'" ' Table C 2 : Transition Probability M a t r i x — W S S L I per Capita, by Province 1-year transitions Upper endpoint (Number) -0.283 -0.116 0.082 0.198 0.425 (62) 0.95 0.05 (59) 0.92 0.08 (58) 0.07 0.88 0.05 (60) 0.02 0.08 0.90 (61) 0.03 0.97 Steady-state distribution 0.364 0.226 0.000 0.296 0.115 3-year transitions Upper endpoint (Number) -0.291 -0.117 0.085 0.201 0.425 (60) 0.85 0.15 (48) 0.04 0.83 0.13 (52) 0.12 . 6.77 0.12 (53) 0.19 0.75 0.06 (57) 0.11 0.89 Steady-state distribution 0.082 0.321 0.196 0.296 0.105 1 C.3 Provincial Employment per Capita Annual employment per capita, relative to the national average, for 1966-1992, is plotted i n Figure C 5 . Cross-section dispersion has changed little over the sample period: the range is largely constant. There is also limited movement within the cross-section distribution: most provinces*maintain a particular rank within the cross-section (Alberta and Ontario at the top, Newfoundland and P E I at the bottom). This is confirmed by results from estimation of transition probability matrices (Table C3). The probability of remaining i n the bottom (top) iquantile is at least 91% (96%) with a three-year transition p e r i o d . The ergodic distributions, where they exist, suggest convergence towards the bottom,half of the distribution. , 86 Ji Even with a six-year transition period, these probabilities remain at 90% and 95% respectively. 86 97 Table C 3 : Transition Probability Matrix—Employment per Capita, Relative to the National Average, by Province 1-year transitions Upper endpoint (Number) -0.098 -0.017 0.049 0.097 0.194 (51) 0.92 0.08 (49) 0.06 0.94 (50) 0.90 0.10 (49) 0.10 0.86 0.04 (51) 0.06 0.94 Steady-state distribution 0.562 0.438 0.000 0.000 0.000 3-year transitions Upper endpoint (Number) -0.100 -0.022 ;0.049 0.099 0.194 (46) 0.91 0.09 0.02 0.98 (42) (41) 0.85 0.15 (46) 0.20 0.80 0.04 0.96 (45) Steady-state distribution n/a n/a n/a n/a n/a C . 4 Industrial G D P per Employee 1 : i ' f ; Annual index data on G D P per person at work and G D P per hour worked are available for seven different industry sectors for 1961-1994. Figures C6 and C7 show that the cross-section distributions of each measure are similar: values for agriculture (No. 1) are highly variable; those for communications (No. 5) show a steady rise; but productivity appears to have fallen i n the service sector (No. 7). Transition probability matrices (Tables C4 and) C5) show similar degrees of mobility within the cross-section. Most persistence is at the ends of the distribution, particularly the top end; mobility is highest in the middle. Mobility is higher than for provincial disaggregation using both W S S L I and employment per capita data. Ergodic distributions are unimodal for both measures. 87 The sectors studied are agriculture, manufacturing, construction, transportation and storage, communications, wholesale and retail trade, and services.' 87 7 Table C4: Transition Probability M a t r i x — G D P per Person at Work, Relative to the Industry Average, by Industry Sector 1-year transitions Upper endpoint (Number) -0.104 -0.027 0.011 0.097 0.499 (44) 0.77 0.16 0.07 (46) 0.15 0.39 0.30 0.15 (44) 0.02 0.32 0.34 0.30 0.02 (45) 0.13 0.49 0.31 0.07 (45) 0.02 0.09 0.89 Steady-state distribution 0.214 0.221 0.222 0.178 0.165 3-year transitions Upper endpoint (Number) -0.101 -0.027 .0.013 0.101 0.499 (41) 0.61 0.20 0.10 0.10 (39) 0.13 0.33 0.23 0.28 0.03 (42) 0.05 0.29 0.33 0.26 0.07 (41) 0.24 0.07 0.27 0.37 0.05 (40) 0.05 0.20 0.75 Steady-state distribution 0.241 0.229 0.138 0.151 0.240 Table C 5 : Transition Probability M a t r i x — G D P per Person Hour at Work, Relative to the Industry Average, by Industry Sector 1-year transitions Upper endpoint (Number) -0.088 -0.028 0,011 0.009 0.422 (44) 0.68 0.02 0.25 ; 0.05 (46) 0.22 0.39 0.28 0.11 i (45) 0.04 0.27 0:38 0.31 (44) 0.02 0.09 0.50 0.07 6.32 (45) ' 0.09 6.91 Steady-state distribution 0.219 0.220 0.168 0.187 0.206 3-year transitions .y j Upper endpoint (Number) -0.087 -0.028 0.009 0.106 0.422 (40) 0.57 0.20 • .0.15 0.08 (40) 0.13 0.40 0.20 0.25 0.03 (41) 0.12 0.29 0-27 0.27 0.05 (42) 0.12 0.10 0.40 0.33 0.05 (40) 0.08 0.17 0.75 Steady-state distribution 0.188 0.226 0.240 0.232 0.114 s •i 1 99 Figure C1 Growth of Provincial Personal Income per Capita (Relative to National Average) 1927-1993 0.8 -, 0.6 1 100 Figure C2 Quantiles for Growth of Provincial Personal Income per Capita (1929-1992) 101 Figure C3 Log of Provincial WSSLI per Capita (Relative to National Average) 1961-1992 102 Figure C4 Quantiles for Log of Provincial WSSLI per Capita (1963-1991) Year quantO — - — - quant3 quantl — - - — quant4 103 quant2 quant5 Figure C5 Provincial Employment per Capita (Relative to National Average) 1966-1993 1.4 i 0.4 H, 0.2 4 ( O N C O O S O T - C M C O m <o >• cooooscnococ ) O) m O) 1 r-. r- N co h» CF) o ^ CM co co co co co Oi Oi O) Cfi IT) CO Oi (O t-~ 0 0 CO O ) CJ) O) CO C31 Year -•—Newfoundland — N e w Brunswick H Manitoba — • — P E I Quebec Saskatchewan British C o l u m b i a 104 — * — N o v a — • Scotia Ontario Alberta O 0) Q T " (V O) O ) Q] C l Figure C6 Index of GDP per Employee 1986=100 (Relative to National Average) 1961-1994 i — i — r —r i i u -i—I—I—I—i—I—I—I—I— CO iCO Ul N O) in CO CD (0 CO CD 05 cn CD CO CD CT) CD CO - i i i*. i N CO i cn i N CO i T- co CO < CO CO cn m CO cn r-co o> i—i—i—i CO OO i - CD i—i—i CO CD CD CD CD Year • agriculture Wr manufacturing *• communications services 105 —if —construction — wholesale and retail trade Figure C 7 Index of G D P per P e r s o n H o u r W o r k e d 1986=100 (Relative to National A v e r a g e ) 1961-1994 U r-1—i—i—i—i—i—i—r i o> CO in N (O CD CO cn CD cn CO cji co cn • i i—i—i—I—I—I—i—i—i—r~i— i— I O1 co co m r-» o> 00 00 hK Nco o> O) Oi Oi o> T- CO OS CO T o>1 oo cn co 1i1 1 cj) 0 ) 0 1 co cj) 1 0) Year —•—— agriculture -B—-manufacturing -communication services 106 — — — construction — wholesale and retail trade D Tests For Robustness Table D I : Transition Probability Matrix—Income by Province, Ontario '_ 1-year transitions Upper endpoint (Number) -0.543 -0.376 -0.258 -0.134 0.94 (118) 0.05 0.01 0.82 (118) 0.03 0.04 o.ii (115) 0.09 0.80 0.10 (116) 0.01 0.03 0.09 0.77 (114) 0.01 0.10 Steady-state distribution 0.114 0.232 0.186 0.223 3-year transitions Upper endpoint (Number) -0.545 -0.381 -0.258 -0.130 (117) 0.86 0.13 0.01 (117) 0.03 0.71 . 0:21 0.04 (107) 0.01 0.13 0.66 0.18 (108) 0.02 0.14 0.69 0.02 (105) 0.03 0.13 Steady-state distribution 0.081 0.156 0.228 0.266 il fill ' 107 Relative to 0.126 0.02 0.09 0.89 0.245 ' 6.126 0.02 : 0.02 0.16 0.82 • 0.270 Table D2: Transition Probability M a t r i x - -Income by Province (subtracting Ontario) 1-year transitions Upper endpoint (Number) -1.561 -1.203 -0.811 -0.519 0.891 (117) 0.14 0.85 i 0:01 (114) 0.12 0.02 0.78 0.08 0.72 (117) 0.10 0.15 0.03 (117) 0.01 0.17 0.74 0.08 (116) 0.03 0.09 0.89 Steady-state distribution 0.184 0.204 0.218 0.197 0.197 3-year transitions Upper endpoint (Number) -1.561 -1.192 ^0:808 -0.513 0.891 0.02 (116) 0.75 0.23 (105) 0.17 0.66 0.16 0.01 (112) 0.14 0.01 0.58 0.22 0.04 (114) 0.21 0.01 0.01 0.66 0.11 (107) 0.01 0.82 0.06 0.11 Steady-state distribution 0.211 0.160 0.220 0.209 0.201 Table D3: Transition Probability Matrix—Income by Province, Aggregating some Provinces _ 1-year transitions Upper endpoint (Number) -0.102 0.045 0.137 0.260 0.453 0.91 0.08 0:01 (80) i(78) 0.06 0.67 0.24 0.03 0.01 0.26 0.05 (77) 6.68 (79) 0.84 0.01 0.06 6.09 (76) 0.09 0.91 Steady-state distribution 0.212 0.241 0.239 0.183 0.126 3-year transitions Upper endpoint (Number) -0.112 0.043 0.137 0.255 0.453 0.92 0.07 (76) .0.01 , 0.04 (74) 0.03 0.61 6.32 (69) 0.03 0.30 0.13 6:54 0.04 (75) 6:i7 0.68 0:11 \.l 1 (78) 0.01 0.15 0.83 ISteady-state distribution 0.178 0.250 0.253 0.194 0.124 1 108 Table D4: Transition Probability Matrix—]GDP per Employee, Manufacturing Industries, Ontario 1-year transitions Upper endpoint '. (Number) -0.355 -0.120 0.015 0.204 0.970 (77) 0.81 0.16 0.01 0.03 0.62 (79) 0.18 0.18 0.01 0.01 (77) 0.18 0.55 0.25 0.03 (80) 0.01 0.01 0.23 0.63 0.13 0.02 (86) 0.03 0.08 0.86 Steady-state distribution 0.202 0.176 .0.181 0.187 0.255 3-year transitions Upper endpoint (Number) -0.350 -0.119 i 0.026 0.211 0.802 0.02 (66) 0.70 0.18 0.03 0.08 (66) 0.23 0.48 0.18 0.06 0.05 0.03 0.25 0.24 (71) 0.45 0.03 0.21 (67) 0.03 0.06 0.57 0.13 0.02 (66) 0.09 0.11 0.79 Steady-state distribution 0.184 0.179 0.186 0.197 0.253 Table D5: Transition Probability M a t r i x — G D P per Employee, Manufacturing Industries, Western Provinces 1-year transitions .; Upper endpoint (Number) -0.601 -0.341 -0,117 0.301 2.000 <0.02 0.02 0.77 0.17 0.03 (60) 0.24 (72) 0.57 017 0.01 0.01 0.04 0.21 0.60 0.13 0.01 (70) (75) 0.03 0.63 0.19 6.16 0.03 0.01 0.01 0.20 0.75 (75) Steady-state distribution 0.259 0.203 0.179 0.187 0,173 3-year transitions ,, ,! , iL Upper endpoint (Number) -0.572 -0.308 -0.082 0.328 1.616 (60) 0.78 0.07 0.05 0.05 0.05 0.16 0.51 0:33 (51) 1 (62) 0.37 0.50 0.08 0.05 (64) 0.02 0.03 0.11 0.56 0.28 0.02 0.02 0.05 0.31 0.60 (58) Steady-state distribution 0.202 0.206 0.205 0.196 0.190 109 It. Table D6: Transition Probability M a t r i x — G D P per Employee, Manufacturing Industries, Eastern Provinces 1-year transitions Upper endpoint (Number) -0.504 -0.203 -0.008 0.208 1.476 (70) 0.87 0.13 0.12 0.62 (69) 0.23 0.01 0.01 (71) 0.23 0.49 0.23 0.06 (71) 0.23 0.59 0.18 0.12 0.01 0.07 0.80 (89) Steady-state distribution 0.169 0.204 0.194 0.246 0.187 3-year transitions Upper endpoint (Number) -0.504 -0.198 -0.001 0.217 1.476 (62) 0.77 0.19 0.03 (61) 0.11 0.49 0.26 0.08 0.05 (61) 0.21 . 0.49 0.21 0.08 0.03 0.23 0.45 0.29 (62) (64) 0.03 0-03 0.19 0.75 Steady-state distribution 0.075 0.148 0.198 0.225 0.355 ii • t' i' • 110 • • i '. 1 Figure D1 Quantiles for Log of Provincial Personal Income per Capita, Divided by Ontario Data (1928-1993) 0.2 1 .1.2-1 Year quanto quantl quant2 — - — - quant3 — - - — quant4 quant5 111 Figure D2 Quantiles for Log of Provincial Personal Income per Capita, Subtracting Ontario Data (1928-1993) 1 i -2.5 J Year quantO quantl quant2 — - — - quant3 — - - — quant4 — — — quant5 112 Figure D3 Quantiles for Log of Aggregated Provincial Personal Income per Capita (1928-1991) 113 Figure D4 Quantiles for Log of Manufacturing Industry GDP per Employee: Ontario (1973-1990) 0.4 + .1.2 J Year quanto quantl quant3 quant4 114 quant2 - quant5 Figure D5 Quantiles for Log of Manufacturing Industry GDP per Employee: Western Provinces (1973-1990) 115 Figure D6 Quantiles for Log of Manufacturing Industry GDP per Employee: Eastern Provinces (1973-1990) 116 E R e s u l t s U s i n g C o u l o m b e a n d Lee D a t a This Appendix summarizes the results arising from studying the mobility dynamics of data used by Coulombe and Lee in their various studies of convergence across Canadian provinces (see [21], [22] and [47])'. E.l Data In their first paper, Coulombe and Lee [21] draw a distinction between provincial product and income data for the period 1961-1990. They use the measures shown in Table E l . In their second paper, [47], Coulombe and Lee examine convergence in productivity, using the variables listed i n Table E2 below for the period 1966-1990. The most recent paper, [22], uses four measures of per capita income over the longer period 1926-1991 (see Table E3). From these papers, Coulombe and Lee conclude that there is evidence of convergence in both product and income measures, with the speed of convergence increasing as disaggregation moves from a per capita, to a per worker, to a per hour basis. Most of the convergence occurred during the period 1950-1977. Regional shocks in the 1930s and 1940s outweighed any evidence for /3-convergence prior to 1950. Nri, • Regional disparities in unemployment rates are the most significant factors in slowing down regional labour market adjustments, which, in turn, hinder convergence i n living standards. Provincial price deflators diminish the speed ; Table E l : Data Used i n Coulombeand Lee [21] Definition Gross provincial product at factor cost per capita, provincial price indeces GPPNAT Gross provincial product at factor cost per capita, national price index EI Earned income per capita i PIT Personal income minus government transfers per capita PI Personal income per capita i r- . t PDI Personal disposable income per capita Variable GPPRO ; > 1. 117 Variable GPPC GPPW GPPH EIC EIW EIH Variable PI PIT PDI EI Table E 2 : Data Used i n Coulombe and Lee [47] Definition Gross provincial product at factor cost per capita Gross provincial product at factor cost per worker Gross provincial product at factor cost per hour worked Earned income per capita Earned income per worker Earned income per hour worked Table E3: Data Used i n Coulombe and Lee [22] Definition Personal income per capita 'i';'';. Personal income minus government transfer payments per capita Personal disposable income per capita Earned income per capita of convergence of regional output but do not appear to affect the rate of convergence of income measures. This difference is attributed to the uneven distribution of economic activity across Canada. Consumption patterns appear to be similar. E.2 M o b i l i t y D y n a m i c s I estimate transition probability matrices for each o f the variables listed above, and present the results for the different variables i n turn. E.2.1 O u t p u t a n d income p e r c a p i t a The first variable I study is gross provincial product ( G P P ) per capita, deflated by (i) a national ( G P P N A T ) and (ii) provincial price indices ( G P P R O ) (see Coulombe and Lee [21]). . Alberta has the highest average G P P per capita (relative to the national average) over the sample period, followed by Ontario, British Columbia and Saskatchewan. The Maritime provinces are at the bottom of the ranking. Provincial price deflators increase cross-section mobility and tend to smooth the output d a t a . There is a slight narrowing of the range between the ends of the cross-section distribution over the sample period. This range is smaller for the nationally deflated data. This supports the Coulombe and Lee [21] finding that provincial price deflators reduce the speed of convergence. Estimation of transition probability matrices for each output measure suggest convergence, both i n terms of a reduction i n cross section dispersion (shape) 1 1 88 r 8 8 This suggests noise in the provincial price data. 118 ' I and within the cross-section distribution (mobility). It appears that convergence is towards state 2. Most of the income measures exhibit high mobility i n state 2 and high persistence at the top of the cross-section distribution. The top is occupied mainly by Ontario; the volatile agricultural economy of Saskatchewan influences mobility further down the distribution. Most quantiles dip i n 1980, except the fifth quantile—which includes Alberta data—which increases in the same period. This reflects the recession that hit Ontario particularly hard and reveals an interesting distinction between the output and income data. A l b e r t a had the highest output per capita (relative to the national average) over the sample period, but for the income data, it drops to third place, behind Ontario and British Columbia. The boost to provincial output from the oil boom was not immediately transmitted into higher income. Results for per capita earned income (EI) andiiper capita personal income minus government transfers (PIT) are very similar^ Provincial rankings are identical. The quantiles exhibit a reduction in dispersion over time. Therefore, the addition of interest and dividend payments does not appear to affect the pattern of convergence—a result consistent with Coulombe and Lee [21]. The ergodic distributions peak in state 4. The addition of government transfer payments (PI) pushes the ergodic peak up the cross-section distribution. In comparison with P I T , values for the top two quantiles are lower, and the bottom three quantiles higher: transfer payments help to reduce the disparity in personal incomes., The Maritime provinces, led by Newfoundland, received the highest value per capita of government transfer payments to persons over the sample period. PI shows more mobility at the bottom end of the cross-section distribution than the other income measures. However, contrary to expectations, this increased mobility is down the distribution rather than up it. Personal disposable income per'capita (PDI) data—which removes direct taxes—covers a narrower range, but exhibits similar patterns of mobility. '• ' 89 ' I 1 E.2.2 P r o d u c t i v i t y measures I now examine the different measures of provincial "productivity—output and income per capita, per worker, and per hour worked—used by Coulombe and Lee in [47]. Comparing first the output measures, gross provincial product per worker ( G P P W ) displays the most mobility, particularly at the top end of the crosssection distribution. Both G P P per capita ( G P P C ) and G P P W converge towards the middle of the cross-section distribution in the long term.' Ranking the provinces by the average value of G P P per hour ( G P P H ) (relative to the national average) over the sample period, places the two western-most provinces at the top, Newfoundland i n the middle , arid the remaining Maritime provinces at the bottom. This suggests that, although there is high unemployment in Newfoundland, those industries that'do employ workers are highly 1 Newfoundland moves up to 8th position—from 9th—using this measure. 89 119 "productive", i n terms of output per hour worked. G P P H exhibits more mobility at the bottom of the cross-section distribution, and is the only variable to display divergence across the cross-section. This divergence could be caused by the size of the distance between quantiles. The end quantiles are further away from the others. Alternatively, it could be attributed to the type and range of industries represented i n the different provinces i' (see Chapter 2 for a fuller discussion). .,' . Rankings of provinces vary across the different income variables: Ontario, A l berta and British Columbia share the top for all three measures, but Saskatchewan appears near the bottom for earned income per worker ( E I W ) and earned income per hour (EIH), whilst Newfoundland moves up the distribution to 7th place for E I W and to 4th place using E I H data. This reflects the impact of unemployment and the type of industries i n each province on the results. The Saskatchewan result can be attributed to the agricultural base of the province. This economy is heavily dependent upon agricultural prices and the success of harvests. A s a consequence, incomes tend to be highly variable and largely seasonal. This volatility is reflected i n the transition probability matrices. A l l the income data exhibit high mobility i n state 2. Probabilities and persistence are high at the top end of the cross-section distribution. Dispersion has declined in the earned income per capita (EIC) and E I H data, but less so i n the E I W series. Further investigation is required to deterrnine whether this reduction arises from a tendency for different industries to locate more widely (rather than concentrating production an particular provinces), or if it reflects a convergence i n productivity across' industries. The E I W results could be attributed (as i n Coulombe and Lee [2.1]), to consistently disparate unemployment rates across provinces. Ergodic distributions for the E I C and E I W series are similar, peaking i n states 3 and 4 respectively. 90 E.2.3 L o n g e r sample p e r i o d „ Personal income per capita (PI) and personal disposable income per capita (PDI) produce very similar transition probability matrices and ergodic distributions. In the long term, observations within the cross-section distribution converge to the middle of the distribution. E.3 Conclusions Mobility dynamics, for both output and income variables—per capita and per worker—suggest convergence both across (shape) and within (mobility) the cross-section distributions. There is high persistence at the top of the distribution, particularly for the income measures of economic growth. In particular, Ontario remains consistently at the peak of the distribution. High mobility I n 1994, Newfoundland's unemployment rate was 20.4%;' while the Canadian average was 10.4%. The food and printing and publishing industries areithe most important manufacturing industries in Newfoundland (in terms of contribution to GDP). GDP per employee in food manufacturing exceeds the average across all manufacturing industries. 90 120 is evident lower down the cross-section distribution, and is attributable to the uncertain, volatile incomes those provinces generate. The different output measures of productivity used i n [47] suggest different long run patterns of convergence. Both G P P C and G P P W exhibit convergence, but G P P H indicates divergence. The income data all suggest convergence i n the long term. These results indicate that there does not appear to be a direct one-toone relationship between output and income measures. Measures of income have more of an opportunity to converge. Government transfers to persons and migration tend to be the dominant sources that reduce dispersion. Measures of output are more heterogeneous, and possess relatively wide variations across provinces. This reflects the diverse natural resources, climates and industries of the different provinces. These findings mainly support Coulombe and Lee's ([21] and [47]) results, but provide additional insight into the driving forces behind economic growth. Where my results differ from Coulombe and Lee is i n finding divergence i n the cross-section distribution for G P P per hour worked. Coulombe and Lee [47] claim that the speed of convergence increases as disaggregation moves from per capita to per worker to per hour. Although my results show that dispersion may be falling over time, activity within the cross-section distribution may be diverging. G P P per hour tends to accumulate at opposite ends of the distribution. : 121 F Quah and Sargent's Dynamic Index Model The dynamic index model developed by Quah and Sargent (QS) [76] uses information from a wide range of cross-section data to explain aggregate fluctuations. F . l Theoretical Model Let {yu, i = 1 , 2 , . . . , JV, t = 1,2,...,T} be an observed segment of a random field where the cross-section (JV) and time series (T) dimensions are of the same order of magnitude. For example, j/,- may be income per capita (relative to the national average) i n province i. Suppose each j/,- is affected by J common factors u = ( « i , «2) • • • > « j ) ' and an idiosyncratic disturbance Z{. This can be written as, J ( G . l ) j/,t = Y <Xij(L)ujt + j=l z. it where u is a vector of orthogonal random walks; and z,- is zero mean, stationary, and has entries uncorrelated across all i as well as with the increments in u at all leads and lags. Observed yu are therefore decomposed into first, factors common across all i and, second, factors specific to each i and orthogonal to all other i. QS [76] assume the following structure for the impact of the common factors on j/,-, (G.2)a (X) = Y aij(rn)L , m 0 M a < oo. The Uj's are integrated processes, with their first differences pairwise orthogonal. They have a finite autoregressive represention, (G.3)7(L)Au, = rtut, • 'it''* ' where j(L) is diagonal, with j - t h entry given by, (G.4) 1 - (l)L gj - (2)L 2 9j - ... - (M )L ; M Bj a M t : g < oo. T) is white noise, with mean 0 and the identity covariance matrix. QS [76] further assume that each Zi is a finite order autogression, u 122 (G.5) A-(£)*,•« = e for alii it with (G.6)#(£) = 1 - - bi(2)L - . . . - & , ( M ) L > , 2 M 6 M < oo. 6 Combining equations ( F l ) and (F5) and denning 4>ij(L) = j3i(L)ctij(L), produces, j (G.7) f3i(L)y = <t>ij{L)u + e . it jt it i=i Exogenous variables, such as a constant or time trend—represented by — can be added (without loss of generality) to (F7). The model given by (F7) is translated into stateispace form, represented by a measurement equation (F8) and a transition equation ( F 9 ) , 91 (G.8) y = ax + dw + z , t t t t and ( G . 9 ) x i = cx t + t +j +it The state vector x has dimension O(N), which is 'computationally intractable' for large 7Y and large T . The state-space representation given by equations (F8) and (F9) describes N time series in terms of an 0(iV)-dimensional vector process. However a large part of x—the idiosyncratic disturbances—is directly observable. The unobservable part—the common factors—has dimension 0(J), which is independent of 7Y. QS [76] exploit this structure to give Kalman smoother projections and moment matrix calculations that are (effectively) i n dependent of N. : 1 The intuition is that increasing the cross-section dimension JV can only help to estimate more precisely the state and its cross moments. This must imply the same property for estimates of the parameters. Deleting entries of y does not change the orthogonality properties on a reduced version of (F8) and (F9). Therefore—if the model is correctly specified—estimators that exploit the orthogonality conditions i n (F8) and (F9) remain consistent independent of N. Smaller TV-systems imply less efficient estimators:' - . When the unknown distribution of (y, u, z, iv) generates a conditional expectation E(u\y, z) that equals the linear projection of u on y and z, then standard Kalman smoother calculations yield the conditional expectations E(x x l\y, w) and E(x \y, w), taking as given a particular setting for the unknown parameters. Iteration on this scheme is the E M a l g o r i t h m . Each iteration of E M requires a Kalman filter and smoother followed by straightforward regression calculations. The algorithm consists of two steps which are iterated to convergence: t t t 92 State space representation is a means of summarizing finite processes. See QS [76] for details of how this is done. ,',, , T h e E M algorithm is a method for maximizing a likelihood function when there are missing observations. 9 1 92 i• u , - 123 <; : I' ''. 1 an estimation and a maximization step. The maximization step calculates the maximum likelihood estimates of all unknown parameters conditional on a full data set. The estimation step constructs estimates of the sufficient statistics of the problem, conditional on the observed data and the parameters. Missing observations are estimated based on the parameter values at one step of the iteration and then the likelihood function is maximized assuming that this is the full observable data set in the other. Under weak regularity conditions the E M algorithm is guaranteed to converge to a point that solves first order conditions. 93 F.2 A p p l i c a t i o n I investigate the extent to which two observable measures, aggregate employment and national G N P , can capture the cross-correlation across provincial G D P data (1967-1990), and compare this with one- and two-index representations. Figures F l and F2 plot the sample standard deviations from second-order autoregressions over the period 1968-1988 (including a constant and time-trend). Figure F l graphs the residual sample standard deviation when employment at lags -2 through 2 are included as additional regressors, versus that without. If points fall below the 45 degree line, comovements are better described by the model represented on the vertical axis than that on the horizontal. Total employment does appear to contain information on provincial G D P fluctuations. A similar result arises with national G N P data (Figure F2). Figure F3 compares these two observable measures—both appear to give similar descriptions of underlying comovements i n the data. A one-index model provides little additional explanatory power over the observable measures, Figures F4 and F5 plot standard deviations bf innovations in the disturbances under a two-index model, against the residuals in province-by-province projections including total employment (Figure F4) and total G N P (Figure F5) respectively. The vertical axis describes the innovations upon removing two common unobservable indexes and imposing extensive orthogonality conditions; the horizontal axis describes the innovations upon removing total employment or total G N P respectively, without requiring the resulting innovations to be orthogonal across provinces. The results suggest that the two-index model provides a better description of the underlying common factors i n provincial G D P than does either aggregate G N P or aggregate employment. The two-index model provides more information about the cross-correlation in the data, but the combined effect of the two indices in the two-index model is very similar to the single-index m o d e l . '. ; 94 Watson and Engle [81] present the methodology for estimation using the E M algorithm in more detail. The estimation step gives parameter values,at step k + 1 based on moment matrices estimated in step k. Each of the moment matrices can be constructed from the output of one pass through the Kalman smoother. " The correlation coefficient between the single index and the summation of the two indices in the two-index model is 0.996508. -,v b; 93 ! 9 4 j Figure F1 Sample Standard Deviations: including Employment (1968-1988) Standard-Deviation 0.1 0.2 0.3 0.4 lnnov'-AR(2) 125 0.5 Figure F2 Sample Standard Deviations: Including G N P (1968-1988) Standard-Deviation 0.1 0.2 0.3 0.4 lnnov-AR(2) 126 0.5 Figure F3 Sample Standard Deviations: Including Employment vs. Including GNP (1968-1988) 127 Figure F4 Sample Standard Deviations: Two-Index Model vs. Total Employment (1968-1988) Standard-Deviation 0.1 0.2 0.3 0.4 lnnov-AR(2)-lncl-Total-Empl 128 0.5 Figure F5 S a m p l e S t a n d a r d Deviations: Two-Index M o d e l v s . Total G N P (1968-1988) Standard-Deviation 0.1 0.2 0.3 0.4 lnnov-AR(2)'-lncl-GNP 129 0.5
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Disaggregate dynamics and economic growth in Canada
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Disaggregate dynamics and economic growth in Canada Wakerly, Elizabeth Clare 1997
pdf
Page Metadata
Item Metadata
Title | Disaggregate dynamics and economic growth in Canada |
Creator |
Wakerly, Elizabeth Clare |
Date Issued | 1997 |
Description | This thesis takes the form of three essays in which I use disaggregate and aggregate information to examine Canadian economic growth. In the first essay, I present evidence that the process of economic growth differs for low income per capita provinces and industries. This contrasts with results from traditional studies of economic convergence. In those papers, estimates of a rate of convergence suggest that poor provinces eventually "catch up" to rich provinces by growing faster. Unfortunately, this approach ignores the pattern of economic growth within the cross-section distribution. Explicitly modelling the evolving distribution, I find little mobility in the cross-sectional ordering and some evidence of divergence. In the long run, the poor stay (relatively) poor and the rich remain (relatively) rich. In the second essay, I examine the dynamic effects of aggregate and disaggregate disturbances on both economic growth and the interaction between disaggregates. The approach is motivated by the class of models which predict two-way interaction between aggregate and disaggregate behaviour, such as Durlauf [28]. The disaggregate disturbance is identified as having no long-run impact on aggregate economic growth. I find that the aggregate shock has a large impact on aggregate income in both the short and long run; and accounts for most of its variation. The disaggregate shock contains some information for aggregate activity at business cycle horizons. Most interaction is explained by the disaggregate disturbance; the aggregate shock contributes little. In the third essay, I present results from a variety of unit root tests on provincial and manufacturing industry panel income data. Standard Dickey- Fuller unit root tests applied to panels require averaging of data across the cross-section. More powerful tests allow pooling of cross-section and time-series information. Using these methods, I find that the null hypothesis of a unit root is rejected—strongly contrasting with results obtained using the standard Dickey-Fuller methodology. |
Extent | 5766666 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2009-04-06 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0088284 |
URI | http://hdl.handle.net/2429/6822 |
Degree |
Doctor of Philosophy - PhD |
Program |
Economics |
Affiliation |
Arts, Faculty of Vancouver School of Economics |
Degree Grantor | University of British Columbia |
Graduation Date | 1997-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_1997-196658.pdf [ 5.5MB ]
- Metadata
- JSON: 831-1.0088284.json
- JSON-LD: 831-1.0088284-ld.json
- RDF/XML (Pretty): 831-1.0088284-rdf.xml
- RDF/JSON: 831-1.0088284-rdf.json
- Turtle: 831-1.0088284-turtle.txt
- N-Triples: 831-1.0088284-rdf-ntriples.txt
- Original Record: 831-1.0088284-source.json
- Full Text
- 831-1.0088284-fulltext.txt
- Citation
- 831-1.0088284.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0088284/manifest