set.seed(123)
x_rand <- matrix(rnorm(50*5), ncol=5)
colnames(x_rand) <- paste0("Var", 1:5)
y_rand <- factor(rep(c("A", "B"), 25))

featurePlot(x = x_rand, y = y_rand, plot = "box")

x_sepal <- iris$Sepal.Length
par(mfrow = c(2, 2))
plot(discretize(x_sepal, method = "interval", breaks = 3), main = "Interval")
plot(discretize(x_sepal, method = "frequency", breaks = 3), main = "Frequency")
plot(discretize(x_sepal, method = "cluster", breaks = 3), main = "Cluster")
plot(discretize(x_sepal, method = "fixed", breaks = c(4, 5.5, 7, 8)), main = "Fixed")

data("Ozone")
ozone_clean <- na.omit(Ozone)
set.seed(123)
boruta_res <- Boruta(V4 ~ ., data = ozone_clean, doTrace = 0)
plot(boruta_res, las = 2, cex.axis = 0.7)