data(inf_M)
#Seasonal dummy only takes quarters 1:3,
#becasue of the presence of common intercept.
OUT<-ipsHEGY(inf_M, deterministic = c(1,1,1), lag.method = "AIC", maxlag = 12,CIPS=TRUE)
OUT$P_HEGY
OUT$U_HEGY
# Simulation of critical values
Run the code above in your browser using DataLab