clear all
set more off
set memory 1000m
cd "c:\ferum\cursos\decon\seniales\datos"
tempfile f1 f2 f3 f4
* view zoom para ampliar, ctrl alt +

import excel "pib_trim.xls", sheet("Hoja1") firstrow

generate fecha = yq(anio,trimestre)
format fecha %tq

tsset fecha, quarterly

gen lpib = ln(pib_ivf)
gen lpib_des=ln(pib_ivf_des)

* tsline lpib
* tsline lpib || tsline lpib_des

* ssc install sax12
* findit sax12

* Desestacionalizacion con el X12-ARIMA

* sax12 lpib ,satype(single)
* sax12 lpib, satype(single) transfunc(auto) ammodel((0,3,2)(0,0,1)) outao(ao2002.2) outauto(ao ls tc) 
 
* importo d11 que es la serie ajustada por estacionalidad 
*sax12im "lpib.out" ,ext(d11)

* tsline lpib_des || tsline lpib_d11 

* Borro los archivos auxiliares que genero
*sax12del lpib. keep(d11)

* Periodograma
* pergram lpib

* densidad espectro parametrica
*arima lpib, arima(3,0,1)

*arima d1.lpib, arima(3,0,0)

*psdensity new_sd new_freq 
*twoway line new_sd new_freq 

tsfilter hp ciclo_pib_hp = lpib_des 

gen ten_pib_hp = lpib_des - ciclo_pib_hp

*tsline ciclo_pib_hp  ten_pib_hp lpib_des
tsline   ten_pib_hp lpib_des


/*

tsfilter bk ciclo_TCR_bk = tcr_2010 , minperiod(18)
tsfilter cf ciclo_TCR_cf = tcr_2010 , minperiod(18)

tsline ciclo_TCR_hp ciclo_TCR_bk ciclo_TCR_cf 

corr ciclo_TCR_hp ciclo_TCR_bk ciclo_TCR_cf