#ESTIMATIONS UNIT ROOT TESTS OF TIME SERIES

#############
library(urca)
#############

# run alg_dataimport.txt first

#bd
adf.bd.nfl1<-ur.df(na.omit(dat.bd[,"bd.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.bd.nfl1)
ers.bd.nfl1<-ur.ers(na.omit(dat.bd[,"bd.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=5); summary(ers.bd.nfl1)
za.bd.nfl<-ur.za(na.omit(dat.bd[,"bd.nfl.gdp"]),model="both",lag=5); summary(za.bd.nfl)

#fr
adf.fr.nfl1<-ur.df(na.omit(dat.fr[,"fr.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.fr.nfl1)
ers.fr.nfl1<-ur.ers(na.omit(dat.fr[,"fr.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=1); summary(ers.fr.nfl1)
za.fr.nfl<-ur.za(na.omit(dat.fr[,"fr.nfl.gdp"]),model="both",lag=5); summary(za.fr.nfl)

#fn
adf.fn.nfl1<-ur.df(na.omit(dat.fn[,"fn.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.fn.nfl1)
ers.fn.nfl1<-ur.ers(na.omit(dat.fn[,"fn.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=2); summary(ers.fn.nfl1)
za.fn.nfl<-ur.za(na.omit(dat.fn[,"fn.nfl.gdp"]),model="both",lag=5); summary(za.fn.nfl)

#bg
adf.bg.nfl1<-ur.df(na.omit(dat.bg[,"bg.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.bg.nfl1)
ers.bg.nfl1<-ur.ers(na.omit(dat.bg[,"bg.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=4); summary(ers.bg.nfl1)
za.bg.nfl<-ur.za(na.omit(dat.bg[,"bg.nfl.gdp"]),model="both",lag=5); summary(za.bg.nfl)

#nl
adf.nl.nfl1<-ur.df(na.omit(dat.nl[,"nl.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.nl.nfl1)
ers.nl.nfl1<-ur.ers(na.omit(dat.nl[,"nl.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=5); summary(ers.nl.nfl1)
za.nl.nfl<-ur.za(na.omit(dat.nl[,"nl.nfl.gdp"]),model="both",lag=5); summary(za.nl.nfl)

#oe
adf.oe.nfl1<-ur.df(na.omit(dat.oe[,"oe.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.oe.nfl1)
ers.oe.nfl1<-ur.ers(na.omit(dat.oe[,"oe.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=1); summary(ers.oe.nfl1)
za.oe.nfl<-ur.za(na.omit(dat.oe[,"oe.nfl.gdp"]),model="both",lag=5); summary(za.oe.nfl)

#it
adf.it.nfl1<-ur.df(na.omit(dat.it[,"it.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.it.nfl1)
ers.it.nfl1<-ur.ers(na.omit(dat.it[,"it.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=2); summary(ers.it.nfl1)
za.it.nfl<-ur.za(na.omit(dat.it[,"it.nfl.gdp"]),model="both",lag=5); summary(za.it.nfl)

#es
adf.es.nfl1<-ur.df(na.omit(dat.es[,"es.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.es.nfl1)
ers.es.nfl1<-ur.ers(na.omit(dat.es[,"es.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=2); summary(ers.es.nfl1)
za.es.nfl<-ur.za(na.omit(dat.es[,"es.nfl.gdp"]),model="both",lag=5); summary(za.es.nfl)

#pt
adf.pt.nfl1<-ur.df(na.omit(dat.pt[,"pt.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.pt.nfl1)
ers.pt.nfl1<-ur.ers(na.omit(dat.pt[,"pt.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=4); summary(ers.pt.nfl1)
za.pt.nfl<-ur.za(na.omit(dat.pt[,"pt.nfl.gdp"]),model="both",lag=5); summary(za.pt.nfl)

#gr
adf.gr.nfl1<-ur.df(na.omit(dat.gr[,"gr.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.gr.nfl1)
ers.gr.nfl1<-ur.ers(na.omit(dat.gr[,"gr.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=5); summary(ers.gr.nfl1)
za.gr.nfl<-ur.za(na.omit(dat.gr[,"gr.nfl.gdp"]),model="both",lag=5); summary(za.gr.nfl)

#uk
adf.uk.nfl1<-ur.df(na.omit(dat.uk[,"uk.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.uk.nfl1)
ers.uk.nfl1<-ur.ers(na.omit(dat.uk[,"uk.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=1); summary(ers.uk.nfl1)
za.uk.nfl<-ur.za(na.omit(dat.uk[,"uk.nfl.gdp"]),model="both",lag=5); summary(za.uk.nfl)

#sd
adf.sd.nfl1<-ur.df(na.omit(dat.sd[,"sd.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.sd.nfl1)
ers.sd.nfl1<-ur.ers(na.omit(dat.sd[,"sd.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=4); summary(ers.sd.nfl1)
za.sd.nfl<-ur.za(na.omit(dat.sd[,"sd.nfl.gdp"]),model="both",lag=5); summary(za.sd.nfl)

#no
adf.no.nfl1<-ur.df(na.omit(dat.no[,"no.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.no.nfl1)
ers.no.nfl1<-ur.ers(na.omit(dat.no[,"no.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=5); summary(ers.no.nfl1)
za.no.nfl<-ur.za(na.omit(dat.no[,"no.nfl.gdp"]),model="both",lag=5); summary(za.no.nfl)

#dk
adf.dk.nfl1<-ur.df(na.omit(dat.dk[,"dk.nfl.gdp"]),type="drift",lags=5, selectlags="AIC"); summary(adf.dk.nfl1)
ers.dk.nfl1<-ur.ers(na.omit(dat.dk[,"dk.nfl.gdp"]),type="DF-GLS",model="constant",lag.max=1); summary(ers.dk.nfl1)
za.dk.nfl<-ur.za(na.omit(dat.dk[,"dk.nfl.gdp"]),model="both",lag=5); summary(za.dk.nfl)




