The pnorm function in R does not work properly. For the standard normal distribution, the pnorm should be less than 1.
> pnorm(1.96) [1] 1.96 # suppose to be 0
