Scribd
Upload
7 views4 pages

Economic Growth Analysis with Lag

The document contains MATLAB code that analyzes economic data from 1970 to 2005 using two different time lags (10 years and 5 years). It plots the relationships between GDP, labor, capital, and technology over time, and compares the plots to steady state values. Multiple figures are generated to show growth contributions and deviations from steady state.

Uploaded by

Lipin Juan
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
7 views4 pages

Economic Growth Analysis with Lag

The document contains MATLAB code that analyzes economic data from 1970 to 2005 using two different time lags (10 years and 5 years). It plots the relationships between GDP, labor, capital, and technology over time, and compares the plots to steady state values. Multiple figures are generated to show growth contributions and deviations from steady state.

Uploaded by

Lipin Juan
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

11/17/10 4:25 PM C:\Documents and Settings\Lee\My Documents\My Dro...\graph1.

m 1 of 4

% =MLPutMatrix("Time", D1:AV1)

% =MLPutMatrix("Y10", D17:AV17)
% =MLPutMatrix("L10", D18:AV18)
% =MLPutMatrix("K10", D19:AV19)
% =MLPutMatrix("S10", D20:AV20)

% =MLPutMatrix("Y5", D40:AV40)
% =MLPutMatrix("L5", D41:AV41)
% =MLPutMatrix("K5", D42:AV42)
% =MLPutMatrix("S5", D43:AV43)

x10 = 0.000447428;
n10 = 0.003890305;
x5 = 0.000773;
n5 = 0.003890305;
sl = 0.61429;
sk = 0.38571;

time = Time;
i = find(time==1970):find(time==2005);
j = datenum(time',1,1);
j = j';
j = j(1,i);
y10 = Y10(1,i);
l10 = L10(1,i);
k10 = K10(1,i);
s10 = S10(1,i);

y5 = Y5(1,i);
l5 = L5(1,i);
k5 = K5(1,i);
s5 = S5(1,i);

y101 = y10*100;

l101 = sl*l10./y10*100;
k101 = sk*k10./y10*100;
s101 = sl*s10./y10*100;

y51 = y5*100;
l51 = sl*l5./y5*100;
k51 = sk*k5./y5*100;
s51 = sl*s5./y5*100;

disp('Steady State');
ya = 0.015238;
Csl10 = (sl*(n10))/ya*100
Csk10 = (sk*(x10+n10))/ya*100
Csa10 = (sl*(x10))/ya*100
devl10 = l101 - Csl10;
devk10 = k101 - Csk10;
deva10 = s101 - Csa10;

figure(1);
py10 = plot(j,y10,'-r*');
11/17/10 4:25 PM C:\Documents and Settings\Lee\My Documents\My Dro...\graph1.m 2 of 4

set(py10,'Color','Black','LineWidth',1,'MarkerSize',5);
hold on;
pl10 = plot(j,l10,'-rs');
set(pl10,'Color','Black','LineWidth',1,'MarkerSize',4);
pk10 = plot(j,k10,'-ro');
set(pk10,'Color','Black','LineWidth',1,'MarkerSize',5);
ps10 = plot(j,s10,'-r^');
set(ps10,'Color','Black','LineWidth',1,'MarkerSize',5);
datetick('x','yyyy');
% ylabel('Contribution to Growth');
title('Method II with 10-year lag');
legend('\DeltaGDP/DGP','\DeltaL/L','\DeltaK/K','\DeltaA/A');
axis([datenum(1970,1,1) datenum(2005,1,1) -0.07 0.12])
grid on;
hold off;

figure(2);
py10 = plot(j,y101,'-r*');
set(py10,'Color','Black','LineWidth',1,'MarkerSize',5);
hold on;
pl10 = plot(j,l101,'-rs');
set(pl10,'Color','Black','LineWidth',1,'MarkerSize',4);
pk10 = plot(j,k101,'-ro');
set(pk10,'Color','Black','LineWidth',1,'MarkerSize',5);
ps10 = plot(j,s101,'-r^');
set(ps10,'Color','Black','LineWidth',1,'MarkerSize',5);
datetick('x','yyyy');
% ylabel('Contribution to Growth');
title('Contribution to Growth');
legend('\DeltaGDP/DGP','\DeltaL/L','\DeltaK/K','\DeltaA/A');
% axis([datenum(1970,1,1) datenum(2005,1,1) -0.07 0.12])
grid on;
hold off;

figure(3);
py10 = plot(j,l101,'-r*');
hold on;
set(py10,'Color','Black','LineWidth',1,'MarkerSize',5);
refline(0,Csl10);
ylabel('percentage');
title('Deviation of Labor force')
datetick('x','yyyy');
hold off;

figure(4);
py10 = plot(j,k101,'-r*');
hold on;
set(py10,'Color','Black','LineWidth',1,'MarkerSize',5);
refline(0,Csk10);
ylabel('percentage');
title('Deviation of capital')
datetick('x','yyyy');
hold off;

figure(5);
py10 = plot(j,s101,'-r*');
11/17/10 4:25 PM C:\Documents and Settings\Lee\My Documents\My Dro...\graph1.m 3 of 4

hold on;
set(py10,'Color','Black','LineWidth',1,'MarkerSize',5);
refline(0,Csa10);
ylabel('percentage');
title('Deviation of technology')
datetick('x','yyyy');
hold off;

figure(6);
pl10 = plot(j,devl10,'-rs');
hold on;
set(pl10,'Color','Black','LineWidth',1,'MarkerSize',4);
pk10 = plot(j,devk10,'-ro');
set(pk10,'Color','Black','LineWidth',1,'MarkerSize',5);
ps10 = plot(j,deva10,'-r^');
set(ps10,'Color','Black','LineWidth',1,'MarkerSize',5);
datetick('x','yyyy');
ylabel('Deviation in Contribution');
title('Deviation from Steady State (10-year lag)');
% legend('\DeltaGDP/DGP','\DeltaL/L','\DeltaK/K','\DeltaA/A');
axis([datenum(1970,1,1) datenum(2005,1,1) -150 150])
grid on;
hold off;

%% five year lag


figure(7);
y5 = Y5(1,i);
l5 = L5(1,i);
k5 = K5(1,i);
s5 = S5(1,i);
py5 = plot(j,y5,'-r*');
set(py5,'Color','Black','LineWidth',1,'MarkerSize',5);
hold on;
pl5 = plot(j,l5,'-rs');
set(pl5,'Color','Black','LineWidth',1,'MarkerSize',4);
pk5 = plot(j,k5,'-ro');
set(pk5,'Color','Black','LineWidth',1,'MarkerSize',5);
ps5 = plot(j,s5,'-r^');
set(ps5,'Color','Black','LineWidth',1,'MarkerSize',5);
datetick('x','yyyy');
% ylabel('Contribution to Growth');
title('method II with 5-year lag');
legend('\DeltaGDP/DGP','\DeltaL/L','\DeltaK/K','\DeltaA/A');
axis([datenum(1970,1,1) datenum(2005,1,1) -0.07 0.12])
grid on;
hold off;

figure(8);
py10 = plot(j,y51,'-r*');
set(py10,'Color','Black','LineWidth',1,'MarkerSize',5);
hold on;
pl10 = plot(j,l51,'-rs');
set(pl10,'Color','Black','LineWidth',1,'MarkerSize',4);
pk10 = plot(j,k51,'-ro');
set(pk10,'Color','Black','LineWidth',1,'MarkerSize',5);
11/17/10 4:25 PM C:\Documents and Settings\Lee\My Documents\My Dro...\graph1.m 4 of 4

ps10 = plot(j,s51,'-r^');
set(ps10,'Color','Black','LineWidth',1,'MarkerSize',5);
datetick('x','yyyy');
% ylabel('Contribution to Growth');
title('Contribution to Growth (5 year lag)');
legend('\DeltaGDP/DGP','\DeltaL/L','\DeltaK/K','\DeltaA/A');
% axis([datenum(1970,1,1) datenum(2005,1,1) -0.07 0.12])
grid on;
hold off;

You might also like