view talk/code/plots.py @ 45:6daaf6a8e431

Add code to plot boys and girls
author Jordi Gutiérrez Hermoso <jordigh@octave.org>
date Sat, 14 May 2016 21:24:47 -0400
parents
children 7887a0b32539
line wrap: on
line source

import matplotlib.pyplot as plt
import matplotlib.cbook as cbook
import numpy as np

with open("../../data/men") as f:
    men = [float(x) for x in f.readlines()]

with open("../../data/women") as f:
    women = [float(x) for x in f.readlines()]

men_stats = cbook.boxplot_stats(men)
women_stats = cbook.boxplot_stats(women)

fig = plt.figure(figsize=(6,4))

# setup the figure and axes
bpAx = fig.add_axes([0.05, 0.7, 0.4, 0.2])   # left, bottom, width, height:
                                            # (adjust as necessary)
histAx = fig.add_axes([0.05, 0.2, 0.4, 0.5]) # left specs should match and
                                            # bottom + height on this line should
                                            # equal bottom on bpAx line
# plot stuff
bpAx.bxp(men_stats, vert=False, flierprops={"marker": 'x'})
histAx.hist(men, bins=100, color=[0.3, 0.3, 1])

# confirm that the axes line up
xlims = np.array([bpAx.get_xlim(), histAx.get_xlim()])
for ax in [bpAx, histAx]:
    ax.set_xlim([xlims.min(), xlims.max()])

bpAx.set_xticklabels([])  # clear out overlapping xlabels
bpAx.set_yticks([])  # don't need that 1 tick mark
bpAx.set_title("Ages of actors")



# setup the figure and axes
bpAx = fig.add_axes([0.55, 0.7, 0.4, 0.2])   # left, bottom, width, height:
                                            # (adjust as necessary)
histAx = fig.add_axes([0.55, 0.2, 0.4, 0.5]) # left specs should match and
                                            # bottom + height on this line should
                                            # equal bottom on bpAx line
# plot stuff
bpAx.bxp(women_stats, vert=False, flierprops={"marker": 'x'})
histAx.hist(women, bins=100, color='pink')

# confirm that the axes line up
xlims = np.array([bpAx.get_xlim(), histAx.get_xlim()])
for ax in [bpAx, histAx]:
    ax.set_xlim([xlims.min(), xlims.max()])

bpAx.set_xticklabels([])  # clear out overlapping xlabels
bpAx.set_yticks([])  # don't need that 1 tick mark
bpAx.set_title("Ages of actresses")

plt.show()