Vcftools r2 Error: Insufficient sites remained after filtering or No data left for analysis!
1
5
Entering edit mode
5.2 years ago
hellbio ▴ 440

I am having the below error while calcuating r2 using vcftools:

vcftools --vcf GATK.QC.MAF5.recode.vcf --hap-r2

VCFtools - v0.1.12b
(C) Adam Auton and Anthony Marcketta 2009

Parameters as interpreted:
--vcf GATK.QC.MAF5.recode.vcf
--max-alleles 2
--min-alleles 2
--hap-r2
--phased

After filtering, kept 133 out of 133 Individuals
Outputting Pairwise LD (phased bi-allelic only)
Error: Insufficient sites remained after filtering

I tried using --geno-r2 --phased and met with the obvious error as the variants are unphased.

vcftools --vcf 133Samples.GATK.QC.MAF5.recode.vcf --geno-r2  --phased

VCFtools - v0.1.12b
(C) Adam Auton and Anthony Marcketta 2009

Parameters as interpreted:
--vcf 133Samples.GATK.QC.MAF5.recode.vcf
--geno-r2
--max-alleles 2
--min-alleles 2
--phased

After filtering, kept 133 out of 133 Individuals
Outputting genotype pairwise LD (bi-allelic only) for a set of SNPs versus all others.
After filtering, kept 0 out of a possible 9792407 Sites
No data left for analysis!
Run Time = 444.00 seconds

Could anyone comment on how to calculate r2 using vcftools?

vcftools r2 LD • 3.6k views
ADD COMMENT
0
Entering edit mode

There was a bug in some of these functions in VCFtools 0.1.12. Perhaps a place to start would be to use 0.1.14 (the latest)?

ADD REPLY
0
Entering edit mode

In my case, I saw that I had a non-phased individual that had been merged with the phased ones and all sites got removed as a result.

ADD REPLY
3
Entering edit mode
4.7 years ago

hello, i have this problem too. Did you solve it yet?

ADD COMMENT
0
Entering edit mode

Did you solve it? I'm facing the same problem.

ADD REPLY
1
Entering edit mode

vcftools --vcf 133Samples.GATK.QC.MAF5.recode.vcf --geno-r2 It's OK --phased is not nessessary

ADD REPLY
0
Entering edit mode

its means that the output of that command show the r2 for each SNPs ? and with this information I can delete SNPs with r2 > 0.7?

ADD REPLY

Login before adding your answer.

Traffic: 841 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6