GWAS Analysis using Plink
1
0
Entering edit mode
4 weeks ago
123 • 0

Hi everyone,

I have only one allele information for each SNP in my ped file. How can I do an association study in plink?

enter image description here

Here 0 = SNP1, 1 = SNP2 and so on.

How do I associate this phenotype to my SNPs?

GWAS PLINK • 140 views
ADD COMMENT
0
Entering edit mode
4 weeks ago

Are you sure PLINK is reading the file correctly and hasn't split the snp across as two columns, so col's 0 and 1 should be merged (Could have been a header issue?). Did you convert your data to PED at some point from vcf, as something could have gone wrong there.

Unless you expect your data to be haploid, which it might be? I don't think PLINK is well suited for such data, but if you do, I would maybe reformat your ped file by simple duplication so A becomes AA (this would require bash awk/sed commands). And then you then might want to check the REF/ALT alleles by generating a frequency file of your dataset

plink --file test --freq 

and then this should recode your ped file

plink --file test --set-hh-missing

There's a bit more info here, https://zzz.bwh.harvard.edu/plink/data.shtml, but this is just a very loose suggestion and I haven't tried it myself (nor do I think most people would have!)

ADD COMMENT

Login before adding your answer.

Traffic: 1994 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