HetzDra / turboGliph

R implementation of GLIPH (Grouping of Lymphocyte Interactions by Paratope Hotspots), an algorithm developed by Glanville et al to identify specificity groups in the T cell receptor repertoire based on local (motif sharing) and global (hamming distance) similarities.
17 stars 4 forks source link

combine problem #7

Open zunyun-Gong opened 6 months ago

zunyun-Gong commented 6 months ago

Thank you for sharing your work. It is a well established R package. When I use it, I notice some problem in gliph2 function. when in some kind of input format, there will be a simpleError .

a519f401509032aa97939278d972655

this might because in the gliph function there is potential hazard showed below

6568a6b8bdb0284e45b429d1bc45eef

if there is only one column in cdr3_sequences which sequences do not have, the sequences will add the new column but with the header cdr3_sequences[, !(base::colnames(cdr3_sequences) %in% base::colnames(sequences))], this may cause the simpleError mentioned above.