#
# Load the non-white men's distribution
# downloaded from
# http://csde.washington.edu/~handcock/RelDist/Data/S/RDBnonwmftfy.RData
#
load("RDBnonwmftfy.RData")
#
# Convert to thousands of dollars
#
ebw88 <- ibw88 / 1000
#
pdf(file = "FigLec6.1.pdf",width=10.5,height=7.5,horiz=T)
par(mfrow=c(1,2))
#
hist(ebw88, type='n', xlim=c(0,80), probability=T,
  main="Earnings of Black Women", ylim=c(0,0.06),
  xlab = "Earnings (thousands of dollars)", ylab = "Probability")
lines(density(ebw88, bw="SJ"))
#
bw.SJ(ebw88)
#
# add a legend
#
legend(x=c(20,20),y=c(0.05,0.05),lty=1,legend=c("default=Sheather-Jones"))
#
hist(ebw88, type='n', xlim=c(0,80), probability=T,
  main="Earnings of Black Women", ylim=c(0,0.06),
  xlab = "Earnings (thousands of dollars)", ylab = "Probability")
lines(density(ebw88, bw=1.541919/2))
lines(density(ebw88, bw=1.541919*2), lty=2)
#
# add a legend
#
legend(x=c(30,30),y=c(0.05,0.05),lty=1:2,legend=c("h=0.77","h=3.08"))