Question: Two Plink map/ped files refuse to merge.
2
gravatar for devenvyas
3.8 years ago by
devenvyas570
Stony Brook
devenvyas570 wrote:

Two map/ped files refused to merge for some reason.

I followed this method here (Plink: Quick Way To Combine Two Datasets Using Only Common Markers by ff.cc.cc) to reduce the datasets to the same list of sites.

I used 

plink --bfile Deven --bmerge Harvard.fam Harvard.bim Harvard.bed --merge-mode 6 --out mismatchB

to identify the sites that had flipped genotypes, and I straightened them out using

plink --bfile Harvard --flip mismatchB.missnp --recode  --make-bed --out HarvardFlipped

I converted them to map/ped from the binary format, and I tried to merge them and got an inexplicable map file error

ERROR: Problem with MAP file line:
HGDP00001    HGDP00001    0    0    1    -9    G T    A A    C C    G G    A G    G G    A A    C T    G G    A G    T T    T T    G G    C C    C C    A A    T C    C C    C C    G G    G G    T T    T T    G G    A G    T T    C C    A A    A G    C T    A G    C C    G G    A A    G G    C C    T C    C T    A G    A G    G G    G G    A A    T C    T T    T T    G G    G G    A G    G G    C C    T C    G A    G G    G G    C C    T

I looked at the map and ped files and everything was fine. Any clue what is going on? I know that the Harvard dataset originally has coordinates using Hg18 and the Deven dataset has coordinates the using Hg19 (but the rs ids come from the same annotation file, which provided the conversion of Affy #s to dbSNP #s)

 

 

 

plink • 2.5k views
ADD COMMENTlink modified 9 weeks ago by zx87547.1k • written 3.8 years ago by devenvyas570
1
gravatar for devenvyas
3.8 years ago by
devenvyas570
Stony Brook
devenvyas570 wrote:

UPDATE: I merged them as binary, and then converted them to map/ped and it worked... You all can probably ignore this.

ADD COMMENTlink written 3.8 years ago by devenvyas570
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1564 users visited in the last hour