This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Looping over plink files

Hi, I have .bed .bim. fam files for various chromosomes of some individuals and I would like merge with other plink files of some other individuals again for each chromosomes I would like to write a loop that does the merging automatically for all the files in just one command line Something like this:

for file in /home/users/chrs/*.{bed,bid,fam};
do plink --bfile filename1.chr$ --bmerge filename2.chr$.bed filename2.chr$.bim filename2.chr$.fam --out filename1_filename2_merged;
done

Here $ means any number from 1 to 22 for each files of the chromosomes.

Any help will be highly appreciated,thanks

bash shell linux

you need another inner loop: for j in seq 1 22

I see, but about the $ sign, it is ok to write like this?will it take any number from 1 to 22 ?

for j in $(seq 1 22); do echo chr$j; done

Hi, sorry, I tried the two loops but it did not work, can you be more specific please?

Could you write what did you try?

I see,I have got it Thanks

I see, yes, I know about that option but I have not tried it. I will do that then. Cheers

0 answers

No answers yet.

Log in to answer this question.