data = csvread('dane.csv');

%oczyszczamy dane ze œredniej
data = data - repmat(mean(data),size(data,1),1);

%wektory ze zmiennymi
yh_aux = data(:,1);
yf_aux = data(:,2);
yn_aux = data(:,3);
ynf_aux = data(:,4);
c_aux = data(:,5);
cf_aux = data(:,6);
pih = data(:,7);
pif = data(:,8);
pin = data(:,9);
pinf = data(:,10);
dw = data(:,11);
dwf = data(:,12);
i_aux = data(:,13);
i_f = data(:,14);
de = data(:,15);

%oczyszczamy z trendu (filtr Hodricka-Prescotta)
yh = yh_aux - hpfilter(yh_aux,1600);
yf = yf_aux - hpfilter(yf_aux,1600);
yn = yn_aux - hpfilter(yn_aux,1600);
ynf = ynf_aux - hpfilter(ynf_aux,1600);
c = c_aux - hpfilter(c_aux,1600);
cf = cf_aux - hpfilter(cf_aux,1600);
i = i_aux - hpfilter(i_aux,1600);

%save prepared data
clear data yh_aux yf_aux yn_aux ynf_aux c_aux cf_aux i_aux
save('polish_data.mat');

dynare NK_open_estimation