\begin{Schunk}
\begin{Sinput}
 par(xpd=T)
 h1 <- hist(dat$nummen, breaks=50, prob=T, main="Male Sexual Partners (nummen)", xlab="Number of Male Sexual Partners", ylim=c(0,0.5)) 
 dennummen <- density(dat$nummen, na.rm=T)
 lines(dennummen, lty=2, lwd=1, col="green")
 text(35, 0.15, label=c("kernel density \n estimate (green)"),pos=3)
 arrows(30, 0.15, 18,0.03, lty=2, length=0.1)
\end{Sinput}
\end{Schunk}
