# creation of the sample alpha <- 0.2; beta <- 0.8; # the following assumes alpha < 0.5 < beta rexemple <- function(n){ res <- list(); res$x <- cbind(runif(n), runif(n)); res$y <- array(n); u <- runif(n); T = 2 * res$x[,1] + res$x[,2] < 1.5; res$y[T] <- u[T]