code: def test_median_biases(n):
      """Test vectorization of median_biases."""
      want_biases = np.concatenate(
          ([1.0, 1.0], [_median_bias(ii) for ii in range(2, n + 1)])
      )
      got_biases = _median_biases(n)
      assert_allclose(want_biases, got_biases)
      assert_allclose(got_biases[n], _median_bias(n))
      assert_allclose(got_biases[:3], 1.0)
  
file_path: mne-python/mne/time_frequency/tests/test_psd.py
imports: import numpy
  import pytest
  from numpy.testing import assert_allclose, assert_array_almost_equal, assert_array_equal
  from scipy.signal import welch
  from mne.time_frequency import psd_array_multitaper, psd_array_welch
  from mne.time_frequency.multitaper import _psd_from_mt
  from mne.time_frequency.psd import _median_biases
  from mne.utils import catch_logging
label_desc: Use image node_id 8 for calling a global function with example usage: test_median_biases(n) without return types