diff --git a/HowBadIsMyBatch.ipynb b/HowBadIsMyBatch.ipynb index d99a9d62502..efd3a7a71fb 100644 --- a/HowBadIsMyBatch.ipynb +++ b/HowBadIsMyBatch.ipynb @@ -1226,6 +1226,73 @@ "\n", "sns.boxplot(x = \"VAX_LOT_PREFIX\", y = \"ADRs\", data = batchCodeTable)" ] + }, + { + "cell_type": "markdown", + "id": "731c27a5", + "metadata": {}, + "source": [ + "#### Moderna Batches" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "b4a9c489", + "metadata": {}, + "outputs": [], + "source": [ + "import re\n", + "\n", + "dataFrame = DataFrameFilter().filterByCovid19(vaers)\n", + "dataFrame = DataFrameFilter().filterBy(dataFrame, manufacturer = \"MODERNA\")\n", + "batchCodeTable = BatchCodeTableFactory._createSummationTableByVAX_LOT(dataFrame)[['ADRs']].reset_index()\n", + "modernaBatchCodePrefix = re.compile(r'^[0-9]{3}[a-zA-Z]')\n", + "batchCodeTable = batchCodeTable[batchCodeTable['VAX_LOT'].apply(lambda vax_lot: bool(modernaBatchCodePrefix.match(vax_lot)))]\n", + "batchCodeTable['CONCENTRATION'] = batchCodeTable['VAX_LOT'].str[3]\n", + "batchCodeTable = batchCodeTable.sort_values(by = 'CONCENTRATION', ascending = True)\n", + "batchCodeTable = batchCodeTable[batchCodeTable['ADRs'] > 400]\n", + "batchCodeTable" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "e26c9d85", + "metadata": {}, + "outputs": [], + "source": [ + "import seaborn as sns\n", + "\n", + "order = ['J', 'K', 'L', 'M', 'A', 'B', 'C', 'D', 'E', 'F', 'G', 'H']\n", + "\n", + "sns.set(rc = {'figure.figsize': (11.7, 8.27)})\n", + "sns.set_theme()\n", + "chart = sns.stripplot(x = \"CONCENTRATION\", y = \"ADRs\", data = batchCodeTable, order = order)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "d1de13c7", + "metadata": {}, + "outputs": [], + "source": [ + "sns.pointplot(x = \"CONCENTRATION\", y = \"ADRs\", data = batchCodeTable, estimator = np.mean, order = order)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "29ae8ca2", + "metadata": {}, + "outputs": [], + "source": [ + "import seaborn as sns\n", + "sns.set_theme(style = \"ticks\", palette = \"pastel\")\n", + "\n", + "sns.boxplot(x = \"CONCENTRATION\", y = \"ADRs\", data = batchCodeTable, order = order)" + ] } ], "metadata": {