#!/usr/bin/env python
# coding: utf-8

# In[1]:


import numpy as np
from sklearn.metrics import median_absolute_error as skmedian_absolute_error


# In[2]:


def median_absolute_error(y_true, y_pred):
    return np.median(np.abs(y_true - y_pred))


# In[3]:


for i in range(10):
    rng = np.random.RandomState(i)
    y_true = rng.rand(10)
    y_pred = rng.rand(10)
    score1 = median_absolute_error(y_true, y_pred)
    score2 = skmedian_absolute_error(y_true, y_pred)
    assert np.isclose(score1, score2)