code: def _get_response(S0s, lambdas):
      S0 = np.mean(S0s) if S0s.size else 0
      # Check if lambdas is empty
      if not lambdas.size:
          response = (np.zeros(3), S0)
          ratio = 0
          return response, ratio
      l01 = np.mean(lambdas, axis=0) if S0s.size else 0
      evals = np.array([l01[0], l01[1], l01[1]])
      response = (evals, S0)
      ratio = evals[1] / evals[0]
      return response, ratio
  
imports: import warnings
  import numpy
  import numbers
  from scipy.integrate import quad
  from scipy.special import lpn, gamma
  import scipy.linalg
  import scipy.linalg.lapack
  from dipy.data import small_sphere, get_sphere, default_sphere
  from dipy.core.geometry import cart2sphere
  from dipy.core.ndindex import ndindex
  from dipy.sims.voxel import single_tensor
  from dipy.reconst.multi_voxel import multi_voxel_fit
  from dipy.reconst.dti import TensorModel, fractional_anisotropy
  from dipy.reconst.shm import sph_harm_ind_list, real_sh_descoteaux_from_index, sph_harm_lookup, lazy_index, SphHarmFit, real_sh_descoteaux, sh_to_rh, forward_sdeconv_mat, SphHarmModel
  from dipy.reconst.utils import _roi_in_volume, _mask_from_roi
  from dipy.direction.peaks import peaks_from_model
  from dipy.core.geometry import vec2vec_rotmat
  from dipy.utils.deprecator import deprecate_with_version, deprecated_params
label_desc: Use image node_id 12 for calling a global function with example usage: _get_response(S0s, lambdas) and returns: response,  ratio, response,  ratio