# Generate cointegrated series with one break
set.seed(123)
n <- 100
e1 <- rnorm(n)
e2 <- rnorm(n)
x <- cumsum(e1)
y <- 0.5 * x + cumsum(e2)
y[51:100] <- y[51:100] + 2 # Add structural break
# Run Maki test
data <- cbind(y, x)
result <- coint_maki(data, m=1, model=0)
print(result)
Run the code above in your browser using DataLab