P <- c(0.0000616649,0.0012931876,0.0014858932,0.0034575074,0.0095432743,0.0112784208,0.0198140078, 0.0260565422,0.0378525090,0.0476971273,0.0665802025,0.1160787054,0.1561110462,0.1741858728, 0.2592136466,0.3849843314,0.3970805883,0.4387950690,0.5686058809,0.5880746208,0.6367807765, 0.7164637107,0.7548314071,0.8594174096,0.8637551603,0.8852179374,0.8854362373,0.8904200780, 0.9319782385,0.9411071229,0.9474470330,0.9605232158,0.9621474910,0.9716238220,0.9750371388, 0.9800862502,0.9856935080,0.9923052342,0.9993104279,0.9994746329,0.9997647547,0.9999417310, 0.9999506389,0.9999650462,0.9999825779,0.9999967088,0.9999994243,0.9999999681) ordinato <- sort(P) n <- length(ordinato) plotpos <- seq(0.5/n, (n - 0.5)/n, by = 1/n) postscript(file="fig1.eps",family="NimbusSan",width=8,height=8, horiz=F,onefile=F,pointsize=16,paper="special") # was 20, changed for sinica paper plot(ordinato, plotpos, xlab="Expected probability", ylab="Observed probability") abline(0,1,lty=3) grid() dev.off()