options crt;  ? 11e1
in mroz;

title 'Exercise 1, part (a) - full sample means + standard deviations';
smpl 1,753;
list vars lfp whrs kl6 k618 wa we ww rpwg hhrs ha he hw faminc mtr
          wmed wfed un cit ax;
msd(terse) vars;

title 'Exercise 1, part (b) - working women:  means + standard deviations';
select lfp;
msd(terse) vars;

title 'Exercise 1, part (c) - property income';
select 1;
prin = faminc - whrs*ww;
msd(terse) prin;

title 'Exercise 1, part (d) - (predicted) log wage';
ax2 = ax*ax;
wa2 = wa*wa;
lww1 = 0;
select lfp;
  lww = log(ww);
  lww1 = lww;
  msd(terse) lww;
  olsq lww c wa we cit ax wa2;
select lfp=0;
  forcst flww;
  lww1 = flww;
  msd(terse) flww;
select 1;
msd(terse) lww1;

title 'Save new variables in databank';
out mroz;
keep prin lww1 ax2 wa2;