A Novel Approach for Testing Fractional Cointegration in Panel Data Models with Fixed Effects

Fractional cointegration in time series data has been explored by several authors, but panel data applications have been largely neglected. A previous study of ours discovered that the Chen and Hurvich fractional cointegration test for time series was fairly robust to a moderate degree of heterogeneity across sections, out of the six tests considered. Therefore, this paper advances a customized version of the Chen and Hurvich methodology to detect cointegrating connections in panels with unobserved fixed effects. Specifically, we develop a test statistic that accommodates variation in the long-term cointegrating vectors and fractional cointegration parameters across observational units. The behaviour of our proposed test is examined through extensive Monte Carlo experiments under various data generating processes and circumstances. The findings reveal that our modified test performs quite well comparatively and can successfully identify fractional cointegrating relationships in panels even in the presence of idiosyncratic disturbances unique to each cross-sectional unit. Furthermore, the proposed modified test procedure established the presence of long-run equilibrium between the exchange rate and labour wage of 36 countries agricultural markets.

Keywords:

Subject: Business, Economics and Management - Econometrics and Statistics

MSC: 94A08; 68U10

1. Introduction

While cointegration techniques have proven useful in economic and financial research for decades in examining equilibrium relationships between non-stationary time series data, the ability of fractional cointegration to accommodate non-integer degrees of integration and additional nonstationary elements addresses limitations of traditional approaches. Fractional cointegration relaxes assumptions regarding integration orders and potential randomness in time series, allowing for more robust modeling. As econometric tools progress to incorporate more complex real-world variations, fractional techniques lend increased accuracy to describing equilibrium conditions over time. Continued exploration of fractional specifications holds promise for deeper insights into dynamics between financial and economic variables [1,2,3].

The complex research surrounding fractional cointegration in multivariate time series has led to an abundance of assessment techniques for detecting such interdependencies. A wealth of studies have created and evaluated various statistical approaches for testing fractional cointegration, ranging from early explorations that provided foundational methods to more recent innovations that have built upon previous work. Among the many contributions to this area of inquiry are tests proposed by Johansen and Nielsen [1], Marmol and Velasco [4], Chen and Hurvich [5], Johansen [6], Robinson [7], Nielsen [8], Avarucci and Velasco [9], asak [10], Johansen and Nielsen [11], Souza et al. [12] - each offering novel ways of quantifying the perplexing dynamics within nonstationary data comprised of multiple correlated components evolving through fractal processes over time.

Leschinski et al. [13] research uncovered substantial shortcomings in how some fractional cointegration tests operate with limited data, as evidenced by their thorough Monte Carlo analysis of these methods with time series. Certain procedures demonstrated weak detection abilities or skewed false positive rates when short-run elements were interrelated. This generalization covered the strategies put forth by Nielsen and Shimotsu [14] and somehow Nielsen and Shimotsu [14] (or Robinson and Yajima [15]), Marmol and Velasco [4], and Hualde and Velasco [16]. While investigating spurious test outcomes with confined samples, their work highlighted necessary improvements for robust evaluation of cointegrated associations over the long term.

In previous work by Olaniran and Ismail [17], the author specifically focused on determining the appropriateness of the existing fractional cointegration test for the fixed effects panel model. The Chen and Hurvich fractional cointegration test [5] originally developed for time series data was observed to withstand a moderate level of cross-sectional heterogeneity. Therefore, in this paper, we modified the test procedure to fully accommodate the fixed effect fractional cointegration structure in panel data.

2. Related Work

Both fractional and standard cointegration were defined concurrently by Engle and Granger [18], but standard cointegration has received significantly more attention. As a result, there has been a need to focus heavily on fractional cointegration in recent years due to the paucity of research, particularly for panel data as opposed to time series data.

Panel data analysis is a means of investigating a specific subject across numerous places over a fixed length of time that is becoming increasingly popular among social and behavioural science academics. Panel data combines time series and cross-sections to improve data quality and quantity in ways that would be impossible if only one of these two dimensions were utilized [19].

Standard cointegration only allows integer values for the memory parameters

δ

and

ν

, and unit root theory is used to test for cointegration’s existence. The fractional cointegration framework, on the other hand, is more general because it enables fractional values for the memory parameter and any positive real number for

ν - δ

. As indicated by the work of [5,7,20], fractional cointegration has gotten a lot of attention recently. For time series data, all of these publications assume that the observed series is either bivariate or that the cointegrating rank is one.

Robinson [21] observed that in the presence of correlation between

X_{t}

and

ϵ_{t}

, the OLS estimator will be inconsistent, and proposed a narrowband least squares estimator (NBLSE) of the long memory parameter in the frequency domain, which was further studied in Marinucci and Robinson [20]. Robinson and Yajima [15] presented methods for finding the cointegrating rank, as in Chen and Hurvich [5], who focused on the estimation of the space of cointegrating vectors without emphasizing the use of time and cross-sections.

While Chen and Hurvich [22] tapered narrowband least squares estimator of the cointegration parameter utililzes a family of tapers, alternative approaches exist. Their methodology difference the data

p - 1

times to account for potential

{(p - 1)}^{t h}

order polynomial trends, rendering the resultant estimator insensitive to deterministic polynomial trends within the series. However, other scholars propose directly modelling such trends to comparable effect. After differencing, they propose multiplying the data by a tapered sequence of constants dependent on the positive integer p to emphasize more recent observations. Yet some argue for applying non-linear tapers or transforming variables prior to estimation to better address non-stationarities. In short, while their proposal offers one valid method, room remains to explore complementary techniques addressing cointegration under non-stationary conditions.

Chen and Hurvich [5] developed an innovative test for fractional cointegration based on their multivariate common-components model allowing memory parameters to vary across factors. This novel approach recognizes the cointegrating rank may exceed one, with the original vectors split into orthogonal fractional subspaces yielding errors with fluctuating memory attributes from vectors drawn across different subspaces. Their fractional cointegration technique using a common-components structure with changeable memory dynamics across factors proved insightful for assessing long-run relationships between nonstationary series. Also, the complex relationships shared between financial variables often emerge through time. Resolving these intricate connections requires parsing cointegrating residuals to uncover the inherent long-term patterns. By applying a Gaussian semiparametric approach using an appropriately scaling bandwidth, one can derive a stable gauge of the persistence inherently linking economic factors within a specified cointegrating framework. With this metric in hand, we can neatly isolate cointegrating subsets and rigorously examine whether fractional processes underlie the observed synchronization. Such techniques help bring order to the initially chaotic web of interdependencies that defines economic landscapes.

Ergemen and Velasco [23] explored a large dataset with numerous variables and time periods using a panel data model containing fixed factors that enable cross-sectional interdependence, persistent information, and fractionally built-in inaccuracies. Their investigation gives a broad analysis for stationary and non-stationary signs, as corroborated through Monte Carlo tests that exhibited effectiveness in implementation. As suggested in a prior study by Robinson and Velasco [24], the strategy was further expanded to alternative estimation methods involving fixed effects and Generalized Method of Moments, a technique applying instruments to acquire consistent and proficient estimates from unobservable variables.

Idiosyncratic shocks can be utilized to highlight spatial dependency in a broader model by permitting nonfactor endogeneity, which results in a cointegrated system analysis in the traditional sense, as in Ergemen [25]. By comparing mean group estimates with pooled estimates produced from a homogeneous version of the model, panel unit-root and related hypothesis testing, as well as slope parameter homogeneity tests, could be constructed.

Leschinski et al. [13] study explored fractional cointegration across various econometric approaches. Nine fractional cointegration tests evaluating spectral density and residual-based methods were compared. When short-term influences contained interdependences, several techniques uncovered noticeably amplified size. Their efficacy diverged as well, with power varying between common-component frameworks opposite triangular systems. Overall, the analysis illuminated how the presence of interconnected transient elements impacts the performance of fractional cointegration testing procedures. It was also found that the effectiveness of the test procedures diverged substantially between the two modelling situations. The test procedure of Chen and Hurvich [5] specifically had considerably higher power for stationary systems under the common component specification, whereas the techniques of Robinson and Yajima [15] and Hualde and Velasco [16] proved to be less able at identifying fractional cointegration. While the test by Chen and Hurvich [5] strongly determined cointegrating associations assuming steady components, Robinson and Yajima [15] methodology and that of Hualde and Velasco [16] saw diminishing prospective for detecting fractional cointegrating connections.

However, the methods presented by Robinson [7], Nielsen [8], Zhang et al. [26] are robust to short-term correlation and appealing due to their simplicity. None have developed a comprehensive test for fractional integration and cointegration that is applicable to all panel data works. Olaniran and Ismail [17] identified this gap by evaluating the existing fractional cointegration tests of time series data for their suitability in fixed effect panel data exhibiting moderate to high heterogeneity. It was found that only the method developed by Chen and Hurvich [5] could withstand the moderate level of heterogeneity seen in panel data. Therefore, in this paper, we expand upon the work of Olaniran and Ismail [17] by modifying Chen and Hurvich [5] test to specifically address fractional cointegration in fixed effect panel data exhibiting varying degrees of heterogeneity between observational units.

3. Modified Chen and Hurvich [5] Test for Fixed Effect Panel Model

Consider a balanced panel model with fixed effects estimating cointegration relationships. The model examines the long-run equilibrium between variables across

(n)

individual units over

(t)

time periods, creating a full sample of

(N = n * t)

observations. The model is defined as follows:

\begin{matrix} u_{i j} & = α_{i} + θ v_{i j} + {(1 - H)}^{- δ} ε_{1 i j} \\ v_{i j} & = {(1 - H)}^{- ν} ε_{2 i j}, \end{matrix}

(1)

This framework assumes a constant cointegrating parameter

(θ)

linking independent variables

(v_{i j})

and dependent variables

(u_{i j})

for all

(i)

individuals. The model controls for individual-specific characteristics through fixed effects parameters

(α_{i})

. Fractional cointegration exists if two conditions are met. First, the independent and dependent variables share a common stochastic trend making them nonstationary yet moving together in the long run. Second, their linear combination is stationary, implying convergence to equilibrium in the error correction form. Proper specification requires testing the orders of integration to validate the theoretical underpinnings of the cointegrating relationship. Overall, this approach explores how divergences from long-run steady states across units are self-correcting over repeated time periods. The null hypothesis of no fractional cointegration holds if

ν = δ

, whereas the alternative of cointegration is supported if

ν > δ

The bivariate system described in (1) can be simplified to Phillips [27] classical cointegration model if

δ

equals 0 and

ν

is equal to 1, known as

C I (1, 1)

. The fractional lag operator

{(1 - H)}^{- ν}

is derived using

{(1 - H)}^{- ν} = \frac{Γ (j + ν)}{Γ (ν) Γ (j + 1)}

, where the Gamma function

(Γ (h))

is defined as

Γ (h) = \int_{0}^{\infty} ω^{h} e^{- ω} d ω

Equation (1) can be represented as:

\begin{matrix} \begin{matrix} u_{i j} - α_{i} & = θ v_{i j} + {(1 - H)}^{- δ} ε_{1 i j} \\ v_{i j} & = {(1 - H)}^{- ν} ε_{2 i j} . \end{matrix} \end{matrix}

(2)

If we let

w_{i j} = u_{i j} - α_{i}

, we have:

\begin{matrix} \begin{matrix} w_{i j} & = θ v_{i j} + {(1 - H)}^{- δ} ε_{1 i j} \\ v_{i j} & = {(1 - H)}^{- ν} ε_{2 i j} . \end{matrix} \end{matrix}

(3)

It is evident that from the demeaning applied to equation (1), equation (3) removes the individual cross-sectional parameters

(α_{i})

from the original model, consequently simplifying the time series representation. As elucidated by Chen and Hurvich [5], the revised formulation within equation (3) allows for estimation using the Tapered Narrow Band Least Square (TNBLS) technique, an approach aptly suited for this adjusted specification. Nonetheless, one must take care that the demeaning process sufficiently accounts for individual effects so as not to confound the simplified structure with omitted nuances across units.

In the paper by Chen and Hurvich [5], a complex-valued taper function

(q_{j})

is defined as follows for analysing time series data. The taper function varies the weight applied to observations across the time period in a smooth manner. It is mathematically expressed as

q_{j} = \frac{1}{2} (1 - {exp}^{- i 2 π (j - 1 / 2) s^{- 1}}), j = 1, 2, \dots, s .

(4)

Correspondingly, the discrete tapered Fourier transform of a series

η_{j}

and the cross-periodogram are provided by equations using the taper weights to calculate a weighted sum and power:

ω_{η, k}^{'} = {(2 π \sum_{j = 1}^{s} {| q_{j}^{p - 1} |}^{2})}^{- 0.5} \sum_{j = 1}^{s} q_{j}^{p - 1} η_{j} {exp}^{- i λ_{k} j},

(5)

I_{η \bar{η}, k}^{'} = ω_{η, k}^{'} {\bar{ω}}_{η, k}^{'}

(6)

respectively.

The averaged tapered periodogram estimated using m number of bandwidths is then determined as the mean of the periodograms across the frequencies. This provides a smoothed view of the variation between the two time series being examined over different cycles in the data.

{\hat{F}}_{η \bar{η}, k}^{'} (m) = 2 π s^{- 1} \sum_{k = 1}^{m} ℜ I_{η \bar{η}, k}^{'}, 1 \leq m \leq \frac{s}{2} .

(7)

Thus, the long-memory parameter

θ

in equation (3) was estimated as shown in equation (8) below:

{\hat{θ}}_{m} = \frac{{\hat{F}}_{v w}^{'} (m)}{{\hat{F}}_{v v}^{'} (m)}

(8)

where the fixed parameter

m \geq 1

. Alternatively, by setting

m = s / 2

and avoiding differencing and tapering techniques, we derive the familiar ordinary least squares (OLS) approach. The OLS estimator provides a simpler means to the same end compared to holding (m) fixed, foregoing pre-processing of the data. In both cases, the goal is to obtain the most accurate approximation of

θ

possible, though the complexity and computational demands vary depending on the chosen estimation technique.

Theorem 1.

In the context of a fixed effect fractional cointegrated panel model defined by equation (3) and adhering to the specified assumptions, the long-memory parameter can be estimated using equation (8). Consequently, the modified test statistic is given by

M_{C H} = \frac{\hat{ν} - \hat{δ}}{\sqrt{V_{C H} / m}}

, where

V_{C H} = 0.5 (\frac{Γ (4 p - 3) Γ^{4} (p)}{Γ^{4} (2 p - 1)})

. Here, p denotes the a priori differencing parameter for

U_{i j} = {u_{i j}, v_{i j}}

. Under the null hypothesis

H_{0}

, the test statistic

M_{C H}

follows a standard normal distribution, i.e.,

M_{C H} \overset{d}{\to} N (0, 1)

Proof.

Required to show that:

(M_{C H} = \frac{\hat{ν} - \hat{δ}}{\sqrt{V_{C H} / m}} \sim N (0, 1),

under

H_{0}

Suppose we redefine

M_{C H}

such that

\hat{a} = \hat{ν} - \hat{δ}

, then:

M_{C H} = \sqrt{m / V_{C H}} \hat{a}

Lemma 1.

The characteristics function of a continuous random variable U is defined as

ψ_{u} (t) = E (e^{i t u})

Thus, the characteristics function of

\hat{a}

ψ_{\hat{a}} (t) = E (e^{i t \hat{a}})

E (e^{i t \hat{a}}) = ψ (0) + i ψ^{'} (0) t + i^{2} ψ^{″} (0) t^{2} + i^{3} O (t^{3})

, where

i^{2} = - 1

. Consequently, the characteristics function of

M_{C H}

is:

\begin{matrix} \begin{matrix} ψ_{M_{C H}} (t) & = ψ_{\hat{a}} (\sqrt{m / V_{C H}} t) \\ ψ_{\hat{a}} (\sqrt{m / V_{C H}} t) & = ψ (0) + i ψ^{'} (0) \sqrt{m / V_{C H}} t + i^{2} ψ^{″} (0) {(\sqrt{m / V_{C H}} t)}^{2} + i^{3} O ({(\sqrt{m / V_{C H}} t)}^{3}), \end{matrix} \end{matrix}

(9)

since

ψ (0) = 1

ψ^{'} (0) = 0 \equiv E (\hat{a}) = 0

under

H_{0}

, similarly,

ψ^{″} (0) = V_{C H} \equiv V a r (\hat{a})

. So,

\begin{matrix} \begin{matrix} ψ_{\hat{a}} (\sqrt{m / V_{C H}} t) & = 1 - V_{C H} {(\sqrt{m / V_{C H}} t)}^{2} - i O ({(\sqrt{m / V_{C H}} t)}^{3}) \\ = 1 - m t^{2} - i O (m^{3 / 2}) \end{matrix} \end{matrix}

(10)

If we assume

m = s / 2

, we have

\begin{matrix} \begin{matrix} ψ_{\hat{a}} (\sqrt{m / V_{C H}} t) & = 1 - s \frac{t^{2}}{2} - i O ({(\frac{s}{2})}^{3 / 2}) \overset{s \to \infty}{\to} e^{- t / 2} \end{matrix} \end{matrix}

(11)

Lemma 2.

The characteristics function of a Gaussian random variable

Q \sim N (0, 1)

with mean 0 and variance 1 is defined as

ψ_{Q} (t) = E (e^{i t Q}) = e^{- t / 2}

Thus,

M_{C H} \overset{d}{\to} N (0, 1)

under

H_{0}

. □

Remark 1.

The theorem establishes that for a fixed effect fractional cointegrated panel model as defined in equation (3), fulfilling assumptions A₁ and A₂, the long-memory parameter can be effectively estimated using equation (8). The modified test statistic

M_{C H} = \frac{\hat{d} - \hat{γ}}{\sqrt{V_{C H} / m}}

converges in distribution to the standardized normal distribution under

H_{0}

(i.e.,

M_{C H} \overset{d}{\to} N (0, 1)

under

H_{0}

4. Simulation Study

To analyze the performance of a balanced fixed effect panel model, we start by considering a scenario where n represents the number of cross-sectional units and t denotes the time periods, leading to a total sample size

N = n \times t

. The model is defined as follows:

\begin{matrix} \begin{matrix} u_{i j} & = α_{i} + θ v_{i t} + {(1 - H)}^{- δ} ϵ_{1 i j}, \\ v_{i j} & = {(1 - H)}^{- ν} ϵ_{2 i j} . \end{matrix} \end{matrix}

(12)

Here,

u_{i t}

and

v_{i t}

are the dependent and independent variables, respectively. The term

μ_{i}

represents the fixed effects specific to each cross-sectional unit, capturing the unobserved heterogeneity. The parameters

θ

δ

, and

ν

are constants that need to be estimated, with

β

being the cointegration parameter. The model includes the lag operator

(1 - H)

, which accounts for the fractional differencing in the error terms

ϵ_{1 i j}

and

ϵ_{2 i j}

In our model, the intercepts

α_{i}

are assigned fixed values of

5, 10, 15, 20,

and 25 for the units

i = 1, 2, \dots, 5

. These intercepts reflect the fixed effects across the different cross-sectional units, highlighting the variation in the baseline levels of the dependent variable across units. To evaluate the performance of the model in finite samples, we conducted Monte Carlo simulations. The data

{(u_{i j}, v_{i j})}^{'}

were generated using the specified model with

θ = 1

. The error terms

ϵ_{i j} = {(ϵ_{1 i j}, ϵ_{2 i j})}^{'}

follow a Gaussian white noise process, where

E (ϵ_{i j}) = 0

V a r (ϵ_{1 i j}) = V a r (ϵ_{2 i j}) = 1

, and

C o v (ϵ_{1 i j}, ϵ_{2 i j}) = ρ

. We examined scenarios with

ρ = 0.0

and

ρ = 0.5

to understand the impact of different levels of correlation between the error terms.

The sample sizes considered in the simulations were

N = 500, 1250,

and 5000, corresponding to time periods

t = 100, 250,

and 1000, respectively. By varying the sample sizes, we aimed to assess how the model performs under different conditions of data availability. Larger sample sizes allow for more robust estimates and better convergence properties, while smaller sample sizes can highlight potential issues related to finite sample biases. This approach helps in understanding the practical applicability of the model in different empirical contexts.

Similar methodologies have been employed in previous studies such as those by [17,28,29,30,31,32]. These studies have explored various aspects of panel data models, including fractional cointegration and the effects of different parameter settings on model performance. By comparing our results with these prior works, we can validate our findings and contribute to the broader literature on fractional cointegration in panel data models.

Table 1 presents the size for the original and modified Chen and Hurvich test. At

ν = δ = 0.3

corresponding to the weakly stationary fractional panel model, for the three chosen significance levels, the original Chen and Hurvich test sizes were highly overestimated. In fact, the null hypothesis was

100 %

incorrectly rejected. On the other hand, the sizes of the modified Chen and Hurvich test were found to be close to the threshold values. These results were equally observed for both low and moderate correlation values. At

ν = δ = 0.6

corresponding to moderate non-stationary fractional panel model, for the three chosen significance levels, the sizes of the original Chen and Hurvich test were highly overestimated, and in fact, the null hypothesis was at least

58 %

incorrectly rejected. On the other hand, the modified Chen and Hurvich test sizes were found to be close to the threshold values. These results were equally observed for both low and moderate correlation values. At

ν = δ = 0.9

corresponding to a high non-stationary fractional panel model, for the three chosen significance levels, the sizes of both the original and modified Chen and Hurvich test were slightly overestimated. These results were equally observed for both low and moderate correlation values. In addition, the sizes of both tests were inconsistent with increasing sample sizes.

Table 2 shows the empirical power for the two tests at varying significance levels. The effect sizes (ES)

E S = ν - δ

varies between

0.3

and

0.6

. At lower ES when both

ν

and

δ

are greater than

0.5

, the empirical power of the original Chen and Hurvich test was found to be extremely lower compared to the modified test, which in most cases achieved

100 %

power. However, when

ν

is greater than

0.5

but

δ < 0.5

, the empirical power of the original Chen and Hurvich test competes with the modified test. Also, when both

ν

and

δ

are less than

0.5

, the empirical power of the original Chen and Hurvich test is slightly higher than the modified test.

Figure 1 and Figure 2 were used to confirm the distribution of the proposed test statistic

M_{C H}

under both

H_{0}

and

H_{1}

. Under

H_{0}

, the proposed density function was found to be approximately close to

N (0, 1)

distribution. However, under

H_{1}

, a Gaussian distribution with different mean and variance values is suspected.

5. Application to the Monthly Price Estimates of 36 Countries

The validity of the proposed

M_{C H}

test procedure was ascertained using the monthly average estimates of the exchange rate and labour wage in agricultural markets for 36 countries. The aim here is to determine if there is a fractional cointegrating relationship between the United States Dollar ($ USD) exchange rate of each country and the labour wage. These 36 countries are classified as having similar agricultural economies according to the World Bank Development Economics Data Group (DECDG). The dataset was extracted from the World Bank Microdata Library. The dataset covers the period from January 2007 to June 2024 for 27 of the 36 countries, from January 2008 to June 2024 for Armenia, Guinea, and Myanmar, from January 2009 to June 2024 for Yemen, from January 2011 to June 2024 for the Republic of Congo and the Syrian Arab Republic, from January 2012 to June 2024 for Iraq and Lebanon, and from January 2017 to June 2024 for Libya.

Suppose we let

E X_{i t}

represent the $ USD exchange rate for

i = 1, \dots, 36

countries and

t = 1, \dots, 210; 198; 186; 162; 150; 90

months, and

L W_{i t}

represent the labour wage. The underlying fractional cointegration model to test the cointegration hypothesis is given by:

\begin{matrix} \begin{matrix} L W_{i j} & = α_{i} + θ E X_{i j} + Δ^{- δ} ϵ_{1 i j} \\ E X_{i j} & = Δ^{- ν} ϵ_{2 i j} . \end{matrix} \end{matrix}

(13)

We estimated the fractional cointegration parameters

ν

and

δ

by applying a bandwidth parameter

η_{m} = 0.75

, which corresponds to

m = t^{0.75}

for each individual country as well as for the combined panel of countries. The hypotheses tested to determine if the two variables exhibit fractional cointegration are as follows:

$H_{0} :$ $E X_{i j}$ and $L W_{i j}$ are not fractionally cointegrated if $ν = δ$ .
$H_{1} :$ $E X_{i j}$ and $L W_{i j}$ are fractionally cointegrated if $(ν > δ)$ . If the null hypothesis is rejected, it implies there is a long-run relationship between the exchange rate and labour wage such that if there is a deviation, it is only temporary.

Figure 3 shows the time plot of exchange rate

E X (t)

and labour wage

L W (t)

for Afghanistan, Armenia, Bangladesh, Burkina Faso, Burundi, Cameroon, Central African Republic, Chad, Democratic Republic of Congo, Republic of Congo, Gambia and Guinea. The plots show that co-movement which is an indication of cointegration between exchange rate and labour wage is suspected across the 12 countries and that it is more evident in Gambia.

Figure 4 shows the time plot of exchange rate

E X (t)

and labour wage

L W (t)

for Guinea-Bisau, Haiti, Indonesia, Iraq, Kenya, Lao PDR, Lebanon, Liberia, Libya, Malawi, Mali and Mauritania. The plots show that co-movement which is an indication of cointegration between exchange rate and labour wage is suspected across the 12 countries and that it is more evident in Indonesia, Haiti, and Malawi.

Figure 5 shows the time plot of exchange rate

E X (t)

and labour wage

L W (t)

for Mozambique, Myanmar, Niger, Nigeria, Philippines, Senegal, Somalia, South Sudan, Sri Lanka, Sudan, Syrian Arab Republic, and Yemen Republic. The plots show that co-movement which is an indication of cointegration between exchange rate and labour wage is suspected across the 12 countries and that it is more evident in Nigeria, South Sudan, Yemen Rep and the Syrian Arab Republic.

Table 3 presents the estimates of the fractional cointegration parameters

ν

δ

, and

θ

, along with the proposed test statistic

M_{C H}

and its corresponding p-value for 36 countries and their panel. These results allow us to assess whether there is evidence of fractional cointegration between the exchange rate and labour wage in these countries. The fractional cointegration relationship is deemed significant if the test statistic

M_{C H}

is sufficiently large and the p-value is below a conventional threshold (e.g., 0.05).

The majority of the countries exhibit significant fractional cointegration, as indicated by p-values less than 0.05. For example, Afghanistan (

M_{C H} = 8.42

P (Q > | M_{C H} |) = 0.0000

), Bangladesh (

M_{C H} = 7.12

P (Q > | M_{C H} |) = 0.0000

), and Burkina Faso (

M_{C H} = 3.16

P (Q > | M_{C H} |) = 0.0010

) all show strong evidence of fractional cointegration. This means that in these countries, the exchange rate and labour wage share a long-term equilibrium relationship, with deviations from this equilibrium expected to revert over time.

Conversely, ten countries do not exhibit significant fractional cointegration, as their p-values exceed the 0.05 threshold. For instance, Cameroon (

M_{C H} = - 0.82

P (Q > | M_{C H} |) = 0.7940

), the Democratic Republic of Congo (

M_{C H} = - 2.30

P (Q > | M_{C H} |) = 0.9890

), and Libya (

M_{C H} = - 1.54

P (Q > | M_{C H} |) = 0.9380

) do not show evidence of a fractional cointegrating relationship between their exchange rates and labour wages. For these countries, the lack of significant cointegration suggests that the two variables do not share a stable long-term relationship.

The panel data results, which aggregate the information across all 36 countries, indicate significant fractional cointegration (

M_{C H} = 2.23

P (Q > | M_{C H} |) = 0.0130

). This aggregated finding underscores the presence of a common long-term relationship between exchange rates and labour wages across the entire sample of countries, reinforcing the robustness of the individual country results most especially the 26 countries that showed evidence of long-run fractional cointegrating relationship.

6. Conclusion

In this paper, we proposed a modified Chen and Hurvich test for fixed-effect fractional cointegrated panels. An illustration is demonstrated using the simulation of several fractional cointegrated panels with varying real-valued order spanning through weak stationarity to high non-stationarity. The test converges to standard normal distribution under the null hypothesis and diverges under the alternative. Simulation results showed that the proposed test maintained the imposed Type I error rate at a considerably low Type II error, which resulted in high power. Although the proposed test exhibited high power for nonstationary panels, the Type I error is not much different from the original Chen and Hurvich test. Thus, there is still room for improvement in the validity of the test for series with

ν

close to 1. Furthermore, the panel results from the monthly estimates from the 36 countries’ markets confirm the presence of a common fractional cointegration relationship across the entire sample. This highlights the generalizability and robustness of the findings, suggesting that the observed cointegration patterns are not isolated incidents but rather indicative of broader economic dynamics in agricultural markets. Overall, the study underscores the importance of considering fractional cointegration in understanding the interplay between exchange rates and labour wages in the global agricultural economy.

Author Contributions

Conceptualization S.F.O. and O.R.O; methodology, S.F.O.; software, S.F.O. and O.R.O.; validation, J.A. and A.A.A.; formal analysis, S.F.O. and O.R.O; investigation, J.A. and A.A.A.; resources, S.F.O.; data curation, S.F.O.; writing—original draft preparation, S.F.O.; writing—review and editing, S.F.O., O.R.O., J.A. and A.A.A.; supervision, O.R.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The authors confirm that the data supporting the findings of this study are available within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

Johansen, S.; Nielsen, M. . Likelihood Inference for a Fractionally Cointegrated Vector Autoregressive Model. Econometrica 2012, 80, 2667–2732. [Google Scholar] [CrossRef]
Khodamoradi, T.; Salahi, M.; Najafi, A.R. A Note on CCMV Portfolio Optimization Model with Short Selling and Risk-neutral Interest Rate. Statistics, Optimization & Information Computing 2020, 8, 740–748. [Google Scholar] [CrossRef]
Mukhodobwane, R.M.; Sigauke, C.; Chagwiza, W.; Garira, W. Volatility Modelling of the BRICS Stock Markets. Statistics, Optimization & Information Computing 2020, 8, 749–772. [Google Scholar] [CrossRef]
Marmol, F.; Velasco, C. Consistent Testing of Cointegrating Relationships. Econometrica 2004, 72, 1809–1844. [Google Scholar] [CrossRef]
Chen, W.W.; Hurvich, C.M. Semiparametric Estimation of Fractional Cointegrating Subspaces. The Annals of Statistics 2006, 34, 2939–2979. [Google Scholar] [CrossRef]
Johansen, S. Representation of Cointegrated Autoregressive Processes with Application to Fractional Processes. Econometric Reviews 2008, 28, 121–145. [Google Scholar] [CrossRef]
Robinson, P.M. Diagnostic Testing for Cointegration. Journal of econometrics 2008, 143, 206–225. [Google Scholar] [CrossRef]
Nielsen, M. . Nonparametric Cointegration Analysis of Fractional Systems with Unknown Integration Orders. Journal of Econometrics 2010, 155, 170–187. [Google Scholar] [CrossRef]
Avarucci, M.; Velasco, C. A Wald Test for the Cointegration Rank in Nonstationary Fractional Systems. Journal of Econometrics 2009, 151, 178–189. [Google Scholar] [CrossRef]
asak, K. Likelihood Based Testing for no Fractional Cointegration. Journal of Econometrics 2010, 158, 67–77. [Google Scholar] [CrossRef]
Johansen, S.; Nielsen, M. . Nonstationary Cointegration in the Fractionally Cointegrated VAR Model. Journal of Time Series Analysis 2019, 40, 519–543. [Google Scholar] [CrossRef]
Souza, I.V.M.; Reisen, V.A.; Franco, G.d.C.; Bondon, P. The Estimation and Testing of the Cointegration Order based on the Frequency Domain. Journal of Business & Economic Statistics 2018, 36, 695–704. [Google Scholar] [CrossRef]
Leschinski, C.; Voges, M.; Sibbertsen, P. A Comparison of Semiparametric Tests for Fractional Cointegration. Statistical Papers, 2020; 1–20. [Google Scholar] [CrossRef]
Nielsen, M. .; Shimotsu, K. Determining the Cointegrating Rank in Nonstationary Fractional Systems by the Exact Local Whittle Approach. Journal of Econometrics 2007, 141, 574–596. [Google Scholar] [CrossRef]
Robinson, P.M.; Yajima, Y. Determination of Cointegrating Rank in Fractional Systems. Journal of Econometrics 2002, 106, 217–241. [Google Scholar] [CrossRef]
Hualde, J.; Velasco, C. Distribution-Free Tests of Fractional Cointegration. Econometric Theory 2008, 24, 216–255. [Google Scholar] [CrossRef]
Olaniran, S.F.; Ismail, M.T. A Comparative Analysis of Semiparametric Tests for Fractional Cointegration in Panel Data Models. Austrian Journal of Statistics 2022, 51, 96–119. [Google Scholar] [CrossRef]
Engle, R.F.; Granger, C.W. Co-integration and Error correction: Representation, Estimation, and Testing. Econometrica: Journal of the Econometric Society, 1987; 251–276. [Google Scholar] [CrossRef]
Gujarati, D.N. Basic Econometrics. Forth Edition. Singapura: McGraw-Hill, 2003. [Google Scholar]
Marinucci, D.; Robinson, P.M. Narrow-band analysis of nonstationary processes. The Annals of Statistics 2001, 29, 947–986. [Google Scholar] [CrossRef]
Robinson, P.M. Efficient tests of nonstationary hypotheses. Journal of the american statistical association 1994, 89, 1420–1437. [Google Scholar] [CrossRef]
Chen, W.W.; Hurvich, C.M. Estimating fractional cointegration in the presence of polynomial trends. Journal of Econometrics 2003, 117, 95–121. [Google Scholar] [CrossRef]
Ergemen, Y.E.; Velasco, C. Estimation of Fractionally Integrated Panels with Fixed Effects and Cross-section Dependence. Journal of econometrics 2017, 196, 248–258. [Google Scholar] [CrossRef]
Robinson, P.M.; Velasco, C. Efficient inference on fractionally integrated panel data models with fixed effects. Journal of Econometrics 2015, 185, 435–452. [Google Scholar] [CrossRef]
Ergemen, Y.E. System estimation of panel data models under long-range dependence. Journal of Business & Economic Statistics 2019, 37, 13–26. [Google Scholar]
Zhang, R.; Robinson, P.; Yao, Q. Identifying cointegration by eigenanalysis. Journal of the American Statistical Association 2019, 114, 916–927. [Google Scholar] [CrossRef]
Phillips, P.C. Optimal Inference in Cointegrated Systems. Econometrica: Journal of the Econometric Society, 1991; 283–306. [Google Scholar] [CrossRef]
Jamil, S.A.M.; Abdullah, M.A.A.; Kek, S.L.; Olaniran, O.R.; Amran, S.E. Simulation of Parametric Model Towards the Fixed Covariate of Right Censored Lung Cancer Data. Journal of Physics: Conference Series. IOP Publishing, 2017, Vol. 890, p. 012172. [CrossRef]
Olaniran, O.R.; Abdullah, M.A.A. Subset Selection in High-Dimensional Genomic Data using Hybrid Variational Bayes and Bootstrap priors. Journal of Physics: Conference Series. IOP Publishing, 2020, Vol. 1489, p. 012030. [CrossRef]
Popoola, J.; Yahya, W.B.; Popoola, O.; Olaniran, O.R. Generalized Self-Similar First Order Autoregressive Generator (GSFO-ARG) for Internet Traffic. Statistics, Optimization & Information Computing 2020, 8, 810–821. [Google Scholar]
Bala, N.M.; bin Safei, S. A Hybrid Harmony Search and Particle Swarm Optimization Algorithm (HSPSO) for Testing Non-functional Properties in Software System. Statistics, Optimization & Information Computing, 2021. [Google Scholar]
Olaniran, S.F.; Olaniran, O.R.; Allohibi, J.; Alharbi, A.A.; Ismail, M.T. A Generalized Residual-Based Test for Fractional Cointegration in Panel Data with Fixed Effects. Mathematics 2024, 12, 1172. [Google Scholar] [CrossRef]

Figure 1. Normal Q-Q plot of the modified [5] test under

H_{0}

Figure 1. Normal Q-Q plot of the modified [5] test under

H_{0}

Figure 2. Normal Q-Q plot of the modified [5] test under

H_{1}

Figure 2. Normal Q-Q plot of the modified [5] test under

H_{1}

Figure 3. Time plot of exchange rate

E X (t)

and labour wage

L W (t)

for Afghanistan, Armenia, Bangladesh, Burkina Faso, Burundi, Cameroon, Central African Republic, Chad, Democratic Republic of Congo, Republic of Congo, Gambia and Guinea.

Figure 3. Time plot of exchange rate

E X (t)

and labour wage

L W (t)

for Afghanistan, Armenia, Bangladesh, Burkina Faso, Burundi, Cameroon, Central African Republic, Chad, Democratic Republic of Congo, Republic of Congo, Gambia and Guinea.

Figure 4. Time plot of exchange rate

E X (t)

and labour wage

L W (t)

for Guinea-Bisau, Haiti, Indonesia, Iraq, Kenya, Lao PDR, Lebanon, Liberia, Libya, Malawi, Mali and Mauritania.

Figure 4. Time plot of exchange rate

E X (t)

and labour wage

L W (t)

for Guinea-Bisau, Haiti, Indonesia, Iraq, Kenya, Lao PDR, Lebanon, Liberia, Libya, Malawi, Mali and Mauritania.

Figure 5. Time plot of exchange rate

E X (t)

and labour wage

L W (t)

for Mozambique, Myanmar, Niger, Nigeria, Philippines, Senegal, Somalia, South Sudan, Sri Lanka, Sudan, Syrian Arab Republic, and Yemen Republic.

Figure 5. Time plot of exchange rate

E X (t)

and labour wage

L W (t)

for Mozambique, Myanmar, Niger, Nigeria, Philippines, Senegal, Somalia, South Sudan, Sri Lanka, Sudan, Syrian Arab Republic, and Yemen Republic.

Table 1. The simulation results of empirical Type I error rates for both the original test proposed by Chen and Hurvich [5], denoted as

T_{C H}

, and the modified version of the same test, referred to as

M_{C H}

, are evaluated across a range of parameters. This evaluation considers varying levels of the parameter

ν

, which represents the fractional differencing parameter, different correlation levels denoted by

ρ

, different level of significance denoted by

γ

and various sample sizes denoted by N.

Table 1. The simulation results of empirical Type I error rates for both the original test proposed by Chen and Hurvich [5], denoted as

T_{C H}

, and the modified version of the same test, referred to as

M_{C H}

, are evaluated across a range of parameters. This evaluation considers varying levels of the parameter

ν

, which represents the fractional differencing parameter, different correlation levels denoted by

ρ

, different level of significance denoted by

γ

and various sample sizes denoted by N.

			$γ = 0.01$			$γ = 0.05$			$γ = 0.10$
	$ν$	Method/N	500	1250	5000	500	1250	5000	500	1250	5000
$ρ = 0$	$0.3$	$T_{C H}$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$
	$0.3$	$M_{C H}$	$0.019$	$0.015$	$0.014$	$0.060$	$0.053$	$0.048$	$0.090$	$0.092$	$0.080$
	$0.6$	$T_{C H}$	$0.896$	$0.917$	$0.871$	$0.967$	$0.976$	$0.946$	$0.980$	$0.989$	$0.966$
	$0.6$	$M_{C H}$	$0.032$	$0.023$	$0.018$	$0.077$	$0.067$	$0.057$	$0.116$	$0.118$	$0.097$
	$0.9$	$T_{C H}$	$0.041$	$0.041$	$0.050$	$0.131$	$0.109$	$0.132$	$0.198$	$0.169$	$0.195$
	$0.9$	$M_{C H}$	$0.031$	$0.027$	$0.045$	$0.095$	$0.082$	$0.110$	$0.152$	$0.151$	$0.173$
$ρ = 0.5$	$0.3$	$T_{C H}$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$
	$0.3$	$M_{C H}$	$0.018$	$0.015$	$0.009$	$0.054$	$0.051$	$0.045$	$0.089$	$0.089$	$0.071$
	$0.6$	$T_{C H}$	$0.747$	$0.724$	$0.577$	$0.891$	$0.870$	$0.778$	$0.937$	$0.934$	$0.846$
	$0.6$	$M_{C H}$	$0.024$	$0.018$	$0.019$	$0.077$	$0.063$	$0.056$	$0.116$	$0.108$	$0.097$
	$0.9$	$T_{C H}$	$0.017$	$0.017$	$0.040$	$0.073$	$0.071$	$0.102$	$0.136$	$0.119$	$0.148$
	$0.9$	$M_{C H}$	$0.027$	$0.021$	$0.041$	$0.079$	$0.075$	$0.102$	$0.132$	$0.129$	$0.154$

Table 2. The simulation results of empirical power for both the original test proposed by Chen and Hurvich [5], denoted as

T_{C H}

, and the modified version of the same test, referred to as

M_{C H}

, are evaluated across a range of parameters. This evaluation considers varying levels of the parameter

ν

and

δ

which represents the fractional differencing parameter for the series and residual, different correlation levels denoted by

ρ

, different level of significance denoted by

γ

and various sample sizes denoted by N.

Table 2. The simulation results of empirical power for both the original test proposed by Chen and Hurvich [5], denoted as

T_{C H}

, and the modified version of the same test, referred to as

M_{C H}

, are evaluated across a range of parameters. This evaluation considers varying levels of the parameter

ν

and

δ

which represents the fractional differencing parameter for the series and residual, different correlation levels denoted by

ρ

, different level of significance denoted by

γ

and various sample sizes denoted by N.

			$γ = 0.01$			$γ = 0.05$			$γ = 0.10$
	$δ$	Method/N	500	1250	5000	500	1250	5000	500	1250	5000
$ν = 0.9$	$0.6$	$T_{C H}$	$0.073$	$0.243$	$0.869$	$0.183$	$0.370$	$0.911$	$0.267$	$0.440$	$0.924$
		$M_{C H}$	$0.881$	$0.992$	$1.000$	$0.946$	$0.999$	$1.000$	$0.965$	$1.000$	$1.000$
	$0.3$	$T_{C H}$	$0.106$	$0.297$	$0.901$	$0.228$	$0.404$	$0.942$	$0.311$	$0.480$	$0.956$
		$M_{C H}$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$
$ν = 0.6$	$0.3$	$T_{C H}$	$0.880$	$0.897$	$0.815$	$0.968$	$0.968$	$0.930$	$0.986$	$0.985$	$0.959$
		$M_{C H}$	$0.877$	$0.995$	$1.000$	$0.954$	$1.000$	$1.000$	$0.978$	$1.000$	$1.000$
	0	$T_{C H}$	$0.866$	$0.900$	$0.890$	$0.966$	$0.969$	$0.963$	$0.983$	$0.985$	$0.979$
		$M_{C H}$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$
$ν = 0.4$	$0.1$	$T_{C H}$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$	$1.000$
		$M_{C H}$	$0.943$	$0.999$	$1.000$	$0.849$	$0.995$	$1.000$	$0.970$	$1.000$	$1.000$

Table 3. Fractional cointegration parameters

ν, δ, θ

estimates and the proposed test statistic and p-value

M_{C H}

and

P (Q > | M_{C H} |)

for the 36 countries and their panel.

Table 3. Fractional cointegration parameters

ν, δ, θ

estimates and the proposed test statistic and p-value

M_{C H}

and

P (Q > | M_{C H} |)

for the 36 countries and their panel.

	N	$\hat{θ}$	$\hat{ν}$	$\hat{δ}$	$M_{C H}$	$P (Q > \| M_{C H} \|)$	Fractional Cointegration ?
Afghanistan	210	-0.60	0.87	-0.25	8.42	0.0000	Yes
Bangladesh	210	-0.66	0.45	-0.50	7.12	0.0000	Yes
Burkina Faso	210	0.59	0.34	-0.08	3.16	0.0010	Yes
Burundi	210	0.00	0.43	0.08	2.65	0.0040	Yes
Cameroon	210	0.63	-0.02	0.09	-0.82	0.7940	No
Central African Republic	210	0.62	0.32	-0.07	2.91	0.0020	Yes
Chad	210	0.64	0.35	0.03	2.40	0.0080	Yes
Congo, Dem. Rep.	210	0.27	0.12	0.43	-2.30	0.9890	No
Gambia, The	210	-0.08	0.46	0.00	3.45	0.0000	Yes
Guinea-Bissau	210	0.62	0.36	0.00	2.67	0.0040	Yes
Haiti	210	0.10	0.62	0.53	0.70	0.2430	No
Indonesia	210	0.00	0.76	0.04	5.40	0.0000	Yes
Kenya	210	-0.28	0.35	-0.21	4.20	0.0000	Yes
Lao PDR	210	0.44	-0.21	-0.05	-1.18	0.8800	No
Liberia	210	0.10	0.73	0.01	5.35	0.0000	Yes
Malawi	210	1.36	0.89	0.01	6.54	0.0000	Yes
Mali	210	0.56	0.33	0.06	2.01	0.0220	Yes
Mauritania	210	-1.31	0.34	-0.02	2.66	0.0040	Yes
Mozambique	210	-0.06	0.83	-0.11	7.04	0.0000	Yes
Niger	210	0.59	0.33	0.11	1.63	0.0510	Yes
Nigeria	210	1.50	0.67	0.00	4.99	0.0000	Yes
Philippines	210	-0.60	0.63	0.34	2.17	0.0150	Yes
Senegal	210	0.60	0.28	0.04	1.81	0.0350	Yes
Somalia	210	0.00	0.05	-0.08	0.92	0.1790	No
South Sudan	210	0.51	0.88	0.02	6.47	0.0000	Yes
Sri Lanka	210	1.37	0.02	0.91	-6.62	1.0000	No
Sudan	210	0.41	0.56	0.36	1.52	0.0640	Yes
Armenia	198	0.18	1.00	0.05	6.93	0.0000	Yes
Guinea	198	0.38	-0.01	0.10	-0.78	0.7830	No
Myanmar	198	0.10	0.78	-0.02	5.85	0.0000	Yes
Yemen, Rep.	186	0.80	0.59	-0.13	5.07	0.0000	Yes
Congo, Rep.	162	0.76	-0.23	-0.05	-1.22	0.8890	No
Syrian Arab Republic	162	0.80	0.43	-0.09	3.48	0.0000	Yes
Iraq	150	0.00	0.01	-0.14	0.92	0.1790	No
Lebanon	150	0.00	0.48	0.01	3.08	0.0010	Yes
Libya	90	0.00	-0.43	-0.15	-1.54	0.9380	No
Panel	6732	1.17	0.85	0.77	2.23	0.0130	Yes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

A Novel Approach for Testing Fractional Cointegration in Panel Data Models with Fixed Effects

Abstract

1. Introduction

2. Related Work

3. Modified Chen and Hurvich [5] Test for Fixed Effect Panel Model

4. Simulation Study

5. Application to the Monthly Price Estimates of 36 Countries

6. Conclusion

Author Contributions

Funding

Data Availability Statement

Conflicts of Interest

References

MDPI Initiatives

Important Links

Subscribe