About the construction of lines in the MS

[nablator]

It’s a simple score.

Score (start, end) = base 2 logarithm of ratio of (start, end) probability to interior probability

Python code:

Code:

for _, r in fam_pos.iterrows():     fam = r["family"]     p_end = r["end"] / tot_end if tot_end else 0.0     p_start = r["start"] / tot_start if tot_start else 0.0     p_int = r["interior"] / tot_int if tot_int else 0.0     # Positive value = enriched at line edge relative to interior.     A_end_score[fam] = np.log2(p_end / p_int) if p_end > 0 and p_int > 0 else np.nan     A_start_score[fam] = np.log2(p_start / p_int) if p_start > 0 and p_int > 0 else np.nan