mhguo1 / TRAPD

Burden testing against public controls
MIT License
50 stars 32 forks source link

MAF and AF #7

Open sureyana956 opened 5 years ago

sureyana956 commented 5 years ago

Hi, thanks for providing this tool. I have a question regarding how MAF cutoff was applied to subset SNP lists. Is alt always considered as the minor allele in your scripts? Take rs7954255 for example, Latino and African actually have a AF greater than 0.5, which means the reference allele becomes the minor allele. Did you take that into consideration in your count_controls.py and make_snp_file.py scripts?

mhguo1 commented 5 years ago

Hi, The code currently only considers the alternate allele as the minor allele. It is a good point and I will try to revise TRAPD to take the minor allele rather than alternate allele. It may be tricky in certain situations (especially if the minor allele is different across populations).

Thanks, Mike

On Jun 6, 2019, at 7:18 PM, sureyana956 notifications@github.com wrote:

Hi, thanks for providing this tool. I have a question regarding how MAF cutoff was applied to subset SNP lists. Is alt always considered as the minor allele in your scripts? Take rs7954255 for example, Latino and African actually have a AF greater than 0.5, which means the reference allele becomes the minor allele. Did you take that into consideration in your count_controls.py and make_snp_file.py scripts?

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/mhguo1/TRAPD/issues/7?email_source=notifications&email_token=AD3L3NISXJBXIVLWRVD2633PZGLNXA5CNFSM4HVL75HKYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4GYEO6YQ, or mute the thread https://github.com/notifications/unsubscribe-auth/AD3L3NNNL24KTGBXBWKWMH3PZGLNXANCNFSM4HVL75HA.