Entering edit mode
7.2 years ago
Shicheng Guo
★
9.6k
Any one have hg19 map for hapmap 3 ped file?
I downloaded the hapmap 3 data from following link. I think some formal download service can be provided since all of us will use it, right?
they only provided hg18 ped and map files. So I transfer hg18 map to hg19 map with lifeOver and trim plink again.
wget http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/liftOver
wget http://hgdownload.soe.ucsc.edu/goldenPath/hg18/liftOver/hg18ToHg19.over.chain.gz
plink --bfile hapmap3_r1_b36_fwd_consensus.qc.poly.recode --recode --tab --out hapmap3_r1_b36_fwd_consensus.qc.poly.recode
awk '{print "chr"$1,"\t",$4,"\t",$4+1,"\t",$2}' hapmap3_r1_b36_fwd_consensus.qc.poly.recode.map > hapmap3.hg18.bed
./liftOver hapmap3.hg18.bed hg18ToHg19.over.chain.gz hapmap3.hg19.bed unmap
awk '{print $1,"\t",$4,"\t",0,"\t",$2}' hapmap3.hg19.bed > hapmap3_r1_b37_fwd_consensus.qc.poly.recode.map
perl -p -i -e 's/chr//g' hapmap3_r1_b37_fwd_consensus.qc.poly.recode.map
mv hapmap3_r1_b36_fwd_consensus.qc.poly.recode.ped hapmap3_r1_b37_fwd_consensus.qc.poly.recode.ped
plink --file hapmap3_r1_b37_fwd_consensus.qc.poly.recode --make-bed --out hapmap3_r1_b37_fwd_consensus.qc.poly.recode
========================================================
Anyone who is interesting in hapmap 3 hg19 based plink files (ped, map and binary form: fam, bim, bed) can download from this link:
https://pan.baidu.com/s/1ElJ-heuW9jZXDeU9R-dZqg
========================================================