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

SNPs 0,1,2 #24

Open
DavidGO464 opened this issue Feb 28, 2023 · 2 comments
Open

SNPs 0,1,2 #24

DavidGO464 opened this issue Feb 28, 2023 · 2 comments

Comments

@DavidGO464
Copy link

Hi!

I am using this snp from illumina exm1615904... whe I use plink it only uses the name of the SNP to get the geenotype (eg, as in this case normal would be CC vs GC or GG).
When usign BEDmatrix, I need to add a '_X' (underscore somthing) to the SNP name and then when I call the SNP it shows me 0,1,2, ... So I wonder what doees it means.

Kind regards

image

@agrueneberg
Copy link
Member

Hi David,

you can omit the minor allele from the column names by creating the BEDMatrix object using X <- BEDMatrix(file, simple_names = TRUE). The output values are allelic dosages (0/1/2/NA), counting the number of minor alleles.

Hope that helps,
Alex

@DavidGO464
Copy link
Author

DavidGO464 commented Feb 28, 2023

Hi David,

you can omit the minor allele from the column names by creating the BEDMatrix object using X <- BEDMatrix(file, simple_names = TRUE). The output values are allelic dosages (0/1/2/NA), counting the number of minor alleles.

Hope that helps, Alex

Thank you for your prompt response. I think that the number is just the coding for a category regarding if it is minor allele homozygus (0: in case the variation in C>T... it would be TT; thats why is minor) and then heterozygus (eg, 1: CT) and major allele homozygus (2 : eg, CC)
i found this:

https://www.cs.cmu.edu/~sssykim/teaching/s13/slides/Lecture1.pdf

So, it would actually count the number of Major alleles (I mean, seems to be that way)

image

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

2 participants