#library(tseries)
library(vars)
library(dynlm)
library(urca)
library(tsDyn)

sample1<-("http://kanggc.iptime.org/book/data/korea(77-93).txt")
sample1_dat<-read.delim(sample1,header=T)

gnp<-ts(sample1_dat$gnp, start=c(1977,1), frequency=4)
m2<-ts(sample1_dat$m2, start=c(1977,1), frequency=4)

n<-length(gnp)

ly<-log(gnp)
lm<-log(m2)

dy=diff(ly)
dm=diff(lm)

lylm<-data.frame(ly, lm)
dydm<-data.frame(dy, dm)

lm.res<-lm(ly~lm-1)
summary(lm.res)

summary(ur.df(lm.res$resid, type="drift", lags = 4))
summary(ur.df(lm.res$resid, type="none", lags = 4))
summary(ur.df(lm.res$resid, type="none", selectlags = c("AIC")))
summary(ur.df(lm.res$resid, type="drift", selectlags = c("AIC")))
summary(ur.df(lm.res$resid, type="none", lags = 1))