
addpath C:\Dynare\5.4\matlab; %C:\dynare\5.4\matlab
addpath("C:\Dynare\5.4\matlab\modules\dseries\src");  % Agregar ruta del módulo dseries de Dynare
addpath("C:\Dynare\5.4\matlab\missing\rows_columns"); % Agregar ruta de manejo de valores faltantes
initialize_dseries_class(); % Inicializar Dseries de Dynare

Data = readtable("C:\Work\1NL_B\Models\U\IN.xlsx","Sheet","Hoja1", "Range", "A1:C159");
U_NSA = dseries(Data.U_NSA, '1984Q1');% Transformar los datos en objetos Dseries

cd(tempdir()); % Carpeta para datos de ajuste estacional X-13ARIMA-SEATS
Aux_X13 = x13(U_NSA).transform('function', 'log').arima('model', '(0 1 1)4'); % Parametros ajuste estacional X-13ARIMA-SEATS 
Aux_X13.x11('save', '(d10 d11)'); % Guardar las series desestacionalizadas y las diferencias
Aux_X13.run(); % Corre el ajuste estacional X-13ARIMA-SEATS
U = Aux_X13.results.d11.data; % Datos desestacionalizados
cd('C:\Work\1NL_B\Models\U'); % Vuelve al directorio de trabajo original

figure;
plot(U.dates, [U.data U_NSA.data]);
legend('U', 'U\_NSA', 'Location', 'best');
title('Comparación de las series U y U\_NSA');

U_SA_table = table(Data.Observation_date, double(U.data), 'VariableNames', {'Date', 'U_SA'});% Crear una tabla con los datos desestacionalizados
writetable(U_SA_table, "C:\Work\1NL_B\Models\U\IN.xlsx", 'Sheet', 'U_SA', 'WriteVariableNames', true);% Exportar la tabla a un archivo de Excel

IN_AUX = table(double(U.data), 'VariableNames', {'U'}) ; 
writetable(IN_AUX, "IN_AUX.csv") ; % Saves the data in a csv file, for use in the model code

dynare UForecast_F.mod ; 

figure;
plot(U.dates,[oo_.SmoothedVariables.U oo_.SmoothedVariables.U_BAR]);
% plot([oo_.forecastMean.U oo_.forecastMean.U_BAR]) ;
