Skip to main content

The Official Journal of the Pan-Pacific Association of Input-Output Studies (PAPAIOS)

Restatement of the I-O Coefficient Stability Problem

Abstract

The capacity of input-output tables to reflect the structural peculiarities of an economy and to forecast, on this basis, its evolution, depends essentially on the characteristics of the matrix A—matrix of I-O (or technical) coefficients. However, the temporal behaviour of these coefficients is yet an open question. In most applications, the stability of matrix A is usually admitted. This is a reasonable assumption only for a short-medium term. In the case of longer intervals, the question is much more complicated.

We shall empirically discuss this problem by using Romanian input-output tables. Our statistical option was motivated inter alia by the existence of official annual data for two decades (1989–2009).

As an introduction, Sect. 1 characterises the general framework of paper. Section 2—The main characteristics of I-O coefficients as statistical time series—examines the variability of technical coefficients expressed in both volume and value terms. The analysis is convergent to other previous works, confirming that the evolution of these coefficients in real and nominal terms is roughly similar. The main finding of this section is that, on one hand, the I-O coefficients are volatile, but on the other, they are serially correlated.

Consequently, Sect. 3—Attractor hypothesis—examines a possible presence of attractors in corresponding statistical series. The paper describes a methodology to approximate these using new indicators obtained by summation—in columns and rows—of the technical coefficients (colsums sca j and rowsums sra i ). The RAS method is involved as a connecting technique between these indicators and sectoral data.

Section 4—Conclusions—presents the main conclusions of the research and outlines several possible future developments. The database and econometric analysis are presented in Statistical and Econometric Appendix.

JEL Classification: C12, C32, C43, C67.

1 Introduction

1. The capacity of input-output tables to reflect the structural peculiarities of the economy and to forecast, on this basis, its evolution, depends essentially on the characteristics of matrix A of I-O (or technical) coefficients. The so-called Leontief matrix [ ( I A ) 1 ] has proven to be a powerful analytical tool in the investigation of propagated effects induced by inter-industry production chains. Our paper utilises the methodological framework developed in [23, 24, 28, 41, 44].

The temporal behaviour of I-O coefficients is yet an open question. In most applications, the stability of matrix A is usually assumed. This comes from both classical and extended interpretations of the Cobb–Douglas production function. According to Sawyer (p. 327 in [38]), “Under the first of these alternative hypotheses, the a ij will be stable in volume terms. Under the second, the a ij will be stable in value terms”. Generally, the relative stability of the technical coefficients can be considered as a reasonable assumption for a short-medium term. In the case of longer intervals, the question is much more complicated.

2. We shall empirically discuss this problem by using Romanian input-output tables. Our statistical option was motivated inter alia by the existence of official annual data for two decades (1989–2009).

These tables are built on an extended classification comprising 105 branches [17]. To simplify computational operations, the present research relates to a more compact version of 10 sectors [11, 33], as described in Table 1.

Table 1 Sectoral structure of the Romanian input-output tables

The correspondence of this collapsed structure to the original extended nomenclature is detailed in [12]. As in any aggregation, the one proposed in Table 1 implies some losses of information.

Nevertheless, the chosen analysis classification remains sufficiently complex and relevant to involve in this discussion some conceptual anchors of chaos theory. Specifically, we investigate whether the I-O coefficients series could contain sets of attractor points. To answer this question, a methodology for their numerical estimation will be applied to the available data.

3. The robustness of structural changes analysis and of the sectoral dynamic general equilibrium models depends mainly on the temporal behaviour of I-O coefficients. These can be estimated:

• in volume terms (at constant prices), denoted as ca ij ; and

• in value terms (at current prices), usually denoted as a ij .

The first estimation concerns the real economy, while the second relates to the nominal one. These determinations are mediated by the relative prices ( reP ij ).

If cx ij represents the part of sector i’s production (at constant prices p 0 i ) used in sector j, and cX j —total output of the sector j (at constant prices p 0 j ), then:

ca ij = cx ij / cX j
(1)

and

a ij = x ij / X j
(2)

in which the same components of the above ratio are expressed in current prices ( p i and p j , respectively).

Introducing the indices P i = p i / p 0 i and P j = p j / p 0 j , we obtain

a ij = x ij / X j = cx ij P i /( cX j P j )=( cx ij / cX j )( P i / P j )= ca ij reP ij
(3)

where reP ij = P i / P j .

The I-O coefficients at constant prices were estimated using formula (3), which is equivalent to ca ij = a ij / reP ij .

Econometric estimations involve several aggregative indicators resulted from the technical coefficients in value terms, namely:

• Colsums ( sca j ), which summarises the I-O coefficients in columns,

sca j = j a ij with j=fixed;i=1,2,,n
(4)

These approximate the weight of intermediary consumption in the total output of every sector.

• Rowsums ( sra i ), which summarises the I-O coefficients in rows,

sra i = i a ij with i=fixed;j=1,2,,n
(5)

These approximate the contribution of each sector to the intermediary consumption of the entire economy.

2 The Main Characteristics of I-O Coefficients as Statistical Time Series

In the evaluation of the temporal features of I-O coefficients, three questions are relevant:

• Do some peculiarities exist in the co-movement of I-O coefficients real-nominal expression?

• Are I-O coefficients really stable?

• Are these coefficients serially correlated?

The following sections attempt to find answers to these problems.

1. Relating to the first question, in principle the dynamics of real and nominal I-O coefficients are interdependent. On the supply side, the modifications in production costs (reflected by ca ij ) influence the current prices of transactions. On the other hand, the changes in relative prices (reflected by a ij ) have an impact on the demand structure and, consequently, on the size of the output and the conditions (technology, human capital, etc.) in which this is achieved. Due to the complexity of economic life, in each historical period this interdependence has some specific features. This is the reason why statistical evaluation becomes important. Given these, the estimation of the synchronisation degree (SDa) of changes in a ij and ca ij can be conclusive.

1.1. Starting from some proposals advanced in the literature about economic structures and cycles, three concrete formulae are considered.

(a) The first could be referred to as the cosine synchronisation degree (SDa1) since it is estimated as a vectorial angle between time series of I-O coefficients in their double expressions:

SDa1= t ( a ij , t ca ij , t )/ [ ( t a ij , t 2 ) 1 / 2 ( t ca ij , t 2 ) 1 / 2 ]
(6)

(b) The well-known correlation coefficient is often applied in statistical comparisons of real-nominal economic time series (see, for instance, [1, 8, 9, 16, 20, 26, 35, 39]). This Galtung–Pearson synchronisation degree (SDa2) is calculated as a ratio of covariance of series a ij and ca ij to the product of their standard deviations, respectively:

SDa 2 = ( n t a ij , t ca ij , t t a ij , t tca ij , t ) / { [ ( n t a ij , t 2 ( t a ij , t ) 2 ) 1 / 2 ] [ ( n t ca ij , t 2 ( t ca ij , t ) 2 ) 1 / 2 ] }
(7)

(c) A third method used in the economic literature for such analysis is worth mentioning [6, 9, 16]. We shall refer to it as the binary synchronisation degree (SDa3), which measures the proportion in which the compared series evolve in the same direction. Technically, a dummy variable is used, its value being 1 when the respective I-O coefficient increases, and 0 when it decreases or stagnates. If such an alternative assignment is denoted as da ij for series a ij , and, correspondingly, as dca ij for series ca ij , then SDa3 is given as

Sda3= { t ( da ijt dca ijt ) + ( 1 da ijt ) ( 1 dca ijt ) } /n
(8)

n being the number of observations in the sample.

1.2. The above described SDa1, SDa2, and SDa3 do not raise special computational problems, and moreover, are easy to interpret. They have been applied in the series of all 100 technical coefficients, and the obtained results are synthesised in Fig. 1. Therefore, 95 % of SDa1 is positioned within 0.75–1 limits, and only 5 % do not exceed 0.75. At the same time, SDa2 is less than 0.5 in only one-fourth of cases; it is between 0.5–0.75 in 12 % of cases, and exceeds 0.75 in the rest (63 %). The last indicator is even more conclusive: SDa3 is within 0.75–1 in 87 % of cases, and less than 0.65 in none of the cases.

Fig. 1
figure 1

Synchronisation degree (SDa) of changes in a ij and ca ij

Summarising, all calculated synchronisation degrees of changes in a ij and ca ij indicate that the I-O coefficients in both their expressions—in volume and value terms—evolve in a similar manner.

1.3. A more nuanced understanding of this interdependence could be obtained by determining the global variability degree of changes in all I-O coefficients, avca for ca ij and ava for a ij :

avca t = j ( wq it ( i ( ca ijt ca ijt 1 ) 2 ) 1 / 2 )
(9)
ava t = j ( wq it ( i ( a ijt a ijt 1 ) 2 ) 1 / 2 )
(10)

where wq i represents the weight of sector i in the total output of economy.

There were applied two unit root tests for ava and avca: ADF—Augmented Dickey–Fuller and PP—Phillips–Perron. All available options concerning the exogenous (no one, constant, constant plus linear trend) have been computed. The results are detailed in Table 2. Indulgently accepting the stationarity assumption, the pairwise Granger test statistically accredits a certain interconnection between the respective series only on a short run, with the causality direction from avca toward ava (probability of null hypothesis = 0.0881) for one lag, and converse, from ava toward avca (probability of null hypothesis = 0.0943) for two lags. More appropriate for non-stationary series, the test Toda–Yamamoto [43] indicates again on a short run (two lags) an influence of ava on avca (according to F-statistic, and Chi-square, the probability for null hypothesis “ava does not cause avca” represents 0.1107 and, respectively, 0.0869).

Table 2 Unit root tests for ava and avca

Except for 4 years (1991, 2002–2003, and 2005), the ratio of ava to avca was <1 in all periods. This means that the changes in relative prices somehow attenuated the shifts in technical coefficients in volume terms.

2. The examination of the co-movement pattern of changes in the real and nominal expressions of I-O coefficients does not clarify if these are relatively stable (small annual changes) or significantly volatile. This is important for our analysis.

In the case of I-O coefficients, we shall adopt a larger interpretation of volatility as an integrating measure of the frequency and size of the changes registered in their evolution. A comprehensive analysis of volatility determinants exceeds the thematic perimeter of this paper. Briefly, we recall the following factors:

• the performance of preponderantly used technologies that redound to most aspects of costs (labour productivity, energy and raw material intensities, quality of goods and services, length of productive cycles, etc.);

• the dimension, and structure of domestic demand, which influence the scale efficiency and relative prices;

• the openness degree of the country, with its impact on firms’ access to external markets, on import substitution effects, and on productive factors migration;

• the institutional reforms that have a great role in both emerging and developed economies; and

• the operational consequences of macroeconomic policies that can facilitate or, on the contrary, hinder the fructification of comparative advantages for the respective economy.

Quantitatively, the volatility of a given indicator will be approximated by its variation coefficient calculated (for the entire available time series) as follows. If q t is the value of this indicator at moment t (t=1,2,,s) and ω q its level admitted as referential, then this coefficient ( C V ) is determined by

C V = [ ( t ( q t / ω q 1 ) 2 ) / s ] 1 / 2
(11)

In principle, ω q can differ depending on the objectives of analysis. As a first choice, we adopt the sample mean, accommodating expression (11) to the standard deviation formula largely used in modern statistics. Such an approach is suitable in forecasting the volatility for different interested horizons by simple extrapolation of its statistically registered level.

The proposed procedure consists of the following steps:

• For each interval two estimations of the respective indicator are determined: an upper and a lower level. The first is obtained by multiplying the mean of the previous series by (1+ C V ), while the other results similarly but using (1 C V ) as a multiplier. We shall designate these values as Y for the upper level and y for the lower one.

• On this basis, two new means are also computed, mixing the corresponding previous series with Y and y: they will be represented by the symbols M, and m, respectively. The statistical volatility is applied again by multiplying the new M by (1+ C V ) and m by (1 C V ). This procedure is continued as much as it is considered useful (the forecast period being denoted by τ=1,2,,n).

• The difference (Yy) can be admitted as an error ( ef V ) attributable to the initially estimated volatility. The interpretation of results would be facilitated by equalising the starting sample mean to unity.

More formally, for the upper level, we have

Y τ 1 =(1+ C V ) M τ 2 ,τ=1,2,,n
(12)
M τ 1 = ( ( s + τ 2 ) M τ 2 + Y τ 1 ) / ( s + τ 1 ) = ( ( s + τ 2 ) M τ 2 + ( 1 + C V ) M τ 2 ) / ( s + τ 1 ) = M τ 2 ( s + τ 2 + 1 + C V ) / ( s + τ 1 ) = M τ 2 ( s + τ 1 + C V ) / ( s + τ 1 ) = M τ 2 ( 1 + C V / ( s + τ 1 ) )
(13)
Y τ =(1+ C V ) M τ 1 =(1+ C V ) M τ 2 ( 1 + C V / ( s + τ 1 ) )
(14)

A simplification can be obtained by passing to indices ( IY τ = Y τ / Y τ 1 ):

IY τ = { ( 1 + C V ) M τ 2 ( 1 + C V / ( s + τ 1 ) ) } / [ ( 1 + C V ) M τ 2 ] = ( 1 + C V / ( s + τ 1 ) )
(15)

This relationship is valid for τ2 since Y 0 = M 0 =1 and Y 1 =(1+ C V ) M 0 =(1+ C V ). Finally, we have

Y n =(1+ C V ) τ ( 1 + C V / ( s + τ 1 ) ) for τ=2,,n
(16)

Symmetrically, the expression of y n is determined as

y n =(1 C V ) τ ( 1 C V / ( s + τ 1 ) ) ,again for τ=2,,n
(16a)

and

ef Vn = [ ( 1 + C V ) τ ( 1 + C V / ( s + τ 1 ) ) ] [ ( 1 C V ) τ ( 1 C V / ( s + τ 1 ) ) ]
(17)

Therefore, ef Vn is influenced mainly by C V , s, and τ. Figures 2(a) and 2(b) illustrate some indifference curves of the initial C V depending on s and m, estimated under the conditions given in Table 3.

Fig. 2
figure 2

(a) Estimation of the initial C V depending on s and the final desirable ef V (Variant A). (b) Estimation of the initial C V depending on s and the final desirable ef V (Variant B)

Table 3 Estimation of the initial C V depending on s and the final desirable ef V

The presented algorithm can be used in establishing a kind of taxonomy scale of I-O coefficients volatility. Toward this aim, it would be necessary to determine the desirable levels of ef V and the length of τ (that is, the value of n). A possible starting point in this sense can be the expectable financial risk induced by economic decisions linked to forecasted I-O coefficients. Addressing this question requires further research. A possible solution to this problem could be adequately extrapolated in other socio-economic fields.

Returning to the Romanian I-O tables, the variation coefficient, based on formula (11), was computed for all statistical series in 1989–2009 (100 ca ij and 100 corresponding a ij ). The results are summarised in Table 4, which shows that there is no I-O coefficient with C V <0.05 and only one with C V <0.1; instead, 85 % of ca ij and 73 % of a ij are characterised by C V >0.3. The hypothesis that the mean of all C V would be between 0.4–0.65 was tested for both series C V ca ij and C V a ij . The results are presented in Fig. 3.

Fig. 3
figure 3

Probability for the mean of C V ca ij and C V a ij to be situated between 0.4–0.65 (tabulated on abscissa)

Table 4 Tabulation of statistical variation coefficients ( C V )

In many cases, the volatility is so high that the calculated ef V becomes abnormal even for very short intervals. As an example, the evolution of the error attributable to the initially estimated volatility ( ef V ) was determined for three cases: for C V =0.1 (variant 1), C V =0.2 (variant 2), and C V =0.3 (variant 3), during τ=1,2,,15. The results of this exercise are denoted as ef V 1, ef V 2, and ef V 3, and are summarised in Fig. 4. We recall that the computed data represent indices comparatively to the mean level of the statistical series (the mean equalised to 1). For C V =0.3, the difference between the forecasted limits of the respective indicator can reach 0.7 in five years and 0.8 in ten. Even for C V =0.1, the potential forecasting error is hardly acceptable. As we have already shown, the levels calculated for Romanian I-O tables are overall much higher than the simulated (in Fig. 4) values of C V .

Fig. 4
figure 4

Simulated ef V for three variants of C V

3. Like other previous studies, the analysis of Romanian I-O tables confirms that the technical coefficients are volatile. What needs to be documented is the nature of this volatility, and the highly questionable factor is the presence of non-linearities in the respective statistical series. Such a possibility has been revealed in many economic indicators [3, 34]. In the case of Romanian I-O tables, we shall also examine whether the data regarding the technical coefficients are independent or, on the contrary, serially correlated.

It is widely accepted that: “The correlation sum in various embeddings can…be used as a measure of determinism in a time series” (p. 313 in [40]). The BDS test is sensitive to a large variety of possible deviations from independence in time series, including linear dependence, non-linear dependence, or chaos. Concerning this technique, our turns to the conceptual and applicative framework developed in [2, 6, 32]. Thus, the null hypothesis of independent and identically distributed (i.i.d.) data is checked against an unspecified alternative.

For the I-O tables examined in this paper, the BDS test was applied to both categories of coefficients—at constant ( ca ij ) and current prices ( a ij ). Concerning the embedding dimension, we sought to cover an extended range of possibilities. Due to the insufficient length of the statistical series, five such variants were adopted: 2, 3, 4, 5, and 6. As a principal guiding mark, the p-value for the tested null hypothesis was retained, computed for the sample data (normal probability) and for their random repetitions (bootstrap probability). Recent software provided both probabilities (normal and bootstrap) for three options related to the distance used for testing: the fraction of pairs, the standard deviations, and the fraction of range. Therefore, 30 p-values were computed for each technical coefficient, resulting in five dimensions, two tested series (original and bootstrap), and three distances.

The characterisation of the global distributions of the obtained p-values for all series of technical coefficients will be discussed. Two classifications are significant.

First, the p-values for all 3000 estimations are classified according to the following thresholds: under 0.05, 0.05–0.1, 0.1–0.25, and 0.25–1, presented in Fig. 5. This shows that in the case of ca ij , over 75 % of p-values (2252) are below 0.05; if the group 0.05–0.1 is added, the proportion reaches 80 %. The picture is similar for a ij : almost 72 % of tests are estimated with p-values of under 0.05, and approximately 76 % have p-values of less than 0.1. This means that, generally, the series of I-O coefficients (either at constant or current prices) are not independent.

Fig. 5
figure 5

Distribution of BDS tests (whole sample) in terms of p-value

The second application sorts I-O coefficients depending on the number of registered BDS p-values under 0.05. Toward this aim, six classes are delimited: up to 5 times, 5–10, 10–15, 15–20, 20–25, and 25–30. Evidently, the sum of classes is equal to 100 (the totality of coefficients). Figure 6 synthesises this distribution, showing that in each of the 86 ca ij , at least 15 tests had p-values of under 0.05. The result is no different in the case of a ij coefficients: among 90 cases, at least 10 p-values were under 0.05. The similarity of the ca ij and a ij series suggests that the volatility of relative prices does not substantially influence the presence of serial correlation in the data.

Fig. 6
figure 6

Classification of the technical coefficients depending on registered p-values under 0.05

Thus, in this section, we can conclude that, on one hand, the I-O coefficients are volatile, but on the other, they are serially correlated. Both statements have statistical support. More simply stated, we acknowledge a paradox because the high volatility indicates rather the presence of a quasi-disorder, while the serial correlation indicates a possible stable pattern in the analysed time series. The following section focuses on this exciting matter.

3 Attractor Hypothesis

The revealed contradictory combination of relatively high volatility of data and their consistent serial correlation generates a legitimate question: Is this contradiction a sign of a possible presence of an attractor in statistical series?

1. Generally, an attractor is considered a point or a closed subset of points (lines, surfaces, volumes), toward which a given system tends to evolve independently of its initial (starting) state [2931, 36, 37]. Three types are frequently mentioned:

• stable steady states,

• different types of cycles, and

• strange attractors.

The first type is relatively usual in Economics (“At best, the notion of equilibrium might, in practice, be identified with the notion of <attractor>”; p. 34 in [14]). The list of such examples is long, from the optimal rates of accumulation to the extended palette of Phillips curves.

Such points or lines need to be regarded rather as historical (that is, contextually determined) phenomena than as permanent, inflexible benchmarks. It is worth mentioning that some authors considered the “natural rate of unemployment” as a rather weak attractor (p. xiii in [4]).

Taking into account the numerous such applications in economics, the following systematisation of types of stable steady states would be useful:

• stable points,

• constant rates of movement (in different expressions, such as indices, elasticities, ratios, spreads, etc.), and

• bands of evolution.

All these are interesting perspectives in researching I-O tables. However, such a target would require many and sustained efforts. Our target is very narrow, namely, to attempt to identify in the studied statistical series some fixed points as possible attractors. This hypothesis will be used in two sub-variants: fixed points as such or slightly variable points with gradually decreasing influence of unknown factors (cumulated over a time parameter). Besides, the econometric analysis will concentrate on the dynamics of each I-O coefficient, considered separately and not in connection with other series.

Therefore, the evolution of I-O coefficients is conceived as an auto-regressive adaptive process, the differences between their actual and long-run levels being influenced by the past deviations. In the simplest form, such an application for Romanian input-output tables was developed in [10]. In a general notation, if y is the time series of interest, we would have the following relationship:

y t = y ˜ α ( y ( 1 ) y ˜ ) = y ˜ (1+α)αy(1)
(18)

where y ˜ represents the long-run levels of y (or the attractor according to this paper’s terminology). It is assumed that 0<|α|<1, which means that y tends asymptotically towards y ˜ . Correspondingly, the first-order difference operator d(y) is defined as

d ( y ) = y y ( 1 ) = y ˜ ( 1 + α ) α y ( 1 ) y ( 1 ) = y ˜ ( 1 + α ) ( 1 + α ) y ( 1 ) = a 0 a 1 y ( 1 )
(19)

The expression (19) contains the equivalencies a 0 = y ˜ (1+α) and a 1 =(1+α).

To be more realistic, this determination will be relaxed by two amendments. On one hand, the last formula will be extended, with gradually diminishing influence of time. On the other, the auto-regressive process may involve lags of higher orders, not only of the first one, as in (19).

2. Even under such modifications, the approximation of possible attractor points requires the presence of at least one non-differentiated observation in the computational formula. Therefore, it would be preferable to use the statistical series stationary in levels (I(0)). Unfortunately, most of the available data do not observe such a restriction. From this point of view, two already mentioned unit root tests were applied: ADF—Augmented Dickey–Fuller and PP—Phillips–Perron test. Each was computed in three versions for the exogenous variables:

• none (denoted as 1),

• individual effects (denoted as 2), and

• individual effects and individual linear trends (denoted 3).

The p-values calculated for all 100 technical coefficients were grouped as follows: 0–0.05, 0.05–0.1, 0.01–0.25, and 0.25–1.

The corresponding distribution for the technical coefficients at constant prices ( ca ij ) is presented in Figs. 7 and 8. Both unit root tests (ADF and PP) show that in around 80 % of the cases, the p-values exceed 0.1. The same result is found for the technical coefficients at current prices (see Figs. 9 and 10).

Fig. 7
figure 7

ADF tests for ca ij

Fig. 8
figure 8

PP tests for ca ij

Fig. 9
figure 9

ADF tests for a ij

Fig. 10
figure 10

PP tests for a ij

At this point, we are confronted with a problem. The BDS test indicated the presence of temporal correlation in the data for technical coefficients (either at constant or at current prices). As previously mentioned, this finding would justify the identification of possible attractor points in their evolution. Since the series are not stationary in levels, in order to avoid the calculation of attractor points (as levels) by first- or second-order differentiation (a difficult computational task), an indirect way to approximate such points will be proposed.

The first step is to determine colsums ( sca j ) and rowsums ( sra i ) for the technical coefficients at current prices. The resulting series are given in Statistical and Econometric Appendix. With respect to these time series, PANEL analysis did not reveal compelling signs of common explicative parameters. For this reason, they were examined separately. Table 5 shows the p-values of the ADF and PP tests for the sca i series. In only three cases ( sca 2 , sca 3 , and sca 4 ) are the corresponding p-values situated in the proximity of 0.25. Consequently, the series sca i will be used as such in regressions.

Table 5 ADF and PP tests for sca i

Table 6 presents the same indicators for sra i . The introduction of econometric estimations for series sra 5 , sra 8 , and sra 10 as such would clearly be too risky. Consequently, the first two were recalculated by the Hodrick–Prescott filter, obtaining for each the sub-series denoted as HP and HPd (difference between filter and primary data), respectively. The third series ( sra 10 ) was replaced with the corresponding logarithms. Table 7 shows the unit root test results, based on which the new series for sra 5 , sra 8 , and sra 10 were used in regressions.

Table 6 ADF and PP tests for sra i
Table 7 ADF and PP tests for derived series sra 5 , sra 8 , and sra 10

The formula (19) with the mentioned amendments was investigated using different specifications. The proposed selection considered, beside the mentioned premises, the results of tests for omitted or redundant variables, and outliers, also. It has also tried to reduce the econometric compromises as much as possible. For the current paper, several types of relationships were retained according to the scheme given in Table 8. Sometimes dummy variables were introduced to decrease the influence of data outliers.

Table 8 Main econometric relationships

3. The OLS-solution of system SyS1scr (Statistical and Econometric Appendix) was submitted to econometric controls from four standpoints: (a) variance inflation factors, (b) Breusch–Pagan–Godfrey heteroskedasticity test, (c) correlogram squared residuals, and (d) stationarity of residuals.

Concerning the variance inflation factors (Table 9), it is conclusive that more than 77 % of the centred VIFs do not exceed 2, and approximately 15 % are situated between 2 and 3; even the rest do not surpass 5.3. Based on these results, we could accept that the specification of the system SyS1scr is not contaminated in an alarming manner by collinearity effects.

Table 9 Variance Inflation Factors—SyS1scr

The test Breusch–Pagan–Godfrey (Table 10) indicates high enough probabilities for the rejection of heteroskedasticity hypothesis.

Table 10 SyS1scr: heteroskedasticity test Breusch–Pagan–Godfrey

The correlogram of squared residuals was computed for five lags (Table 11). In most cases, Q-statistics are associated with relatively large p-values, which attest a weak serial correlation in the residuals.

Table 11 Correlogram of residuals squared—SyS1scr

Concerning the stationarity of residuals, both unit root tests ADF and PP were applied again, in all available options for exogenous (Table 12). There were thus generated 132 values of the probability the respective residual has a unit root. Out of these, 76.52 % are placed under 0.05, and 10.61 % between 0.05–0.1.

Table 12 ADF and PP unit root tests of residuals SyS1scr

The above presented tests (for collinearity, heteroskedasticity, serial correlation, and stationarity of residuals) show that OLS could be acceptable to estimate the system SyS1scr.

4. The system SyS1scr has been solved using other four techniques: Weighted Least Squares (WLS), Seemingly Unrelated Regression (SUR), Generalised linear models (GLM), and Generalised Method of Moments (GMM). The obtained results are detailed in Statistical and Econometric Appendix.

The solution induced by Weighted Least Squares slightly ameliorates the standard errors, maintaining, however, the parameters of equations practically at the same level as OLS. The differences between Seemingly Unrelated Regression and OLS regarding estimators and coefficients of determination are also insignificant. The same conclusion is valid for the Generalised Linear Models (applied with bootstrap).

The Generalised Method of Moments was involved in variant HAC for the time series (Bartlett and Variable Newey–West). Despite the large number enough of trials, the results were inconclusive. First, in order to obtain a plausible solution, it was necessary to break SyS1scr into three sub-systems—SyS1scaG, SyS1sraG, and SyS1sra8G—which have been separately computed. Secondly, the algorithm did not work with dummies, or these were not introduced casually, but according to the specification test about outliers.

Briefly, the comparative analysis of different techniques suggests as acceptable OLS method. Nevertheless, a problem persists. According to Statistical and Econometric Appendix (System Residual Cross-Correlations—OLS), the disturbances of some relationships represented in SyS1scr are correlated. They reflect, at great extent, the indubitable fact of inter-industry linkages. Obviously, there must be a consistent solution of the question hereby discussed. It could result from a re-specification of the entire system by explicit inclusion in the equations of the factors inducing cross-correlations among input-output technical coefficients, and subsequently applying computational methods that avoid simultaneity effects. But such an approach should need further interdisciplinary research. Until then, I am reluctant to involve techniques which somehow mechanically constrict the cross-correlations of I-O coefficients. Consequently, for the present OLS will keep being involved in the succeeding steps of our approach.

5. Based on the previous system, the fitted sca j f and sra i f can be obtained, but not a ij f as such. To approximate these, the RAS technique was applied. During its half-century existence [42], this method has registered extended applications, including in recent researches [7, 18, 19, 21, 22, 25, 27]. Usually, the starting matrix for every t is the statistical matrix A t 1 , which is adjusted by successive bi-proportional corrections in dependence on exogenously given sectoral outputs. The applicability of such a method for an emergent economy such as in Romania has already been documented [13].

The present paper slightly modifies this procedure, using sca j f and sra i f as column and row restrictions in a RAS algorithm. The resulting technical coefficients (denoted as ra ij ) are relevant from the present research perspective. Notably, ra ij are calculated using the fitted sca j f and sra i f. The formulae, however, are based on the hypothesis that the respective original statistical series contain attractor points. Consequently, the analysis of the differences resra ij = a ij ra ij can be informative. Given the independency of these differences, the assumption that sca j f and sra i f include attractor points and that the derived ra ij contain such compatible points becomes plausible since both sca j and sra i represent simple summations of the corresponding a ij .

Consequently, we return to the BDS test. As in the previous application, the test was applied to both probabilities (normal and bootstrap) in three options related to the distance (fraction of pairs, standard deviations, and fraction of range) and in five dimensions (2, 3, 4, 5, and 6). For each resra ij , 30 p-values were again computed (as before). The distribution of all 3000 p-values is described in Fig. 11. Only one fifth of the p-values do not exceed 0.05. This proportion falls to 8 % in the case of the bootstrap method, which is more relevant for relatively short series.

Fig. 11
figure 11

Distribution of the BDS tests for resra ij

For this reason, as a general approximation, the serial independence of resra ij differences was assumed. Consequently, the probability of attractor points in the data for a ij cannot be neglected.

6. Further on, the attractor points will be estimated based on the following additional assumptions:

• It is admitted that in the proximity of an attractor point, the values of the respective technical coefficients are relatively stable. In other words, first- and higher-order differences tend to disappear.

• In terms of level, the value of the technical coefficient coincides or is close to that of the attractor point. The importance of the presence of observations in level (I(0) problem) in econometric formulae has already been outlined.

• The attractor points are conceived at long-run levels. For large values of t, it is admitted that t 1 0 and t/(t+1)1.

The scheme containing the main econometric relationships will be adapted to these assumptions, the result being the algebraical expressions of attractors in the 9 types of specifications (Table 13) included in SyS1scr. Their symbols are given the prefix a: asca j and asra i .

Table 13 Algebraical attractor definitions

Table 14 presents the approximated attractors for colsums ( asca j ) and rowsums ( asra i ) of the I-O coefficients. These estimations were included as column–row restrictions in a new RAS application concerning all a ij . This algorithm was applied on a matrix compounded by the average levels of the respective statistical coefficients (for the entire interval 1989–2009). Table 15 presents the so-obtained attractor points ( aa ij ).

Table 14 Attractor-points for the colsums and rowsums of technical coefficients
Table 15 Attractor-points for individual technical coefficients ( aa ij )

4 Conclusions

The analysis of Romanian I-O tables (based on surveys for 21 consecutive years) reveals new evidence in favour of the statement that the technical coefficients are volatile (illustrated by the relatively high standard deviation of corresponding series). This affects both determinations of I-O coefficients, either in volume ( ca ij ) or in value terms ( a ij ); the first is referred to as real volatility and the second as nominal volatility. Their dynamic pattern is similar, as confirmed by three measures: (a) the vectorial angle between the series a ij and ca ij , (b) the Galtung–Pearson correlation (also a cosine of the vectorial angle but between their deviations against the mean) and (c) the binary synchronisation degree.

To verify whether or not the I-O coefficients are serially correlated, the BDS procedure was used as a test covering a large variety of possible deviations from independence in the time data. Again, both forms of technical coefficients were studied. Generally, the serial correlation could not be statistically rejected. It is important to mention that this conclusion resulted from a relatively extended database.

Due to these two circumstances—high volatility and serial correlation—the possible presence of attractors in the technical coefficients series was taken into consideration. Such points would be flexibly interpreted not as unchangeable levels but rather as historical (contextually determined) phenomena. This approach is similar to the manner in which other authors regarded the natural rate of unemployment, for instance, as a weak attractor. Consequently, the evolution of I-O coefficients was conceived as an auto-regressive adaptive process, the differences between the actual coefficients and their long-run levels being influenced by the precedent deviations. Since the available series for sectoral coefficients are, as a rule, non-stationary, more aggregate indicators were employed in econometric analysis (column and row sums of I-O coefficients). The RAS technique was used to transform these into sectoral estimations.

The paper’s approach can be considered as an attempt to conciliate the assumption of I-O coefficients’ stability with their undisputable volatility.

Further research could improve on the econometric estimations through structural specifications of the technical coefficients, including their stable co-movements. Thus, more complex econometric specifications must be cautiously adopted, but based on a solid economic motivation.

The possible presence of attractors in the series of I-O coefficients also opens a large research space. A deeper investigation of their determinants—technologies, inter-industry linkages, institutional factors—would be interesting from both the theoretical and the applicative perspective. In addition, it would be relevant to clarify the temporal stability of the attractors themselves.

Statistical and Econometric Appendix

Table 16 Column-sums of the technical coefficients at current prices
Table 17 Row-sums of the technical coefficients at current prices
Table 18 System SYS1scr: Specification
Table 19 SYS1scr estimated by different methods—sample 1990–2009: OLS—ordinary least squares
Table 20 SYS1scr estimated by different methods—sample 1990–2009: WLS—weighted least squares
Table 21 SYS1scr estimated by different methods—sample 1990–2009: SUR—seemingly unrelated regression
Table 22 SYS1scr estimated by different methods—sample 1990–2009: GLM—generalized linear models with bootstrap
Table 23 Comparative estimation output OLS–SUR
Table 24 Generalized method of moments—time series (HAC): Kernel: Bartlett, bandwidth: Variable Newey–West (5), no prewhitening
Table 25 System residual cross-correlations—OLS: ordered by variables, 5 lags

References

  1. Afonso A, Furceri D (2007) Business cycle synchronization and insurance mechanisms in the EU. Working paper no 844/December, European Central Bank, Eurosystem. http://www.ecb.int/pub/pdf/scpwps/ecbwp844.pdf. Accessed 10 June 2012

    Google Scholar 

  2. Belaire-Franch J, Contreras-Bayarri D (2002) The BDS test: a practitioner’s guide. DT 02–01, University of Valencia. http://centros.uv.es/web/departamentos/D10/data/investigacion/PDF132.pdf. Accessed 14 March 2012

    Google Scholar 

  3. Bevilacqua F, van Zon A (2001) Random walks and non-linear paths in macroeconomic time series: some evidence and implications. Maastricht, November. http://www.wu.ac.at/inst/vw1/gee/papers/gee!wp22.pdf. Accessed 12 April 2012

    Google Scholar 

  4. Blanchard O: Preface. In The natural rate of unemployment—reflections on 25 years of the hypothesis. Edited by: Cross R. Cambridge University Press, New York; 1995.

    Google Scholar 

  5. Bordo MD, Helbling T (2003) Have national business cycles become more synchronized? Working paper 10130, NBER. http://www.nber.org/papers/w10130. Accessed 13 July 2011 Bordo MD, Helbling T (2003) Have national business cycles become more synchronized? Working paper 10130, NBER. http://www.nber.org/papers/w10130. Accessed 13 July 2011

  6. Brock W, Dechert D, Scheickman J, LeBaron B: A test for independence based on the correlation dimension. Econom Rev 1996, 15: 197–235. 10.1080/07474939608800353

    Article  Google Scholar 

  7. Capilit E (2009) RAS method. Asian Development Bank, RETA 6483—adopting the supply and use framework towards 1993. System of national accounts compliance in selected DMCs, 2nd data review workshop, 7–12 December. http://www2.adb.org/statistics/icp/files/6483–2DRW/RAS-Method-Capilit.pdf. Accessed 13 July 2011

    Google Scholar 

  8. Cashin P, McDermott CJ, Scott A (1999) The myth of co-moving commodity prices. Discussion paper G99/9, Reserve Bank of New Zealand, December. http://www.rbnz.govt.nz/research/discusspapers/g99_9.pdf. Accessed 8 June 2012

    Google Scholar 

  9. Darvas Z, Szapáry G (2004) Business cycle synchronization in the enlarged EU. Paper revised in October 2004 for the joint ECB-IMF workshop on global financial integration, stability and business cycles: exploring the links, Frankfurt. http://www.imf.org/external/np/seminars/eng/2004/ecbimf/pdf/szapar.pdf. Accessed 11 June 2012

    Google Scholar 

  10. Dobrescu E (2006) Integration of macroeconomic behavioural relationships and the input–output block (Romanian modelling experience). Paper presented at the international conference on policy modelling (Ecomod 2006), Hong Kong, June 28–30, 2006

    Google Scholar 

  11. Dobrescu E (2009) Measuring the interaction of structural changes with inflation. Rom J Econ Forecast, Supplement

    Google Scholar 

  12. Dobrescu E, Gaftea V, Scutaru C: Using the Leontief matrix to estimate the impact of investments upon the global output. Rom J Econ Forecast 2010, 13: 176–187.

    Google Scholar 

  13. Dobrescu E, Gaftea V: On the accuracy of RAS method in an emergent economy. Amfiteatru Econ 2012, 14: 1–20.

    Google Scholar 

  14. Giarini O: Science and economics: the case of uncertainty & disequilibrium. Cadmus 2011,1(2):25–34.

    Google Scholar 

  15. Harding D, Pagan A (1999) Dissecting the cycle. Working paper 13/99, Melbourne Institute. http://fmwww.bc.edu/RePEc/es2000/1164.pdf. Accessed 14 June 2012

    Google Scholar 

  16. Inklaar R, Jong-A-Pin R, De Haan J (2005) Trade and business cycle synchronization in OECD countries. A re-examination. Working paper 1546, category 6: monetary policy and international finance, CESIFO. http://www.cesifo.de/pls/guestci/download/CESifo%20Working%20Papers%202005/CESifo%20Working%20Papers%20September%202005/cesifo1_wp1546.pdf. Accessed 11 June 2012

    Google Scholar 

  17. INSSE—TEMPO (2012) Input–output tables. http://statistici.insse.ro/shop/?lang=en

    Google Scholar 

  18. Jackson RW, Murray AT (2003) Alternate input–output matrix updating formulations. Research paper 2003–6. http://rri.wvu.edu/pdffiles/jacksonmurraywp2003–6.pdf. Accessed 10 Jan 2012

    Google Scholar 

  19. Jian X, Jing H, Yanyun Z (2007) Assessment study for the RAS method based on China’s input–output tables. Management School of Graduate University, Chinese Academy of Sciences. http://www.iariw.org/papers/2007/jian.pdf. Accessed 13 July 2011

    Google Scholar 

  20. Kalemli-Ozcan S, Papaioannou E, Peydró J-L (2010) Financial regulation, financial globalization and the synchronization of economic activity. Working paper 1221/July, European Central Bank, Eurosystem. http://www.ecb.int/pub/pdf/scpwps/ecbwp1221.pdf. Accessed 10 June 2012

    Google Scholar 

  21. Kratena K, Zakarias G (2002) Technical coefficients change by bi-proportional econometric adjustment functions. Paper presented at the 14th international conference on input–output techniques, Montréal, Canada, October 10–15, 2002. http://www.iioa.org/pdf/14th%20conf/kratena_zakarias.pdf. Accessed 7 January 2012

    Google Scholar 

  22. Lahr ML, de Mesnard L: Biproportional techniques in input–output analysis: table updating and structural analysis. Econ Syst Res 2004, 16: 115–134. 10.1080/0953531042000219259

    Article  Google Scholar 

  23. Leontief W: Analiza input–output. Editura Stiintifica, Bucharest; 1970.

    Google Scholar 

  24. Leontief W: Input–output economics. 2nd edition. Oxford University Press, New York; 1986.

    Google Scholar 

  25. Linden JA, Dietzenbacher E: The determinants of structural change in the European Union: a new application of RAS. Environ Plan A 2000, 32: 2205–2229. 10.1068/a32105

    Article  Google Scholar 

  26. Li L, Zhang N, Willett TD (2011) A survey of measurements of interdependence. Claremont Institute for Economic Policy Studies. http://www.cgu.edu/PDFFiles/SPE/Willett/Papers/A%20Survey%20of%20Measurements%20of%20Interdependence.pdf. Accessed 11 June 2012

    Google Scholar 

  27. Mesnard L (2000) Failure of the normalization of the RAS method: absorption and fabrication effects are still incorrect. LATEC—document de travail—economie (1991–2003) 2000–01, LATEC, Laboratoire d’Analyse et des Techniques Economiques, CNRS UMR 5118, Université de Bourgogne. http://leg.u-bourgogne.fr/documents-de-travail/e2000–01.pdf. Accessed 20 July 2011

    Google Scholar 

  28. Miller RE, Blair PD: Input–output analysis: foundations and extensions. 2nd edition. Cambridge University Press, Cambridge; 2009.

    Book  Google Scholar 

  29. Milnor J: On the concept of attractor. Commun Math Phys 1985, 99: 177–195. 10.1007/BF01212280

    Article  Google Scholar 

  30. Milnor J: Comments on the concept of attractor: correction and remarks. Commun Math Phys 1985, 102: 517–519. 10.1007/BF01209298

    Article  Google Scholar 

  31. Milnor JW (2006) Attractor. Scholarpedia 1(11):1815. http://www.scholarpedia.org/article/Attractor. Accessed 5 June 2012

    Article  Google Scholar 

  32. Moloney K, Raghavendra S: Testing for nonlinear dependence in the credit default swap market. Econ Res Int 2011., 2011: Article ID 708704. doi:10.1155/2011/708704 Article ID 708704. doi:10.1155/2011/708704

    Google Scholar 

  33. National Commission for Prognosis, Romania (2012) Data base of the macromodel-Dobrescu from the project strengthening the institutional capacity in evaluating and formulating of macroeconomic policies for economic convergence with EU within the National Commission for Prognosis, SMIS code 27153

    Google Scholar 

  34. Nelson CR, Plosser CI: Trends and random walks in macroeconomic time series: some evidence and implications. J Monet Econ 1982, 10: 139–162. 10.1016/0304-3932(82)90012-5

    Article  Google Scholar 

  35. Nguyen T (2007) Determinants of business cycle synchronization in East Asia: an extreme bound analysis. Working paper 2007/14, DEPOCEN. http://www.depocenwp.org/upload/pubs/NguyenNgocToan/Determinants_of_Business_Cycle_Synchronization_in_East_Asia-rv2_DEPOCENWP.pdf. Accessed 10 June 2012

  36. Rosser JB Jr, Kramer KL Jr: Integrating the complexity vision into mathematical economics. In Complexity and the teaching of economics. Edited by: Colander D. Edward Elgar, Cheltenham Glos; 2000:209–230. http://cob.jmu.edu/rosserjb/mathecon.doc. Accessed 5 June 2012 http://cob.jmu.edu/rosserjb/mathecon.doc. Accessed 5 June 2012

    Google Scholar 

  37. Ruelle D: What is…a strange attractor? Not Am Math Soc 2006,53(7):764–765.

    Google Scholar 

  38. Sawyer JA: Forecasting with input–output matrices: are the coefficients stationary? Econ Syst Res 1992, 4: 325–348. doi:10.1080/09535319200000030 doi:10.1080/09535319200000030 10.1080/09535319200000030

    Article  Google Scholar 

  39. Schirwitz B, Wälde K (2004) Synchronization of business cycles in G7 and EU14 countries. European Commission Directorate General Economic and Financial Affairs, October. http://www.waelde.com/pdf/SchirwitzWaeldeTIERsynchronization.pdf. Accessed 10 June 2012 Schirwitz B, Wälde K (2004) Synchronization of business cycles in G7 and EU14 countries. European Commission Directorate General Economic and Financial Affairs, October. http://www.waelde.com/pdf/SchirwitzWaeldeTIERsynchronization.pdf. Accessed 10 June 2012

  40. Sprott JC: Chaos and time-series analysis. Oxford University Press, New York; 2003.

    Google Scholar 

  41. Stone R (1961) Input–output and national accounts organisation for European economic co-operation, Paris Stone R (1961) Input–output and national accounts organisation for European economic co-operation, Paris

  42. Stone R, Brown JAC: A long-term growth model for the British economy. In Europe’s future in figures. Edited by: Geary RC. North-Holland, Amsterdam; 1962.

    Google Scholar 

  43. Toda HY, Yamamoto T: Statistical inference in vector autoregressions with possibly integrated processes. J Econom 1995, 66: 225–250. 10.1016/0304-4076(94)01616-8

    Article  Google Scholar 

  44. United Nations (1999) Handbook of input–output table: compilation and analysis, studies in methods. Handbook of national accounting, series F, No. 74, Department for Economic and Social Affairs, Statistics Division, New York. http://unstats.un.org/unsd/publication/SeriesF/SeriesF_74E.pdf. Accessed 13 July 2011 United Nations (1999) Handbook of input–output table: compilation and analysis, studies in methods. Handbook of national accounting, series F, No. 74, Department for Economic and Social Affairs, Statistics Division, New York. http://unstats.un.org/unsd/publication/SeriesF/SeriesF_74E.pdf. Accessed 13 July 2011

Download references

Acknowledgements

The author thanks to V. Gaftea, D. Jula, M. Matei, B. Pauna, and C. Saman for their computational assistance. He is also highly grateful to the anonymous referent of the “Journal of Economic Structures” for the suggested analytical extensions.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Emilian Dobrescu.

Additional information

Competing Interests

The author declares that he has no competing interests.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Dobrescu, E. Restatement of the I-O Coefficient Stability Problem. Economic Structures 2, 2 (2013). https://doi.org/10.1186/2193-2409-2-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/2193-2409-2-2

Keywords