%
% This file displays Figures 1-5 in
%
%      McGrattan and Rogerson
%      Changes in the Distribution of Family Hours Worked Since 1950
%

%      Tetyana Dubovyk, 4-5-07
%      Revised, ERM, 8-30-07
%      Revised, Tetyana Dubovyk, 9-16-07

iprint   = input('Would you like to print the figures? (1=yes) ');

%
% Figures 1-2B. Aggregate Hours Per Person, 1950--2005
%
yr1      = 1950:10:1990;
yr2      = [1980:10:2000,2005];
yr       = [1950:10:2000,2005];
fig1hrs2 = [25.2731
            25.4848
            26.3655
            27.9675
            30.5731];
fig1uhrs = [27.3031
            29.8749
            29.7818
            30.2901];
fig1     = [fig1hrs2;fig1uhrs(3:4)*fig1hrs2(5)/fig1uhrs(2)];

figure(1)
plot(yr,fig1,'LineWidth',2,'Color',[1,0,0], ...
     'Marker','s','MarkerFaceColor',[1,0,0])
title('Figure 1. Aggregate Hours Per Person, Ages 25-64',  ...
      'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Average Weekly Hours')
grid
axis([1950,2005,0,45]);
if iprint==1; print -depsc figure1; end;
 
figure(2)
plot(yr1,fig1hrs2,'LineWidth',2,'Color',[1,0,0], ...
     'Marker','s','MarkerFaceColor',[1,0,0])
line(yr2,fig1uhrs,'LineWidth',2,'Color',[1,0,0], ...
     'LineStyle','--' ,'Marker','o','MarkerFaceColor',[1,0,0])
title(['Appendix Figure 1. Aggregate Hours Per Person, Ages 25-64'], ...
       'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Average Weekly Hours')
legend('Based on HRSWORK2', ...
       'Based on UHRSWRK','Location','Southeast')
grid
axis([1950,2005,0,45]);

 
fig2Ahrs2 = [24.8428             31.4252             28.8010
             25.4954             29.1782             28.0608
             26.6718             28.0432             27.0204
             28.3452             29.9582             28.2560
             31.4300             31.1863             29.5931];
fig2Auhrs = [27.7245             29.2301             27.4554
             30.7583             30.4547             28.6576
             31.0708             30.1850             28.3644
             30.9783             31.5366             28.2482];
fig2A     = [fig2Ahrs2;
             fig2Auhrs(3:4,:).*(ones(2,1)*(fig2Ahrs2(5,:)./fig2Auhrs(2,:)))];


figure(3)
plot(yr,fig2A(:,1),'LineWidth',2,'Color',[1,0,0], ...
     'Marker','s','MarkerFaceColor',[1,0,0])
line(yr,fig2A(:,2),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','o','MarkerFaceColor',[0,1,0])
line(yr,fig2A(:,3),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','d','MarkerFaceColor',[0,0,1])
title('Figure 2A. Average Weekly Hours Per Person of Families and Singles', ...
      'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Average Weekly Hours')
legend('Families', ...
       'Single men', ...
       'Single women','Location','Southeast');
grid
axis([1950,2005,0,45]);
if iprint==1; print -depsc figure2A; end;

figure(4)
plot(yr1,fig2Ahrs2(:,1),'LineWidth',2,'Color',[1,0,0], ...
     'Marker','s','MarkerFaceColor',[1,0,0])
line(yr2,fig2Auhrs(:,1),'LineWidth',2,'Color',[1,0,0], ...
     'LineStyle','--','Marker','s','MarkerFaceColor',[1,0,0])
line(yr1,fig2Ahrs2(:,2),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','o','MarkerFaceColor',[0,1,0])
line(yr2,fig2Auhrs(:,2),'LineWidth',2,'Color',[0,1,0], ...
     'LineStyle','--','Marker','o','MarkerFaceColor',[0,1,0])
line(yr1,fig2Ahrs2(:,3),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','d','MarkerFaceColor',[0,0,1])
line(yr2,fig2Auhrs(:,3),'LineWidth',2,'Color',[0,0,1], ...
     'LineStyle','--','Marker','d','MarkerFaceColor',[0,0,1])
title(['Appendix Figure 2A. Average Weekly Hours Per Person of Families ', ...
       'and Singles'],'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Average Weekly Hours')
legend('Families, Based on HRSWORK2', ...
       'Families, Based on UHRSWORK', ...
       'Single men, Based on HRSWORK2', ...
       'Single men, Based on UHRSWORK', ...
       'Single women, Based on HRSWORK2', ...
       'Single women, Based on UHRSWORK','Location','Southeast');
grid
axis([1950,2005,0,45]);

 
fig2Bhrs2 = [41.6833              8.0024
             40.1376             10.8532
             39.5603             13.7832
             38.6713             18.0191
             39.5396             23.3182];
fig2Buhrs = [38.0975             17.3514
             38.8413             22.6766
             38.0883             24.0529
             37.9264             24.0581];
fig2B     = [fig2Bhrs2;
             fig2Buhrs(3:4,:).*(ones(2,1)*(fig2Bhrs2(5,:)./fig2Buhrs(2,:)))];


figure(5)
plot(yr,fig2B(:,1),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','o','MarkerFaceColor',[0,1,0])
line(yr,fig2B(:,2),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','d','MarkerFaceColor',[0,0,1])
title(['Figure 2B. Average Weekly Hours Per Person of Married ', ...
       'Men and Women'],'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Average Weekly Hours')
grid
legend('Married men', ...
       'Married women','Location','Southeast')
axis([1950,2005,0,45]);
if iprint==1; print -depsc figure2B; end;
 
figure(6)
plot(yr1,fig2Bhrs2(:,1),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','o','MarkerFaceColor',[0,1,0])
line(yr2,fig2Buhrs(:,1),'LineWidth',2,'Color',[0,1,0], ...
     'LineStyle','--','Marker','o','MarkerFaceColor',[0,1,0])
line(yr1,fig2Bhrs2(:,2),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','d','MarkerFaceColor',[0,0,1])
line(yr2,fig2Buhrs(:,2),'LineWidth',2,'Color',[0,0,1], ...
     'LineStyle','--','Marker','d','MarkerFaceColor',[0,0,1])
title(['Appendix Figure 2B. Average Weekly Hours Per Person of Married ', ...
       'Men and Women'],'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Average Weekly Hours')
legend('Married men, Based on HRSWORK2', ...
       'Married men: Based on UHRSWORK', ...
       'Married women: Based on HRSWORK2', ...
       'Married women: Based on UHRSWORK','Location','Southeast')
grid
axis([1950,2005,0,45]);


%
% Figures 3A-3B. Dispersion in hours, 1950--2005
%
 
fig3Ahrs2 = [0.4486              0.4841              0.4753
             0.4706              0.5155              0.4835
             0.4604              0.5172              0.4944
             0.4722              0.4851              0.4762
             0.4485              0.4695              0.4723];
fig3Auhrs = [0.4896              0.5063              0.4926
             0.4634              0.4888              0.4909
             0.4804              0.5031              0.5061
             0.4678              0.4614              0.4983];
fig3A     = [fig3Ahrs2;
             fig3Auhrs(3:4,:).*(ones(2,1)*(fig3Ahrs2(5,:)./fig3Auhrs(2,:)))];

figure(7)
plot(yr,fig3A(:,1),'LineWidth',2,'Color',[1,0,0], ...
     'Marker','s','MarkerFaceColor',[1,0,0])
line(yr,fig3A(:,2),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','o','MarkerFaceColor',[0,1,0])
line(yr,fig3A(:,3),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','d','MarkerFaceColor',[0,0,1])
title('Figure 3A. Dispersion in Hours for Families and Singles','FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Standard Deviation of Hours/ Mean of Hours')
legend('Families', ...
       'Single men', ...
       'Single women','Location','Northeast');
grid
axis([1950,2005,.3,.8]);
if iprint==1; print -depsc figure3A; end;

figure(8)
plot(yr1,fig3Ahrs2(:,1),'LineWidth',2,'Color',[1,0,0], ... 
     'Marker','s','MarkerFaceColor',[1,0,0])
line(yr2,fig3Auhrs(:,1),'LineWidth',2,'Color',[1,0,0], ...
     'LineStyle','--','Marker','s','MarkerFaceColor',[1,0,0])
line(yr1,fig3Ahrs2(:,2),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','o','MarkerFaceColor',[0,1,0])
line(yr2,fig3Auhrs(:,2),'LineWidth',2,'Color',[0,1,0], ...
     'LineStyle','--','Marker','o','MarkerFaceColor',[0,1,0])
line(yr1,fig3Ahrs2(:,3),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','d','MarkerFaceColor',[0,0,1])
line(yr2,fig3Auhrs(:,3),'LineWidth',2,'Color',[0,0,1], ...
     'LineStyle','--','Marker','d','MarkerFaceColor',[0,0,1])
title('Appendix Figure 3A. Dispersion in Hours for Families and Singles', ...
      'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Standard Deviation of Hours/ Mean of Hours')
legend('Families, Based on HRSWORK2', ...
       'Families, Based on UHRSWORK', ...
       'Single men, Based on HRSWORK2', ...
       'Single men, Based on UHRSWORK', ...
       'Single women, Based on HRSWORK2', ...
       'Single women, Based on UHRSWORK','Location','Northeast');
grid
axis([1950,2005,.3,.8]);

 
fig3Bhrs2 = [0.3713              2.0324
             0.4004              1.6250
             0.4015              1.3366
             0.4435              1.0703
             0.4506              0.8490];
fig3Buhrs = [0.4660              1.1102
             0.4683              0.8789
             0.5098              0.8365
             0.5046              0.8339];
fig3B     = [fig3Bhrs2;
             fig3Buhrs(3:4,:).*(ones(2,1)*(fig3Bhrs2(5,:)./fig3Buhrs(2,:)))];
 
figure(9)
plot(yr,fig3B(:,1),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','o','MarkerFaceColor',[0,1,0])
line(yr,fig3B(:,2),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','d','MarkerFaceColor',[0,0,1])
title('Figure 3B. Dispersion in Hours for Married Men and Women','FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Standard Deviation of Hours/ Mean of Hours')
legend('Married men', ...
       'Married women','Location','Northeast')
grid
axis([1950,2005,.3,2.2]);
if iprint==1; print -depsc figure3B; end;
 
figure(10)
plot(yr1,fig3Bhrs2(:,1),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','o','MarkerFaceColor',[0,1,0])
line(yr2,fig3Buhrs(:,1),'LineWidth',2,'Color',[0,1,0], ...
     'LineStyle','--' ,'Marker','o','MarkerFaceColor',[0,1,0])
line(yr1,fig3Bhrs2(:,2),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','d','MarkerFaceColor',[0,0,1])
line(yr2,fig3Buhrs(:,2),'LineWidth',2,'Color',[0,0,1], ...
     'LineStyle','--','Marker','d','MarkerFaceColor',[0,0,1])
title('Appendix Figure 3B. Dispersion in Hours for Married Men and Women', ...
      'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Standard Deviation of Hours/ Mean of Hours')
legend('Married men, Based on HRSWORK2', ...
       'Married men: Based on UHRSWORK', ...
       'Married women: Based on HRSWORK2', ...
       'Married women: Based on UHRSWORK','Location','Northeast')
grid
axis([1950,2005,.3,2.2]);

%
% Figure 4.  Correlation of hours between couples, 1950--2005
%
 
fig4hrs2 = [-0.0144
             0.0115
             0.0195
             0.0763
             0.1210]*100;
fig4uhrs = [ 0.0743
             0.1165
             0.1395
             0.0926]*100;
fig4     = [fig4hrs2;fig4uhrs(3:4)*fig4hrs2(5)/fig4uhrs(2)];

figure(11)
plot(yr,fig4(:,1),'LineWidth',2,'Color',[1,0,0], ...
     'Marker','s','MarkerFaceColor',[1,0,0])
title('Figure 4. Correlation of Hours Between Husbands and Wives', ...
      'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Correlation of Hours (%)')
grid
axis([1950,2005,-100,100]);
if iprint==1; print -depsc figure4; end;
 
figure(12)
plot(yr1,fig4hrs2(:,1),'LineWidth',2,'Color',[1,0,0], ...
     'Marker','s','MarkerFaceColor',[1,0,0])
line(yr2,fig4uhrs(:,1),'LineWidth',2,'Color',[1,0,0], ...
     'LineStyle','--','Marker','s','MarkerFaceColor',[1,0,0])
title('Appendix Figure 4. Correlation of Hours Between Husbands and Wives', ...
      'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Correlation of Hours (%)')
legend('Based on HRSWORK2', ...
       'Based on UHRSWRK','Location','Southeast')
grid
axis([1950,2005,-100,100]);


%
% Figures 5A-5C. Hours and Employment Conditional on Husband's Wage
%

deciles  = 0:10;
fig5Ahrs = [12.0397     18.4989
            18.4346     22.9408
            20.5240     25.0089
            20.7368     25.8725
            19.5855     26.2959
            18.4826     26.3792
            17.5080     26.3672
            17.5486     26.0225
            16.5527     24.6857
            15.0379     23.6242
            12.6971     18.9563];

fig5Aemp = [34.0985     48.2091
            52.3355     60.2370
            57.7017     65.2594
            58.8993     68.1919
            56.8508     69.2540
            54.6173     70.3857
            51.7661     70.3407
            52.7344     69.8999
            50.2487     67.5509
            47.2819     65.2027
            41.6221     54.0741];

 
figure(13)
bar(deciles,fig5Ahrs)
title('Figure 5A. Wife''s Hours Conditional on Husband''s Wage', ...
      'FontSize',13)
set(gca,'XTick',0:10,'FontSize',13)
xlabel('Wage Decile of Husband (0=Not in Labor Force)')
ylabel('Average Weekly Hours of Wives')
legend('1980','2005',2)
grid
axis([-1,11,0,50]);
if iprint==1; print -depsc figure5A; end;
 
figure(14)
bar(deciles,fig5Aemp)
title(['Appendix Figure 5A. Wife''s Employment Conditional ', ...
       'on Husband''s Wage'],'FontSize',13)
set(gca,'XTick',0:10,'FontSize',13)
xlabel('Wage Decile of Husband (0=Not in Labor Force)')
ylabel('Employment Rate of Wives')
legend('1980','2005',2)
grid
axis([-1,11,0,100]);

fig5Bhrs = [ 0.0000      0.0000
            40.3919     39.7128
            41.7008     41.3048
            41.9847     41.5255
            41.2826     42.1092
            40.9988     42.0469
            41.2128     42.2020
            41.2714     41.8648
            40.6526     41.7818
            40.5827     41.0130
            36.6732     43.0595];
fig5Bemp = [0.0000      0.0000
            85.7941     87.9240
            92.4024     91.9397
            93.8585     92.0937
            94.8720     92.8612
            94.7310     93.6129
            94.8300     93.3448
            95.8100     93.1092
            95.4976     94.3212
            95.1933     93.4484
            92.9187     92.7777];

figure(15)
bar(deciles,fig5Bhrs)
title('Figure 5B. Husband''s Hours Conditional on Husband''s Wage', ...
      'FontSize',13)
set(gca,'XTick',0:10,'FontSize',13)
xlabel('Wage Decile of Husband (0=Not in Labor Force)')
ylabel('Average Weekly Hours of Husband')
legend('1980','2005',2)
grid
axis([-1,11,0,50]);
if iprint==1; print -depsc figure5B; end;
 
figure(16)
bar(deciles,fig5Bemp)
title(['Appendix Figure 5B. Husband''s Employment Conditional ', ...
       'on Husband''s Wage'],'FontSize',13)
set(gca,'XTick',0:10,'FontSize',13)
xlabel('Wage Decile of Husband (0=Not in Labor Force)')
ylabel('Employment Rate of Husbands')
legend('1980','2005',2)
grid
axis([-1,11,0,110]);
        
 
fig5Chrs = [ 6.0198      9.1695
            29.4132     31.5204
            31.1124     33.3269
            31.3608     33.8380
            30.4340     34.3347
            29.7407     34.3259
            29.3604     34.3827
            29.4100     34.0420
            28.6026     33.3368
            27.8103     32.4156
            24.6852     31.1401];
fig5Cemp = [17.0492     23.8963
            69.0648     74.4002
            75.0520     78.8780
            76.3789     80.3551
            75.8614     81.2549
            74.6742     82.1667
            73.2981     81.9853
            74.2722     81.6487
            72.8731     81.0974
            71.2376     79.4831
            67.2704     73.6382];
        
figure(17)
bar(deciles,fig5Chrs)
title('Figure 5C. Family''s Hours Conditional on Husband''s Wage', ...
      'FontSize',13)
set(gca,'XTick',0:10,'FontSize',13)
xlabel('Wage Decile of Husband (0=Not in Labor Force)')
ylabel('Average Weekly Hours of Each Family Member')
legend('1980','2005',2)
grid
axis([-1,11,0,50]);
if iprint==1; print -depsc figure5C; end;
 
figure(18)
bar(deciles,fig5Cemp)
title(['Appendix Figure 5C. Family''s Employment Conditional ', ...
       'on Husband''s Wage'],'FontSize',13)
set(gca,'XTick',0:10,'FontSize',13)
xlabel('Wage Decile of Husband (0=Not in Labor Force)')
ylabel('Employment Rate of Family Members')
legend('1980','2005',2)
grid
axis([-1,11,0,100]);
        
 
%
% Figures 6A-6C. Hours and Employment Conditional on Husband's Education
%

quartiles = 1:4;
str1(1)   = {'Group Legend:'};
str1(2)   = {'  1 = HS Dropout'};
str1(3)   = {'  2 = High School'};
str1(4)   = {'  3 = Some College'};
str1(5)   = {'  4 = College Plus'};
fig6Ahrs = [ 11.0872             15.9616             18.7828
             11.0546             18.2116             25.0056
             11.2778             19.6400             26.4687
              8.8639             18.8844             24.7954];
fig6Aemp  = [31.0380             44.5770             48.8217
             31.3277             51.6969             64.6840
             31.7658             55.5974             68.4792
             26.1533             55.0214             65.4733];
 
figure(19)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig6Ahrs)
title('Figure 6A. Wife''s Hours Conditional on Husband''s Education', ...
      'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Average Weekly Hours of Wives')
legend('1960','1980','2005','Location','Northwest')
grid
axis([0,5,0,50]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);
if iprint==1; print -depsc figure6A; end;

figure(20)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig6Aemp)
title(['Appendix Figure 6A. Wife''s Employment Rate Conditional on,' ...
       ' Husband''s Education'], ...
       'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Employment Rate of Wives')
legend('1960','1980','2005','Location','Northwest')
grid
axis([0,5,0,100]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);

fig6Bhrs = [ 38.4538             32.8175             31.6576
             41.9867             39.1616             36.7856
             41.6613             40.8126             39.0658
             43.1510             43.1503             42.0628];
fig6Bemp = [ 89.3192             77.2042             72.5050
             96.0357             88.9324             81.4950
             95.3471             91.4779             85.1931
             97.2595             95.0510             89.5538];
figure(21)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig6Bhrs)
title('Figure 6B. Husband''s Hours Conditional on Husband''s Education', ...
      'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Average Weekly Hours of Husbands')
legend('1960','1980','2005','Location','Northwest')
grid
axis([0,5,0,50]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);
if iprint==1; print -depsc figure6B; end;

figure(22)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig6Bemp)
title(['Appendix Figure 6B. Husband''s Employment Rate Conditional on', ...
       'Husband''s Education'],'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Employment Rate of Husbands')
legend('1960','1980','2005','Location','Northwest')
grid
axis([0,5,0,100]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);

fig6Chrs = [ 24.7705             24.3896             25.3981
             26.5207             28.6866             30.9686
             26.4696             30.2263             32.8476
             26.0075             31.0174             33.5306];
fig6Cemp = [ 60.1786             60.8906             60.9785
             63.6817             70.3147             73.2116
             63.5565             73.5377             76.9492
             61.7064             75.0362             77.6559];

figure(23)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig6Chrs)
title('Figure 6C. Family''s Hours Conditional on Husband''s Education', ...
      'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Average Weekly Hours of Families')
legend('1960','1980','2005','Location','Northwest')
grid
axis([0,5,0,50]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);
if iprint==1; print -depsc figure6C; end;

figure(24)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig6Cemp)
title(['Appendix Figure 6C. Family''s Employment Rate Conditional on', ...
       ' Husband''s Education'],'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Employment Rate of Families')
legend('1960','1980','2005','Location','Northwest')
grid
axis([0,5,0,100]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);

%
% Figure 7A-7B.
%
yr1=[1960:10:2000,2005];
fig7Ahrs = [11.0219  11.9466  42.4906  43.6502
            14.1959  14.6327  41.6807  41.8260
            18.3271  21.8588  39.5339  43.6836
            22.2581  27.1342  38.4352  44.4630
            23.0064  27.7607  35.9320  44.0980
            23.5284  26.5041  36.1631  42.6147];
fig7Bemp = [31.3488  34.9157  96.6793  97.2337
            41.1952  43.8557  95.8287  96.0296
            52.2376  62.4237  89.7065  95.6530
            61.7551  73.1546  86.6105  94.8679
            61.6767  72.1863  80.1661  92.7486
            61.5255  68.7477  80.3262  90.4592];
str2(1)  = {'Family Type with Both:'};
str2(2)  = {'  1 = HS Educated'};
str2(3)  = {'  2 = College+ Educated'};
 
figure(25)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
plot(yr1,fig7Ahrs(:,3),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','s','MarkerFaceColor',[0,1,0])
line(yr1,fig7Ahrs(:,4),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','o','MarkerFaceColor',[0,1,0])
line(yr1,fig7Ahrs(:,1),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','s','MarkerFaceColor',[0,0,1])
line(yr1,fig7Ahrs(:,2),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','o','MarkerFaceColor',[0,0,1])
title('Figure 7A. Average Weekly Hours for Two Family Types', 'FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Average Weekly Hours')
legend('Married men, Family type 1', ...
       'Married men, Family type 2', ...
       'Married women, Family type 1', ...
       'Married women, Family type 2','Location','Southeast')
grid
axis([1960,2005,0,45]);
set(gcf,'currentaxes',h);
text(0.01,.6,str2);
if iprint==1; print -depsc figure7A; end;
 
figure(26)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
plot(yr1,fig7Bemp(:,3),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','s','MarkerFaceColor',[0,1,0])
line(yr1,fig7Bemp(:,4),'LineWidth',2,'Color',[0,1,0], ...
     'Marker','o','MarkerFaceColor',[0,1,0])
line(yr1,fig7Bemp(:,1),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','s','MarkerFaceColor',[0,0,1])
line(yr1,fig7Bemp(:,2),'LineWidth',2,'Color',[0,0,1], ...
     'Marker','o','MarkerFaceColor',[0,0,1])
title('Figure 7B. Employment Rate for Two Family Types','FontSize',13)
set(gca,'XTick',[1950:10:2000,2005],'FontSize',13)
xlabel('Year')
ylabel('Employment Rate')
legend('Married men, Family type 1', ...
       'Married men, Family type 2', ...
       'Married women, Family type 1', ...
       'Married women, Family type 2','Location','Southeast')
grid
axis([1960,2005,0,100]);
set(gcf,'currentaxes',h);
text(0.01,.6,str2);
if iprint==1; print -depsc figure7B; end;


%
% Figures 8A-8D. Hours and Employment Conditional on Husband's Education
%                and on Wife's Education
%
fig8Ahrs  = [10.1125             13.4226             14.2746
              9.3755             13.7393             15.2191
              9.8032             14.5389             15.1301
              7.6989             13.0502             15.5987];
fig8Aemp  = [28.4703             37.6750             37.6537
             26.4332             39.1167             40.3845
             27.5123             41.0729             40.7074
             22.0548             37.5788             41.5021];
fig8Bhrs  = [12.6701             18.6612             21.1075
             11.0219             18.3271             23.5284
             10.4959             18.3698             23.1843
              7.3838             16.1378             20.6467];
fig8Bemp  = [35.1788             52.0529             55.1606
             31.3488             52.2376             61.5255
             29.4613             52.5696             61.0542
             21.6226             47.5372             56.3202];
fig8Chrs  = [15.9134             22.0876             25.4632
             13.5336             21.7479             27.8502
             11.8334             20.8605             26.6891
              7.2848             17.1292             22.3434];
fig8Cemp  = [43.9739             60.9922             64.9999
             38.3879             61.5402             72.1730
             33.5734             59.2218             69.8458
             22.2479             51.5370             61.6897];
fig8Dhrs  = [23.4535             25.4909             28.8698
             19.3251             26.0646             32.2483
             16.8618             25.7695             31.3897
             11.9466             21.8588             26.5041];
fig8Demp  = [63.4399             68.5112             70.3715
             53.9329             70.7740             80.0978
             47.5671             70.3940             78.5246
             34.9157             62.4237             68.7477];

         
figure(27)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig8Ahrs)
title(['Figure 8A. Hours of Wives with no High School Degree,';
       '            Conditional on Husband''s Education       '], ...
       'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Average Weekly Hours of Wives with no HS Degree')
legend('1960','1980','2005','Location','Northeast')
grid
axis([0,5,0,35]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);
if iprint==1; print -depsc figure8A; end;
 
figure(28)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig8Aemp)
title(['Appendix Figure 8A. Employment Rate of Wives with no High';
       '    School Degree, Conditional on Husband''s Education    '], ...
       'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Employment Rate of Wives with no HS Degree')
legend('1960','1980','2005','Location','Northeast')
grid
axis([0,5,0,100]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);

figure(29)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig8Bhrs)
title(['Figure 8B. Hours of Wives with High School Degree,';
       '          Conditional on Husband''s Education      '], ...
       'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband') 
ylabel('Average Weekly Hours of Wives with HS Degree')
legend('1960','1980','2005','Location','Northeast')
grid
axis([0,5,0,35]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);
if iprint==1; print -depsc figure8B; end;
 
 
figure(30)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig8Bemp)
title(['Appendix Figure 8B. Employment Rate of Wives with High';
       '   School Degree, Conditional on Husband''s Education  '], ...
       'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Employment Rate of Wives with HS Degree')
legend('1960','1980','2005','Location','Northeast')
grid
axis([0,5,0,100]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);

figure(31)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig8Chrs)
title(['Figure 8C. Hours of Wives with Some College Education,';
       '      Conditional on Husband''s Education              '], ...
       'FontSize',13) 
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Average Weekly Hours of Wives with Some College')
legend('1960','1980','2005','Location','Northeast')
grid
axis([0,5,0,35]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);
if iprint==1; print -depsc figure8C; end;
 
 
figure(32)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig8Cemp)
title(['Appendix Figure 8C. Employment Rate of Wives with Some College';
       '        Education, Conditional on Husband''s Education         '], ...
       'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Employment Rate of Wives with Some College')
legend('1960','1980','2005','Location','Northeast')
grid
axis([0,5,0,100]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);

figure(33)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig8Dhrs)
title(['Figure 8D. Hours of Wives with College+ Education,';
       '          Conditional on Husband''s Education      '], ...
       'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Average Weekly Hours of Wives with College+')
legend('1960','1980','2005','Location','Northeast')
grid
axis([0,5,0,35]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);
if iprint==1; print -depsc figure8D; end;
 
figure(34)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.3 .1 .68 .8]);
bar(quartiles,fig8Demp)
title(['Appendix Figure 8D. Employment Rate of Wives with College+';
       '       Education, Conditional on Husband''s Education      '], ...
       'FontSize',13)
set(gca,'XTick',1:4,'FontSize',13)
xlabel('Education Group of Husband')
ylabel('Employment Rate of Wives with College+')
legend('1960','1980','2005','Location','Northeast')
grid
axis([0,5,0,100]);
set(gcf,'currentaxes',h);
text(0.05,.6,str1);
 
 
