rcauchy(n) rnorm(n) * runif(n, 0, 1) x <- rnorm(n, mean = 0, sd = 1); x + x^3 contaminate(n, prob.bad = 0.1, badmu = 2) contaminate(n, prob.bad = 0.2, badstd = 2.5)