diff --git a/interpret_qc.py b/interpret_qc.py index 1967ee166e529524802befa26f2ad59b22aaf218..83d45e465aff9e4d36dd3592adb30860962d6a21 100644 --- a/interpret_qc.py +++ b/interpret_qc.py @@ -33,7 +33,7 @@ def qc_day(qc_path): frame = frame.combine_first(qc_frame) frame = frame.query('missingDataFlag == 0') - qc_frame_sum = qc_frame.sum(axis=0).to_string() + qc_frame_sum = (qc_frame > .95).sum(axis=0).to_string() plots = [] @@ -69,10 +69,10 @@ def plot_variable_qc(frame, loc, filename=None): return mpld3.fig_to_html(fig) def plot_outliers(frame, qc_mask, loc): - if (~np.isnan(qc_mask) & qc_mask > 0).any(): - frame.ix[(np.isnan(qc_mask) | (qc_mask == 0)) & (frame['qc_percent'] == 0), loc].plot(style='b.') - frame.ix[(np.isnan(qc_mask) | (qc_mask == 0)) & (frame['qc_percent'] > 0), loc].plot(style='k.', alpha=.2) - frame.ix[~np.isnan(qc_mask) & (qc_mask > 0) & (frame['qc_percent'] > 0), loc].plot(style='r.') + if (~np.isnan(qc_mask) & qc_mask > .95).any(): + frame.ix[(np.isnan(qc_mask) | (qc_mask < .95)) & (frame['qc_percent'] < .95), loc].plot(style='b.') + frame.ix[(np.isnan(qc_mask) | (qc_mask < .95)) & (frame['qc_percent'] > .95), loc].plot(style='k.', alpha=.2) + frame.ix[~np.isnan(qc_mask) & (qc_mask > .95) & (frame['qc_percent'] > .95), loc].plot(style='r.') plt.xlabel('Time') plt.title(loc)