powered by
pmax(ts)
where T is number of observations.
Note: [1.6] = 1
data(enrollment) p<-pmax(enrollment) library(urca) summary(ur.df(enrollment,type="trend",lag=p,selectlag="BIC"))
Run the code above in your browser using DataLab