\begin{Schunk}
\begin{Sinput}
 par(xpd=T)
 dat$nummen[dat$nummen > 99] <- NA
 mnummen <- round(mean(dat$nummen, na.rm=T),3) 
 sdnummen <- round(sd(dat$nummen, na.rm=T),3) 
 mednummen <- round(median(dat$nummen, na.rm=T),3) 
 cvnummen <- round(sdnummen/mnummen,3)
 h1 <- hist(dat$nummen, breaks=100, prob=T, main="Male Sexual Partners (nummen)", xlab="Number of Male Sexual Partners", ylim=c(0,0.7)) 
 text(max(h1$mids), max(h1$density)-.2, pos=2, label=paste("Mean=", mnummen,"\nStd.Dev=", sdnummen,"\nMedian=", mednummen,"\nC.V.=", cvnummen))  
\end{Sinput}
\end{Schunk}
