code: def range_cum_s(magnitude):
      sigma = np.std(magnitude)
      N = len(magnitude)
      m = np.mean(magnitude)
      s = np.cumsum(magnitude - m) * 1.0 / (N * sigma)
      R = np.max(s) - np.min(s)
      return {"Rcs": R}
  
imports: import pandas
  import numpy
  from statsmodels.tsa.stattools import acf, adfuller, pacf
  import itertools
  from statsmodels.tsa.ar_model import AR
  from scipy.signal import cwt, find_peaks_cwt, ricker, welch
  from scipy.stats import linregress
  import scipy.stats
  from scipy.stats import kurtosis
  from scipy.stats import skew
  import numpy
  from scipy import signal, integrate
  from scipy.interpolate import interp1d
  import math
  from statsmodels.tools.sm_exceptions import MissingDataError
  from numpy.linalg import LinAlgError
label_desc: Use image node_id 50 for calling a global function with example usage: range_cum_s(magnitude) and returns: dict