Question: ReadVCF not reading samples
gravatar for Alexander Skates
4.9 years ago by
United Kingdom
Alexander Skates360 wrote:

Hi, I'm trying to read a VCF file I have using readVcf, however the following occurs:

> fl <- file.path('tmp', 'vcf', 'rs11780156.genotypes.vcf.gz')
> vcf <- readVcf(fl, 'hg19')
> scanVcfHeader(fl)
class: VCFHeader
samples(190): HG00096 HG00097 ... NA12889 NA12890
meta(2): META contig
fixed(2): FILTER ALT
> dim(vcf)
[1] 13161     0
> geno(vcf)
List of length 0

I've looked at the file itself, and it definitely has the sample genotypes in it, so I'm not sure why geno(vcf) is empty.

The file in question is at however it seems to apply to any that I try and load via readVcf. These are just VCF I files I downloaded from 1000 Genomes using tabix.

variantannotation readvcf R • 1.6k views
ADD COMMENTlink modified 2.2 years ago by Biostar ♦♦ 20 • written 4.9 years ago by Alexander Skates360

Sorry, broken link:

ADD REPLYlink written 4.9 years ago by Alexander Skates360

Please paste output of sessionInfo(), gonna guess that you are missing a parameter in the read

ADD REPLYlink modified 4.8 years ago • written 4.8 years ago by Ying W4.0k

Hi Ying, it turned out that the issue was 1000 Genomes (the original source for the data) were missing the FORMAT lines in the header. I got in touch with them and they fixed it, so this is no longer an issue.

ADD REPLYlink written 4.8 years ago by Alexander Skates360
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: 1949 users visited in the last hour