Question: Merging multiple .ped/.map files for PLINK analisys purposes
0
gravatar for giusdalt95
9 weeks ago by
giusdalt950
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
1
gravatar for chrchang523
8 weeks ago by
chrchang5237.6k
United States
chrchang5237.6k wrote:

See the discussion of merge failures at https://www.cog-genomics.org/plink/1.9/data#merge3 .

ADD COMMENTlink written 8 weeks ago by chrchang5237.6k
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: 2453 users visited in the last hour