QTLtools pca detects no varaints
0
0
Entering edit mode
5.8 years ago
GK1610 ▴ 120

Hi I have two vcf files from different sets of individuals

1.vcf file is genotype data and 2.vcf files is wgs data

when I run

QTLtools pca --vcf 1.vcf --scale --center --maf 0.05 --distance 50000 --out 1.genotypes.output
QTLtools pca --vcf 2.vcf --scale --center --maf 0.05 --distance 50000 --out 2.genotypes.output

separately on each file I am able get pcas but when I run QTLtools pca on merged vcf file. I dont get any pcs. I get this message. Can anyone help me with this?

bcftools merge 1.vcf 2.vcf > 1.2.vcf

QTLtools pca --vcf 1.2.vcf --scale --center --maf 0.05 --distance 50000 --out 1.2.genotypes.output

Reading genotype data in [test.vcf.gz]
  * File format detected: VCF
  * Read 100000 lines
  * 0 variants included
  * 126480 variants excluded by user

EXITED: Cannot find genotypes in target region!
snp • 1.4k views
ADD COMMENT

Login before adding your answer.

Traffic: 2451 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