h47r.py

# 0 or 1 indicates absence/presence of variant
data = [
    [0, 1, 0, 1, 0, 1],
    [0, 0, 0, 1, 1, 1],
    [0, 0, 0, 1, 1, 1],
    [0, 0, 0, 1, 1, 1],
    [0, 0, 0, 1, 1, 1],
    [0, 0, 0, 0, 1, 1],
    [0, 0, 0, 0, 0, 0],
    [1, 1, 1, 0, 0, 0],
    ]

from matplotlib import pyplot as plt
fig, ax = plt.subplots()

import numpy as np

alpha = 3 * 1.0 / len(data)

for trace in data:
    trace = np.cumsum(trace)

    plt.step(range(1, len(trace) + 1), trace, color='k', alpha=alpha, lw=15)

plt.xlabel('copy number')
plt.ylabel('cumulative number of H47R variants')
plt.show()