Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Inference questions #129

Open
solcos opened this issue Dec 1, 2023 · 0 comments
Open

Inference questions #129

solcos opened this issue Dec 1, 2023 · 0 comments

Comments

@solcos
Copy link

solcos commented Dec 1, 2023

Hello,

I've been testing somalier for relatedness and also sex inference.

I have some questions regarding sex and relatedness.

  1. Sex inference:
    I've noticed that in some samples, the results for sex inference follow the formulas in (Why this is a male? #125)
    but the result is not 'translated' in the sex column (inferred sex) of the somalier.samples file.

e.g.
sex phenotype X_depth_mean X_het X_hom_alt
-1 unknown 17.45 1 162

and the formula is --> if stats[i].x_het / stats[i].x_hom_alt < 0.05 and stats[i].x_dp.n > 10:

The results should be male, am I correct?
The 'x_dp.n' variable is the 'X_depth_mean' in the samples.tsv?
Or do you take into account other parameters?


  1. Relatedness

Regarding relatedness, I've noticed that in "https://github.com/brentp/somalier/wiki/pedigree-inference" to infer the relationships of the samples, you need a minimum of 3 samples "a trio", have relatedness value between 0.4 and 0.6 to $sample, have IBS0 / IBS2 < 0.005 to $sample, and also that the parents can't have a relatedness higher than 0.06.

Why the 0.06 is needed?

Thanks in advance!

Arnau

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant