Question: Merging multiple .ped/.map files for PLINK analisys purposes
gravatar for giusdalt95
9 weeks ago by
giusdalt950 wrote:

Good moring everyone! I am trying to merge 5 couple of flat-format .ped/.map files with PLINK command:

plink --noweb --file 53519559 --merge-list allfiles.txt --recode --out merged

but the terminal gives me back this error:

Found 54 SNPs that do not match in terms of allele codes
Might include strand flips, although flipped A/T and C/G SNPs will be undetected)
Writing problem SNPs to [ merged.missnp ]

ERROR: Stopping due to mis-matching SNPs -- check +/- strand?

N.B: I am due to specify that I started from 5 VCF files, from which I filtered for SNPs, removed duplicate positions and filtered out SNPs with missing rsID; subsequently I used vcftools to convert the filtered vcf files in plink format, through the code:

vcftools --gzvcf my.vcf.gz --plink --out myfileconverted

If someone could help me I would be very grateful! P.S: obviously don't hesitate to ask if you need more infos about.

vcftools snp plink • 156 views
ADD COMMENTlink modified 8 weeks ago by chrchang5237.6k • written 9 weeks ago by giusdalt950
gravatar for chrchang523
8 weeks ago by
United States
chrchang5237.6k wrote:

See the discussion of merge failures at .

ADD COMMENTlink written 8 weeks ago by chrchang5237.6k
Please log in to add an answer.


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