library(tigerstats) set.seed(1234) n<-100000 z<-rnorm(n,0,1) (mean(z)) (sd(z)) par(mfrow=c(2,1)) hist(z, freq=F,xlim=c(-3,3),breaks=1000) par(new=T) plot(density(z), axes=F, main="", xlim=c(-4,4), lwd=2, col="blue") qnormGC(0.975,region="below",m=0,s=1, graph=T) z00<-rep(NA,10) z01<-rep(NA,10) z02<-rep(NA,10) z03<-rep(NA,10) z04<-rep(NA,10) for(i in 1:10) { z00[i]<-pnorm((i-1)/100, 0, 1)-0.5 } (z00<-round(z00, digits=4)) for(i in 10:19) { z01[i]<-pnorm(i/100, 0, 1)-0.5 } (z01<-round(z01[10:19], digits=4)) for(i in 20:29) { z02[i]<-pnorm(i/100, 0, 1)-0.5 } (z02<-round(z02[20:29], digits=4)) for(i in 30:39) { z03[i]<-pnorm(i/100, 0, 1)-0.5 } (z03<-round(z03[30:39], digits=4)) for(i in 40:49) { z04[i]<-pnorm(i/100, 0, 1)-0.5 } (z04<-round(z04[40:49], digits=4)) zdist<-rbind(z00,z01,z02,z03,z04) (zdist<-round(zdist, digits=4))