e2 <- function(x) {
out <- 0*x
for(i in 1:length(x))
out[i] <-integrate(function(y) qnorm(y),lower=0,upper=x[i])$value
out }
integrate(e2,lower=0, upper=a)$value
above is my code , when a is small , say a<0.45 the result is right .
however , when a>0.5
the result is incorrect . why ? thank you .
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
