# NOT RUN {
y = c(rep(1,200),rep(0,200))
x = rnorm(400)
#Inducing an aberration enrichment signal by perturbing some of the cases
x[1:20]=x[1:20]-3;
res2 = aziz.test(y,x,rep=20000)
print_summary(res2)
which_aberrant(x,x,res2)
which_aberrant(c(-5,1.5,-2.5,-0.5,2),x,res2)#testing if new values are within the aberrant interval
# }
Run the code above in your browser using DataLab