GTOOL and Plink: Single Dosage Probabilities
1
1
Entering edit mode
5.0 years ago

I have SNP data where each SNP is expressed as a single dosage (0...2). According to the methodology, which I am trying to replicate, gtool was used to make hard-calls and convert it to plink's .PED file format.

The schema of the SNP data is as followed....

SNP Subject 1   Subject 2   Subject 3   …
rs0001  0.001      2          1.490 
rs0002  1.002      0.3900        0  
…

As far as I can tell, gtool utilizes 3-probabilities per SNPs. Additionally, I can't think of anyway to convert single probabilities to 3-probabilities, as specified in the .GEN format.

Thus, I am wondering where I am going wrong. Perhaps there were different versions of gtool (this analysis was done 5+ years ago) or other tools/steps that I should be considering.

plink gtool • 1.3k views
ADD COMMENT
2
Entering edit mode
5.0 years ago
zx8754 11k

As far as I know you can't go back from dosage to 3 probabilities. In other words, can you solve below?

x*0 + y*1 + z*2 = 1.49
x + y + z = 1

Plink can take --dosage as input.

ADD COMMENT

Login before adding your answer.

Traffic: 2336 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6