\begin{Schunk}
\begin{Sinput}
 par(xpd=T) 
 wh <- hist(women$sei, xlab="Socio-economic Index", ylab="Proportion", xlim=c(0,100), ylim=c(0, 0.04), breaks=25, prob = TRUE, main="")
 lines(womend, lty=2)
 text(10, 0.01, label=c("kernel \n density \n estimate"),pos=3)
 arrows(10, 0.01, 18,0.003, lty=2, length=0.1)
 text( 0.8*max(women$sei, na.rm = TRUE), 0.03, pos = 3, 
  label = paste("mean=",round(mean(women$sei,na.rm = TRUE),2), "\n std.dev.=", round(sd(women$sei, na.rm = T),2)))
\end{Sinput}
\end{Schunk}
