load ipi.txt

fechasIPI = timefmt ([ 12, 1975, 1, length(ipi) ]);
ipi_fechas=[fechasIPI, ipi];
save ipi_fechas.txt ipi_fechas
plot(ipi)

log_ipi=log(ipi);

plot(log_ipi)

logipi_fechas=[fechasIPI,log_ipi]; 
save logipi_fechas.txt logipi_fechas

help autodhr
PaP=12./(0:6)
TVPaP=[1 1 1 1 1 1 1;1 0 0 0 0 0 0]

[VAR,P,TVP,oar]=autodhr(log_ipi,12,[],[],PaP,TVPaP)

NVR=VAR(2:8)./VAR(1)

filt=0;
[trend,season,cycle,irreg]=dhrfilt(log_ipi,P,TVP,VAR,12,filt,0);

DHR_estacional=season(:,1);

DHR_tendencia=trend(:,1);

logipi_COMP_fechas=...
   [timefmt2(12,1975,1,size(log_ipi,1)),log_ipi,trend(:,1),season(:,1),irreg];
save logipi_COMP_fechas.txt logipi_COMP_fechas

plot([log_ipi,trend(:,1)])

dtrend=diff(trend(:,1));

logipi_DT_fechas=[timefmt2(12,1975,2,size(dtrend,1)),dtrend];
save logipi_DT_fechas.txt logipi_DT_fechas

plot(dtrend)

plot(season(:,1))

plot(irreg)

[VAR2,P2,TVP2]=autodhr(log_ipi,12,[],[34],PaP,TVPaP,1)

NVR2=VAR2(2:8)./VAR2(1)

[trend2,season2,cycle2,irreg2]=dhrfilt(log_ipi,P2,TVP2,VAR2,12,filt,1);
logipi_COMP2_fechas=...
   [timefmt2(12,1975,1,size(log_ipi,1)),log_ipi,trend2(:,1),season2(:,1),irreg2];
save logipi_COMP2_fechas.txt logipi_COMP2_fechas

% plot([log_ipi,trend2(:,1)])
% plot(season2(:,1))
% plot(irreg2)

dtrend2=diff(trend2(:,1));

logipi_DT2_fechas=[timefmt2(12,1975,2,size(dtrend,1)),dtrend,dtrend2];
save logipi_DT2_fechas.txt logipi_DT2_fechas

plot([dtrend,dtrend2])

[VAR3,P3,TVP3,oar3]=autodhr(log_ipi,12,[.1,.45,36,500,2],[],PaP,TVPaP)
NVR3=VAR3(2:8)./VAR3(1)

[VAR4,P4,TVP4,oar4]=autodhr(log_ipi,12)
NVR4=VAR4(2:8)./VAR4(1)

load vabindus.txt

log_vabindus=log(vabindus);

plot(log_vabindus)

logvab_fechas=[timefmt2(4,1980,1,size(log_vabindus,1)),log_vabindus];
save logvab_fechas.txt logvab_fechas

[logipi_trim] = acumula(log_ipi,3,0); %trimestralizacion

logipi_trim_fechas=[timefmt2(4,1975,1,size(logipi_trim,1)),logipi_trim];
save logipi_trim_fechas.txt logipi_trim_fechas
plot(logipi_trim)

[VAR5,P5,TVP5,oar5]=autodhr(logipi_trim,4,[],[],[],[1 1;1 0],1)

[trend5,season5,cycle5,irreg5]=dhrfilt(logipi_trim,P5,TVP5,VAR5,4,filt,1);
logipi_COMP5_fechas=...
    [timefmt2(4,1975,1,size(logipi_trim,1)),logipi_trim,trend5(:,1),season5(:,1),irreg5];
save logipi_COMP5_fechas.txt logipi_COMP5_fechas

% plot([logipi_trim,trend5(:,1)])
% plot(season5(:,1))
% plot(irreg5)

dtrend5=diff(trend5(:,1));

logipi_DT5_fechas=[timefmt2(4,1975,2,size(dtrend5,1)),dtrend5];
save logipi_DT5_fechas.txt logipi_DT5_fechas

plot(dtrend5)

plot([dtrend,dtrend2])

[trend6,season6,cycle6,irreg6]=dhrfilt(log_vabindus,[Inf],[1;1],VAR5(1:2),4,filt,1);

logvab_COMP6_fechas=...
    [timefmt2(4,1980,1,size(log_vabindus,1)),log_vabindus,trend6(:,1),irreg6];
save logvab_COMP6_fechas.txt logvab_COMP6_fechas
dtrend6=diff(trend6(:,1));
logvab_DT6_fechas=[timefmt2(4,1980,2,size(dtrend6,1)),dtrend6];
save logvab_DT6_fechas.txt logvab_DT6_fechas

%plot([log_vabindus,trend6(:,1)])
%plot([irreg5(21:119), irreg6])

plot([dtrend5(21:118) dtrend6])


[VAR7,P7,TVP7,oar7]=autodhr(logipi_trim,4,[],25)

[trend7,season7,cycle7,irreg7]=dhrfilt(logipi_trim,P7,TVP7,VAR7,4,filt,1);
logipi_COMP7_fechas=...
   [timefmt2(4,1975,1,size(logipi_trim,1)),logipi_trim,trend7(:,1),season7(:,1),irreg7];
save logipi_COMP7_fechas.txt logipi_COMP7_fechas

% plot([logipi_trim,trend7(:,1)])
% plot(season7(:,1))
% plot(irreg7)

dtrend7=diff(trend7(:,1));

logipi_DT7_fechas=[timefmt2(4,1975,2,size(dtrend7,1)),dtrend7];
save logipi_DT7_fechas.txt logipi_DT7_fechas

plot([dtrend5,dtrend7])

[trend8,season8,cycle8,irreg8]=dhrfilt(log_vabindus,P7(1),TVP7(:,1),VAR7(1:2),4,filt,1);

logvab_COMP8_fechas=...
   [timefmt2(4,1980,1,size(log_vabindus,1)),log_vabindus,trend8(:,1),irreg8];
save logvab_COMP8_fechas.txt logvab_COMP8_fechas
dtrend8=diff(trend8(:,1));
logvab_DT8_fechas=[timefmt2(4,1980,2,size(dtrend8,1)),dtrend8];
save logvab_DT8_fechas.txt logvab_DT8_fechas

% plot([vabindus,trend8(:,1)])
% plot([irreg7(21:119), irreg8])

plot([dtrend7(21:118) dtrend8])