?rnorm
dnorm(0)
exp(-0^2/2)/sqrt(2*pi)
pnorm(0)
pnorm(1)-pnorm(-1)
qnorm(0.5)
qnorm(0.025)
sims=rnorm(1000)
hist(sims)
summary(sims)
pnorm(1,mean=2,sd=1)
qnorm(0.5,1100207,23)
sims=rnorm(1000,mean=100,sd=20)
hist(sims)
hist(sims, prob=TRUE)
curve(dnorm(x, mean=100, sd=20), col="red", 
      lwd=2, add=TRUE)
summary(sims)
bp=boxplot(sims)
bp$out
