It can happen that propLess is equal to 1. Then, z0 = quantile("normal", propLess) does not work. What should be done then?