#### Abstract

This paper investigates the synchronization of fixed income markets within Eurozone countries using the new wavelet based methodology. Conventional wavelet methods that use multivariate set of variables to calculate pairwise correlation and cross correlation lead to spurious correlation due to possible relationships with other variables, amplification of type-1 errors, and results, in the form of large set of erroneous graphs. Given these disadvantages of conventional wavelet based pairwise correlation and cross-correlation method, we avoid these limitations by using wavelet multiple correlation and multiple cross correlations to analyze the relationships in Eurozone fixed income markets. Our results based on this methodology indicate that Eurozone fixed income markets are highly integrated and this integration grows with timescales, and hence there is almost no scope for independent monetary policy and bond diversification in these countries.

#### 1. Introduction

On January 1, 1999, a monetary union Known as European Monetary Union (EMU) was formed by group of eleven European Union countries. This was followed by the introduction of euro as the common currency. The European Central Bank began operating and carrying out the common monetary policy, and all EU-11 government bills and bonds began to be denominated in Euro. The introduction of the euro has been one of the most important events for global financial markets in the last decade and it has become one of the three major currencies in the world after its introduction, taking its place alongside the US dollar and the Japanese yen. An immediate consequence of the adoption of the common currency has been the convergence of Eurozone interest rates and the integration of fixed income markets (Hartmann et al. [1]).

A voluminous body of literature has focused on the integration of equity markets (Bartram and Wang [2], Hardouvelis et al. [3], and Diamandis [4]); however, there is relatively scant literature which investigates integration between the bond markets. Synchronization of bond markets may be induced by diverse types of reasons like similarity in monetary policy reactions functions, increased monetary policy coordination, and globalization induced closer alignment of economic fundamentals (Hunter and Simson [5]). Under such synchronized and integrated system of international markets, conduct of monetary policy by central banks becomes even much more challenging. It is generally argued that the financial market globalization is responsible for reducing the scope for independent monetary policy action by strengthening the relationship between national fixed income markets. However, if the long-term interest rates are determined by national rather than international factors, monetary policy makers are still independent in framing the domestic monetary policy (Barassi et al. [6]). This is true not only for large open rather small open economies also. In essence, long-term interest rates are synonymous to other assets prices or policy instruments. They may also be expected to move together when financial markets are substantially deregulated. It is therefore important to determine the extent to which interest rates are synchronized in different countries for a number of reasons. Fujihara and Mougoué [7] underlined several reasons as to why researchers should know the integration in financial markets in general and fixed income markets in particular. These include the following. (i) Knowledge of linkages in international financial markets could be useful in knowing how economic disturbances are transmitted between countries. (ii) The integration of international financial markets could affect the conduct of monetary policy. This is because of the fact that the degree of independence of domestic stabilization policies from foreign policies is inversely related to the extent to which domestic real interest rates are linked with foreign rates.

There is a scant but significant literature available on the comovement between long term interest rates and much of this research mainly focuses on cointegration and error correction between the long-term interest rates (see, e.g., Barassi et al. [6], Degennaro et al. [8], Kiran [9], Smith [10], and Yang [11]). However, these econometric methodologies deal with just two time frames: the short and the long run with restrictive assumption of variables being integrated of order one (nonstationary).

An alternative methodology of correlation is based on wavelets that considers nonstationarity as an inherent tendency of data and has ability to analyze the time series over different time scales. The technique of wavelets has frequently been applied in number of recent works (Dar et al. [12, 13] and Tiwari et al. [14]). Some of the recent work has also focused on the comovement of equity prices in different countries. Rua and Nunes [15] and Raghavan et al. [16], for example, use conventional wavelets to calculate pairwise multiscale correlations of stock returns for different countries. Conventional wavelet based correlation and cross correlation, however, use simple pairwise correlations to analyze the relationship between different time series. This comes with large number of disadvantages; for example, one needs to calculate, plot, and compare a large number of wavelet correlation and cross-correlation graphs (In this study we are using wavelets to decompose the eleven Eurozone stock returns into four detail components. Using the pairwise wavelet correlation and cross correlation, we would have end up with 55 wavelet correlation graphs and 220 cross-correlation graphs.). Moreover, within the multivariate context, pairwise correlation coefficient could be spurious due to possible relationships with the stock prices of other countries. Finally, pairwise correlations in multiscale are prone to amplify the type-1 errors due to* experimentwise error rate* (see Cohen et al. [17]). For instance, at a given wavelet scale, 5% level of significance doing pairwise wavelet correlation significance tests with nine unrelated series of stock returns (number of stock returns we are analyzing in our study) will inflate the chance of type-1 error to . This increases the chance of finding significant correlation to 94% at given wavelet scale somewhere among 55 tests instead of mere 5%.

Given these disadvantages of conventional wavelet correlation and cross correlation, this study uses the wavelet multiple correlation and cross correlation recently proposed by Fernández-Macho [18] to analyze the integration in the Eurozone bond market over different timescales. The proposed methodology is useful at least in three respects. Overall correlation within the multivariate set of different timescales in stock markets can be viewed just in two plots, spurious correlation obtained from the pairwise correlations within the multivariate set of stock returns could be avoided, and finally, the proposed methodology could be useful in providing the protection against type-1 errors. Our results based on this methodology indicate that the Eurozone bond markets are highly integrated at studied timescales and this integration grows stronger with lower frequencies (higher time scales).

The remainder of the paper is organized as follows. The motivation and methodology used in the paper is described in Section 2. Section 3 gives the data description and discussion of results and finally, Section 4 draws the conclusions and policy implications.

#### 2. Motivation and Methodology

Much of the empirical research on integration in cross-country bond markets focuses on the short run and long run relationship between long term interest rates within cointegration and error correction frame work. One of the limitations of these models is that the definition of long run and short run is narrow. Moreover, these models have been designed to deal with just two time frames: the short and the long run with restrictive assumption of variables being integrated of order one . However, central banks do not operate at two timescales only; rather, they operate at more than two timescales separately (see Aguiar-Conraria et al. [19]). This operation of central banks at different timescales also leads to heterogeneous behaviour of long term interest rates across different timescales. Therefore, it becomes imperative to analyze the integration of long term interest rates at different timescales (frequencies).

##### 2.1. The Wavelet Transforms

To overcome the problems pointed above, wavelet analysis has been proposed. Wavelet analysis shares several features in common with Fourier analysis but has the advantage of capturing features in the underlying series that vary across both time and frequency. Wavelets are mathematical functions that decompose data into different frequency components, after which each component is studied with a resolution matched to its scale. These functions are generated by dilations and translations of a single function called a* mother wavelet*. More explicitly, wavelets are defined as
where is a normalization factor to ensure that wavelets are comparable across scales and time series. Here, the dilation factor “” controls the width of the wavelet: it stretches it into a long wavelet function to measure the low frequency movements (long-run cycles), and it is compressed into a short wavelet function to measure the high frequency movements whereas the translation factor “” determines the time location of the wavelet. To be a mother wavelet, must obey a basic rule, known as the* wavelet admissibility condition*:
where is the Fourier transform of . To guarantee that , the wavelet function must satisfy the conditions and . The first condition means that must be an oscillatory function with* zero mean* and the second condition ensures that the wavelet function has* unit energy*.

The* Continuous wavelet transform* (CWT) of a time series with respect to is given by the following convolution:
The inverse wavelet transform can be defined so that can be reconstructed by means of the formula
Provided satisfies the admissibility condition (2), the* discrete wavelet transform* (DWT) (Daubechies [20]), in particular, allows the decomposition of a signal into components of different frequency. There are two basic wavelet functions, the* father wavelet * and the* mother wavelet *. The formal definitions of the father and mother wavelets are the functions
where in a level decomposition and denotes the maximum scale sustainable by the number of data points and the two types of wavelets stated above satisfiy
Thus, the scaling function integrates to 1 and reconstructs the smooth and low-frequency parts of a signal, whereas the wavelet function integrates to 0 and describes the detailed and high-frequency parts of a signal. In this way, by applying a -level multiresolution decomposition analysis we can provide a complete reconstruction of the signal partitioned into a set of -frequency components so that each component corresponds to a particular range of frequencies.

The wavelet representation of a discrete signal in is given by
where is the number of multiresolution levels (or scales) and ranges from 1 to the number of coefficients in the specified level. The coefficients , in (7) are the wavelet transform coefficients and can be approximated by the following relations:
The coefficients , known as the* smooth coefficients* represent the underlying smooth behaviour of the time series at the coarse scale , while , known as the* detailed coefficient,* describes the coarse scale deviations from the smooth behaviour and provides progressively finer scale deviations from the smooth behaviour.

The* Maximal overlap discrete wavelet transform* (MODWT) (Percival and Walden [21]) is similar to DWT in that high-pass and low-pass filters are applied to the input signal at each level. However, in the MODWT, the output signal is not subsampled (not decimated) but the filters are upsampled at each level. MODWT is preferred over the DWT, since the MODWT can handle any sample size , whilst the DWT of level restricts the sample size to a multiple of . Also MODWT (wavelet and scaling coefficients) is invariant to circularly shifting the time series under study and its multiresolution detail and smooth coefficients are associated with zero phase filters, two properties that the DWT does not hold and excepting the last few coefficients, the MODWT is not affected by the arrival of new information. Finally, the MODWT wavelet variance estimator is asymptotically more efficient than the same estimator based on DWT, which in turn makes it more suitable when calculating wavelet correlations. More details on MODWT can be found in Percival and Walden [21].

Mathematically, decomposing a time series using the MODWT to -levels involves the application of pairs of filters. The filtering operation at the th level consists of applying a rescaled father wavelet to yield a set of detail coefficients
and a rescaled mother wavelet to yield a set of scaling coefficients
for all times , where is the function to be decomposed. The* rescaled mother * and* father wavelets * for the th level are a set of scale dependent localized differencing and averaging operators and can be regarded as rescaled versions of the originals. The th level equivalent filter coefficients have a width , where is the width of the base filter.

##### 2.2. The Wavelet Correlation and Cross Correlation

The* wavelet variance* is defined as the expected value of if we consider only the nonboundary coefficients. An* unbiased estimator* of the wavelet variance is formed by removing all the coefficients that are affected by the boundary conditions and is given by
where is the number of nonboundary coefficients at the th level. Similarly, the* wavelet covariance* between functions and is defined to be the covariance of the wavelet coefficients at a given scale. The* unbiased estimator* of the wavelet covariance at the th level is given by
where all the wavelet coefficients affected by the boundary are removed. Since the* wavelet correlation estimator* is usually constructed by the MODWT and is simply made up of the wavelet covariance and the wavelet variance for and , The MODWT estimator of the* wavelet correlation* can be expressed as
The* wavelet cross-correlation estimator* is a natural extension of the wavelet correlation estimator and in fact has similar properties (see Percival and Walden [21]).

##### 2.3. The Wavelet Multiple Correlation and Cross Correlation

In the existing wavelet literature, the comovement of macroeconomic variables within the multivariate set of variables has been tested mostly by using pairwise correlations and cross correlations. However, owing to its several limitations the wavelet multiple correlation and cross-correlation methodology has been recently proposed by Fernández-Macho [18] to study Eurozone stock market returns.

In this novel methodology, for any multivariate stochastic process , the* wavelet* coefficients are obtained by applying the MODWT to each process at scale . The* wavelet multiple correlation* (WMC) for the stochastic process , at each scale is defined as the square root of the regression coefficient of determination in that linear combination of variables , , for which such coefficient of determination is maximum. As we know, the coefficient of determination corresponding to the regression of a variable on a set of regressors, , is obtained by the relation , where is the th diagonal element of the inverse of the complete correlation matrix . Therefore, the equation for WMC can be formulated as
where refers to the correlation matrix of , and the max. diag(·) operator provides selection for the largest element in the diagonal of the argument. Alternatively, one can also define the WMC equation in terms of wavelet variances and covariance as
where is chosen so as to maximize and are the fitted values in the regression of on the rest of wavelet coefficients at scale .

Similarly, we may also define the* wavelet multiple cross correlation* (WMCC) as
where is a lag between observed and fitted values of the variable selected as the criterion variable at each scale . Further, it is clear from (15) and (16) that, for , the WMC and WMCC will coincide with the standard wavelet correlation and cross correlation because at the scale , and , where is the estimated coefficient in the regression of on .

Let be the realization of the multivariate stochastic process , . Then, after applying the MODWT of order to each of the univariate time series , , we would obtain length- vectors of MODWT coefficients
Therefore, based on either (14) or (15), a consistent* wavelet multiple correlation estimator* is given by
where the wavelet variances and covariances similar to (11), (12), and (13) are given
respectively, with . Similarly, the consistent wavelet multiple cross correlation estimator may be defined as
The* confidence intervals* are constructed by assuming that is a realization of with as the vectors of wavelet coefficients given by (17). Also, if represents the wavelet correlation obtained by (18), then
where , , and stands for the folded normal distribution (see Fernández-Macho [18]). Therefore, the* confidence intervals* (CI) associated with the sample wavelet correlation coefficient can be constructed by the following relation:

#### 3. Data Description, Results, and Discussion

In this section we use the wavelet multiple correlation and wavelet multiple cross correlation with data from the eleven main Eurozone fixed income markets as follows: Germany, France, Italy, Spain, Netherlands, Belgium, Austria, Greece, Finland, Portugal, and Ireland. The monthly data of long term interest rates from January 2001 to September 2012 have been collected, considering that this period coincides with the formation of European Monetary Union. The data have been taken from the website of European central bank. Selected descriptive statistics of long term interest rates for all the stock market indices are presented in Table 1. The sample mean of long term interest rate is positive for all the eleven stock market indices under consideration. The measure of* skewness* indicates that the long term interest rates for Germany, France, Finland, Austria, and Netherlands are* negative* and for the rest of stock indices it* skewed positively.* The distribution of long term interest rate exhibits* excess kurtosis* for stock indices of Germany, Spain, Finland, Greece, Ireland, Portugal, and Italy, which indicate that data series are* leptokurtic* relative to normal distribution. However, for rest of the stock indices the distributions are* platykurtic* in character. Excluding France, Belgium, and Austria, the* Jarque-Bera* normality test was rejected for all the countries.

In order to calculate the wavelet multiple correlation we begin by decomposing the long-term bond rates into different time scales using MODWT to each of the long term interest rate series. Considering the smaller sample size, we decomposed the monthly long-term interest rates by applying the MODWT with a* Daubechies least asymmetric* (LA) wavelet filter of length *,* commonly denoted by* LA *. This filter is favored because of our small sample size.

Given the maximum decomposition level , we could have attained the maximum level of seven decompositions. However, the number of feasible wavelet coefficients gets small for higher levels; we therefore chose to carry out the wavelet analysis with so that four wavelet coefficients and one scaling coefficient were produced for each country, that is, , , , , and , respectively. The time dynamics of each wavelet coefficient in Table 2 is given by [, ] (Since frequency and time are inversely related, the wavelet coefficients are associated with the frequency interval given by [, ].). The first wavelet coefficient therefore represents the dynamics of time series within two to four months, the second one within four to eight months, and so on. Wavelet multiple correlation within 95% confidence intervals obtained from all the bond rates is shown in Figure 1 and Table 3. It is clear from the correlation plots and the ninth column of Table 3 that multiple correlations are significant and high at all the timescales (Results should, however, be interpreted with caution. Confidence intervals are based on fisher’s result for usual bivariate correlation.). We find that correlation is strong at all the frequencies and this correlation grows stronger with lower frequencies. In particular correlation at the highest frequency (16~32 months) is nearly perfect. The discrepancies between Eurozone bond markets are therefore small and negligible. The wavelet multiple cross correlations for the different wavelet scales with leads and lags are presented in Figure 2 and Table 3. The country that maximizes the multiple correlation against the linear combination of other countries is shown in the upper right corner of Figure 2. We find that multiple cross correlations get stronger with lower frequencies and France maximizes the multiple correlations against the linear combination of other countries. This indicates that France has a potential to lead or lag other markets at all the studied frequencies. However, given the symmetry (zero skewness) in Figure 2, it is difficult to say whether French fixed income market leads or lags the linear combination of other markets. Overall, it is found that Eurozone fixed income markets are nearly perfectly integrated especially at higher time scales.

The investigation of fixed income market linkages in international financial markets is of utmost importance. It has important policy implications for monetary policy independence, modeling long term interest rates, bond portfolio diversification, and fiscal deficit financing. While the conventional correlation has often been used to study the market integration for different countries, it does not reveal the correlation at different time horizons. We study fixed income market integration at different time horizons by exploiting the methodology of wavelets. Our results show that the fixed income markets in Eurozone are highly integrated at all the studied frequencies. The Eurozone fixed markets behave as a single coherent market. The synchronization in these countries therefore provides very little scope of independent monetary policy by member countries. Monetary policy, in essence, enables the central bank to react to the specific problems of the economy. For example, to curb the internal and external demand and supply shocks and subsequently create an ambience of price stability and overall economic stability, monetary policy is important. However, the formation of Eurozone contributes a significant loss to the independence of monetary policy making by member countries. Further, given the very high correlation at different frequencies, there is almost no scope for bond diversification in Eurozone countries. This study could be extended to analyze the dynamic comovement of bond markets during the tranquil as well as recession periods.

#### 4. Conclusion and Policy Implications

We applied a new methodology based on wavelet multiple correlation and wavelet multiple cross correlation to study the integration in the Eurozone fixed income markets. Our results based on wavelet multiple correlation and cross correlation indicated that Eurozone fixed income markets are highly integrated with integration becoming stronger with lower frequencies. At the lowest frequency (long run) corresponding to time dynamics of 16~32 months, Eurozone fixed income markets were found to be almost perfectly integrated. These results are the unique contribution of this study which we have been able to calculate with single shot. Further, we allowed certain number of lags between observed and fitted values from the same linear combination as before at each of the wavelet scales. Our results in particular showed that French fixed income market maximizes the multiple correlation against the linear combination of other countries, hence indicating that France has a potential to lead or lag the linear combination of other markets. Over all, we found strong to very strong linkages in the Eurozone fixed income markets, especially at lower frequencies. From the perspective of national central banks in Eurozone area, we conclude that there is less scope for independent monetary policy. Also, a very high correlation at different frequencies indicates that there is almost no scope for bond diversification in Eurozone countries.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.