Skip to content

Commit 8a7640e

Browse files
committed
Update tutorial to mention tcrblosum usage
1 parent 6c4c5e6 commit 8a7640e

1 file changed

Lines changed: 5 additions & 3 deletions

File tree

docs/tutorials/tutorial_3k_tcr.ipynb

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1423,10 +1423,11 @@
14231423
"\n",
14241424
"To this end, we need to set `metric=\"tcrdist\"` and specify a `cutoff` parameter.\n",
14251425
"The distance is based on the [BLOSUM62](https://en.wikipedia.org/wiki/BLOSUM) matrix.\n",
1426-
"For instance, a distance of `10` is equivalent to 2 Rs mutating into N.\n",
1427-
"This appoach was initially proposed by Dash et al. {cite}`TCRdist` and is based on the [tcrdist3](https://github.com/kmayerb/tcrdist3) implementation.\n",
1426+
"This approach was initially proposed by Dash et al. {cite}`TCRdist` and is based on the [tcrdist3](https://github.com/kmayerb/tcrdist3) implementation.\n",
1427+
"For instance, a distance of `12` is equivalent to two `R`s mutating into `K`s with the default parameters.\n",
1428+
"Alternatively, [TCRBLOSUM](https://doi.org/10.1093/bib/bbae602) alpha/beta substitution matrices can be selected with `base_matrix=\"tcrblosum\"`.\n",
14281429
"\n",
1429-
"All cells with a distance between their CDR3 sequences lower than `cutoff` will be connected in the network.\n"
1430+
"Cells are connected in the network when their CDR3 sequence distances are at or below `cutoff`. In this case, this must hold for both receptor arms (`receptor_arms=\"all\"`, i.e. both VJ and VDJ chains), considering any possible dual immune receptor chain pairing (`dual_ir=\"any\"`).\n"
14301431
]
14311432
},
14321433
{
@@ -1464,6 +1465,7 @@
14641465
" metric=\"tcrdist\",\n",
14651466
" sequence=\"aa\",\n",
14661467
" cutoff=15,\n",
1468+
" # base_matrix=\"tcrblosum\",\n",
14671469
")"
14681470
]
14691471
},

0 commit comments

Comments
 (0)