The existing code would deduplicate LR pairs for humans but not mice.
For humans:
lrs_full_arr = np.unique(np.array(lrs_full))
For mice it reused lrs_full, not lrs_full_arr:
genes1 = [lr_.split("_")[0] for lr_ in lrs_full] genes2 = [lr_.split("_")[1] for lr_ in lrs_full]
lrs_full_arr = np.array(
[
genes1[i][0]
+ genes1[i][1:].lower()
+ "_"
+ genes2[i][0]
+ genes2[i][1:].lower()
for i in range(len(lrs_full))
],
)
Change - lrs is now a set and mouse output for the multi-database calls goes from 6364 entries (2293 duplicated) to 4071 unique, sorted. Existing mouse results won't reproduce exactly.
The existing code would deduplicate LR pairs for humans but not mice.
For humans:
For mice it reused lrs_full, not lrs_full_arr:
Change - lrs is now a set and mouse output for the multi-database calls goes from 6364 entries (2293 duplicated) to 4071 unique, sorted. Existing mouse results won't reproduce exactly.