library(tidyverse)
Boston = MASS::Boston
length(Boston)

Boston %>%
  keep(is.numeric) %>% 
  gather() %>% 
  ggplot(aes(value)) +
  facet_wrap(~ key, scales = "free") +
  geom_histogram()

Boston %>%
  keep(is.numeric) %>% 
  gather() %>% 
  ggplot(aes(value)) +
  facet_wrap(~ key, scales = "free") +
  geom_boxplot()

Boston %>%
  keep(is.numeric) %>% 
  gather() %>% 
  ggplot(aes(value)) +
  facet_wrap(~ key, scales = "free") +
  geom_histogram() + scale_x_log10()

library(GGally)
Boston %>%
select(medv, age, crim, rad, indus, nox) %>% ggpairs()


library(corrplot)
corrplot(cor(Boston), tl.col = "black", order = "hclust", hclust.method = "average", addrect = 4, tl.cex = 0.7)