def downstddv(x): if x.size < 4: return np.nan median = np.nanpercentile(x, 50) return np.nanstd(x[x < median])