% MAKEFIGS    produces Figures 1 through 6 in
%
%                 "Why Are Married Women Working So Much?"
%                  Staff Report 317
%                  by Larry Jones, Rody Manuelli, and Ellen McGrattan
%

iprint = 1;     % Results saved in eps format if iprint =1;
load hours.dat
load blau.dat
load nipaq2.dat
load qshare.dat
load figqh.dat
load figqw.dat
load figqs.dat
load figsh.dat
load figsw.dat
load figss.dat
load figdh.dat
load figdw.dat
load figds.dat
load figdhc.dat
load figdwp.dat
load figdap.dat

%
% Figure 1. U.S. Hours and Wage Ratios, 1950-2010
%  A. Market Hours of Work Per Week

figure(1)
plot(hours(:,1),hours(:,2),'LineWidth',2,'Color',[1,0,0], ...
     'Marker','s','MarkerFaceColor',[1,0,0])
line(hours(:,1),hours(:,3),'LineWidth',2,'Color',[0,0,0], ....
     'LineStyle','-.' ,'Marker','s','MarkerFaceColor',[0,0,0])
line(hours(:,1),hours(:,4),'LineWidth',2,'Color',[0,1,0], ...
     'LineStyle','--' ,'Marker','s','MarkerFaceColor',[0,1,0])
line(hours(:,1),hours(:,5),'LineWidth',2,'Color',[0,0,1], ...
     'LineStyle',':' ,'Marker','s','MarkerFaceColor',[0,0,1])
set(gca,'XTick',[1950:10:2010],'FontSize',14)
axis([1950,2010,0,50]);
text(1955,43.5,'Married men',  'Fontsize',14)
text(1955, 6.0,'Married women','Fontsize',14)
text(1955,33.0,'Single men',   'Fontsize',14)
text(1955,25.0,'Single women', 'Fontsize',14)
if (iprint==1);
  print -depsc fig1a
end;

%
% Figure 1. U.S. Hours and Wage Ratios
%  B. Ratio of Women's to Men's Wage Rates (%)
%

figure(2)
plot(blau(:,1),blau(:,2),'LineWidth',2,'Color',[1,0,0], ...
     'Marker','s','MarkerFaceColor',[1,0,0])
set(gca,'XTick',[1950:10:2010],'FontSize',14)
axis([1950,2010,40,80]);
if (iprint==1);
  print -depsc fig1b
end;

%
% Figure 2. U.S. Deflators and Expenditure Shares, 1929-2010
%  A. Ratio of Home Investment Deflators to GDP Deflator (1950==100)

figure(3)
plot(nipaq2(:,1),nipaq2(:,4),'LineWidth',1.5,'Color',[1,0,0], ...
     'MarkerFaceColor',[1,0,0])
line(nipaq2(:,1),nipaq2(:,5),'LineWidth',1.5,'Color',[0,0,0], ....
     'LineStyle','--' ,'MarkerFaceColor',[0,0,0])
set(gca,'XTick',[1930:10:2010],'FontSize',14)
axis([1929,2010,0,130]);
text(1959,100.5,'Durable Consumption and','Fontsize',14)
text(1959, 95,'Residential Investment','Fontsize',14)
text(1942, 50,'Household Appliances','Fontsize',14)
if (iprint==1);
  print -depsc fig2a
end;


%
% Figure 2. U.S. Deflators and Expenditure Shares, 1929-2010
%  B. Home Investment Expenditure Shares (1950=100)

figure(4)
plot(nipaq2(:,1),nipaq2(:,2),'LineWidth',1.5,'Color',[1,0,0], ...
     'MarkerFaceColor',[1,0,0])
line(nipaq2(:,1),nipaq2(:,3),'LineWidth',1.5,'Color',[0,0,0], ....
     'LineStyle','--' ,'MarkerFaceColor',[0,0,0])
set(gca,'XTick',[1930:10:2010],'FontSize',14)
axis([1929,2010,0,120]);
text(1958, 95.5,'Durable Consumption and','Fontsize',14)
text(1958, 90,'Residential Investment','Fontsize',14)
text(1950, 46,'Household Appliances','Fontsize',14)
if (iprint==1);
  print -depsc fig2b
end;


%
% Figure 3. Predictions of Model with Falling Home Capital Price
%  A. Market Hours of Work Per Week

figure(5)
plot(figqh(:,1),figqh(:,2),'LineWidth',2,'Color',[1,0,0])
line(figqh(:,1),figqh(:,3),'LineWidth',2,'Color',[0,0,0],'LineStyle','-.');
line(figqh(:,1),figqh(:,4),'LineWidth',2,'Color',[0,1,0],'LineStyle','--');
line(figqh(:,1),figqh(:,5),'LineWidth',2,'Color',[0,0,1],'LineStyle',':');
hold
scatter(hours(:,1),hours(:,2),30,[1,0,0],'Marker','s','Markerfacecolor',[1,0,0])
scatter(hours(:,1),hours(:,3),30,[0,0,0],'Marker','s','Markerfacecolor',[0,0,0])
scatter(hours(:,1),hours(:,4),30,[0,1,0],'Marker','s','Markerfacecolor',[0,1,0])
scatter(hours(:,1),hours(:,5),30,[0,0,1],'Marker','s','Markerfacecolor',[0,0,1])
hold
axis([1950,2010,0,50]);
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
text(1955,43.8,'Married men',  'Fontsize',14)
text(1955, 9.0,'Married women','Fontsize',14)
text(1955,36.0,'Single men',   'Fontsize',14)
text(1955,26.0,'Single women', 'Fontsize',14)
if (iprint==1);
  print -depsc fig3a
end;

%
% Figure 3. Predictions of Model with Falling Home Capital Price
%  B. Ratio of Women's to Men's Wage Rates (%)
%

figure(6)
plot(figqw(:,1),figqw(:,2),'LineWidth',2,'Color',[1,0,0]);
hold
scatter(blau(:,1),blau(:,2),30,[1,0,0],'Marker','s','Markerfacecolor',[1,0,0])
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
axis([1950,2010,40,80]);
hold
if (iprint==1);
  print -depsc fig3b
end;

%
% Figure 3. Predictions of Model with Falling Home Capital Price
%  C. Ratio of Home Investment to Market Output (1950=100)
%

figure(7)
plot(figqs(:,1),figqs(:,2),'LineWidth',2,'Color',[1,0,0]);
hold
scatter(qshare(:,1),qshare(:,2),30,[1,0,0],'Marker','s','Markerfacecolor',[1,0,0])
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
axis([1950,2010,30,140]);
hold
if (iprint==1);
  print -depsc fig3c
end;



%
% Figure 4. Predictions of Model with Inferior Home Goods
%  A. Market Hours of Work Per Week

figure(8)
plot(figsh(:,1),figsh(:,2),'LineWidth',2,'Color',[1,0,0])
line(figsh(:,1),figsh(:,3),'LineWidth',2,'Color',[0,0,0],'LineStyle','-.');
line(figsh(:,1),figsh(:,4),'LineWidth',2,'Color',[0,1,0],'LineStyle','--');
line(figsh(:,1),figsh(:,5),'LineWidth',2,'Color',[0,0,1],'LineStyle',':');
hold
scatter(hours(:,1),hours(:,2),30,[1,0,0],'Marker','s','Markerfacecolor',[1,0,0])
scatter(hours(:,1),hours(:,3),30,[0,0,0],'Marker','s','Markerfacecolor',[0,0,0])
scatter(hours(:,1),hours(:,4),30,[0,1,0],'Marker','s','Markerfacecolor',[0,1,0])
scatter(hours(:,1),hours(:,5),30,[0,0,1],'Marker','s','Markerfacecolor',[0,0,1])
hold
axis([1950,2010,0,50]);
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
text(1955,43.8,'Married men',  'Fontsize',14)
text(1955, 6.0,'Married women','Fontsize',14)
text(1955,36.0,'Single men',   'Fontsize',14)
text(1955,26.0,'Single women', 'Fontsize',14)
if (iprint==1);
  print -depsc fig4a
end;

%
% Figure 4. Predictions of Model with Inferior Home Goods
%  B. Ratio of Women's to Men's Wage Rates (%)
%

figure(9)
plot(figsw(:,1),figsw(:,2),'LineWidth',2,'Color',[1,0,0]);
hold
scatter(blau(:,1),blau(:,2),30,[1,0,0],'Marker','s','Markerfacecolor',[1,0,0])
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
axis([1950,2010,40,80]);
hold
if (iprint==1);
  print -depsc fig4b
end;

%
% Figure 4. Predictions of Model with Inferior Home Goods
%  C. Ratio of Home Investment to Market Output (1950=100)
%

figure(10)
plot(figss(:,1),figss(:,2),'LineWidth',2,'Color',[1,0,0]);
hold
scatter(qshare(:,1),qshare(:,2),30,[1,0,0],'Marker','s','Markerfacecolor',[1,0,0])
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
axis([1950,2010,30,140]);
hold
if (iprint==1);
  print -depsc fig4c
end;


%
% Figure 5. Predictions of Model with Changes in Discrimination
%  A. Market Hours of Work Per Week

figure(11)
plot(figdh(:,1),figdh(:,2),'LineWidth',2,'Color',[1,0,0])
line(figdh(:,1),figdh(:,3),'LineWidth',2,'Color',[0,0,0],'LineStyle','-.');
line(figdh(:,1),figdh(:,4),'LineWidth',2,'Color',[0,1,0],'LineStyle','--');
line(figdh(:,1),figdh(:,5),'LineWidth',2,'Color',[0,0,1],'LineStyle',':');
hold
scatter(hours(:,1),hours(:,2),30,[1,0,0],'Marker','s','Markerfacecolor',[1,0,0])
scatter(hours(:,1),hours(:,3),30,[0,0,0],'Marker','s','Markerfacecolor',[0,0,0])
scatter(hours(:,1),hours(:,4),30,[0,1,0],'Marker','s','Markerfacecolor',[0,1,0])
scatter(hours(:,1),hours(:,5),30,[0,0,1],'Marker','s','Markerfacecolor',[0,0,1])
hold
axis([1950,2010,0,50]);
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
text(1955,43.8,'Married men',  'Fontsize',14)
text(1955, 9.0,'Married women','Fontsize',14)
text(1955,36.0,'Single men',   'Fontsize',14)
text(1955,25.0,'Single women', 'Fontsize',14)
if (iprint==1);
  print -depsc fig5a
end;

%
% Figure 5. Predictions of Model with Changes in Discrimination
%  B. Ratio of Women's to Men's Wage Rates (%)
%

figure(12)
plot(figdw(:,1),figdw(:,2),'LineWidth',2,'Color',[1,0,0]);
hold
scatter(blau(:,1),blau(:,2),30,[1,0,0],'Marker','s','Markerfacecolor',[1,0,0])
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
axis([1950,2010,40,80]);
hold
if (iprint==1);
  print -depsc fig5b
end;

%
% Figure 5. Predictions of Model with Changes in Discrimination
%  C. Ratio of Home Investment to Market Output (1950=100)
%

figure(13)
plot(figds(:,1),figds(:,2),'LineWidth',2,'Color',[1,0,0]);
hold
scatter(qshare(:,1),qshare(:,2),30,[1,0,0],'Marker','s','Markerfacecolor',[1,0,0])
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
axis([1950,2010,30,140]);
hold
if (iprint==1);
  print -depsc fig5c
end;


%
% Figure 6. Auxiliary Predictions of Model with Discrimination
%  A. Human Capital Relative to Married Men in 1950

figure(14)
plot(figdhc(:,1),figdhc(:,2),'LineWidth',2,'Color',[1,0,0])
line(figdhc(:,1),figdhc(:,3),'LineWidth',2,'Color',[0,0,0],'LineStyle','-.');
line(figdhc(:,1),figdhc(:,4),'LineWidth',2,'Color',[0,1,0],'LineStyle','--');
line(figdhc(:,1),figdhc(:,5),'LineWidth',2,'Color',[0,0,1],'LineStyle',':');
axis([1950,2010,.2,1.1]);
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
text(1955,1.02,'Married men',  'Fontsize',14)
text(1955,.35, 'Married women','Fontsize',14)
text(1955,.90, 'Single men',   'Fontsize',14)
text(1955,.50, 'Single women', 'Fontsize',14)
if (iprint==1);
  print -depsc fig6a
end;

%
% Figure 6. Auxiliary Predictions of Model with Discrimination
%  B. Ratio of Single to Married Wage Rate by Sex (%)
%

figure(15)
plot(figdwp(:,1),figdwp(:,2),'LineWidth',2,'Color',[1,0,0]);
line(figdwp(:,1),figdwp(:,3),'LineWidth',2,'Color',[0,0,0],'LineStyle','-.');
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
axis([1950,2010,90,120]);
text(1955,112,'Women','Fontsize',14)
text(1955, 97,'Men',  'Fontsize',14)
if (iprint==1);
  print -depsc fig6b
end;

%
% Figure 6. Auxiliary Predictions of Model with Discrimination
%  C. Average Market Productivity Relative to Trend (1950=100)
%

figure(16)
plot(figdap(:,1),figdap(:,2),'LineWidth',2,'Color',[1,0,0]);
set(gca,'XTick',[1950:10:2010],'FontSize',14)
set(gcf,'PaperPosition',[.25,2.5,11,6])
axis([1950,2010,85,105]);
if (iprint==1);
  print -depsc fig6c
end;
