diff --git a/jupyter/random_angles_vs_generated.ipynb b/jupyter/random_angles_vs_generated.ipynb index ddc9509..58b927d 100644 --- a/jupyter/random_angles_vs_generated.ipynb +++ b/jupyter/random_angles_vs_generated.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": 19, + "execution_count": 25, "metadata": {}, "outputs": [ { @@ -11,7 +11,7 @@ "PosixPath('/home/wukevin/projects/protdiff/plots')" ] }, - "execution_count": 19, + "execution_count": 25, "metadata": {}, "output_type": "execute_result" } @@ -22,6 +22,7 @@ "import json\n", "\n", "import numpy as np\n", + "from scipy import stats\n", "from matplotlib import pyplot as plt\n", "import seaborn as sns\n", "\n", @@ -45,6 +46,34 @@ " generated_sctm = json.load(source)" ] }, + { + "cell_type": "code", + "execution_count": 31, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "MannwhitneyuResult(statistic=417832.5, pvalue=4.374068292433547e-60)" + ] + }, + "execution_count": 31, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "def to_arr(x):\n", + " x = np.array(list(x))\n", + " return x[~np.isnan(x)]\n", + "\n", + "\n", + "stats.mannwhitneyu(\n", + " to_arr(generated_sctm.values()),\n", + " to_arr(baseline_sctm.values()),\n", + ")" + ] + }, { "cell_type": "code", "execution_count": 24, diff --git a/plots/sctm_scores_gen_vs_baseline.pdf b/plots/sctm_scores_gen_vs_baseline.pdf index adf9f35..555741b 100644 Binary files a/plots/sctm_scores_gen_vs_baseline.pdf and b/plots/sctm_scores_gen_vs_baseline.pdf differ