Skip to content

Commit d34faa3

Browse files
committed
Update evo opp threshold
1 parent 0c33196 commit d34faa3

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

notebooks/compute_rates.ipynb

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -470,12 +470,12 @@
470470
"cell_type": "markdown",
471471
"metadata": {},
472472
"source": [
473-
"Compute genome-wide rates, grouping data across all subtypes. Only analyze mutations with an evolutionary opportunity >1, corresponding to 1 mutation count per site in a segment across all branches where a mutation was possible. This focuses on the mutations with the most data (though results do not really change if all mutations are included)."
473+
"Compute genome-wide rates, grouping data across all subtypes. Only analyze mutations with an evolutionary opportunity >0.1, corresponding to 0.1 synonymous mutation count per site in a segment across all branches where a mutation was possible. This focuses on the mutations with the most data (though results do not really change if all mutations are included)."
474474
]
475475
},
476476
{
477477
"cell_type": "code",
478-
"execution_count": 5,
478+
"execution_count": null,
479479
"metadata": {
480480
"execution": {
481481
"iopub.execute_input": "2026-03-09T18:48:57.969616Z",
@@ -605,7 +605,7 @@
605605
],
606606
"source": [
607607
"# Compute genome-wide rates for each host category (all, human, avian)\n",
608-
"genome_wide_evo_opp_threshold = 1.0\n",
608+
"genome_wide_evo_opp_threshold = 0.1\n",
609609
"genome_wide_rates = (\n",
610610
" counts_df[counts_df['evo_opp'] > genome_wide_evo_opp_threshold]\n",
611611
" .groupby(['mut_type', 'mut_class', 'host'])[['actual_count', 'evo_opp']].sum()\n",
@@ -1613,7 +1613,7 @@
16131613
"name": "python",
16141614
"nbconvert_exporter": "python",
16151615
"pygments_lexer": "ipython3",
1616-
"version": "3.10.18"
1616+
"version": "3.13.12"
16171617
}
16181618
},
16191619
"nbformat": 4,

0 commit comments

Comments
 (0)