Learn R Programming

kardl (version 0.1.1)

narayan: NARAYAN test

Description

This function performs the Narayan test, which is designed to assess cointegration using critical values specifically tailored for small sample sizes. Unlike traditional cointegration tests that may rely on asymptotic distributions, the Narayan test adjusts for the limitations of small samples, providing more accurate results in such contexts. This makes the test particularly useful for studies with fewer observations, as it accounts for sample size constraints when determining the presence of a long-term equilibrium relationship between variables.

Usage

narayan(model, case = 3, signif_level = "auto")

Value

A list with class "kardl" containing the following components:

  • type: The type of the test, which is "cointegration".

  • case: The case of the test, which is a numeric value (1, 2, 3, 4, or 5).

  • statistic: The calculated F-statistic for the test.

  • k: The number of long-run variables in the model.

  • Cont: The conclusion of the test, indicating whether cointegration is found, inconclusive, or not found.

  • BoundNum: A numeric value indicating the result of the test: 1 for cointegration, 0 for inconclusive, and -1 for no cointegration.

  • siglvl: The significance level used in the test, which can be "auto", "0.10", "0.1", "0.05", "0.025", or "0.01".

  • criticalValues: A numeric vector containing the critical values for the test based on the specified case and significance level.

  • parameter: A character vector containing the names of the long-run parameters in the model.

  • FH0: The null hypothesis of the test, which is a character string describing the hypothesis being tested.

  • Fmodel: The detailed F test results, including the F-statistic, p-value, degrees of freedom, and residual degrees of freedom.

  • warnings: A character vector containing any warnings generated during the test, such as issues with the model specification or data length.

  • method: The method used for the test.

Arguments

model

The kardl obejct

case

Numeric or character. Specifies the case of the test to be used in the function. Acceptable values are 1, 2, 3, 4, 5, and "auto". If "auto" is chosen, the function determines the case automatically based on the model's characteristics. Invalid values will result in an error.

  • 1: No intercept and no trend. This case is not supported by the Narayan test.

  • 2: Restricted intercept and no trend.

  • 3: Unrestricted intercept and no trend.

  • 4: Unrestricted intercept and restricted trend.

  • 5: Unrestricted intercept and unrestricted trend.

signif_level

Character or numeric. Specifies the significance level to be used in the function. Acceptable values are "auto", "0.10", "0.1", "0.05", "0.025", and "0.01". If a numeric value is provided, it will be converted to a character string. If "auto" is chosen, the function determines the significance level automatically. Invalid values will result in an error.

Hypothesis testing

The null hypothesis (H0) of the F Bound test is that there is no cointegration among the variables in the model. In other words, it tests whether the long-term relationship between the variables is statistically significant. If the calculated F-statistic exceeds the upper critical value, we reject the null hypothesis and conclude that there is cointegration. Conversely, if the F-statistic falls below the lower critical value, we fail to reject the null hypothesis, indicating no evidence of cointegration. If the F-statistic lies between the two critical values, the result is inconclusive.

$$ \Delta {y}_t = \psi + \varphi t + \eta _0 {y}_{t-1} + \sum_{i=1}^{k} { \eta _i {x}_{i,t-1} } + \sum_{j=1}^{p} { \gamma_{j} \Delta {y}_{t-j} }+ \sum_{i=1}^{k} {\sum_{j=0}^{q_i} { \beta_{ij} \Delta {x}_{i,t-j} } }+ e_t $$

Cases 1, 3, 5:

$$\mathbf{H_{0}:} \eta_0 = \eta_1 = \dots = \eta_k = 0$$ $$\mathbf{H_{1}:} \eta_{0} \neq \eta_{1} \neq \dots \neq \eta_{k} \neq 0$$

Case 2:

$$\mathbf{H_{0}:} \eta_0 = \eta_1 = \dots = \eta_k = \psi = 0$$ $$\mathbf{H_{1}:} \eta_{0} \neq \eta_{1} \neq \dots \neq \eta_{k} \neq \psi \neq 0$$

Case 4:

$$\mathbf{H_{0}:} \eta_0 = \eta_1 = \dots = \eta_k = \varphi = 0$$ $$\mathbf{H_{1}:} \eta_{0} \neq \eta_{1} \neq \dots \neq \eta_{k} \neq \varphi \neq 0$$

References

Narayan, P. K. (2005). The saving and investment nexus for China: evidence from cointegration tests. Applied economics, 37(17), 1979-1990.

See Also

pssf psst banerjee recmt

Examples

Run this code

kardl_model<-kardl(imf_example_data,
                   CPI~ER+PPI+asym(ER)+deterministic(covid)+trend,
                   mode=c(1,2,3,0))
A<-narayan(kardl_model)
A
cat(paste0("The F statistics=",A$statistic," where k=",A$k,"."))
cat(paste0("\nWe found '",A$Cont, "' at ",A$siglvl,"."))
cat(paste0("\nThe proboblity is = ",sprintf("%.10f",A$Fmodel[["Pr(>F)"]][2])))

# Utilising magrittr:

library(magrittr)
imf_example_data %>% kardl(CPI~ER+PPI+asym(ER)+deterministic(covid)+trend,
                   mode=c(1,2,3,0)) %>% narayan()


# Getting details of the test.
mySummary<-summary(A)
mySummary

# Utilising magrittr:
imf_example_data %>% kardl(CPI~ER+PPI+asym(ER)+deterministic(covid)+trend,
                   mode=c(1,2,3,0)) %>% narayan()  %>% summary()


# Critical Values are
A$criticalValues
# The null hypothesis :
A

# getting H0
mySummary$H0

# Detailed F test
A$Fmodel

Run the code above in your browser using DataLab