notebooks/plot_old_versus_new_phenodigm_resnik_and_jaccard.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "id": "6gxb5C_USH_B"
   },
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "!python -m pip install oaklib"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# node.js needs to be at least version 12 or so for OAK\n",
    "!node --version"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "id": "jTF-gDH-UMb_"
   },
   "outputs": [],
   "source": [
    "# Files\n",
    "old_scores = 'https://kg-hub.berkeleybop.io/frozen_incoming_data/hp-mp-phenodigm-cache.txt.gz'\n",
    "new_scores = 'https://kg-hub.berkeleybop.io/kg-phenio/hp-mp_similarities.tar.gz'\n",
    "\n",
    "old_tar_file = 'hp-mp-phenodigm-cache.txt.gz'\n",
    "new_tar_file = 'hp-mp_similarities.tar.gz'\n",
    "\n",
    "old_sim_file = 'hp-mp-phenodigm-cache.txt'\n",
    "new_sim_file = 'KGPhenio_similarities'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "id": "WK1CujpFSq5I"
   },
   "outputs": [],
   "source": [
    "!wget {old_scores}\n",
    "!gunzip {old_tar_file}\n",
    "\n",
    "!wget {new_scores}\n",
    "!tar -xvzf {new_tar_file}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/"
    },
    "id": "EISuj4jUSwnH",
    "outputId": "e4728013-63d9-4234-c4ed-aeeac74837cc"
   },
   "outputs": [],
   "source": [
    "!ls"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 424
    },
    "id": "zu6BCrj5UPnv",
    "outputId": "77ceb655-5f68-44ad-dc34-6eee7c95b1b1"
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>source_old</th>\n",
       "      <th>destination_old</th>\n",
       "      <th>jaccard_old</th>\n",
       "      <th>resnik_score_old</th>\n",
       "      <th>subsumer_old</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0004014</td>\n",
       "      <td>0.470588</td>\n",
       "      <td>2.659134</td>\n",
       "      <td>MP_0002160;</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0020998</td>\n",
       "      <td>0.434783</td>\n",
       "      <td>3.300796</td>\n",
       "      <td>MP_0009198;</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0020957</td>\n",
       "      <td>0.363636</td>\n",
       "      <td>2.659134</td>\n",
       "      <td>MP_0002160;</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0020869</td>\n",
       "      <td>0.212766</td>\n",
       "      <td>3.300796</td>\n",
       "      <td>MP_0009198;</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0020824</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>3.300796</td>\n",
       "      <td>MP_0009198;</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13033730</th>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>MP:0004176</td>\n",
       "      <td>0.243243</td>\n",
       "      <td>5.040392</td>\n",
       "      <td>HP_0031703;</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13033731</th>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>MP:0030154</td>\n",
       "      <td>0.928571</td>\n",
       "      <td>7.560589</td>\n",
       "      <td>MP_0000049;</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13033732</th>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>MP:0030155</td>\n",
       "      <td>0.866667</td>\n",
       "      <td>7.560589</td>\n",
       "      <td>MP_0000049;</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13033733</th>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>MP:0030156</td>\n",
       "      <td>0.650000</td>\n",
       "      <td>7.560589</td>\n",
       "      <td>MP_0000049;</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13033734</th>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>MP:0030158</td>\n",
       "      <td>0.433333</td>\n",
       "      <td>7.560589</td>\n",
       "      <td>MP_0000049;</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>13033735 rows × 5 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          source_old destination_old  jaccard_old  resnik_score_old  \\\n",
       "0         HP:0000033      MP:0004014     0.470588          2.659134   \n",
       "1         HP:0000033      MP:0020998     0.434783          3.300796   \n",
       "2         HP:0000033      MP:0020957     0.363636          2.659134   \n",
       "3         HP:0000033      MP:0020869     0.212766          3.300796   \n",
       "4         HP:0000033      MP:0020824     0.500000          3.300796   \n",
       "...              ...             ...          ...               ...   \n",
       "13033730  HP:0008609      MP:0004176     0.243243          5.040392   \n",
       "13033731  HP:0008609      MP:0030154     0.928571          7.560589   \n",
       "13033732  HP:0008609      MP:0030155     0.866667          7.560589   \n",
       "13033733  HP:0008609      MP:0030156     0.650000          7.560589   \n",
       "13033734  HP:0008609      MP:0030158     0.433333          7.560589   \n",
       "\n",
       "         subsumer_old  \n",
       "0         MP_0002160;  \n",
       "1         MP_0009198;  \n",
       "2         MP_0002160;  \n",
       "3         MP_0009198;  \n",
       "4         MP_0009198;  \n",
       "...               ...  \n",
       "13033730  HP_0031703;  \n",
       "13033731  MP_0000049;  \n",
       "13033732  MP_0000049;  \n",
       "13033733  MP_0000049;  \n",
       "13033734  MP_0000049;  \n",
       "\n",
       "[13033735 rows x 5 columns]"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "old_sim_pd = pd.read_csv(old_sim_file, sep=\"\\t\", header=None)\n",
    "mapping = {old_sim_pd.columns[0]:'source_old',\n",
    "           old_sim_pd.columns[1]: 'destination_old', \n",
    "           old_sim_pd.columns[2]: 'jaccard_old',            \n",
    "           old_sim_pd.columns[3]: 'resnik_score_old',\n",
    "           old_sim_pd.columns[4]: 'subsumer_old',\n",
    "           }\n",
    "old_sim_pd = old_sim_pd.rename(columns=mapping)\n",
    "\n",
    "old_sim_pd['source_old'] = old_sim_pd['source_old'].str.replace('_',':')\n",
    "old_sim_pd['destination_old'] = old_sim_pd['destination_old'].str.replace('_',':')\n",
    "\n",
    "old_sim_pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 294
    },
    "id": "xeSa7hT8Wxwn",
    "outputId": "e9c245fd-4350-4382-f2c6-310159334fce"
   },
   "outputs": [],
   "source": [
    "old_sim_pd['jaccard_old'].plot(kind='hist', xlim=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 294
    },
    "id": "2-HXLoHVW6WN",
    "outputId": "2a137acf-c8c9-49f5-a30c-a3c20d1822fa"
   },
   "outputs": [],
   "source": [
    "old_sim_pd['resnik_score_old'].plot(kind='hist', bins=15, xlim=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 424
    },
    "id": "5WYrbeyiUrZX",
    "outputId": "5a4c03c4-14c1-4521-9f55-f31c44bb8e3e"
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>source</th>\n",
       "      <th>destination</th>\n",
       "      <th>resnik_score</th>\n",
       "      <th>jaccard</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>MP:3000005</td>\n",
       "      <td>MP:3000005</td>\n",
       "      <td>14.526168</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>HP:0030350</td>\n",
       "      <td>HP:0030350</td>\n",
       "      <td>14.526168</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>MP:0008925</td>\n",
       "      <td>MP:0008925</td>\n",
       "      <td>14.526168</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>MP:0008926</td>\n",
       "      <td>MP:0008926</td>\n",
       "      <td>14.526168</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>MP:0008926</td>\n",
       "      <td>MP:0013185</td>\n",
       "      <td>14.526168</td>\n",
       "      <td>0.857143</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51383633</th>\n",
       "      <td>HP:0001802</td>\n",
       "      <td>HP:0200054</td>\n",
       "      <td>5.041991</td>\n",
       "      <td>0.214286</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51383634</th>\n",
       "      <td>HP:0100142</td>\n",
       "      <td>MP:0009001</td>\n",
       "      <td>5.041991</td>\n",
       "      <td>0.117647</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51383635</th>\n",
       "      <td>HP:0010388</td>\n",
       "      <td>MP:0009001</td>\n",
       "      <td>5.041991</td>\n",
       "      <td>0.133333</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51383636</th>\n",
       "      <td>HP:0010323</td>\n",
       "      <td>HP:0001802</td>\n",
       "      <td>5.041991</td>\n",
       "      <td>0.214286</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51383637</th>\n",
       "      <td>HP:0010374</td>\n",
       "      <td>MP:0009001</td>\n",
       "      <td>5.041991</td>\n",
       "      <td>0.214286</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>51383638 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "              source destination  resnik_score   jaccard\n",
       "0         MP:3000005  MP:3000005     14.526168  1.000000\n",
       "1         HP:0030350  HP:0030350     14.526168  1.000000\n",
       "2         MP:0008925  MP:0008925     14.526168  1.000000\n",
       "3         MP:0008926  MP:0008926     14.526168  1.000000\n",
       "4         MP:0008926  MP:0013185     14.526168  0.857143\n",
       "...              ...         ...           ...       ...\n",
       "51383633  HP:0001802  HP:0200054      5.041991  0.214286\n",
       "51383634  HP:0100142  MP:0009001      5.041991  0.117647\n",
       "51383635  HP:0010388  MP:0009001      5.041991  0.133333\n",
       "51383636  HP:0010323  HP:0001802      5.041991  0.214286\n",
       "51383637  HP:0010374  MP:0009001      5.041991  0.214286\n",
       "\n",
       "[51383638 rows x 4 columns]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "new_sim_pd = pd.read_csv(new_sim_file, sep=\",\")\n",
    "new_sim_pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 294
    },
    "id": "D1goCywhXLnE",
    "outputId": "aa54b8a3-e60e-414b-ce73-bed93208ca46"
   },
   "outputs": [],
   "source": [
    "new_sim_pd['jaccard'].plot(kind='hist', xlim=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 294
    },
    "id": "ToXYHHHNXUyn",
    "outputId": "81d198d3-5263-4a5b-dfc8-af8d5fd29b69"
   },
   "outputs": [],
   "source": [
    "new_sim_pd['resnik_score'].plot(kind='hist', bins=15, xlim=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "id": "KZWylyZ-XvBs",
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>source_old</th>\n",
       "      <th>destination_old</th>\n",
       "      <th>jaccard_old</th>\n",
       "      <th>resnik_score_old</th>\n",
       "      <th>subsumer_old</th>\n",
       "      <th>source</th>\n",
       "      <th>destination</th>\n",
       "      <th>resnik_score</th>\n",
       "      <th>jaccard</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0004014</td>\n",
       "      <td>0.470588</td>\n",
       "      <td>2.659134</td>\n",
       "      <td>MP_0002160;</td>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0004014</td>\n",
       "      <td>5.358735</td>\n",
       "      <td>0.153846</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0020998</td>\n",
       "      <td>0.434783</td>\n",
       "      <td>3.300796</td>\n",
       "      <td>MP_0009198;</td>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0020998</td>\n",
       "      <td>7.041238</td>\n",
       "      <td>0.166667</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0020957</td>\n",
       "      <td>0.363636</td>\n",
       "      <td>2.659134</td>\n",
       "      <td>MP_0002160;</td>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0020957</td>\n",
       "      <td>7.038994</td>\n",
       "      <td>0.153846</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0020824</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>3.300796</td>\n",
       "      <td>MP_0009198;</td>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0020824</td>\n",
       "      <td>6.481542</td>\n",
       "      <td>0.363636</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0002771</td>\n",
       "      <td>0.400000</td>\n",
       "      <td>3.300796</td>\n",
       "      <td>MP_0009198;</td>\n",
       "      <td>HP:0000033</td>\n",
       "      <td>MP:0002771</td>\n",
       "      <td>7.041238</td>\n",
       "      <td>0.153846</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2361470</th>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>MP:0004204</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>7.560589</td>\n",
       "      <td>MP_0000049;</td>\n",
       "      <td>MP:0004204</td>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>10.195435</td>\n",
       "      <td>0.230769</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2361471</th>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>MP:0030154</td>\n",
       "      <td>0.928571</td>\n",
       "      <td>7.560589</td>\n",
       "      <td>MP_0000049;</td>\n",
       "      <td>MP:0030154</td>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>10.195435</td>\n",
       "      <td>0.272727</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2361472</th>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>MP:0030155</td>\n",
       "      <td>0.866667</td>\n",
       "      <td>7.560589</td>\n",
       "      <td>MP_0000049;</td>\n",
       "      <td>MP:0030155</td>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>10.195435</td>\n",
       "      <td>0.250000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2361473</th>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>MP:0030156</td>\n",
       "      <td>0.650000</td>\n",
       "      <td>7.560589</td>\n",
       "      <td>MP_0000049;</td>\n",
       "      <td>MP:0030156</td>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>10.195435</td>\n",
       "      <td>0.272727</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2361474</th>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>MP:0030158</td>\n",
       "      <td>0.433333</td>\n",
       "      <td>7.560589</td>\n",
       "      <td>MP_0000049;</td>\n",
       "      <td>MP:0030158</td>\n",
       "      <td>HP:0008609</td>\n",
       "      <td>8.446235</td>\n",
       "      <td>0.230769</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5833600 rows × 9 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         source_old destination_old  jaccard_old  resnik_score_old  \\\n",
       "0        HP:0000033      MP:0004014     0.470588          2.659134   \n",
       "1        HP:0000033      MP:0020998     0.434783          3.300796   \n",
       "2        HP:0000033      MP:0020957     0.363636          2.659134   \n",
       "3        HP:0000033      MP:0020824     0.500000          3.300796   \n",
       "4        HP:0000033      MP:0002771     0.400000          3.300796   \n",
       "...             ...             ...          ...               ...   \n",
       "2361470  HP:0008609      MP:0004204     0.500000          7.560589   \n",
       "2361471  HP:0008609      MP:0030154     0.928571          7.560589   \n",
       "2361472  HP:0008609      MP:0030155     0.866667          7.560589   \n",
       "2361473  HP:0008609      MP:0030156     0.650000          7.560589   \n",
       "2361474  HP:0008609      MP:0030158     0.433333          7.560589   \n",
       "\n",
       "        subsumer_old      source destination  resnik_score   jaccard  \n",
       "0        MP_0002160;  HP:0000033  MP:0004014      5.358735  0.153846  \n",
       "1        MP_0009198;  HP:0000033  MP:0020998      7.041238  0.166667  \n",
       "2        MP_0002160;  HP:0000033  MP:0020957      7.038994  0.153846  \n",
       "3        MP_0009198;  HP:0000033  MP:0020824      6.481542  0.363636  \n",
       "4        MP_0009198;  HP:0000033  MP:0002771      7.041238  0.153846  \n",
       "...              ...         ...         ...           ...       ...  \n",
       "2361470  MP_0000049;  MP:0004204  HP:0008609     10.195435  0.230769  \n",
       "2361471  MP_0000049;  MP:0030154  HP:0008609     10.195435  0.272727  \n",
       "2361472  MP_0000049;  MP:0030155  HP:0008609     10.195435  0.250000  \n",
       "2361473  MP_0000049;  MP:0030156  HP:0008609     10.195435  0.272727  \n",
       "2361474  MP_0000049;  MP:0030158  HP:0008609      8.446235  0.230769  \n",
       "\n",
       "[5833600 rows x 9 columns]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "intersected_scores1 = old_sim_pd.merge(new_sim_pd, how='inner',\n",
    "                                      left_on=['source_old', 'destination_old'], right_on=['source', 'destination'])\n",
    "intersected_scores2 = old_sim_pd.merge(new_sim_pd, how='inner',\n",
    "                                      left_on=['source_old', 'destination_old'], right_on=['destination', 'source'])\n",
    "intersected_scores = pd.concat([intersected_scores1, intersected_scores2])\n",
    "intersected_scores"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "ax1 = intersected_scores.plot.scatter(x='resnik_score',\n",
    "                                      y='resnik_score_old',\n",
    "                                      c='DarkBlue',\n",
    "                                      s=0.01\n",
    "                                     )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "ax1 = intersected_scores.plot.scatter(x='jaccard',\n",
    "                                      y='jaccard_old',\n",
    "                                      c='DarkBlue',\n",
    "                                      s=0.01\n",
    "                                     )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Which pairs have high scores in the old set and low scores in the new set?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>source_old</th>\n",
       "      <th>destination_old</th>\n",
       "      <th>jaccard_old</th>\n",
       "      <th>resnik_score_old</th>\n",
       "      <th>subsumer_old</th>\n",
       "      <th>source</th>\n",
       "      <th>destination</th>\n",
       "      <th>resnik_score</th>\n",
       "      <th>jaccard</th>\n",
       "      <th>resnik_diff</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>962688</th>\n",
       "      <td>HP:0002607</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>14.359005</td>\n",
       "      <td>MP_0003281;</td>\n",
       "      <td>HP:0002607</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>5.159422</td>\n",
       "      <td>0.300000</td>\n",
       "      <td>9.199583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>632818</th>\n",
       "      <td>HP:0040183</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>0.944444</td>\n",
       "      <td>14.359005</td>\n",
       "      <td>MP_0003281;</td>\n",
       "      <td>HP:0040183</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>5.159422</td>\n",
       "      <td>0.272727</td>\n",
       "      <td>9.199583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1327215</th>\n",
       "      <td>HP:0012700</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>0.529412</td>\n",
       "      <td>14.351005</td>\n",
       "      <td>HP_0012700;</td>\n",
       "      <td>HP:0012700</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>5.159422</td>\n",
       "      <td>0.333333</td>\n",
       "      <td>9.191583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1326917</th>\n",
       "      <td>HP:0012701</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>14.351005</td>\n",
       "      <td>HP_0012700;</td>\n",
       "      <td>HP:0012701</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>5.159422</td>\n",
       "      <td>0.300000</td>\n",
       "      <td>9.191583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>878712</th>\n",
       "      <td>HP:0012700</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>0.529412</td>\n",
       "      <td>14.351005</td>\n",
       "      <td>HP_0012700;</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>HP:0012700</td>\n",
       "      <td>5.159422</td>\n",
       "      <td>0.333333</td>\n",
       "      <td>9.191583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1071039</th>\n",
       "      <td>HP:0410188</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>0.282051</td>\n",
       "      <td>2.508834</td>\n",
       "      <td>MP_0008469;</td>\n",
       "      <td>HP:0410188</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>13.139874</td>\n",
       "      <td>0.300000</td>\n",
       "      <td>-10.631039</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>728809</th>\n",
       "      <td>HP:0410198</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>0.321429</td>\n",
       "      <td>2.508834</td>\n",
       "      <td>MP_0008469;</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>HP:0410198</td>\n",
       "      <td>13.139874</td>\n",
       "      <td>0.300000</td>\n",
       "      <td>-10.631039</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1069220</th>\n",
       "      <td>HP:0410198</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>0.321429</td>\n",
       "      <td>2.508834</td>\n",
       "      <td>MP_0008469;</td>\n",
       "      <td>HP:0410198</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>13.139874</td>\n",
       "      <td>0.300000</td>\n",
       "      <td>-10.631039</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2020913</th>\n",
       "      <td>HP:0002553</td>\n",
       "      <td>MP:0003761</td>\n",
       "      <td>0.222222</td>\n",
       "      <td>2.512377</td>\n",
       "      <td>HP_0000271;</td>\n",
       "      <td>HP:0002553</td>\n",
       "      <td>MP:0003761</td>\n",
       "      <td>13.833021</td>\n",
       "      <td>0.153846</td>\n",
       "      <td>-11.320644</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1361777</th>\n",
       "      <td>HP:0002553</td>\n",
       "      <td>MP:0003761</td>\n",
       "      <td>0.222222</td>\n",
       "      <td>2.512377</td>\n",
       "      <td>HP_0000271;</td>\n",
       "      <td>MP:0003761</td>\n",
       "      <td>HP:0002553</td>\n",
       "      <td>13.833021</td>\n",
       "      <td>0.153846</td>\n",
       "      <td>-11.320644</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5833600 rows × 10 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         source_old destination_old  jaccard_old  resnik_score_old  \\\n",
       "962688   HP:0002607      MP:0003281     1.000000         14.359005   \n",
       "632818   HP:0040183      MP:0003281     0.944444         14.359005   \n",
       "1327215  HP:0012700      MP:0003281     0.529412         14.351005   \n",
       "1326917  HP:0012701      MP:0003281     0.500000         14.351005   \n",
       "878712   HP:0012700      MP:0003281     0.529412         14.351005   \n",
       "...             ...             ...          ...               ...   \n",
       "1071039  HP:0410188      MP:0005562     0.282051          2.508834   \n",
       "728809   HP:0410198      MP:0005562     0.321429          2.508834   \n",
       "1069220  HP:0410198      MP:0005562     0.321429          2.508834   \n",
       "2020913  HP:0002553      MP:0003761     0.222222          2.512377   \n",
       "1361777  HP:0002553      MP:0003761     0.222222          2.512377   \n",
       "\n",
       "        subsumer_old      source destination  resnik_score   jaccard  \\\n",
       "962688   MP_0003281;  HP:0002607  MP:0003281      5.159422  0.300000   \n",
       "632818   MP_0003281;  HP:0040183  MP:0003281      5.159422  0.272727   \n",
       "1327215  HP_0012700;  HP:0012700  MP:0003281      5.159422  0.333333   \n",
       "1326917  HP_0012700;  HP:0012701  MP:0003281      5.159422  0.300000   \n",
       "878712   HP_0012700;  MP:0003281  HP:0012700      5.159422  0.333333   \n",
       "...              ...         ...         ...           ...       ...   \n",
       "1071039  MP_0008469;  HP:0410188  MP:0005562     13.139874  0.300000   \n",
       "728809   MP_0008469;  MP:0005562  HP:0410198     13.139874  0.300000   \n",
       "1069220  MP_0008469;  HP:0410198  MP:0005562     13.139874  0.300000   \n",
       "2020913  HP_0000271;  HP:0002553  MP:0003761     13.833021  0.153846   \n",
       "1361777  HP_0000271;  MP:0003761  HP:0002553     13.833021  0.153846   \n",
       "\n",
       "         resnik_diff  \n",
       "962688      9.199583  \n",
       "632818      9.199583  \n",
       "1327215     9.191583  \n",
       "1326917     9.191583  \n",
       "878712      9.191583  \n",
       "...              ...  \n",
       "1071039   -10.631039  \n",
       "728809    -10.631039  \n",
       "1069220   -10.631039  \n",
       "2020913   -11.320644  \n",
       "1361777   -11.320644  \n",
       "\n",
       "[5833600 rows x 10 columns]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "intersected_scores['resnik_diff'] = intersected_scores['resnik_score_old'] - intersected_scores['resnik_score']\n",
    "intersected_scores.sort_values(by='resnik_diff', ascending=False, inplace=True)\n",
    "intersected_scores"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "intersected_scores['resnik_diff'].plot(kind='hist')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "intersected_scores.boxplot(column='resnik_diff')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(-6.1562406229018185, -0.16051698226410993)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Find pairs +/- 1.5 iqr\n",
    "q3 = np.quantile(intersected_scores['resnik_diff'], 0.75)\n",
    "q1 = np.quantile(intersected_scores['resnik_diff'], 0.25)\n",
    "iqr = q3 - q1\n",
    "lower = q1 - 1.5 * iqr\n",
    "upper = q3 + 1.5 * iqr\n",
    "(lower, upper)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>source_old</th>\n",
       "      <th>destination_old</th>\n",
       "      <th>jaccard_old</th>\n",
       "      <th>resnik_score_old</th>\n",
       "      <th>subsumer_old</th>\n",
       "      <th>source</th>\n",
       "      <th>destination</th>\n",
       "      <th>resnik_score</th>\n",
       "      <th>jaccard</th>\n",
       "      <th>resnik_diff</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>962688</th>\n",
       "      <td>HP:0002607</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>14.359005</td>\n",
       "      <td>MP_0003281;</td>\n",
       "      <td>HP:0002607</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>5.159422</td>\n",
       "      <td>0.300000</td>\n",
       "      <td>9.199583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>632818</th>\n",
       "      <td>HP:0040183</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>0.944444</td>\n",
       "      <td>14.359005</td>\n",
       "      <td>MP_0003281;</td>\n",
       "      <td>HP:0040183</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>5.159422</td>\n",
       "      <td>0.272727</td>\n",
       "      <td>9.199583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1327215</th>\n",
       "      <td>HP:0012700</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>0.529412</td>\n",
       "      <td>14.351005</td>\n",
       "      <td>HP_0012700;</td>\n",
       "      <td>HP:0012700</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>5.159422</td>\n",
       "      <td>0.333333</td>\n",
       "      <td>9.191583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1326917</th>\n",
       "      <td>HP:0012701</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>14.351005</td>\n",
       "      <td>HP_0012700;</td>\n",
       "      <td>HP:0012701</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>5.159422</td>\n",
       "      <td>0.300000</td>\n",
       "      <td>9.191583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>878712</th>\n",
       "      <td>HP:0012700</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>0.529412</td>\n",
       "      <td>14.351005</td>\n",
       "      <td>HP_0012700;</td>\n",
       "      <td>MP:0003281</td>\n",
       "      <td>HP:0012700</td>\n",
       "      <td>5.159422</td>\n",
       "      <td>0.333333</td>\n",
       "      <td>9.191583</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1071039</th>\n",
       "      <td>HP:0410188</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>0.282051</td>\n",
       "      <td>2.508834</td>\n",
       "      <td>MP_0008469;</td>\n",
       "      <td>HP:0410188</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>13.139874</td>\n",
       "      <td>0.300000</td>\n",
       "      <td>-10.631039</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>728809</th>\n",
       "      <td>HP:0410198</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>0.321429</td>\n",
       "      <td>2.508834</td>\n",
       "      <td>MP_0008469;</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>HP:0410198</td>\n",
       "      <td>13.139874</td>\n",
       "      <td>0.300000</td>\n",
       "      <td>-10.631039</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1069220</th>\n",
       "      <td>HP:0410198</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>0.321429</td>\n",
       "      <td>2.508834</td>\n",
       "      <td>MP_0008469;</td>\n",
       "      <td>HP:0410198</td>\n",
       "      <td>MP:0005562</td>\n",
       "      <td>13.139874</td>\n",
       "      <td>0.300000</td>\n",
       "      <td>-10.631039</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2020913</th>\n",
       "      <td>HP:0002553</td>\n",
       "      <td>MP:0003761</td>\n",
       "      <td>0.222222</td>\n",
       "      <td>2.512377</td>\n",
       "      <td>HP_0000271;</td>\n",
       "      <td>HP:0002553</td>\n",
       "      <td>MP:0003761</td>\n",
       "      <td>13.833021</td>\n",
       "      <td>0.153846</td>\n",
       "      <td>-11.320644</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1361777</th>\n",
       "      <td>HP:0002553</td>\n",
       "      <td>MP:0003761</td>\n",
       "      <td>0.222222</td>\n",
       "      <td>2.512377</td>\n",
       "      <td>HP_0000271;</td>\n",
       "      <td>MP:0003761</td>\n",
       "      <td>HP:0002553</td>\n",
       "      <td>13.833021</td>\n",
       "      <td>0.153846</td>\n",
       "      <td>-11.320644</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>211443 rows × 10 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         source_old destination_old  jaccard_old  resnik_score_old  \\\n",
       "962688   HP:0002607      MP:0003281     1.000000         14.359005   \n",
       "632818   HP:0040183      MP:0003281     0.944444         14.359005   \n",
       "1327215  HP:0012700      MP:0003281     0.529412         14.351005   \n",
       "1326917  HP:0012701      MP:0003281     0.500000         14.351005   \n",
       "878712   HP:0012700      MP:0003281     0.529412         14.351005   \n",
       "...             ...             ...          ...               ...   \n",
       "1071039  HP:0410188      MP:0005562     0.282051          2.508834   \n",
       "728809   HP:0410198      MP:0005562     0.321429          2.508834   \n",
       "1069220  HP:0410198      MP:0005562     0.321429          2.508834   \n",
       "2020913  HP:0002553      MP:0003761     0.222222          2.512377   \n",
       "1361777  HP:0002553      MP:0003761     0.222222          2.512377   \n",
       "\n",
       "        subsumer_old      source destination  resnik_score   jaccard  \\\n",
       "962688   MP_0003281;  HP:0002607  MP:0003281      5.159422  0.300000   \n",
       "632818   MP_0003281;  HP:0040183  MP:0003281      5.159422  0.272727   \n",
       "1327215  HP_0012700;  HP:0012700  MP:0003281      5.159422  0.333333   \n",
       "1326917  HP_0012700;  HP:0012701  MP:0003281      5.159422  0.300000   \n",
       "878712   HP_0012700;  MP:0003281  HP:0012700      5.159422  0.333333   \n",
       "...              ...         ...         ...           ...       ...   \n",
       "1071039  MP_0008469;  HP:0410188  MP:0005562     13.139874  0.300000   \n",
       "728809   MP_0008469;  MP:0005562  HP:0410198     13.139874  0.300000   \n",
       "1069220  MP_0008469;  HP:0410198  MP:0005562     13.139874  0.300000   \n",
       "2020913  HP_0000271;  HP:0002553  MP:0003761     13.833021  0.153846   \n",
       "1361777  HP_0000271;  MP:0003761  HP:0002553     13.833021  0.153846   \n",
       "\n",
       "         resnik_diff  \n",
       "962688      9.199583  \n",
       "632818      9.199583  \n",
       "1327215     9.191583  \n",
       "1326917     9.191583  \n",
       "878712      9.191583  \n",
       "...              ...  \n",
       "1071039   -10.631039  \n",
       "728809    -10.631039  \n",
       "1069220   -10.631039  \n",
       "2020913   -11.320644  \n",
       "1361777   -11.320644  \n",
       "\n",
       "[211443 rows x 10 columns]"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Get pairs outside that range (lower than the lower bound, higher than the higher bound)\n",
    "divergent_scores = intersected_scores[(intersected_scores['resnik_diff'] <= lower) | (intersected_scores['resnik_diff'] >= upper)]\n",
    "divergent_scores"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "divergent_scores[(divergent_scores['resnik_diff'] <= lower)]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Overall, most pairs have higher Resnik similarity with new methods than with older ones. More than 211 thousand pairs meet our cutoff."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The top 5 differences in Resnik score are for pairs involving MP:0003281 (fecal incontinence). The top HP hit in the old scores is for HP:0002607 (bowel incontinence) and even includes the synonym \"Fecal incontinence\", so that appears to be a very good match. The low score in the new set may reflect how the ancestors of these terms are otherwise dissimilar: \n",
    "In HP, it looks like\n",
    "```\n",
    "All\n",
    " Phenotypic abnormality\n",
    "  Abnormality of the digestive system\n",
    "   Abnormality of digestive system physiology\n",
    "    Abnormal large intestine physiology\n",
    "     Bowel incontinence\n",
    "  Constitutional symptom\n",
    "   Impairment of activities of daily living\n",
    "    Impaired continence\n",
    "     Bowel incontinence\n",
    "```\n",
    "While in MP, we have\n",
    "```\n",
    "mammalian phenotype\n",
    " digestive/alimentary phenotype\n",
    "  abnormal digestive system physiology\n",
    "   abnormal defecation\n",
    "    fecal incontinence\n",
    "```\n",
    "\n",
    "The other HP terms in this top 5 are HP:0040183 (Encopresis; a child of HP:0002607), HP:0012700 (Abnormal large intestine physiology; see above), HP:0012701 (Bowel urgency; a child of HP:0012700), and HP:0012700 (that's a duplicate - shouldn't be there)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The *lowest* differences (i.e., the pairs where Resnik was previously very low and is now very high) involve HP:0410188 (Decreased glucose-6-phosphate dehydrogenase level in red blood cells), MP:0005562 (decreased mean corpuscular hemoglobin),  HP:0002553 (Highly arched eyebrow), and MP:0003761 (arched palate). Several of these are duplicates. In the latter pair, the lexical similarity of \"arched\" is a clue, as the path between these terms involves UPHENO:0041110 (arched anatomical entity). Related to PATO:0001594 (arched), this term is useful for describing a type of similarity but isn't meaningful for phenotype similarity, at least in this case. For the former pair of HP:0410188 vs MP:0005562, the phenotypes appear similar on the basis of involving blood. The path between them crosses through UPHENO:0034498 (decreased level of chemical entity in erythrocyte)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Let's see a few more of those pairs.\n",
    "divergent_scores.tail(20)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "HP:0002227 (White eyelashes) vs MP:0006242 (white pupil) is a pair with physical similarity in both anatomy and attribute. It has a relatively high score of 13.427556 in the new method, as does HP:0002226 (White eyebrow)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "HP:0006611 (Decreased number of sternal ossification centers) vs MP:0010869 (decreased bone trabecula number): both involve bone, certainly, though the high similarity score may be misrepresenting what is otherwise a dissimilar phenotype."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## OAK and Graph viz"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Time to visualize, courtesy of OAK, obographviz, and KG-Phenio."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Set up the alias first for notebook purposes\n",
    "%alias phenio runoak -i sqlite:obo:phenio"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "phenio info \"HP:0006611\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "phenio viz -p i,p \"HP:0006611\" \"MP:0010869\" -o output2.png --no-view"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Retrieve all UPHENO MRCAs for divergent HP-MP pairs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('HP:0002607', 'MP:0003281')\n",
      "('HP:0040183', 'MP:0003281')\n",
      "('HP:0012700', 'MP:0003281')\n",
      "('HP:0012701', 'MP:0003281')\n",
      "('MP:0003281', 'HP:0012700')\n",
      "('HP:0012702', 'MP:0003281')\n",
      "('HP:0012423', 'MP:0003281')\n",
      "('HP:0100819', 'MP:0003319')\n",
      "('HP:0033279', 'MP:0003319')\n",
      "('HP:0100744', 'MP:0030883')\n"
     ]
    }
   ],
   "source": [
    "div_pairs = list(zip(divergent_scores[\"source\"], divergent_scores[\"destination\"]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "test1 = \"HP:0006611\"\n",
    "test2 = \"MP:0010869\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The following will retrieve all shortest paths for all the divergent pairs. Each requires processing once complete to:\n",
    "* Remove header lines\n",
    "* Select only the UPHENO part of the path\n",
    "* Find the UPHENO midpoint for each\n",
    "\n",
    "Doing the whole set this way takes intractably long, so let's try the last 1000."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [],
   "source": [
    "all_spaths = {}\n",
    "count = 1000"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 1.4 s, sys: 6.63 s, total: 8.03 s\n",
      "Wall time: 47min 44s\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "{('HP:0002607', 'MP:0003281'): 'UPHENO:0002320',\n",
       " ('HP:0040183', 'MP:0003281'): 'UPHENO:0002320',\n",
       " ('HP:0012700', 'MP:0003281'): 'UPHENO:0002320',\n",
       " ('HP:0012701', 'MP:0003281'): 'UPHENO:0002320',\n",
       " ('MP:0003281', 'HP:0012700'): 'UPHENO:0002332',\n",
       " ('HP:0012702', 'MP:0003281'): 'UPHENO:0002320',\n",
       " ('HP:0012423', 'MP:0003281'): 'UPHENO:0002320',\n",
       " ('HP:0100819', 'MP:0003319'): 'UPHENO:0019528',\n",
       " ('HP:0033279', 'MP:0003319'): 'UPHENO:0019528',\n",
       " ('HP:0100744', 'MP:0030883'): 'UPHENO:0020736',\n",
       " ('MP:0012399', 'HP:0011937'): 'UPHENO:0081318',\n",
       " ('HP:0008755', 'MP:0031308'): 'UPHENO:0004536',\n",
       " ('HP:0002786', 'MP:0013753'): 'UPHENO:0004536',\n",
       " ('HP:0033990', 'MP:0013753'): 'UPHENO:0004536',\n",
       " ('HP:0033979', 'MP:0013753'): 'UPHENO:0004536',\n",
       " ('HP:0033979', 'MP:0031308'): 'UPHENO:0004536',\n",
       " ('HP:0002786', 'MP:0031308'): 'UPHENO:0004536',\n",
       " ('HP:0008755', 'MP:0013753'): 'UPHENO:0004536',\n",
       " ('HP:0002779', 'MP:0013753'): 'UPHENO:0004536',\n",
       " ('HP:0033990', 'MP:0031308'): 'UPHENO:0004536',\n",
       " ('HP:0002779', 'MP:0031308'): 'UPHENO:0004536',\n",
       " ('HP:0005939', 'MP:0010822'): 'UPHENO:0004536',\n",
       " ('HP:0006522', 'MP:0010822'): 'UPHENO:0004536',\n",
       " ('HP:0002108', 'MP:0010822'): 'UPHENO:0004536',\n",
       " ('HP:0002107', 'MP:0010822'): 'UPHENO:0004536',\n",
       " ('HP:0004876', 'MP:0010822'): 'UPHENO:0004536',\n",
       " ('HP:0500223', 'MP:0030724'): 'UPHENO:0048724',\n",
       " ('MP:0030724', 'HP:0500227'): 'UPHENO:0048729',\n",
       " ('MP:0030724', 'HP:0500197'): 'UPHENO:0048694',\n",
       " ('HP:0500191', 'MP:0030724'): 'UPHENO:0048723',\n",
       " ('MP:0030724', 'HP:0500191'): 'UPHENO:0048723',\n",
       " ('MP:0030724', 'HP:0500230'): 'UPHENO:0048696',\n",
       " ('MP:0030724', 'HP:0500210'): 'UPHENO:0048728',\n",
       " ('MP:0030724', 'HP:0500208'): 'UPHENO:0048727',\n",
       " ('HP:0500193', 'MP:0030724'): 'UPHENO:0048735',\n",
       " ('MP:0030724', 'HP:0500203'): 'UPHENO:0048726',\n",
       " ('MP:0030724', 'HP:0500223'): 'UPHENO:0048724',\n",
       " ('MP:0030724', 'HP:0500200'): 'UPHENO:0048732',\n",
       " ('MP:0030724', 'HP:0500244'): 'UPHENO:0048702',\n",
       " ('HP:0500200', 'MP:0030724'): 'UPHENO:0048732',\n",
       " ('MP:0030724', 'HP:0500193'): 'UPHENO:0048735',\n",
       " ('HP:0500227', 'MP:0030724'): 'UPHENO:0048729',\n",
       " ('HP:0500212', 'MP:0030724'): 'UPHENO:0048722',\n",
       " ('HP:0500210', 'MP:0030724'): 'UPHENO:0048728',\n",
       " ('HP:0500244', 'MP:0030724'): 'UPHENO:0048702',\n",
       " ('HP:0500230', 'MP:0030724'): 'UPHENO:0048696',\n",
       " ('MP:0030724', 'HP:0500212'): 'UPHENO:0048722',\n",
       " ('HP:0500208', 'MP:0030724'): 'UPHENO:0048727',\n",
       " ('HP:0500203', 'MP:0030724'): 'UPHENO:0048726',\n",
       " ('HP:0500197', 'MP:0030724'): 'UPHENO:0048694',\n",
       " ('MP:0030308', 'HP:0032045'): 'UPHENO:0081566',\n",
       " ('HP:0032045', 'MP:0030308'): 'UPHENO:0019224',\n",
       " ('HP:0410246', 'MP:0002496'): 'UPHENO:0002263',\n",
       " ('MP:0003749', 'HP:0002714'): 'UPHENO:0076800',\n",
       " ('HP:0002714', 'MP:0003749'): 'UPHENO:0020734',\n",
       " ('HP:0032097', 'MP:0030600'): 'UPHENO:0051804',\n",
       " ('HP:0032096', 'MP:0030600'): 'UPHENO:0051804',\n",
       " ('HP:0032098', 'MP:0030600'): 'UPHENO:0051804',\n",
       " ('MP:0010531', 'HP:0004296'): 'UPHENO:0001005',\n",
       " ('HP:0007902', 'MP:0006202'): 'UPHENO:0075679',\n",
       " ('HP:0030848', 'MP:0006136'): 'UPHENO:0020351',\n",
       " ('HP:0030970', 'MP:0006136'): 'UPHENO:0020351',\n",
       " ('HP:0030847', 'MP:0006136'): 'UPHENO:0020351',\n",
       " ('HP:0030849', 'MP:0006136'): 'UPHENO:0020351',\n",
       " ('HP:0002639', 'MP:0006136'): 'UPHENO:0020351',\n",
       " ('HP:0030846', 'MP:0006136'): 'UPHENO:0020351',\n",
       " ('HP:0030969', 'MP:0006136'): 'UPHENO:0020351',\n",
       " ('HP:0025322', 'MP:0006136'): 'UPHENO:0020351',\n",
       " ('MP:0006136', 'HP:0030969'): 'UPHENO:0003964',\n",
       " ('MP:0006136', 'HP:0030846'): 'UPHENO:0002343',\n",
       " ('MP:0009860', 'HP:0009741'): 'UPHENO:0076779',\n",
       " ('HP:0009741', 'MP:0009860'): 'UPHENO:0019487',\n",
       " ('HP:0033367', 'MP:0012549'): 'UPHENO:0076286',\n",
       " ('HP:0012418', 'MP:0012549'): 'UPHENO:0076286',\n",
       " ('MP:0012549', 'HP:0012418'): 'UPHENO:0079822',\n",
       " ('HP:0030874', 'MP:0012549'): 'UPHENO:0076286',\n",
       " ('HP:0034038', 'MP:0012549'): 'UPHENO:0076286',\n",
       " ('MP:0012549', 'HP:0410213'): 'UPHENO:0077824',\n",
       " ('HP:0410213', 'MP:0012549'): 'UPHENO:0076286',\n",
       " ('HP:0001126', 'MP:0005242'): 'UPHENO:0021749',\n",
       " ('HP:0033738', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('MP:0001948', 'HP:0000076'): 'UPHENO:0002332',\n",
       " ('HP:0000076', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('HP:0033741', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('HP:0033733', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('HP:0033735', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('HP:0033740', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('HP:0033739', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('HP:0033737', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('HP:0033734', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('HP:0033736', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('HP:0033742', 'MP:0001948'): 'UPHENO:0002332',\n",
       " ('MP:0021008', 'HP:0007132'): 'UPHENO:0020914',\n",
       " ('HP:0007132', 'MP:0021008'): 'UPHENO:0081435',\n",
       " ('MP:0009792', 'HP:0000964'): 'UPHENO:0001005',\n",
       " ('HP:0010728', 'MP:0009773'): 'UPHENO:0003020',\n",
       " ('MP:0009773', 'HP:0010728'): 'UPHENO:0003020',\n",
       " ('HP:0006693', 'MP:0010535'): 'UPHENO:0076781',\n",
       " ('HP:0011556', 'MP:0010434'): 'UPHENO:0021062',\n",
       " ('HP:0000539', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('HP:0025573', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('HP:0000540', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('HP:0012803', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('MP:0006151', 'HP:0012803'): 'UPHENO:0003020',\n",
       " ('HP:0031728', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('HP:0000545', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('MP:0006151', 'HP:0000539'): 'UPHENO:0003020',\n",
       " ('HP:0008499', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('HP:0500066', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('MP:0006151', 'HP:0000540'): 'UPHENO:0003020',\n",
       " ('HP:0031729', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('HP:0011003', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('HP:0500042', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('MP:0006151', 'HP:0000545'): 'UPHENO:0003020',\n",
       " ('HP:0031624', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('HP:0031730', 'MP:0006151'): 'UPHENO:0020998',\n",
       " ('HP:0006089', 'MP:0013371'): 'UPHENO:0002405',\n",
       " ('HP:0006089', 'MP:0013372'): 'UPHENO:0002405',\n",
       " ('MP:0030843', 'HP:0030310'): 'UPHENO:0076767',\n",
       " ('MP:0030860', 'HP:0030311'): 'UPHENO:0076767',\n",
       " ('MP:0030883', 'HP:0030310'): 'UPHENO:0076767',\n",
       " ('HP:0030310', 'MP:0030843'): 'UPHENO:0020968',\n",
       " ('HP:0031610', 'MP:0030883'): 'UPHENO:0020736',\n",
       " ('HP:0030311', 'MP:0030860'): 'UPHENO:0020968',\n",
       " ('HP:0031610', 'MP:0030843'): 'UPHENO:0020968',\n",
       " ('HP:0003835', 'MP:0030883'): 'UPHENO:0020736',\n",
       " ('HP:0003994', 'MP:0030883'): 'UPHENO:0020736',\n",
       " ('HP:0003835', 'MP:0030843'): 'UPHENO:0020968',\n",
       " ('HP:0002827', 'MP:0030860'): 'UPHENO:0076767',\n",
       " ('MP:0030883', 'HP:0003994'): 'UPHENO:0076767',\n",
       " ('HP:0030310', 'MP:0030883'): 'UPHENO:0020736',\n",
       " ('MP:0030860', 'HP:0002827'): 'UPHENO:0076767',\n",
       " ('HP:0003834', 'MP:0030883'): 'UPHENO:0079872',\n",
       " ('HP:0001374', 'MP:0030860'): 'UPHENO:0076767',\n",
       " ('MP:0030883', 'HP:0031610'): 'UPHENO:0076767',\n",
       " ('MP:0030883', 'HP:0003835'): 'UPHENO:0076767',\n",
       " ('HP:0002283', 'MP:0021008'): 'UPHENO:0012541',\n",
       " ('HP:0002283', 'MP:0021011'): 'UPHENO:0012541',\n",
       " ('MP:0021008', 'HP:0002059'): 'UPHENO:0012541',\n",
       " ('MP:0021011', 'HP:0002283'): 'UPHENO:0012541',\n",
       " ('MP:0021008', 'HP:0002283'): 'UPHENO:0012541',\n",
       " ('MP:0021008', 'HP:0006892'): 'UPHENO:0012541',\n",
       " ('HP:0002506', 'MP:0021008'): 'UPHENO:0012541',\n",
       " ('HP:0002506', 'MP:0021011'): 'UPHENO:0012541',\n",
       " ('MP:0021011', 'HP:0002506'): 'UPHENO:0012541',\n",
       " ('MP:0021008', 'HP:0002506'): 'UPHENO:0012541',\n",
       " ('MP:0021011', 'HP:0002059'): 'UPHENO:0012541',\n",
       " ('MP:0021008', 'HP:0007058'): 'UPHENO:0012541',\n",
       " ('HP:0002059', 'MP:0021011'): 'UPHENO:0012541',\n",
       " ('HP:0006892', 'MP:0021008'): 'UPHENO:0012541',\n",
       " ('MP:0021011', 'HP:0006892'): 'UPHENO:0012541',\n",
       " ('HP:0006892', 'MP:0021011'): 'UPHENO:0012541',\n",
       " ('HP:0002059', 'MP:0021008'): 'UPHENO:0012541',\n",
       " ('HP:0007058', 'MP:0021008'): 'UPHENO:0012541',\n",
       " ('HP:0007058', 'MP:0021011'): 'UPHENO:0012541',\n",
       " ('MP:0021011', 'HP:0007058'): 'UPHENO:0012541',\n",
       " ('HP:0100098', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('HP:0010595', 'MP:0030792'): 'UPHENO:0020967',\n",
       " ('HP:0100091', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('HP:0100095', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('HP:0100097', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('HP:0100097', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('HP:0100098', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('HP:0010593', 'MP:0030792'): 'UPHENO:0020967',\n",
       " ('HP:0010594', 'MP:0030788'): 'UPHENO:0020967',\n",
       " ('HP:0100094', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('HP:0010593', 'MP:0030788'): 'UPHENO:0020967',\n",
       " ('MP:0030788', 'HP:0006400'): 'UPHENO:0003070',\n",
       " ('HP:0010595', 'MP:0030788'): 'UPHENO:0020967',\n",
       " ('HP:0006400', 'MP:0030788'): 'UPHENO:0003070',\n",
       " ('HP:0010594', 'MP:0030792'): 'UPHENO:0020967',\n",
       " ('MP:0030792', 'HP:0006500'): 'UPHENO:0003070',\n",
       " ('HP:0010056', 'MP:0030788'): 'UPHENO:0020967',\n",
       " ('HP:0100091', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('HP:0010329', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('HP:0010323', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0006400'): 'UPHENO:0003070',\n",
       " ('MP:0030788', 'HP:0005715'): 'UPHENO:0003070',\n",
       " ('MP:0030792', 'HP:0005715'): 'UPHENO:0003070',\n",
       " ('HP:0005715', 'MP:0030788'): 'UPHENO:0003070',\n",
       " ('HP:0010631', 'MP:0030792'): 'UPHENO:0020967',\n",
       " ('HP:0010323', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('HP:0010631', 'MP:0030788'): 'UPHENO:0020967',\n",
       " ('HP:0100092', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('HP:0100094', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('HP:0100092', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('MP:0030788', 'HP:0006500'): 'UPHENO:0003070',\n",
       " ('HP:0010335', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('HP:0010329', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('HP:0010335', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('MP:0030788', 'HP:0010631'): 'UPHENO:0003070',\n",
       " ('HP:0010341', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('HP:0006500', 'MP:0030788'): 'UPHENO:0020967',\n",
       " ('HP:0006500', 'MP:0030792'): 'UPHENO:0020967',\n",
       " ('HP:0010160', 'MP:0030792'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0010631'): 'UPHENO:0003070',\n",
       " ('HP:0010056', 'MP:0030792'): 'UPHENO:0020967',\n",
       " ('HP:0100095', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('HP:0010341', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('HP:0010160', 'MP:0030788'): 'UPHENO:0076727',\n",
       " ('HP:0009737', 'MP:0010305'): 'UPHENO:0021474',\n",
       " ('HP:0009737', 'MP:0006232'): 'UPHENO:0021474',\n",
       " ('MP:0003129', 'HP:0012621'): 'UPHENO:0002832',\n",
       " ('HP:0012621', 'MP:0003129'): 'UPHENO:0020738',\n",
       " ('HP:0040050', 'MP:0000425'): 'UPHENO:0026385',\n",
       " ('MP:0000425', 'HP:0040050'): 'UPHENO:0075997',\n",
       " ('MP:0000425', 'HP:0007776'): 'UPHENO:0075997',\n",
       " ('HP:0007776', 'MP:0000425'): 'UPHENO:0026385',\n",
       " ('HP:0100515', 'MP:0011535'): 'UPHENO:0002332',\n",
       " ('MP:0013371', 'HP:0025276'): 'UPHENO:0001005',\n",
       " ('MP:0013372', 'HP:0025276'): 'UPHENO:0001005',\n",
       " ('MP:0013370', 'HP:0025276'): 'UPHENO:0001005',\n",
       " ('MP:0010535', 'HP:0006693'): 'UPHENO:0015324',\n",
       " ('HP:0006976', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('HP:0006929', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('HP:0006846', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('HP:0007030', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('HP:0007069', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('HP:0002480', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0007030'): 'UPHENO:0002433',\n",
       " ('HP:0007335', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0200134'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0006846'): 'UPHENO:0002433',\n",
       " ('HP:0001298', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0007335'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0007239'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0006929'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0006976'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0006789'): 'UPHENO:0002433',\n",
       " ('HP:0006965', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0002480'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0002448'): 'UPHENO:0002433',\n",
       " ('HP:0006789', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0007105'): 'UPHENO:0002433',\n",
       " ('HP:0007105', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('MP:0013806', 'HP:0001298'): 'UPHENO:0002433',\n",
       " ('HP:0007111', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('HP:0007239', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('HP:0002448', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('HP:0200134', 'MP:0013806'): 'UPHENO:0002433',\n",
       " ('HP:0000643', 'MP:0030143'): 'UPHENO:0002320',\n",
       " ('MP:0030143', 'HP:0000643'): 'UPHENO:0002332',\n",
       " ('MP:0006244', 'HP:0000643'): 'UPHENO:0002332',\n",
       " ('HP:0000643', 'MP:0006244'): 'UPHENO:0002332',\n",
       " ('HP:0006490', 'MP:0030790'): 'UPHENO:0020967',\n",
       " ('HP:0030293', 'MP:0030790'): 'UPHENO:0020967',\n",
       " ('MP:0030790', 'HP:0030291'): 'UPHENO:0003070',\n",
       " ('HP:0030293', 'MP:0030794'): 'UPHENO:0020967',\n",
       " ('HP:0030291', 'MP:0030789'): 'UPHENO:0020967',\n",
       " ('MP:0030789', 'HP:0006490'): 'UPHENO:0003070',\n",
       " ('HP:0030293', 'MP:0030789'): 'UPHENO:0020967',\n",
       " ('HP:0006490', 'MP:0030793'): 'UPHENO:0020967',\n",
       " ('HP:0030291', 'MP:0030790'): 'UPHENO:0020967',\n",
       " ('HP:0030291', 'MP:0030793'): 'UPHENO:0020967',\n",
       " ('HP:0006490', 'MP:0030789'): 'UPHENO:0020967',\n",
       " ('HP:0030291', 'MP:0030794'): 'UPHENO:0020967',\n",
       " ('MP:0030790', 'HP:0006490'): 'UPHENO:0003070',\n",
       " ('MP:0030794', 'HP:0030291'): 'UPHENO:0003070',\n",
       " ('MP:0030793', 'HP:0030291'): 'UPHENO:0003070',\n",
       " ('MP:0030793', 'HP:0006490'): 'UPHENO:0003070',\n",
       " ('MP:0030794', 'HP:0006490'): 'UPHENO:0003070',\n",
       " ('HP:0006490', 'MP:0030794'): 'UPHENO:0020967',\n",
       " ('MP:0030789', 'HP:0030291'): 'UPHENO:0003070',\n",
       " ('HP:0030293', 'MP:0030793'): 'UPHENO:0020967',\n",
       " ('HP:0000643', 'MP:0030144'): 'UPHENO:0002320',\n",
       " ('HP:0000643', 'MP:0030147'): 'UPHENO:0002320',\n",
       " ('MP:0030144', 'HP:0000643'): 'UPHENO:0002332',\n",
       " ('MP:0030147', 'HP:0000643'): 'UPHENO:0002332',\n",
       " ('MP:0009274', 'HP:0000501'): 'UPHENO:0001005',\n",
       " ('HP:0006089', 'MP:0013369'): 'UPHENO:0002405',\n",
       " ('HP:0006089', 'MP:0013370'): 'UPHENO:0002405',\n",
       " ('HP:0007410', 'MP:0013369'): 'UPHENO:0002405',\n",
       " ('HP:0007410', 'MP:0013368'): 'UPHENO:0002405',\n",
       " ('HP:0006089', 'MP:0013368'): 'UPHENO:0002405',\n",
       " ('HP:0011581', 'MP:0010427'): 'UPHENO:0021062',\n",
       " ('HP:0010126', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100116'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100121'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100156'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100117'): 'UPHENO:0076806',\n",
       " ('HP:0100100', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100117'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100116'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100233'): 'UPHENO:0076806',\n",
       " ('HP:0010163', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010163', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100226', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100115'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100232'): 'UPHENO:0076806',\n",
       " ('HP:0100226', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100053'): 'UPHENO:0076806',\n",
       " ('HP:0100046', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100046', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100115'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100113'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100087'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100087'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100156'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100170'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100121'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010133'): 'UPHENO:0076806',\n",
       " ('HP:0100115', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100115', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0010141'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100089'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100053'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010134'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010134'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010133'): 'UPHENO:0076806',\n",
       " ('HP:0100111', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100048'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100068'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100048'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100047'): 'UPHENO:0076806',\n",
       " ('HP:0100107', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0010132'): 'UPHENO:0076806',\n",
       " ('HP:0100104', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100104', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100118', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100170'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010132'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100113'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100112'): 'UPHENO:0076806',\n",
       " ('HP:0100116', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100057'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010137'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100195'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100195'): 'UPHENO:0076806',\n",
       " ('HP:0010166', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0010140'): 'UPHENO:0076806',\n",
       " ('HP:0010166', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100104'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100086'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100200'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010166'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010166'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010165'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100086'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100200'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010121'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010121'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100104'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100201'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100201'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100215'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010131'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100107'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100194'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100107'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010136'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100057'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100179'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100179'): 'UPHENO:0076806',\n",
       " ('HP:0010123', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010123', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100068'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100159'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100159'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100112'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100111'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100110'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100232'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100111'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100046'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010137'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100046'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100199'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100199'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010131'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100110'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100165'): 'UPHENO:0076806',\n",
       " ('HP:0100158', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100068', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0010147'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100105'): 'UPHENO:0076806',\n",
       " ('HP:0100195', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100150', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100150', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100195', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100124'): 'UPHENO:0083639',\n",
       " ('MP:0030788', 'HP:0100124'): 'UPHENO:0083639',\n",
       " ('MP:0030788', 'HP:0100106'): 'UPHENO:0076806',\n",
       " ('HP:0100159', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100126'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100192'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100192'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100139'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100126'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100139'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100106'): 'UPHENO:0076806',\n",
       " ('HP:0100201', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100201', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0010147'): 'UPHENO:0076806',\n",
       " ('HP:0100059', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100059', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100105'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100079'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100079'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100148'): 'UPHENO:0076806',\n",
       " ('HP:0010143', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100132'): 'UPHENO:0083639',\n",
       " ('MP:0030788', 'HP:0100132'): 'UPHENO:0083639',\n",
       " ('MP:0030788', 'HP:0100148'): 'UPHENO:0076806',\n",
       " ('HP:0100056', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100052'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010126'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010163'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010163'): 'UPHENO:0076806',\n",
       " ('HP:0100159', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100052'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010135'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100054'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100137'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100137'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100054'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010126'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010136'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010128'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010128'): 'UPHENO:0076806',\n",
       " ('HP:0100206', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100139', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100117', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100139', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100126', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100147'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100070'): 'UPHENO:0076806',\n",
       " ('HP:0100048', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100070'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100233'): 'UPHENO:0076806',\n",
       " ('HP:0100048', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100158', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100088'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100088'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100118'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100069'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100118'): 'UPHENO:0076806',\n",
       " ('HP:0100081', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100056'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100056'): 'UPHENO:0076806',\n",
       " ('HP:0100081', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100150'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100150'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100089'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100165'): 'UPHENO:0076806',\n",
       " ('HP:0100126', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100147'): 'UPHENO:0076806',\n",
       " ('HP:0100217', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100067'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100045'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100045'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010146'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010146'): 'UPHENO:0076806',\n",
       " ('HP:0010139', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0010139', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0010135'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100074'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100123'): 'UPHENO:0083639',\n",
       " ('MP:0030792', 'HP:0100067'): 'UPHENO:0076806',\n",
       " ('HP:0100217', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100123'): 'UPHENO:0083639',\n",
       " ('MP:0030788', 'HP:0100074'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100059'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100059'): 'UPHENO:0076806',\n",
       " ('HP:0100206', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100058'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100227'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100145'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100145'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010140'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100203'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100136'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100085'): 'UPHENO:0076806',\n",
       " ('HP:0100045', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0010170'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010117'): 'UPHENO:0076806',\n",
       " ('HP:0100192', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0010164'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010164'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010117'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010138'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010138'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010116'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100085'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010115'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010115'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010122'): 'UPHENO:0076806',\n",
       " ('HP:0100161', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100161', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010164', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010164', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0010122'): 'UPHENO:0076806',\n",
       " ('HP:0100192', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100045', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0010170'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010123'): 'UPHENO:0076806',\n",
       " ('HP:0100088', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100102'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100081'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100223'): 'UPHENO:0076806',\n",
       " ('HP:0100105', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100173', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100054', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100054', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100078'): 'UPHENO:0076806',\n",
       " ('HP:0100088', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100167'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100182'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100167'): 'UPHENO:0076806',\n",
       " ('HP:0100234', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100173', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100078'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100182'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010172'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010172'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010171'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010171'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010123'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010114'): 'UPHENO:0076806',\n",
       " ('HP:0100228', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100225', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100181'): 'UPHENO:0076806',\n",
       " ('HP:0100184', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100193', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100075'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100075'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100076'): 'UPHENO:0076806',\n",
       " ('HP:0010145', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010128', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100076'): 'UPHENO:0076806',\n",
       " ('HP:0100225', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100217'): 'UPHENO:0076806',\n",
       " ('HP:0010126', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100173'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100173'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100171'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100171'): 'UPHENO:0076806',\n",
       " ('HP:0100056', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100161'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100161'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100160'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010114'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100217'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010143'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100206'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010143'): 'UPHENO:0076806',\n",
       " ('HP:0010172', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010138', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010172', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100225'): 'UPHENO:0076806',\n",
       " ('HP:0100234', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100064'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100190'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100190'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100206'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100181'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100225'): 'UPHENO:0076806',\n",
       " ('HP:0010128', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100065'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100063'): 'UPHENO:0076806',\n",
       " ('HP:0100184', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100065'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100063'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100064'): 'UPHENO:0076806',\n",
       " ('HP:0100193', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0100212'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100212'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100081'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100184'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010139'): 'UPHENO:0076806',\n",
       " ('HP:0010114', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010114', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0010115', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010115', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100193'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100193'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100166'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100166'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100158'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100204'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100204'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100158'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100215'): 'UPHENO:0076806',\n",
       " ('HP:0100167', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0010117', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100167', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100057', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0010117', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0010139'): 'UPHENO:0076806',\n",
       " ('HP:0100057', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0010129'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010145'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010145'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100228'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100203'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100102'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010130'): 'UPHENO:0076806',\n",
       " ('HP:0100068', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0010142'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010142'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010127'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100234'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100234'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100228'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100101'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010127'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010124'): 'UPHENO:0076806',\n",
       " ('HP:0100106', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030788', 'HP:0010144'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010129'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100214'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100214'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0010124'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0100101'): 'UPHENO:0076806',\n",
       " ('HP:0010148', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100067', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100067', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010142', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100079', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100184'): 'UPHENO:0076806',\n",
       " ('HP:0100078', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100078', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010140', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100076', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100076', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010140', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0100223'): 'UPHENO:0076806',\n",
       " ('HP:0010129', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0010137', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0010129', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100070', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100070', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010142', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100080'): 'UPHENO:0076806',\n",
       " ('HP:0010144', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100065', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100065', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010148', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100079', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100228', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0100226'): 'UPHENO:0076806',\n",
       " ('HP:0100087', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100100'): 'UPHENO:0076806',\n",
       " ('HP:0010124', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010124', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030792', 'HP:0010144'): 'UPHENO:0076806',\n",
       " ('MP:0030788', 'HP:0100100'): 'UPHENO:0076806',\n",
       " ('HP:0100089', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0100089', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010127', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010137', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0100087', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('HP:0010131', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010131', 'MP:0030792'): 'UPHENO:0020575',\n",
       " ('MP:0030788', 'HP:0010148'): 'UPHENO:0076806',\n",
       " ('MP:0030792', 'HP:0010148'): 'UPHENO:0076806',\n",
       " ('HP:0010132', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('HP:0010133', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100226'): 'UPHENO:0076806',\n",
       " ('HP:0010134', 'MP:0030788'): 'UPHENO:0081558',\n",
       " ('MP:0030792', 'HP:0100136'): 'UPHENO:0076806',\n",
       " ('MP:0003518', 'HP:0000060'): 'UPHENO:0081790',\n",
       " ('HP:0012815', 'MP:0003518'): 'UPHENO:0020427',\n",
       " ('HP:0000060', 'MP:0003518'): 'UPHENO:0069435',\n",
       " ('HP:0009119', 'MP:0030382'): 'UPHENO:0004536',\n",
       " ('MP:0030382', 'HP:0009119'): 'UPHENO:0080178',\n",
       " ('MP:0030093', 'HP:0009119'): 'UPHENO:0080178',\n",
       " ('MP:0030093', 'HP:0005453'): 'UPHENO:0004536',\n",
       " ('MP:0030382', 'HP:0005453'): 'UPHENO:0004536',\n",
       " ('HP:0005453', 'MP:0030382'): 'UPHENO:0004536',\n",
       " ('HP:0009120', 'MP:0030382'): 'UPHENO:0020974',\n",
       " ('MP:0003591', 'HP:0000809'): 'UPHENO:0002832',\n",
       " ('MP:0011489', 'HP:0000809'): 'UPHENO:0002832',\n",
       " ('HP:0000809', 'MP:0003591'): 'UPHENO:0020738',\n",
       " ('HP:0000809', 'MP:0011489'): 'UPHENO:0019880',\n",
       " ('HP:0033363', 'MP:0010923'): 'UPHENO:0021276',\n",
       " ('HP:0033376', 'MP:0010923'): 'UPHENO:0021276',\n",
       " ('MP:0010923', 'HP:0033376'): 'UPHENO:0004536',\n",
       " ('HP:0033006', 'MP:0010923'): 'UPHENO:0021276',\n",
       " ('HP:4000060', 'MP:0010923'): 'UPHENO:0021276',\n",
       " ('MP:0010923', 'HP:0033208'): 'UPHENO:0004536',\n",
       " ('MP:0010923', 'HP:0033137'): 'UPHENO:0004536',\n",
       " ('HP:0033328', 'MP:0010923'): 'UPHENO:0021276',\n",
       " ('MP:0010923', 'HP:4000060'): 'UPHENO:0004536',\n",
       " ('MP:0010923', 'HP:0033241'): 'UPHENO:0004536',\n",
       " ('HP:0033137', 'MP:0010923'): 'UPHENO:0021276',\n",
       " ('HP:0033246', 'MP:0010923'): 'UPHENO:0021276',\n",
       " ('MP:0010923', 'HP:0033245'): 'UPHENO:0004536',\n",
       " ('MP:0010923', 'HP:0033246'): 'UPHENO:0004536',\n",
       " ('HP:0033245', 'MP:0010923'): 'UPHENO:0021276',\n",
       " ('HP:0033241', 'MP:0010923'): 'UPHENO:0021276',\n",
       " ('MP:0010923', 'HP:0033363'): 'UPHENO:0004536',\n",
       " ('MP:0010923', 'HP:0033006'): 'UPHENO:0004536',\n",
       " ('MP:0010923', 'HP:0033328'): 'UPHENO:0004536',\n",
       " ('HP:0033208', 'MP:0010923'): 'UPHENO:0021276',\n",
       " ('MP:0008159', 'HP:0001783'): 'UPHENO:0020817',\n",
       " ('MP:0008159', 'HP:0010068'): 'UPHENO:0006276',\n",
       " ('MP:0008159', 'HP:0006200'): 'UPHENO:0020651',\n",
       " ('MP:0008159', 'HP:0010027'): 'UPHENO:0006432',\n",
       " ('MP:0008159', 'HP:0003981'): 'UPHENO:0020869',\n",
       " ('MP:0008159', 'HP:0001230'): 'UPHENO:0020648',\n",
       " ('MP:0008159', 'HP:0003886'): 'UPHENO:0020980',\n",
       " ('MP:0008159', 'HP:0003993'): 'UPHENO:0020868',\n",
       " ('MP:0008159', 'HP:0006095'): 'UPHENO:0020651',\n",
       " ('MP:0008159', 'HP:0006009'): 'UPHENO:0020651',\n",
       " ('MP:0004114', 'HP:0006681'): 'UPHENO:0002816',\n",
       " ('HP:0006681', 'MP:0004114'): 'UPHENO:0020158',\n",
       " ('MP:0010533', 'HP:0006681'): 'UPHENO:0002816',\n",
       " ('HP:0006681', 'MP:0010533'): 'UPHENO:0020706',\n",
       " ('MP:0030855', 'HP:0005004'): 'UPHENO:0041360',\n",
       " ('MP:0030855', 'HP:0003998'): 'UPHENO:0079872',\n",
       " ('HP:0003998', 'MP:0030855'): 'UPHENO:0079872',\n",
       " ('HP:0003996', 'MP:0030855'): 'UPHENO:0079872',\n",
       " ('MP:0030855', 'HP:0003996'): 'UPHENO:0079872',\n",
       " ('HP:0005004', 'MP:0030855'): 'UPHENO:0079872',\n",
       " ('HP:0010223', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0010225', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0009154', 'MP:0030882'): 'UPHENO:0076727',\n",
       " ('HP:0010225', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('MP:0030902', 'HP:0003902'): 'UPHENO:0079872',\n",
       " ('HP:0009665', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0010225', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0010223', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0010223', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0010223', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0010222', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0010225', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0010225', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('HP:0010224', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0010223', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('HP:0005013', 'MP:0030855'): 'UPHENO:0079872',\n",
       " ('HP:0010223', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0010224', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('HP:0010224', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0010224', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0010269', 'MP:0030857'): 'UPHENO:0076754',\n",
       " ('HP:0009665', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0010224', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0010269', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009665', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009665', 'MP:0030857'): 'UPHENO:0076754',\n",
       " ('HP:0009665', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0010266', 'MP:0030903'): 'UPHENO:0076754',\n",
       " ('HP:0010224', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0010225', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0010226', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0010276', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0010267', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('MP:0030858', 'HP:0003903'): 'UPHENO:0081802',\n",
       " ('HP:0009154', 'MP:0030856'): 'UPHENO:0076727',\n",
       " ('HP:0009532', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0010018', 'MP:0030903'): 'UPHENO:0076754',\n",
       " ('HP:0010018', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0010018', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0009694', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0010267', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0010267', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0010018', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('MP:0030857', 'HP:0003903'): 'UPHENO:0081802',\n",
       " ('MP:0030882', 'HP:0003903'): 'UPHENO:0081802',\n",
       " ('HP:0010267', 'MP:0030857'): 'UPHENO:0076754',\n",
       " ('HP:0010267', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0010018', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('HP:0010018', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0010265', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0010265', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('MP:0030883', 'HP:0003903'): 'UPHENO:0081802',\n",
       " ('HP:0009521', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0009694', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('MP:0030856', 'HP:0003903'): 'UPHENO:0081802',\n",
       " ('HP:0010276', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0005913', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0010226', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0010226', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0010226', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0010226', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0010226', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('MP:0030882', 'HP:0003902'): 'UPHENO:0079872',\n",
       " ('HP:0010276', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0005913', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0010276', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0009527', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009527', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0010276', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('MP:0030903', 'HP:0003902'): 'UPHENO:0079872',\n",
       " ('HP:0010266', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0010266', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0010266', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0009694', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009527', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0009154', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('HP:0009419', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('MP:0030883', 'HP:0003902'): 'UPHENO:0079872',\n",
       " ('HP:0009512', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0010222', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('HP:0010269', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0005924', 'MP:0030883'): 'UPHENO:0079872',\n",
       " ('HP:0009413', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('HP:0009413', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('MP:0030855', 'HP:0004014'): 'UPHENO:0079872',\n",
       " ('HP:0010269', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0009396', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0005924', 'MP:0030902'): 'UPHENO:0079872',\n",
       " ('HP:0009413', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009418', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0005924', 'MP:0030903'): 'UPHENO:0079872',\n",
       " ('HP:0003898', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0009418', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0009418', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0009418', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0009396', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0009396', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0009187', 'MP:0030856'): 'UPHENO:0076727',\n",
       " ('HP:0009396', 'MP:0030903'): 'UPHENO:0076754',\n",
       " ('HP:0009155', 'MP:0030902'): 'UPHENO:0076727',\n",
       " ('HP:0005913', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009396', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('HP:0010277', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009413', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0005924', 'MP:0030882'): 'UPHENO:0079872',\n",
       " ('HP:0010278', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0009523', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0004000', 'MP:0030855'): 'UPHENO:0079872',\n",
       " ('HP:0010278', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0010278', 'MP:0030857'): 'UPHENO:0076754',\n",
       " ('HP:0010278', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0010277', 'MP:0030903'): 'UPHENO:0076754',\n",
       " ('HP:0010277', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0005924', 'MP:0030855'): 'UPHENO:0079872',\n",
       " ('HP:0005924', 'MP:0030856'): 'UPHENO:0079872',\n",
       " ('HP:0005924', 'MP:0030857'): 'UPHENO:0079872',\n",
       " ('HP:0010277', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0005924', 'MP:0030858'): 'UPHENO:0079872',\n",
       " ('HP:0010277', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0010277', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('HP:0010277', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0009413', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0009396', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009396', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0005913', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('HP:0009695', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('MP:0030858', 'HP:0003902'): 'UPHENO:0079872',\n",
       " ('HP:0010220', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009155', 'MP:0030856'): 'UPHENO:0076727',\n",
       " ('HP:0010220', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0009512', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009187', 'MP:0030902'): 'UPHENO:0076727',\n",
       " ('HP:0009154', 'MP:0030903'): 'UPHENO:0076754',\n",
       " ('HP:0009154', 'MP:0030902'): 'UPHENO:0076727',\n",
       " ('HP:0009521', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0010220', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('HP:0009154', 'MP:0030883'): 'UPHENO:0076727',\n",
       " ('HP:0010222', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0009695', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0010265', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0009695', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0009695', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('MP:0030856', 'HP:0003902'): 'UPHENO:0079872',\n",
       " ('HP:0009695', 'MP:0030903'): 'UPHENO:0076754',\n",
       " ('HP:0010222', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0010222', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('MP:0030857', 'HP:0003902'): 'UPHENO:0079872',\n",
       " ('HP:0010222', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0010220', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0009521', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0009521', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009523', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0009695', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('HP:0009187', 'MP:0030882'): 'UPHENO:0076727',\n",
       " ('HP:0009187', 'MP:0030883'): 'UPHENO:0076727',\n",
       " ('HP:0005913', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0009695', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0009155', 'MP:0030883'): 'UPHENO:0076727',\n",
       " ('HP:0009514', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0005913', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0009418', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0010269', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0010220', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0009514', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0010220', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0009419', 'MP:0030903'): 'UPHENO:0076754',\n",
       " ('HP:0003898', 'MP:0030902'): 'UPHENO:0079872',\n",
       " ('HP:0009155', 'MP:0030882'): 'UPHENO:0076727',\n",
       " ('HP:0003898', 'MP:0030883'): 'UPHENO:0079872',\n",
       " ('HP:0003898', 'MP:0030882'): 'UPHENO:0079872',\n",
       " ('HP:0003898', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('HP:0010266', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('HP:0010247', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0010266', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0009516', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0009526', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0009420', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0009516', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0009420', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0009420', 'MP:0030857'): 'UPHENO:0076754',\n",
       " ('MP:0030902', 'HP:0003903'): 'UPHENO:0081802',\n",
       " ('HP:0009420', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('MP:0030856', 'HP:0003898'): 'UPHENO:0079872',\n",
       " ('MP:0030857', 'HP:0003898'): 'UPHENO:0079872',\n",
       " ('MP:0030858', 'HP:0003898'): 'UPHENO:0079872',\n",
       " ('MP:0030882', 'HP:0003898'): 'UPHENO:0079872',\n",
       " ('HP:0009533', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('MP:0030855', 'HP:0004006'): 'UPHENO:0041360',\n",
       " ('HP:0003898', 'MP:0030857'): 'UPHENO:0079872',\n",
       " ('HP:0009515', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0003903', 'MP:0030883'): 'UPHENO:0079872',\n",
       " ('HP:0009533', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0003903', 'MP:0030882'): 'UPHENO:0079872',\n",
       " ('HP:0003903', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('MP:0030883', 'HP:0003898'): 'UPHENO:0079872',\n",
       " ('MP:0030902', 'HP:0003898'): 'UPHENO:0079872',\n",
       " ('HP:0003903', 'MP:0030902'): 'UPHENO:0079872',\n",
       " ('HP:0009516', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0010266', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0003903', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0010271', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0010271', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0010271', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('HP:0010271', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0010271', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009505', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0009514', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0009419', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0009419', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009505', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0009504', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009504', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0009504', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0009504', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0009503', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009503', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0009503', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0009420', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0009503', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0009526', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0009685', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0009533', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0003898', 'MP:0030856'): 'UPHENO:0079872',\n",
       " ('HP:0003903', 'MP:0030857'): 'UPHENO:0079872',\n",
       " ('HP:0003903', 'MP:0030856'): 'UPHENO:0079872',\n",
       " ('HP:0009194', 'MP:0030902'): 'UPHENO:0084447',\n",
       " ('HP:0200050', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0009525', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0200050', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0009194', 'MP:0030883'): 'UPHENO:0084447',\n",
       " ('HP:0200050', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0009194', 'MP:0030882'): 'UPHENO:0084447',\n",
       " ('HP:0009515', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0200050', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0009683', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009525', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0200050', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('HP:0200050', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009683', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0009525', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0009194', 'MP:0030856'): 'UPHENO:0084447',\n",
       " ('HP:0009515', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0009515', 'MP:0030883'): 'UPHENO:0076754',\n",
       " ('HP:0009523', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009523', 'MP:0030882'): 'UPHENO:0076754',\n",
       " ('HP:0009683', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0009514', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009525', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0009258', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009684', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('HP:0003902', 'MP:0030903'): 'UPHENO:0021677',\n",
       " ('HP:0003902', 'MP:0030902'): 'UPHENO:0079872',\n",
       " ('HP:0003902', 'MP:0030883'): 'UPHENO:0079872',\n",
       " ('HP:0003902', 'MP:0030882'): 'UPHENO:0079872',\n",
       " ('HP:0003899', 'MP:0030856'): 'UPHENO:0079872',\n",
       " ('MP:0030903', 'HP:0003903'): 'UPHENO:0081802',\n",
       " ('HP:0009685', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0009684', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009258', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0009258', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0009258', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('MP:0030903', 'HP:0003898'): 'UPHENO:0079872',\n",
       " ('HP:0003902', 'MP:0030858'): 'UPHENO:0020178',\n",
       " ('HP:0003902', 'MP:0030857'): 'UPHENO:0079872',\n",
       " ('HP:0003902', 'MP:0030856'): 'UPHENO:0079872',\n",
       " ('HP:0009526', 'MP:0030902'): 'UPHENO:0076754',\n",
       " ('HP:0009684', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0009684', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0009684', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0009684', 'MP:0030903'): 'UPHENO:0076754',\n",
       " ('HP:0010271', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0009526', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0010024', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0010256', 'MP:0030856'): 'UPHENO:0076754',\n",
       " ('HP:0010249', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0009673', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0009673', 'MP:0030858'): 'UPHENO:0076754',\n",
       " ('HP:0010254', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0009673', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0010254', 'MP:0030883'): 'UPHENO:0076718',\n",
       " ('HP:0010254', 'MP:0030882'): 'UPHENO:0076718',\n",
       " ('HP:0010254', 'MP:0030857'): 'UPHENO:0076718',\n",
       " ('HP:0010254', 'MP:0030856'): 'UPHENO:0076718',\n",
       " ('HP:0010255', 'MP:0030903'): 'UPHENO:0076754',\n",
       " ('HP:0010255', 'MP:0030902'): 'UPHENO:0076718',\n",
       " ('HP:0009694', 'MP:0030902'): 'UPHENO:0076718'}"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%%time\n",
    "for pair in div_pairs[0:count]:\n",
    "    spath = !runoak -i sqlite:obo:phenio paths -p i,p {pair[0]} --target {pair[1]} --no-autolabel\n",
    "    spath = ((spath[1:])[0]).replace(\"[\", \"\").replace(\"]\", \"\").replace(\"\\t\", \",\").replace(\"\\'\",\"\").replace(\" \",\"\")\n",
    "    spath = spath.split(\",\")\n",
    "\n",
    "    spath_upheno = []\n",
    "    for node in spath:\n",
    "        if node.startswith(\"UPHENO\"):\n",
    "            spath_upheno.append(node)\n",
    "\n",
    "    spath_upheno_mrca = spath_upheno[len(spath_upheno)//2]        \n",
    "\n",
    "    all_spaths[pair] = spath_upheno_mrca\n",
    "\n",
    "all_spaths"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [],
   "source": [
    "with open(\"upheno_paths.tsv\", 'w') as outfile:\n",
    "    for spath in all_spaths:\n",
    "        outfile.write(f'{spath[0]}\\t{spath[1]}\\t{all_spaths[spath]}\\n')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "colab": {
   "provenance": []
  },
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}