code: def glover_time_derivative(
      tr, oversampling=16, time_length=32.0, onset=0.0
  ):
      """Implementation of the flover time derivative hrf (dhrf) model
  
      Parameters
      ----------
      tr: float, scan repeat time, in seconds
      oversampling: int, temporal oversampling factor, optional
      time_length: float, hrf kernel length, in seconds
      onset: float, onset of the response
  
      Returns
      -------
      dhrf: array of shape(length / tr, float),
            dhrf sampling on the provided grid
      """
      do = 0.1
      dhrf = (
          1.0
          / do
          * (
              glover_hrf(tr, oversampling, time_length, onset + do)
              - glover_hrf(tr, oversampling, time_length, onset)
          )
      )
      return dhrf
  
file_path: nipy/nipy/modalities/fmri/hemodynamic_models.py
imports: import warnings
  import numpy
  from scipy.stats import gamma
label_desc: Use image node_id 5 for calling a global function with example usage: glover_time_derivative(tr, oversampling, time_length, onset) and returns: dhrf