timothyfrasier / related

8 stars 7 forks source link

How to differentiate fullsib from parent/offspring diads? #5

Open ronahuel opened 4 years ago

ronahuel commented 4 years ago

Hello everyone. Currently i'm trying to estimate the relatedness among 120+- individuals. I've performed simulations based on my genotype data and there's a high correlation between simulated observed and expected values for almost all of the estimators. So I started to work with trioML. My question is, how can I, or how should I stablish the difference between fullsib-paren/offsprping when my trioML estimate is near 0.5? I made this question because I cannot find a way to perform a boxplot with the assignment, as the simulation values can be plotted according to the tutorial.

Another question I would like to ask, is how to interpret when a particular diad estimates of relatedness is >1 (1.9). These individuals are highly homocigous and had a high F coefficient. How can I estimate the relationship when there's such a high value? Can they be fullsib - parent/offspring?

If you can help me with this I would really appreciate it.

Thanks.

timothyfrasier commented 4 years ago

Hi:

A good way to differentiate between potential relationships with similar r-values (e.g., full-sibs and parent-offspring) is to look at the other moments of relatedness, specifically delta7 and delta8, which are provided in the outfiles. For parent-offspring, the expected delta7 is 0 whereas for full-sibs it is 0.25. For parent-offspring the expected delta8 is 1 whereas for full-sibs it is 0.5. See Wang (2011) Molecular Ecology Resources 11: 141-145 for details.

Again, comparing these values, in addition to the actual relatedness value, should allow better resolution when deciphering relationships, include those with r-values > 1.

-Tim


Timothy R. Frasier Coordinator - Forensic Sciences Program Associate Professor - Department of Biology Saint Mary's University 923 Robie Street Halifax, NS B3H 3C3, Canada Email: timothy.frasier@smu.ca Tel: (902) 491-6382 www.frasierlab.ca http://frasierlab.ca/

On Jul 24, 2020, at 2:46 PM, ronahuel notifications@github.com wrote:

CAUTION: This email is from an external sender. Do not click on links or open attachments unless you trust the sender. Please forward suspected phishing emails to reportphishing@smu.ca

Hello everyone. Currently i'm trying to estimate the relatedness among 120+- individuals. I've performed simulations based on my genotype data and there's a high correlation between simulated observed and expected values for almost all of the estimators. So I started to work with trioML. My question is, how can I, or how should I stablish the difference between fullsib-paren/offsprping when my trioML estimate is near 0.5? I made this question because I cannot find a way to perform a boxplot with the assignment, as the simulation values can be plotted according to the tutorial.

Another question I would like to ask, is how to interpret when a particular diad estimates of relatedness is >1 (1.9). These individuals are highly homocigous and had a high F coefficient. How can I estimate the relationship when there's such a high value? Can they be fullsib - parent/offspring?

If you can help me with this I would really appreciate it.

Thanks.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/timothyfrasier/related/issues/5, or unsubscribe https://github.com/notifications/unsubscribe-auth/AA3NYWKMEOYLINE2CVB4CLLR5HCGNANCNFSM4PG5XCQQ.

ronahuel commented 4 years ago

Thanks for you early answer !! Thank you so much!

ronahuel commented 4 years ago

Hi: A good way to differentiate between potential relationships with similar r-values (e.g., full-sibs and parent-offspring) is to look at the other moments of relatedness, specifically delta7 and delta8, which are provided in the outfiles. For parent-offspring, the expected delta7 is 0 whereas for full-sibs it is 0.25. For parent-offspring the expected delta8 is 1 whereas for full-sibs it is 0.5. See Wang (2011) Molecular Ecology Resources 11: 141-145 for details. Again, comparing these values, in addition to the actual relatedness value, should allow better resolution when deciphering relationships, include those with r-values > 1. -Tim __ Timothy R. Frasier Coordinator - Forensic Sciences Program Associate Professor - Department of Biology Saint Mary's University 923 Robie Street Halifax, NS B3H 3C3, Canada Email: timothy.frasier@smu.ca Tel: (902) 491-6382 www.frasierlab.ca http://frasierlab.ca/ On Jul 24, 2020, at 2:46 PM, ronahuel @.> wrote: CAUTION: This email is from an external sender. Do not click on links or open attachments unless you trust the sender. Please forward suspected phishing emails to @. Hello everyone. Currently i'm trying to estimate the relatedness among 120+- individuals. I've performed simulations based on my genotype data and there's a high correlation between simulated observed and expected values for almost all of the estimators. So I started to work with trioML. My question is, how can I, or how should I stablish the difference between fullsib-paren/offsprping when my trioML estimate is near 0.5? I made this question because I cannot find a way to perform a boxplot with the assignment, as the simulation values can be plotted according to the tutorial. Another question I would like to ask, is how to interpret when a particular diad estimates of relatedness is >1 (1.9). These individuals are highly homocigous and had a high F coefficient. How can I estimate the relationship when there's such a high value? Can they be fullsib - parent/offspring? If you can help me with this I would really appreciate it. Thanks. — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#5>, or unsubscribe https://github.com/notifications/unsubscribe-auth/AA3NYWKMEOYLINE2CVB4CLLR5HCGNANCNFSM4PG5XCQQ.

Thanks again Timo. I've read Wang 2011. I've noted that Table 1 in that paper is not the same than the Coancestry Guideline's Table 1.

There are differences in the interpretation of Rxy for almost all kind of dyads. Have you any idea why?