This is a test version of Biostars. For the public version, visit https://www.biostars.org.
extract SNPs (rs.....) along with subject ID and pedigree ID in PLINK

Hi all,

I am trying to extract SNPs for a specific location with Pedigree ID information in there, however every time I use parameters, I get SNPs but without the pedigree or family ID. I want to merge SNPs with phenotype information using the common ID.

plink --bfile FULL_MERGE_chr17 --chr 17 --from-bp 39921041 --to-bp 39927601 --write-snplist --out FULL_MERGE_chr17.res

It produces a file which only has one column with all SNPs (FULL_MERGE_chr17.res.snpslist), however misses ID.

Could anyone suggest which parameter could be used to add ID information in the same SNP file?

I have tried adding --recode12 and --transpose, --recode-rlist but nothing seems to work! At some point, I got a different number of SNPs, but not all of them.

I would appreciate any help!

Thanks, Best, Priya

plink
• 2,390 views • link
written by pkpdkachroo • pkpdkachroo

Maybe just use --recode this gives you a .ped file and here the ID is next to the SNPs that are in your region

• 0 views • link
written by Floris Brenk Floris Brenk

1 answer

Thanks, but I used that and as before my SNP file has 107 SNPs, but wc -l on the ped file gives me 2183. It does not look like I got the right ID.

• 0 views • link
written by pkpdkachroo • pkpdkachroo

Log in to answer this question.

More posts like this