ma import masked_invalid if x is None: from numpy import arange x = arange(len(rc)) if ax is None: ax = gca() from scipy.stats import scoreatpercentile upper95 = array([scoreatpercentile(rc_controls[:,i], 97.5) for i in range(len(x))]) lower95 = array([scoreatpercentile(rc_...