# delimit ;
set memory 5000;
*************************************************;
* Administrative Commands;
*************************************************;
set more off;
clear;
log using ch7_ee.log,replace;
*************************************************;
* Read in the Dataset;
*************************************************;
use cps92_98.dta;
desc;
**************************************************************;
* Analysis for 1998 Data;
*;
keep if year == 1998;
generate age2 = age*age;
generate lnahe = ln(ahe);
generate fem_bac=female*bachelor;
generate lnage = ln(age);
generate fem_age=female*age;
generate fem_age2=female*age2;
generate bac_age=bachelor*age;
generate bac_age2=bachelor*age2;
generate fem_bac_age = fem_bac*age;
generate fem_bac_age2 = fem_bac*age2;
regress lnahe age female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
regress lnahe age age2 female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
test age age2;
regress lnahe lnage female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
regress lnahe age age2 fem_age fem_age2 female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
test age age2 fem_age fem_age2;
test fem_age fem_age2;
regress lnahe age age2 bac_age bac_age2 female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
test age age2 bac_age bac_age2;
test bac_age bac_age2;
regress lnahe age age2 fem_age fem_age2 bac_age bac_age2 female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
test age age2 fem_age fem_age2 bac_age bac_age2;
test fem_age fem_age2 bac_age bac_age2;
test fem_age fem_age2;
test bac_age bac_age2;
regress lnahe age age2 fem_age fem_age2 bac_age bac_age2 fem_bac_age fem_bac_age2
 female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
test age age2 fem_age fem_age2 bac_age bac_age2 fem_bac_age fem_bac_age2;
test fem_bac_age fem_bac_age2;
test fem_age fem_age2 bac_age bac_age2 fem_bac_age fem_bac_age2;;
test fem_age fem_age2;
test bac_age bac_age2;
clear;
***********************;
* Analysis for 1992;
***********************;
use cps92_98.dta; 
keep if year == 1992;
replace ahe = ahe*(163.0/140.3);
generate age2 = age*age;
generate lnahe = ln(ahe);
generate fem_bac=female*bachelor;
generate lnage = ln(age);
generate fem_age=female*age;
generate fem_age2=female*age2;
generate bac_age=bachelor*age;
generate bac_age2=bachelor*age2;
generate fem_bac_age = fem_bac*age;
generate fem_bac_age2 = fem_bac*age2;
regress lnahe age female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
regress lnahe age age2 female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
test age age2;
regress lnahe lnage female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
regress lnahe age age2 fem_age fem_age2 female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
test age age2 fem_age fem_age2;
test fem_age fem_age2;
regress lnahe age age2 bac_age bac_age2 female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
test age age2 bac_age bac_age2;
test bac_age bac_age2;
regress lnahe age age2 fem_age fem_age2 bac_age bac_age2 female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
test age age2 fem_age fem_age2 bac_age bac_age2;
test fem_age fem_age2 bac_age bac_age2;
test fem_age fem_age2;
test bac_age bac_age2;
regress lnahe age age2 fem_age fem_age2 bac_age bac_age2 fem_bac_age fem_bac_age2
 female bachelor fem_bac, robust;
dis "Adjusted Rsquared = " _result(8);
test age age2 fem_age fem_age2 bac_age bac_age2 fem_bac_age fem_bac_age2;
test fem_bac_age fem_bac_age2;
test fem_age fem_age2 bac_age bac_age2 fem_bac_age fem_bac_age2;;
test fem_age fem_age2;
test bac_age bac_age2;
log close;
clear;
exit;