title: Adding a scatter of points to a boxplot using matplotlib
body: <p>I have seen this wonderful boxplot in <a href="http://www.nmr.mgh.harvard.edu/~iglesias/pdf/n3_bias_field_correction_explained_as_a_bayesian_modeling_method.pdf" rel="noreferrer">this article</a> (Fig.2).</p>
  
  <p><img src="https://i.stack.imgur.com/yO8F6.png" alt="A wonderful boxplot"></p>
  
  <p>As you can see, this is a boxplot on which are superimposed a scatter of black points: x indexes the black points (in a random order), y is the variable of interest. I would like to do something similar using Matplotlib, but I have no idea where to start. So far, the boxplots which I have found online are way less cool and look like this:</p>
  
  <p><img src="https://i.stack.imgur.com/mzql7.png" alt="Usual boxplots"></p>
  
  <p>Documentation of matplotlib:
  <a href="http://matplotlib.org/api/pyplot_api.html#matplotlib.pyplot.boxplot" rel="noreferrer">http://matplotlib.org/api/pyplot_api.html#matplotlib.pyplot.boxplot</a></p>
  
  <p>Ways to colorize boxplots:
  <a href="https://github.com/jbmouret/matplotlib_for_papers#colored-boxes" rel="noreferrer">https://github.com/jbmouret/matplotlib_for_papers#colored-boxes</a></p>