/*Stata Lab 2 */
clear
capture log close
version 11
log using statalab2, replace

cd "m:\ec114 statalab"
use bhps_lab


*Explore the data
describe

tab sex
tab education
tab education sex, column

sum baserate
sum baserate, detail
display as text "Average Basic Hourly Wage = " as result r(mean)
display as text "Standard Deviation = " as result r(sd)
display as text "Median Basic Hourly Wage = " as result r(p50)
display as text "Inter Quartile Range = " as result (r(p75)-r(p25))
display as text "Range = " as result (r(max)-r(min))

sum jbhrs
sum jbhrs, detail
display as text "Average Weekly Working Hours = " as result r(mean)
display as text "Standard Deviation = " as result r(sd)
display as text "Median Weekly Working Hours = " as result r(p50)
display as text "Inter Quartile Range = " as result (r(p75)-r(p25))
display as text "Range = " as result (r(max)-r(min))

bysort sex: sum baserate
bysort sex: sum baserate, detail

tab sex
tab sex, nolabel

sum baserate if sex==1
sum baserate if sex==1, detail
display as text "Average Basic Hourly Wage (Men) = " as result r(mean)
display as text "Standard Deviation = " as result r(sd)
display as text "Median Basic Hourly Wage (Men) = " as result r(p50)
display as text "Inter Quartile Range = " as result (r(p75)-r(p25))
display as text "Range = " as result (r(max)-r(min))

sum baserate if sex==2
sum baserate if sex==2, detail
display as text "Average Basic Hourly Wage (Women) = " as result r(mean)
display as text "Standard Deviation = " as result r(sd)
display as text "Median Basic Hourly Wage (Women) = " as result r(p50)
display as text "Inter Quartile Range = " as result (r(p75)-r(p25))
display as text "Range = " as result (r(max)-r(min))


bysort sex: sum jbhrs
bysort sex: sum jbhrs, detail

sum jbhrs if sex==1
sum jbhrs if sex==1, detail
display as text "Average Weekly Working Hours (Men) = " as result r(mean)
display as text "Standard Deviation = " as result r(sd)
display as text "Median Weekly Working Hours (Men) = " as result r(p50)
display as text "Inter Quartile Range = " as result (r(p75)-r(p25))
display as text "Range = " as result (r(max)-r(min))

sum jbhrs if sex==2
sum jbhrs if sex==2, detail
display as text "Average Weekly Working Hours (Women) = " as result r(mean)
display as text "Standard Deviation = " as result r(sd)
display as text "Median Weekly Working Hours (Women) = " as result r(p50)
display as text "Inter Quartile Range = " as result (r(p75)-r(p25))
display as text "Range = " as result (r(max)-r(min))


log close