# 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