search

sgkit-dev / sgkit

Scalable genetics toolkit
https://sgkit-dev.github.io/sgkit
Apache License 2.0
231 stars 32 forks source link

Genealogical relatedness greater than 1 #1188

Closed fouerghi closed 7 months ago

fouerghi commented 8 months ago

Hi everyone, I was wondering if there were any explanations of how coefficient of relationship could be greater than 1. I am getting values of 1.5 genealogical relatedness between an individual and itself and I am not sure where this could be coming from.

timothymillar commented 8 months ago

Hi @fouerghi, thanks for the question. I can't be certain of what the cause is without seeing your data, but a self-relatedness > 1 is expected if that individual's parents are related to one another. I assume that you are calculating this from a genealogy/pedigree? If so, then you should be able to check that the parents are related. If they are not meant to be related, then this result may be a symptom of a misspecified pedigree.

fouerghi commented 8 months ago

Thank you @timothymillar!

jeromekelleher commented 7 months ago

Closing this as addressed. Feel free to reopen if anyone feels there's something we should follow up on.