#install.packages("urca")
library(urca)

set.seed(123)
e=rnorm(500)
(x=cumsum(e))
plot(x,type="l")


# ADF Unit Root Test 

(dx=diff(x))
n=(length(dx))
n
(lagx=x[2:n-1])
tr=1:498
ndx=dx[2:n]
lagdx=dx[1:n-1]

summary(lm(ndx~0+lagx+lagdx))

df<-ur.df(x, type="none", lags=1)
summary(df)

summary(lm(ndx~lagx+lagdx))

df.d<-ur.df(x, type="drift", lags=1)
summary(df.d)

summary(lm(ndx~tr+lagx+lagdx))

df.t<-ur.df(x, type="trend", lags=1)
summary(df.t)