'   1.- Genera la serie Xt=mu I(t<t0) + Zt  donde Zt = rho Zt-1 + WNt
'   2.- Realiza un ADF test con constante y tendencia, eligiendo el número de rezagos utilizando
'      el criterio de información de Schwartz



create u 1 100

' Parametros : 
scalar rho = 0.5
scalar t0 = 50
scalar mu = 5
scalar m = 10 

vector(100) xv = 0
smpl @first @first +1
series z = 0
smpl @first+1 @last
series z = rho*z(-1)+nrnd
smpl @first @last
stom(z,vz)


for !t = 1 to t0
xv(!t)  = mu + vz(!t)
next
for !t = (t0+1) to 100
xv(!t) = vz(!t)
next
mtos(xv,x)

x.uroot(adf, trend,maxlag=10)