Split plink files by a number of SNPs
0
0
Entering edit mode
5.9 years ago
kakukeshi ▴ 80

Hi guys,

Is there any way to split the bed, bim and fam files by a specific number of SNPs. Let say in the files for chr1 I have 86000 SNPs. I would like to generate 86 files with 1000 SNPs each.

Many thanks

GWAS SNP plink • 3.3k views
ADD COMMENT
1
Entering edit mode

You can extract SNPs of your interest using --extract flag. Loop for every 1,000 SNPs.

ADD REPLY
0
Entering edit mode

Thanks, that's exactly what I'm doing. I created x files with different subsets of snps and the extract them using --extract flag

ADD REPLY
0
Entering edit mode

Please use tags appropriately. As such experts can quickly find your question. In this case "plink" would have been relevant, so I have edited your post to add it. Please keep this in mind for future posts.

ADD REPLY
0
Entering edit mode

if you have one SNP per line, you can use split command in linux. you can split the file either by number of lines per output file or number of output files you want.

ADD REPLY
1
Entering edit mode

If I'm not mistaken some of these plink files are binary, would that work with split?

ADD REPLY
0
Entering edit mode

i am not sure if binary files can be split with split. Got confused bed with binary ped. @ WouterDeCoster

ADD REPLY
0
Entering edit mode

What is the next step, what is the purpose? I worked with millions of SNPs per chromosome and never really hit the wall using plink.

ADD REPLY
0
Entering edit mode

The reason is that I'm running a software called "Trinculo" to perform GWAS on an ordinal variable with multiple categories and the calculations are slow, so I need to split the files into chunks to parallelize the calculation.

ADD REPLY
0
Entering edit mode

there are tools/script to split chromosome wise. Can you not use them? kakukeshi

ADD REPLY
0
Entering edit mode

Out of curiosity, why would someone ask how to do something if they already knew how to do it or which tools to use?

Also, the OP is not asking chromosome-wise they're asking snp-wise.

ADD REPLY

Login before adding your answer.

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