Various code improvements for EDA:

- use log10 Y-axis scale for continuous variables' histogram plots;
- clear R environment at start;
- replace qplot() calls with ggplot() calls for consistency and more flexible control;
- used 'identity' to prevent warnings;
- added and tested Gamma distribution function for QQ plots (currently disabled, may need to use different value of the 'shape' parameter).