Do you mean those dots above and below the two black lines?
from pylab import * spread= rand(50) * 100 center = ones(25) * 50 flier_high = rand(10) * 100 + 100 flier_low = rand(10) * -100 data =concatenate((spread, center, flier_high, flier_low), 0) r = boxplot(data)

Save the returned dict from boxplot, and you can get all the information from it, for example:
top_points = r["fliers"][0].get_data()[1] bottom_points = r["fliers"][2].get_data()[1] plot(np.ones(len(top_points)), top_points, "+") plot(np.ones(len(bottom_points)), bottom_points, "+")

source share