I have recently started reading ISLR and am finding the plots in the
book very useful.
A visualization aid one often uses for exploratory data analysis is a scatter
plot of the response variable against a potential predictor. Overlaying the
ordinary least squares fit line on this scatter provides a readily accessible
visual representation of the effect of the predictor on the response (if any).
Following is a simple snippet that I wrote in R to plot such graphs for any
arbitrary dataset with some numeric response variable. Note that the function
only attempts the plots for predictors which are numeric (or integer). It also
attempts a crude adjustment of the layout of the plot according to the number
of predictors.