% Tetyana Dubovyk 
% June, 15, 2004
% Graphs that are included in the final version of
% McGrattan, Ellen and Richard Rogerson, Changes in Hours Worked: 1950-2000,
% QR summer 2004, Federal Reserve Bank of Minneapolis.

a=20:10:80;
ave99=[ ... 
  23.78  19.93  23.37  23.37  22.27  14.29   3.90
  22.70  20.94  24.33  24.47  22.53  11.99   2.99
  21.64  22.07  25.12  25.45  22.16   8.43   2.21
  20.79  22.96  26.00  26.31  22.60   6.95   1.45
  19.47  23.97  27.24  27.76  23.30   5.11   1.24
  18.44  25.48  27.50  28.40  21.06   5.18   1.55
  18.23  25.60  28.61  28.90  21.03   5.72   1.33
  15.50  27.00  30.98  32.15  22.53   5.32   1.29
  15.37  30.03  33.64  32.53  22.27   5.39   1.31
  17.62  32.18  32.62  33.90  22.67   5.45   1.23
  17.35  31.64  33.95  35.57  23.22   5.36   1.19
  16.33  33.19  34.94  36.71  23.46   5.38   1.16
];
figure( 1)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.2 .1 .67 .67],'FontSize',12); %[left bottom width height]
xlabel('Age (in Years)','FontSize',12)
ylabel('Hours','FontSize',12)
line(a,ave99(1,:),'LineWidth',2,'Color',[0,0,0],'LineStyle',':')
line(a,ave99(2,:),'LineWidth',2,'Color',[0,0,0])
line(a,ave99(3,:),'LineWidth',2,'Color',[0.5,0.5,0.5])
line(a,ave99(4,:),'LineWidth',2,'Color',[0,0,0],'LineStyle','--')
line(a,ave99(5,:),'LineWidth',1,'Color',[0,0,0])
line(a,ave99(6,:),'LineWidth',1,'Color',[0,0,0],'LineStyle','--')
line(a,ave99(7,:),'LineWidth',1,'Color',[0,0,1],'LineStyle','--')
line(a,ave99(8,:),'LineWidth',1,'Color',[0,0,1])
line(a,ave99(9,:),'LineWidth',2,'Color',[0,0,1],'LineStyle','--')
line(a,ave99(10,:),'LineWidth',2,'Color',[1,0,0],'LineStyle','--')
line(a,ave99(11,:),'LineWidth',2,'Color',[1,0,0])
line(a,ave99(12,:),'LineWidth',2,'Color',[1,0,0],'LineStyle',':')
legend(    '1866-75','1876-85', '1886-95','1896-1905','1906-15','1916-25',...
          '1926-35','1936-45', '1946-55','1956-65','1966-75','1976-85')
axis([15,85,0,50]);
grid

str(1) = {'Charts 1-3' };
str(2) = {'Possible Shifts in Hours Worked'};
str(3) = {'     '};
str(4) = {'Extrapolated Average Weekly Hours Worked per Person'};
str(5) = {'by Cohorts at Various Ages in the United States'};
str(6) = {'     '};
str(7) = {'Chart 1 Total Population'};
set(gcf,'currentaxes',h);
text(.2,.90,str,'FontSize',12);
set(gcf,'PaperOrientation','portrait');
set(gcf,  'PaperPositionMode','manual');
set(gcf,         'PaperUnits','inches');
set(gcf,'PaperPosition',[0.25 2.5 8 6]); % [left bottom width height]
print -dpsc fig_1
 
ave10=[ ... 
  37.73  41.75  44.05  42.12  39.06  25.10   7.39
  35.34  41.43  43.49  41.44  37.49  20.74   5.15
  32.95  41.03  42.87  40.88  35.35  13.32   3.65
  30.79  40.78  42.37  40.17  33.65  10.77   2.62
  28.19  40.47  41.82  39.40  32.78   7.65   2.23
  25.76  39.84  41.01  39.19  28.94   7.30   2.70
  24.32  40.01  40.86  38.04  27.06   7.69   2.38
  20.39  39.55  40.18  38.60  27.20   6.66   2.30
  18.46  37.93  39.83  37.19  25.33   6.33   2.33
  20.00  37.96  37.93  36.52  24.13   6.01   2.20
  19.01  36.39  36.96  36.01  23.16   5.45   2.14
  17.74  35.34  35.88  35.15  21.81   5.05   2.09
];
figure( 2)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.2 .1 .67 .67],'FontSize',12); %[left bottom width height]
xlabel('Age (in Years)','FontSize',12)
ylabel('Hours','FontSize',12)
line(a,ave10(1,:),'LineWidth',2,'Color',[0,0,0],'LineStyle',':')
line(a,ave10(2,:),'LineWidth',2,'Color',[0,0,0])
line(a,ave10(3,:),'LineWidth',2,'Color',[0.5,0.5,0.5])
line(a,ave10(4,:),'LineWidth',2,'Color',[0,0,0],'LineStyle','--')
line(a,ave10(5,:),'LineWidth',1,'Color',[0,0,0])
line(a,ave10(6,:),'LineWidth',1,'Color',[0,0,0],'LineStyle','--')
line(a,ave10(7,:),'LineWidth',1,'Color',[0,0,1],'LineStyle','--')
line(a,ave10(8,:),'LineWidth',1,'Color',[0,0,1])
line(a,ave10(9,:),'LineWidth',2,'Color',[0,0,1],'LineStyle','--')
line(a,ave10(10,:),'LineWidth',2,'Color',[1,0,0],'LineStyle','--')
line(a,ave10(11,:),'LineWidth',2,'Color',[1,0,0])
line(a,ave10(12,:),'LineWidth',2,'Color',[1,0,0],'LineStyle',':')
legend(    '1866-75','1876-85', '1886-95','1896-1905','1906-15','1916-25',...
          '1926-35','1936-45', '1946-55','1956-65','1966-75','1976-85')
axis([15,85,0,50]);
grid
str(1) = {'   ' };
str(2) = {'   '};
str(3) = {'     '};
str(4) = {'   '};
str(5) = {'   '};
str(6) = {'     '};
str(7) = {'Chart 2 Males'};
set(gcf,'currentaxes',h);
text(.2,.90,str,'FontSize',12);
set(gcf,'PaperOrientation','portrait');
set(gcf,  'PaperPositionMode','manual');
set(gcf,         'PaperUnits','inches');
set(gcf,'PaperPosition',[0.25 2.5 8 6]);
  print -dpsc fig_2
 
ave20=[ ... 
  10.42   0.00   2.78   3.76   4.38   4.04   1.00
  10.59   1.38   5.35   6.95   6.98   3.82   1.37
  10.79   3.97   7.69   9.82   8.77   4.23   1.27
  11.18   5.95  10.04  12.47  12.31   3.95   0.81
  11.08   8.25  13.07  16.49  14.78   3.17   0.70
  11.41  11.76  14.71  18.46  14.16   3.54   0.88
  12.36  11.90  17.10  20.41  15.67   4.09   0.76
  10.77  15.13  22.13  26.06  18.21   4.14   0.74
  12.39  22.29  27.62  28.00  19.35   4.47   0.75
  15.20  26.41  27.33  31.18  21.08   4.78   0.70
  15.63  26.76  30.74  34.77  22.89   5.00   0.69
  14.84  30.63  33.61  37.67  24.45   5.29   0.67
];
figure( 3)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.2 .1 .67 .67],'FontSize',12); %[left bottom width height]
xlabel('Age (in Years)','FontSize',12)
ylabel('Hours','FontSize',12)
line(a,ave20(1,:),'LineWidth',2,'Color',[0,0,0],'LineStyle',':')
line(a,ave20(2,:),'LineWidth',2,'Color',[0,0,0])
line(a,ave20(3,:),'LineWidth',2,'Color',[0.5,0.5,0.5])
line(a,ave20(4,:),'LineWidth',2,'Color',[0,0,0],'LineStyle','--')
line(a,ave20(5,:),'LineWidth',1,'Color',[0,0,0])
line(a,ave20(6,:),'LineWidth',1,'Color',[0,0,0],'LineStyle','--')
line(a,ave20(7,:),'LineWidth',1,'Color',[0,0,1],'LineStyle','--')
line(a,ave20(8,:),'LineWidth',1,'Color',[0,0,1])
line(a,ave20(9,:),'LineWidth',2,'Color',[0,0,1],'LineStyle','--')
line(a,ave20(10,:),'LineWidth',2,'Color',[1,0,0],'LineStyle','--')
line(a,ave20(11,:),'LineWidth',2,'Color',[1,0,0])
line(a,ave20(12,:),'LineWidth',2,'Color',[1,0,0],'LineStyle',':')
legend(    '1866-75','1876-85', '1886-95','1896-1905','1906-15','1916-25',...
          '1926-35','1936-45', '1946-55','1956-65','1966-75','1976-85')
axis([15,85,0,50]);
grid
str(1) = {'   ' };
str(2) = {'   '};
str(3) = {'     '};
str(4) = {'   '};
str(5) = {'   '};
str(6) = {'     '};
str(7) = {'Chart 3 Females'};
set(gcf,'currentaxes',h);
text(.2,.90,str,'FontSize',12);
set(gcf,'PaperOrientation','portrait');
set(gcf,  'PaperPositionMode','manual');
set(gcf,         'PaperUnits','inches');
set(gcf,'PaperPosition',[0.25 2.5 8 6]);
  print -dpsc fig_3


% Chart 4 in update of QRw98
load goldin.dat
figure(4)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.2 .1 .67 .67],'FontSize',12); %[left bottom width height]
set(gca,'XTick',[1820 1840 1860 1880 1900 1920 1940 1960 1980 2000]);
set(gca,'YTick',[.3 .4 .5 .6 .7 .8]);
T = length(goldin);
line(goldin(1:6,1),goldin(1:6,2),'LineWidth',2,'Color',[1,0,0],...
    'LineStyle','-' ,'Marker','s','MarkerFaceColor',[1,0,0]);
text(1830,.55,'Manufacturing','FontSize',12);
line(goldin(7:T,1),goldin(7:T,2),'LineWidth',2,'Color',[0,0,1],...
    'LineStyle','-' ,'Marker','o','MarkerFaceColor',[0,0,1]);
text(1945,.68,'All Sectors','FontSize',12);
axis([1820,2000,0.3,0.8]);
%title('Chart 4. Ratio of Female to Male Earnings')
str(1) = {'Chart 4' };
str(2) = {'A Narrowing of the Gender Wage Gap'};
str(3) = {'     '};
str(4) = {'The Ratio of Female to Male Before-Tax Earnings'};
str(5) = {'for Full-Time Workers in Manufacturing Between'};
str(6) = {'1820 and 1930 and in All Sectors Between 1900 and 1992,'};
str(7) = {'in the United States'};
set(gcf,'currentaxes',h);
text(.2,.90,str,'FontSize',12);
str1(1) = {'Source: Goldin 1994'};
text(.2,.02,str1,'FontSize',12);
set(gcf,'PaperOrientation','portrait');
set(gcf,  'PaperPositionMode','manual');
set(gcf,         'PaperUnits','inches');
set(gcf,'PaperPosition',[0.25 2.5 8 6]);
print -dpsc fig_4


% Chart 5 in update of QRw98
load gordon.dat
load deflator.dat
p = gordon(:,2)./deflator(19:56,2);
figure(5)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.2 .1 .67 .67],'FontSize',12); %[left bottom width height]
set(gca,'XTick',[1945 1950 1955 1960 1965 1970 1975 1980 1985]);
set(gca,'YTick',[0 0.2 0.4 0.6 0.8 1]);
line(1947:1984, p/p(1),'LineWidth',1,'Color',[1,0,0],...
    'LineStyle','-' ,'Marker','s','MarkerFaceColor',[1,0,0]);
% title('Chart 5. Relative price of appliances')
axis([1945,1985,0,1]);
str(1) = {'Chart 5' };
str(2) = {'Innovating in Home Technology'};
str(3) = {'     '};
str(4) = {'The Price Index for Eight Home Appliances'};
str(5) = {'Relative to the GDP Deflator, 1947-84'};
str(6) = {'   '};
str(7) = {'   '};
set(gcf,'currentaxes',h);
text(.2,.90,str,'FontSize',12);
str1(1) = {'Source: Gordon 1990'};
text(.2,.02,str1,'FontSize',12);
set(gcf,'PaperOrientation','portrait');
set(gcf,  'PaperPositionMode','manual');
set(gcf,         'PaperUnits','inches');
set(gcf,'PaperPosition',[0.25 2.5 8 6]);
print -dpsc fig_5


% Chart 6 in update of QRw98
load benefits.dat
load cpi.dat
b = benefits(:,2)./ (cpi(28:89,2)/cpi(88,2));
figure(6)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.2 .1 .67 .67],'FontSize',12); %[left bottom width height]
set(gca,'XTick',[1940 1950 1960 1970 1980 1990 2000]);
set(gca,'YTick',[0 200 400 600 800 1000]);
line(1940:2001, b,'LineWidth',1,'Color',[1,0,0],...
    'LineStyle','-' ,'Marker','s','MarkerFaceColor',[1,0,0])
% title('Chart 6. Average monthly benefits for retired workers ($2000)')
axis([1940,2000,0,1000]);
str(1) = {'Chart 6' };
str(2) = {'Increasing Retirement Benefits'};
str(3) = {'     '};
str(4) = {'Average Monthly Social Security Benefits to U.S.'};
str(5) = {'Retired Workers in Constant 2000 Dollars, 1940-2000'};
str(6) = {'   '};
str(7) = {'   '};
set(gcf,'currentaxes',h);
text(.2,.90,str,'FontSize',12);
str1(1) = {'Source: U.S. Social Security Administration'};
text(.2,.02,str1,'FontSize',12);
set(gcf,'PaperOrientation','portrait');
set(gcf,  'PaperPositionMode','manual');
set(gcf,         'PaperUnits','inches');
set(gcf,'PaperPosition',[0.25 2.5 8 6]);
print -dpsc fig_6
  

yr1=1950:10:1990;
tab1=[ ... 
             22.3395
             21.5478
             21.1541
             22.0734
             23.8576
];
 
yr2=1980:10:2000;
tab2=[ ... 
             21.6710
             23.4770
             23.4990
];
 
figure(7)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.2 .1 .67 .67],'FontSize',12); %[left bottom width height]
set(gca,'XTick',[1950 1960 1970 1980  1990 2000])
xlabel( 'year', 'FontSize',12)
ylabel( 'average hours', 'FontSize',12)
line(yr1,tab1,'LineWidth',2,'Color',[0,0,1]           ,'Marker','s','MarkerFaceColor',[0,0,1])
line(yr2,tab2,'LineWidth',2,'Color',[1,0,0]           ,'Marker','s','MarkerFaceColor',[1,0,0])
legend(                               'based on hrswork2','based on uhrswork',2)
axis([1950,2000,20,25]);
grid
str(1) = {'Charts B1-B3' };
str(2) = {'Two Measures of Average Hours of Work per Person'};
str(3) = {'     '};
str(4) = {'   '};
str(5) = {'   '};
str(6) = {'     '};
str(7) = {'Chart B1 Total Population'};
set(gcf,'currentaxes',h);
text(.2,.90,str,'FontSize',12);
set(gcf,'PaperOrientation','portrait');
set(gcf,  'PaperPositionMode','manual');
set(gcf,         'PaperUnits','inches');
set(gcf,'PaperPosition',[0.25 2.5 8 6]);
  print -dpsc fig_7
 
yr1=1950:10:1990;
tab1=[ ... 
             34.1804
             31.9268
             29.7506
             28.7002
             29.1089
];
 
yr2=1980:10:2000;
tab2=[ ... 
             28.1930
             28.6220
             27.7840
];
 
figure(8)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.2 .1 .67 .67],'FontSize',12); %[left bottom width height]
set(gca,'XTick',[1950 1960 1970 1980  1990 2000])
xlabel( 'year', 'FontSize',12)
ylabel( 'average hours', 'FontSize',12)
line(yr1,tab1,'LineWidth',2,'Color',[0,0,1]           ,'Marker','s','MarkerFaceColor',[0,0,1])
line(yr2,tab2,'LineWidth',2,'Color',[1,0,0]           ,'Marker','s','MarkerFaceColor',[1,0,0])
legend(                               'based on hrswork2','based on uhrswork',1)
axis([1950,2000,27,35]);
grid
str(1) = {'   ' };
str(2) = {'   '};
str(3) = {'     '};
str(4) = {'   '};
str(5) = {'   '};
str(6) = {'     '};
str(7) = {'Chart B2 Males'};
set(gcf,'currentaxes',h);
text(.2,.90,str,'FontSize',12);
set(gcf,'PaperOrientation','portrait');
set(gcf,  'PaperPositionMode','manual');
set(gcf,         'PaperUnits','inches');
set(gcf,'PaperPosition',[0.25 2.5 8 6]);
  print -dpsc fig_8
 
yr1=1950:10:1990;
tab1=[ ... 
             10.8700
             11.8416
             13.3776
             16.0203
             19.0291
];
 
yr2=1980:10:2000;
tab2=[ ... 
             15.8160
             18.7720
             19.5090
];
 
figure(9)
h = axes('Position',[0 0 1 1],'Visible','off');
axes('Position',[.2 .1 .67 .67],'FontSize',12); %[left bottom width height]
set(gca,'XTick',[1950 1960 1970 1980  1990 2000])
xlabel( 'year', 'FontSize',12)
ylabel( 'average hours', 'FontSize',12)
line(yr1,tab1,'LineWidth',2,'Color',[0,0,1]           ,'Marker','s','MarkerFaceColor',[0,0,1])
line(yr2,tab2,'LineWidth',2,'Color',[1,0,0]           ,'Marker','s','MarkerFaceColor',[1,0,0])
legend(                               'based on hrswork2','based on uhrswork',2)
axis([1950,2000,10,20]);
grid
str(1) = {'   ' };
str(2) = {'   '};
str(3) = {'     '};
str(4) = {'   '};
str(5) = {'   '};
str(6) = {'     '};
str(7) = {'Chart B3 Females'};
set(gcf,'currentaxes',h);
text(.2,.90,str,'FontSize',12);
set(gcf,'PaperOrientation','portrait');
set(gcf,  'PaperPositionMode','manual');
set(gcf,         'PaperUnits','inches');
set(gcf,'PaperPosition',[0.25 2.5 8 6]);
  print -dpsc fig_9
 
