code: def _compute_quantile_wrapper(
      x, q, axis=None, keepdims=False, interpolation="linear", out=None
  ):
      if not _validate_quantile(q):
          raise ValueError("Quantiles must be in the range [0, 1]")
      if interpolation in [
          "linear",
          "lower",
          "higher",
          "midpoint",
          "nearest",
          "nearest_jax",
      ]:
          if interpolation == "nearest_jax":
              return _handle_axis(
                  x, q, _quantile, keepdims=keepdims, axis=axis
              )
          else:
              axis = tuple(axis) if isinstance(axis, list) else axis
  
              return np.quantile(
                  x,
                  q,
                  axis=axis,
                  method=interpolation,
                  keepdims=keepdims,
                  out=out,
              ).astype(x.dtype)
      else:
          raise ValueError(
              "Interpolation must be 'linear', 'lower', 'higher', 'midpoint' or 'nearest'"
          )
  
imports: from typing import Optional, Union, Tuple, Sequence
  import numpy
  import math
  import ivy
  from ivy.func_wrapper import with_unsupported_dtypes
  from .None import backend_version
  from ..statistical import _infer_dtype
  from copy import deepcopy
label_desc: Use image node_id 10 for calling a global function with example usage: _compute_quantile_wrapper(x, q, axis, keepdims, interpolation, out) and returns: _handle_axis, np