DATA=data.frame(x1 = c(sample(c(letters[1:5], NA), 1000, r = T)), x2 = runif(1000), x3 = runif(1000),