data(macro)
y=macro[,"INF"]
KPSS1=kpss_1br(y,model=c("intercept","both")[2],use=c("nw","ba"))
KPSS1$teststat
KPSS1$cval
y[KPSS1$bpoint,]
#Plot
plot.ts(KPSS1$tstats,ylim=range(c(KPSS1$tstats,KPSS1$cval)));grid()
abline(h=KPSS1$cval[1],col="red")
abline(h=KPSS1$cval[2],col="blue")
abline(h=KPSS1$cval[3],col="green")
abline(v=KPSS1$bpoint,col="red",lty=2)
Run the code above in your browser using DataLab