How to integrate lods from different markers and in different families in the same gene?

[changebio]

SEQLinkage/source/Runner.py

Line 376 in aa3f900

elif fam not in lods[theta] or lod > lods[theta][fam]:

Based on the code, it uses max lod in markers and add all lods among families.
lod = log10(P(theta))-log10(P(0.5)), theta is seq(0,0.45,0.05)
About integrating among families.
I think it is fine if the max lods are bigger than 0, we add them up.
If the max lods are less than 0, what is the meaning? should we set them to 0?

[gaow]

LOD is Log odds and can certainly be smaller than 0: it means the odds is smaller than 1. We should keep them as is.