Question: Problems to extract a SNP list from .ped and .map files
0
gravatar for saual.maria
4.5 years ago by
saual.maria10
Spain
saual.maria10 wrote:

Dear all. I want to remove a list of SNPs from my .ped and .map files. I have seen that plink has these instructions --exclude and --extract, that remove SNPs in the first case or create a new file with the SNPs you want to keep in the second case. Both would help me. 

I execute the instruction:

plink --file myfile --extract mySNPlist.txt --recode --out newfile

(or the same in the opposite way)

plink --file myfile --exclude mySNPlist2.txt --recode --out newfile

and it works in both cases BUT when I open the new map file, although the number of SNPs is correct, it set many SNPs to 0 in the column corresponding to chromosome and I don´t understand why, as the SNP list I provide only includes mapped SNPs.

Could you please help me with this? Thank you in advance. María

plink snp subset forum • 2.2k views
ADD COMMENTlink modified 4.4 years ago by Biostar ♦♦ 20 • written 4.5 years ago by saual.maria10
1
gravatar for saual.maria
4.5 years ago by
saual.maria10
Spain
saual.maria10 wrote:

I have already fixed the problem... I was working with plink version 1.07 but the version 1.9 allows to set the number of chromosomes (that was the problem, by default plink assumes data are from humans). Therefore by adding --chr-set numberofchromosomes the problem is solved

ADD COMMENTlink written 4.5 years ago by saual.maria10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 512 users visited in the last hour