title: Overlay Linear Regression Line on Scatter Plot (iPython Notebook)
body: <pre><code>gh_data = ascii.read('http://dept.astro.lsa.umich.edu/~ericbell/data/GHOSTS/M81/ngc3031-    field15.newphoto_radec')
  ra = gh_data['col5'][:]
  dec = gh_data['col6'][:]
  f606 = gh_data['col3'][:] 
  f814 = gh_data['col4'][:]
  plot(f6062-f8142,f8142, 'bo', alpha=0.15)  
  axis([-1,2.5,27,23]) 
  xlabel('F606W-F814W')
  ylabel('F814W')
  title('Field 14')
  </code></pre>
  
  <p>The data set is imported and organized into different columns, I am trying to overlay a line of best fit, or linear regression over the scatterplot created, but I cannot figure out how. Thanks in advance.</p>