## Not run:
# Y.obs <- 11:20
# res <- matrix(0, 10, 50)
# for(i in 1:50){
# Y.rep <- matrix(rpois(10*5000, 15), 10, )
# res[, i] <- cdfU(Y.obs, Y.rep)
# }
# matplot(t(res), type="l")
# abline(h = ppois(11:20, 15))
# ## End(Not run)
Run the code above in your browser using DataLab