NeuroKit: Calculating time of inspiration and expiration using RSP

Hi, I had the idea to include average inspiration and expiration times in the output of rsp_intervalrelated, but my solution right now seems a bit unpythonic (is that a word?). Any ideas for making this cleaner?

# Process it
signals, info = nk.rsp_process(rsp, sampling_rate=250)

# Get indexes of phase changes
phase_peaks = signals[signals.RSP_Phase_Completion == 1]

# Calc diff for each phase change
phase_peaks['phase_len'] = np.diff(phase_peaks.index, prepend=0) 
phase_peaks['phase_time'] = phase_peaks.phase_len/250

# make df for expirations and inspirations
insp = phase_peaks[phase_peaks.RSP_Phase == 0]
exp = phase_peaks[phase_peaks.RSP_Phase == 1]

print("Mean insp_time:", insp.phase_time.mean())
print("Mean exp_time:", exp.phase_time.mean())