goptions reset=all; 
data a; do p = .005 to .9995 by .0005; 
YC = log(-1*log(1-p)); * Complementary log-log;
Y= log(p/(1-p));  
linear = -1 + 0.5*Y;  *linear model in Y; 
PC = 1-exp(-1*exp(linear)); * Inverse of CLL; 
PL = exp(linear)/(1+exp(linear)); *Inverse Logisitc Link; 
PN = probnorm(linear);  *Probit; 
Elliott = linear/(1+abs(linear)); 
HypTan = 2*PL-1; * Hyperbolic Tangent; 
output; end; 

**** Graphics *****;  
symbol1 v=none i=join w=4;  ** W is multiplier for width of line **;
                            ** V is plot symbol **; 
                            ** With no color specified, cycle through 
                               device color list **; 
proc gplot; 
plot (PC PL PN Elliott Hyptan)*Y/overlay legend href=2; 
proc corr; var PC PL PN Elliott Hyptan;
run; quit;