% DRAWZ   Use results from GR_IR to plot GR's Technology, Output, and Hours
%

%         Ellen McGrattan, 3-13-04
%         Revised, ERM, 6-16-04

ly   = log(prd.*hpc);
lh   = log(hpc);
lz   = cumsum(epsz/100);
t0   = (length(ly)-length(lz))+1;
t1   = length(ly);
lh   = lh(t0:t1);
ly   = ly(t0:t1);
dy   = ly-hptrend(ly,1600);
dz   = lz-hptrend(lz,1600);
dh   = lh-hptrend(lh,1600);

plot(t(t0:t1),dz,t(t0:t1),dy,t(t0:t1),dh)
title('HP-filtered GR Technology, Per-capita Output, and Per-capita Hours')

disp('Correlations of GR Technology, Output, and Hours:')
disp(corrcoef([dz,dy,dh]))

