Remove Some Snps In Plink
5
1
Entering edit mode
8.9 years ago
mary ▴ 210

Hello all

I want remove some snp from map and also genotype from ped file. I used this command in plink:

plink --file FILE  --exclude SNP-LIST-FILE

I think the some SNP in map file should be deleted and also genotype in ped, but when I check the number of raw in map and the column of ped file its same as before exclude

is there any one know my snp in snplist excluded or not?

thanks

plink • 14k views
ADD COMMENT
1
Entering edit mode

thank you very much. it worked

ADD REPLY
7
Entering edit mode
8.9 years ago

Your command is correct, but you should add two options: --out and --recode. If you don't specify these, plink executes the commands and prints a summary to the screen (and to the .log file), but does not actually write the results to a ped file nor overwrites your original data.

The --recode option specifies that you want to write a .ped file as output. The --out option specifies the suffix of the output file:

plink --file mydata --exclude snp_list --recode --out filter_data

Also, check that the the SNP-LIST-FILE is formatted correctly. It must contain the ids of the SNPs that you want to remove.

ADD COMMENT
0
Entering edit mode

I also used the above command but it only removes one occurrence and the remaining stays the same. For example, if rs12345 occurs 5 times, it only removes one of them.

Is there a way of removing all occurrences in one go?

ADD REPLY
1
Entering edit mode
6.9 years ago
emblascz ▴ 10

Does this completely remove records or just recode them so PLINK doesn't use the data?  I want the new PED file to actually have fewer records. I am the data manager for a project using using PLINK and we've discovered consent issues with a handful of cases. I want to completely eliminate their data, not just recode data so PLINK doesn't use it.  Is this the best way to do that? 

Thanks,

Emily

ADD COMMENT
1
Entering edit mode

With --recode (or better, --make-bed), the newly written file will not contain the excluded records at all.

ADD REPLY
0
Entering edit mode
8.8 years ago
roll ▴ 330

I also used the above command but it only removes one occurrence and the remaining stays the same. For example, if rs12345 occurs 5 times, it only removes one of them.

Is there a way of removing all occurrences in one go?

ADD COMMENT
0
Entering edit mode
5.9 years ago
chenhr • 0

hello

I want to used .ped and .map file to make population structure by frappe with a maximum likelihood method.first ,I used this command in plink: plink --ped DataResult.ped --recode12 ,then used this command in frappe: frappe MyParamSim.txt, but it aways failed,the error as follow:

marker 664146 has only 2 non-missing alleles (8 required)
marker 664153 has only 2 non-missing alleles (8 required)
marker 664155 has only 4 non-missing alleles (8 required)
marker 664171 has only 2 non-missing alleles (8 required)
marker 664201 has only 2 non-missing alleles (8 required)
marker 664203 has only 6 non-missing alleles (8 required)
marker 664231 has only 2 non-missing alleles (8 required)
marker 664255 has only 2 non-missing alleles (8 required)
marker 664265 has only 6 non-missing alleles (8 required)
marker 664270 has only 6 non-missing alleles (8 required)
marker 664278 has only 4 non-missing alleles (8 required)
marker 664295 has only 6 non-missing alleles (8 required)
marker 664297 has only 6 non-missing alleles (8 required)
marker 664309 has only 6 non-missing alleles (8 required)
marker 664362 has only 6 non-missing alleles (8 required)----

according to the error, I used this command in plink: plink --file DataResult.ped --geno 0.1, but it fail.

please help me, please, thank you!

ADD COMMENT

Login before adding your answer.

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