\begin{Schunk}
\begin{Sinput}
 myx2 <- rnorm(1000, 20, 20)
 par(mfcol=c(1,2))
 hist(myx, breaks=30, prob=TRUE, xlab="A Beautiful Variable", ylab="Density", main="", xlim=c(-40, 120), ylim=c(0, 0.025))
 text(10, 0.020, label=paste("Mean = ", round(mx,2)), cex=1.5) 
 #text(82, 0.018, label=paste("Std. Dev. = ", round(sdx,2)))
 hist(myx2, breaks=30, prob=T, xlab="Another Variable", ylab="Density", main="", xlim=c(-40,120), ylim=c(0,0.025))
 mx2 <- mean(myx2, na.rm=T)
 sdx2 <- sd(myx2, na.rm=T)
 text(45, 0.021, label=paste("Mean = ", round(mx2,2)), cex=1.5) 
\end{Sinput}
\end{Schunk}
