x = xlsread('norges2.xls')
% este arquivo usa o spread da taxa preferencial brasileira como medida de 
% prêmio de financiamento externo. Em função da disponibilidade dos dados
% a amostra vai de mar/05 a dez13

output_gap = x(:,1);
interest_rate = x(:,2);
inflation = x(:,3)*4;
%multipliquei a taxa de inflação por 4 para transformar ela em anual
labor = x(:,4);
consumption = x(:,5);
investment = x(:,6);
fx = x(:,7);
wr = x(:,8);


% os dados têm que ir do 3o tri de 96 até o 1o tri de 2013
% O vetor t abaixo da valores para o eixo x do gráfico
% Assim, o intervalo de 0.25 corresponde a 1 trimestre
%1996.75 corresponde ao 3o tri de 1996 e 2013 corresponde ao 1o tri de 
%2013. Seguindo esta lógica, para montar o vetor de datas, teríamos:
% 25: 1o tri
% 50: 2o tri
% 75: 3o tri
% 00: 4o tri (por exemplo, 2014:00 se refere ao 4o tri de 2013)


t = (2005.25:0.25:2014.00)';

figure
set(gcf,'Name','Data for estimation (hp filtered)');
subplot(421)
plot(t,output_gap,'b.-')
title('Hiato do produto')
axis([min(t)-0.25, max(t)+0.25, min(output_gap)-.015 max(output_gap)+.015])


subplot(422)
plot(t,interest_rate,'b.-')
title('Selic')
axis([min(t)-0.25, max(t)+0.25, min(interest_rate)-0.015 max(interest_rate)+0.015])

subplot(423)
plot(t,inflation,'b.-')
title('Inflação')
axis([min(t)-0.25, max(t)+0.25, min(inflation)-0.015 max(inflation)+0.015])

subplot(424)
plot(t,labor,'b.-')
title('Horas')
axis([min(t)-0.25, max(t)+0.25, min(labor)-0.015 max(labor)+0.015])

subplot(4,2,5)
plot(t,consumption,'b.-')
title('Consumo')
axis([min(t)-0.25, max(t)+0.25, min(consumption)-0.015 max(consumption)+0.015])

subplot(4,2,6)
plot(t,investment,'b.-')
title('Investimento')
axis([min(t)-0.25, max(t)+0.25, min(investment)-0.015 max(investment)+0.015])


subplot(4,2,7)
plot(t,fx,'b.-')
title('Câmbio real')
axis([min(t)-0.25, max(t)+0.25, min(fx)-0.015 max(fx)+0.015])

subplot(4,2,8)
plot(t,wr,'b.-')
title('prêmio de financiamento')
axis([min(t)-0.25, max(t)+0.25, min(wr)-0.015 max(wr)+0.015])

% Correspondência entre os nomes das séries na base de dados e no modelo

y = output_gap;
r = interest_rate;
pi = inflation;
N = labor;
C = consumption;
inv= investment;
Q = fx;
rp = wr;