Question: Problem in merging 1000G and cases in Plink
0
gravatar for fatima
3.0 years ago by
fatima10
fatima10 wrote:

I have some problems in stating the imputation. I am PhD student and I work on imputation in case-only design. My reference panel is build37 of 1000G (phase 3) in VCF format and a set of data including cases (in Plink format). before starting imputation in "Minimac3" prephasing with "shapeit" is necessary. "merging cases and controls(reference panel)" and "flip" them is necessary before prephasin.

But I have some problems in the first steps. I converted the ref. panel from VCF to Plink and tried to merge them with cases. I have this Error:

Warning: Multiple positions seen for variant 'rs201556956'.

.

.

Warning: Multiple positions seen for variant 'rs200991502'.

17838 markers loaded from CD_GermanyKielchr2_mod.bim.

7047141 markers to be merged from ref_b37_ph3.bim.

Of these, 7029359 are new, while 17782 are present in the base dataset.

Error: 7932 variants with 3+ alleles present.

  • If you believe this is due to strand inconsistency, try --flip with merge.missnp.

    (Warning: if this seems to work, strand errors involving SNPs with A/T or C/G alleles probably remain in your data. If LD between nearby SNPs is high, --flip-scan should detect them.)

  • If you are dealing with genuine multiallelic variants, we recommend exporting that subset of the data to VCF (via e.g. '--recode vcf'), merging with another tool/script, and then importing the result; PLINK is not yet suited to handling them.

Then I tried to remove them(but generally I do not want to remove them) and merge them but I have the other error:

Error: 7916 variants with 3+ alleles present.

  • If you believe this is due to strand inconsistency, try --flip with merge1.missnp.

    (Warning: if this seems to work, strand errors involving SNPs with A/T or C/G alleles probably remain in your data. If LD between nearby SNPs is high, --flip-scan should detect them.)

  • If you are dealing with genuine multiallelic variants, we recommend exporting that subset of the data to VCF (via e.g. '--recode vcf'), merging with another tool/script, and then importing the result; PLINK is not yet suited to handling them.

Now I would like to know what is the reason and in which file there is a problem as I have some errors which shows I have duplicate data! So I decided to check all the "multiple_position" SNP and "merge1.missnp" to explore if they are in the cases or in the reference panel.

Could you give me some advice how I can solve this level to go to next step(prephasing)? Is there some algorithms to check or better view points? As I am new in Bioinformatics I guess maybe there are better ways.

ADD COMMENTlink written 3.0 years ago by fatima10

These are common problems when merging PLINK files. SNP coordinates may vary between files due to different genome builds, and the allele coding may vary at some locations. Here is a previous post which may help;

Accounting for problem SNPS when merging multiple plink files

ADD REPLYlink written 3.0 years ago by christopher medway440

Hello fghaderinezhad!

It appears that your post has been cross-posted to another site: http://seqanswers.com/forums/showthread.php?p=199433

This is typically not recommended as it runs the risk of annoying people in both communities.

ADD REPLYlink written 3.0 years ago by WouterDeCoster41k
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: 909 users visited in the last hour