Linear regression plot for magnitude and earthquake stations

An attempt at the bonus. Results do not align with publication that does not read from Github. Would like to discuss issues.

Bonus not included