class BolusArrivalError(Exception):
    pass


def bolus_arrival_index(time_series):
    """
    Determines the DCE bolus arrival series index. The bolus arrival is
    the first series with a difference in average signal larger than
    double the difference from first two points.
 
    :param time_series: the 4D NiFTI scan image file path
    :return: the bolus arrival series index
    :raise BolusArrivalError: if the bolus arrival could not be determined
    """
    import nibabel as nb
    import numpy as np

    nii = nb.load(time_series)
    data = nii.get_data()
    n_vols = data.shape[-1]
    signal_means = np.array([np.mean(data[:,:,:, idx])
                             for idx in xrange(n_vols)])
    signal_diffs = np.diff(signal_means)

    # If we see a difference in average signal larger than double the
    # difference from first two points, take that as bolus arrival.
    base_diff = np.abs(signal_diffs[0])
    for idx, diff_val in enumerate(signal_diffs[1:]):
        if diff_val > 2 * base_diff:
            return idx + 1
    else:
        raise BolusArrivalError("Unable to determine bolus arrival")

@@ -119,30 +119,34 @@
 lus arrival 
-series
+time point
  index. The 
@@ -162,16 +162,16 @@
 ival
- is
 %0A    
+is 
 the 
@@ -180,19 +180,20 @@
 rst 
-series with
+occurence of
  a d
@@ -372,22 +372,26 @@
 arrival 
-series
+time point
  index%0A 

Change series to time point.