if (FALSE) {
n <- 1e6
## rho > 0
mean(rnorm(n) + 1.5 * rexp(n, rate = 2) <= 0.1)
pnormexp(0.1, rate = 2, rho = 1.5)
mean(rnorm(n) + 1.5 * rexp(n, rate = 2) <= 0.9)
pnormexp(0.9, rate = 2, rho = 1.5)
## rho < 0
mean(rnorm(n) - 1.5 * rexp(n, rate = 2) <= 0.1)
pnormexp(0.1, rate = 2, rho = -1.5)
mean(rnorm(n) - 1.5 * rexp(n, rate = 2) <= 0.9)
pnormexp(0.9, rate = 2, rho = -1.5)
}
Run the code above in your browser using DataLab