#We can make a scatterplot matrix by:
pairs(sahd)

#We can make a correlation matrix and visualize it by:
library(corrplot)
library(RColorBrewer)

M <-cor(sahd[,c(-5,-10)])
corrplot(M, method="circle")
corrplot(M, type="upper", order="hclust",
         col=brewer.pal(n=8, name="RdYlBu"))

#We can	 make boxplots of the data by:
boxplot(sahd)
boxplot(M)