if (FALSE) {
n <- 1e6
pnormweibull(0.1, scale = 1/2, shape = 1, rho = 1.5)
pnormweibull(0.8, scale = 1/2, shape = 1, rho = 1.5)
mean(rnorm(n) + 1.5 * rweibull(n, scale = 1/2, shape = 1) <= 0.1)
mean(rnorm(n) + 1.5 * rweibull(n, scale = 1/2, shape = 1) <= 0.8)
pnormweibull(0.1, scale = 1/2, shape = 1, rho = -1.5)
pnormweibull(0.8, scale = 1/2, shape = 1, rho = -1.5)
mean(rnorm(n) - 1.5 * rweibull(n, scale = 1/2, shape = 1) <= 0.1)
mean(rnorm(n) - 1.5 * rweibull(n, scale = 1/2, shape = 1) <= 0.8)
pnormweibull(0.1, scale = 1/2, shape = 2, rho = -1.5)
pnormweibull(0.8, scale = 1/2, shape = 2, rho = -1.5)
mean(rnorm(n) - 1.5 * rweibull(n, scale = 1/2, shape = 2) <= 0.1)
mean(rnorm(n) - 1.5 * rweibull(n, scale = 1/2, shape = 2) <= 0.8)
}
Run the code above in your browser using DataLab