\begin{Schunk}
\begin{Sinput}
 par(xpd=T)
 h1 <- hist(dat$nummen, breaks=50, prob=T, main="Male Sexual Partners", 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")
 legend("topright",legend=c("Histogram", "Kernel Density"), lty=c(1, 2), col=c("black","green"))
\end{Sinput}
\end{Schunk}
