# NOT RUN {
#========================================================================================
#            1)             validation of this function
#========================================================================================
#'
  x<-0.2
 Phi(x)==stats::pnorm(x)
#========================================================================================
#            1)             Build the data
#========================================================================================
#'
 a  <-  0.1;
 NX <-  222;
 x  <-  runif(100,-11,11)
 y  <-  Phi_inv(exp(a/NX) *Phi(x))-x
 plot(x,y)
 a  <-  0.1;
 NX <-  222;
 x  <-  runif(100,0,11)
 y  <-  Phi_inv(exp(a/NX) *Phi(x))-x
 plot(x,y)
 a  <-  0.1;
 NX <-  222;
 x  <-  runif(100,2,4)
 y  <-  Phi_inv(exp(a/NX) *Phi(x))-x
 plot(x,y)
 a  <-  0.01;
 NX <-  222;
 x  <-  runif(100,2,4);
 y  <-  Phi_inv(exp(a/NX) *Phi(x))-x
 plot(x,y)
 a  <-  0.01;
 NX <-  222;
 x  <-  runif(100,3.5,4);
 y  <-  Phi_inv(exp(a/NX) *Phi(x))-x
 plot(x,y)
# }
Run the code above in your browser using DataLab