set.seed(123)
x <- rnorm(100)
UPM.ratio(0, mean(x), x)
## Joint Upper CDF
if (FALSE) {
x <- rnorm(5000) ; y <- rnorm(5000)
plot3d(x, y, Co.UPM(0, sort(x), sort(y), x, y), col = "blue", xlab = "X", ylab = "Y",
zlab = "Probability", box = FALSE)
}
Run the code above in your browser using DataLab