The first figure was plotted using these instructions to R:

moist=c(32.200,32.300,33.100,33.200,33.300,34.500, 35.200,35.300,36.500,36.800,37.000,37.600) obsf=c(1/12,2/12,3/12,4/12,5/12,6/12,7/12,8/12,9/12,10/12,11/12,12/12) plot(moist,obsf,type='s',ylim=c(0,1)) points(qnorm(p,35,2),p,type='l',col='red') p1=pnorm(moist,35,2) p2=rank(moist)/length(moist) p3=c(0,p2[-length(p2)]) points(moist,p1,col='red',pch=20) points(moist,p2,col='blue',pch=20) points(moist,p3,col='green',pch=20) d1=p2-p1 d2=p3-p1 d=max(abs(c(d1,d2))) d

The second figure was plotted using these instructions to R:

moist=c(32.200,32.300,33.100,33.200,33.300,34.500,35.200,35.300,36.500,36.800,37.000,37.600) obsf=c(.5/12,1.5/12,2.5/12,3.5/12,4.5/12,5.5/12,6.5/12,7.5/12,8.5/12,9.5/12,10.5/12,11.5/12) plot(moist,obsf,type='s',ylim=c(0,1)) points(qnorm(p,35,2),p,type='l',col='red') p1=pnorm(moist,35,2) p2=(rank(moist)-0.5)/length(moist) points(moist,p1,col='red',pch=20) points(moist,p2,col='blue',pch=20) d1=p2-p1 d=max(abs(d1)) + (1/(2*length(moist))) d

Giving this: