{ "cells": [ { "cell_type": "code", "execution_count": null, "id": "ad6c7885", "metadata": { "ExecuteTime": { "end_time": "2026-05-11T23:09:56.928583Z", "start_time": "2026-05-11T23:09:54.461627Z" } }, "outputs": [], "source": [ "import glob\n", "import pandas as pd" ] }, { "cell_type": "code", "execution_count": null, "id": "5cb011e5", "metadata": { "ExecuteTime": { "end_time": "2026-05-11T23:10:37.275332Z", "start_time": "2026-05-11T23:10:37.268530Z" } }, "outputs": [], "source": [ "def path_to_image_html(path):\n", " \"\"\"\n", " Converts a file path or URL to an HTML image tag for display in a DataFrame.\n", " \"\"\"\n", " return f''" ] }, { "cell_type": "code", "execution_count": null, "id": "04207506", "metadata": { "ExecuteTime": { "end_time": "2026-05-11T23:21:10.590816Z", "start_time": "2026-05-11T23:21:10.587169Z" } }, "outputs": [], "source": [ "def plot_motifs_table(models, head='counts'):\n", " motif_df = pd.DataFrame(index=range(5), columns=models)\n", " motif_df = motif_df.fillna('')\n", " for model in models:\n", " modisco_path = f'modisco/{model}/{head}/trimmed_logos/pos*.cwm.fwd.png'\n", " num_motifs = len(glob.glob(modisco_path))\n", " num_motifs_to_show = min(5, num_motifs) # show at most top 5 motifs\n", " motif_image_paths = []\n", " for i in range(num_motifs_to_show):\n", " motif_image_path = f'modisco/{model}/{head}/trimmed_logos/pos_patterns.pattern_{i}.cwm.fwd.png'\n", " motif_df[model].iloc[i] = motif_image_path\n", "\n", " motif_df = motif_df.style.format(path_to_image_html)\n", " return motif_df" ] }, { "cell_type": "code", "execution_count": 22, "id": "68e1a965", "metadata": { "ExecuteTime": { "end_time": "2026-05-11T23:21:11.108255Z", "start_time": "2026-05-11T23:21:11.096697Z" } }, "outputs": [ { "data": { "text/html": [ "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
 ALL1_KMT2A_CALL1_KMT2A_N
0
1
2
3
4
\n" ], "text/plain": [ "" ] }, "execution_count": 22, "metadata": {}, "output_type": "execute_result" } ], "source": [ "models = ['ALL1_KMT2A_C', 'ALL1_KMT2A_N']\n", "plot_motifs_table(models, 'counts')" ] }, { "cell_type": "code", "execution_count": 23, "id": "53e6b9d7", "metadata": { "ExecuteTime": { "end_time": "2026-05-11T23:21:13.100845Z", "start_time": "2026-05-11T23:21:13.091248Z" } }, "outputs": [ { "data": { "text/html": [ "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
 ALL1_KMT2A_CALL1_KMT2A_N
0
1
2
3
4
\n" ], "text/plain": [ "" ] }, "execution_count": 23, "metadata": {}, "output_type": "execute_result" } ], "source": [ "plot_motifs_table(models, 'profile')" ] }, { "cell_type": "code", "execution_count": 24, "id": "c86f02a2", "metadata": { "ExecuteTime": { "end_time": "2026-05-11T23:21:14.187194Z", "start_time": "2026-05-11T23:21:14.175569Z" } }, "outputs": [ { "data": { "text/html": [ "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
 RS411_KMT2A_CRS411_KMT2A_N
0
1
2
3
4
\n" ], "text/plain": [ "" ] }, "execution_count": 24, "metadata": {}, "output_type": "execute_result" } ], "source": [ "models = ['RS411_KMT2A_C', 'RS411_KMT2A_N']\n", "plot_motifs_table(models, 'counts')" ] }, { "cell_type": "code", "execution_count": 25, "id": "85345341", "metadata": { "ExecuteTime": { "end_time": "2026-05-11T23:21:17.898083Z", "start_time": "2026-05-11T23:21:17.888465Z" } }, "outputs": [ { "data": { "text/html": [ "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
 RS411_KMT2A_CRS411_KMT2A_N
0
1
2
3
4
\n" ], "text/plain": [ "" ] }, "execution_count": 25, "metadata": {}, "output_type": "execute_result" } ], "source": [ "plot_motifs_table(models, 'profile')" ] }, { "cell_type": "code", "execution_count": 28, "id": "3142ca1a", "metadata": { "ExecuteTime": { "end_time": "2026-05-11T23:21:56.283438Z", "start_time": "2026-05-11T23:21:56.260738Z" } }, "outputs": [ { "data": { "text/html": [ "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
 CD34_KMT2A_CCD34_KMT2A_N
0
1
2
3
4
\n" ], "text/plain": [ "" ] }, "execution_count": 28, "metadata": {}, "output_type": "execute_result" } ], "source": [ "models = ['CD34_KMT2A_C', 'CD34_KMT2A_N']\n", "plot_motifs_table(models, 'counts')" ] }, { "cell_type": "code", "execution_count": 29, "id": "38ce1930", "metadata": { "ExecuteTime": { "end_time": "2026-05-11T23:21:56.416303Z", "start_time": "2026-05-11T23:21:56.395162Z" } }, "outputs": [ { "data": { "text/html": [ "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
 CD34_KMT2A_CCD34_KMT2A_N
0
1
2
3
4
\n" ], "text/plain": [ "" ] }, "execution_count": 29, "metadata": {}, "output_type": "execute_result" } ], "source": [ "plot_motifs_table(models, 'profile')" ] }, { "cell_type": "code", "execution_count": null, "id": "21073e26", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "bpnet", "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.7.12" } }, "nbformat": 4, "nbformat_minor": 5 }