n=100 media=170 dt=10 muestra=rnorm(n,media,dt) xlim=c(media-3*dt,media+3*dt) ylim=c(-3,3) m=mean(muestra) s=sd(muestra) et=s/sqrt(n) #ERROR TIPICO hist(muestra,xlim=xlim) rug(muestra) #plot(c(muestra),c(rep(0.5,n)),xlim=xlim,ylim=ylim) abline(v=media) points(m,0.5,col="#FF0000",cex=1,pch=20) lines(c(m-2*et,m+2*et),c(1,1),col="#0000FF",lty=2)