It would be of more interest to analyze the distribution of the integer results of the decimation, at various "precision" settings.
The analogReadXXbit method returns a float, so will need modify it to return unsigned long, or put debug prints in the library.
float avg_reading = (float)reading_sum/(float)num_samples_to_avg;
return avg_reading;