Change rsID to chromosome:location in bim files for plink
2
1
Entering edit mode
2.6 years ago

Some of the SNPs in my dataset are without rs ids, so in the bim file, they have dots instead of rs ids for them. So for errorless analysis how should replace the rs id column with that of the chromosome: location column instead?

genome snp plink • 4.2k views
ADD COMMENT
5
Entering edit mode
2.6 years ago
zx8754 10k

Using awk, (not tested):

awk '{if($2 == ".") {print $1,$1"_"$3,$4,$5,$6} else {print $0}}' myfile.bim > myfileClean.bim
ADD COMMENT
5
Entering edit mode
2.6 years ago

This can also be done with plink 2.0's --set-all-var-ids flag.

But it's more valuable to understand the awk solution, since that opens the door to a broadly useful set of text-manipulation techniques.

ADD COMMENT

Login before adding your answer.

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