Search
Question: Vranges (somatic signature) VariantAnnotation
0
4.2 years ago by
Brazil
humaasif790 wrote:

Hi all,
I tried to make VRanges and am getting this error

> library(VariantAnnotation)
> fl <- system.file('UVAcl1_haplotype_snps_only.vcf', package="VariantAnnotation")
Error in .io_check_exists(path(con)) : file(s) do not exist:
''

Can any body help me with this error.
i am using this manual Introduction to VariantAnnotation ,Valerie Obenchain, March 18, 2014
is there any you tube video, any workshop video or tutorial how to construct VRanges

Regards

modified 4.2 years ago by Steve Lianoglou12k • written 4.2 years ago by humaasif790
0
4.2 years ago by
Switzerland
federico.comoglio100 wrote:

Hi humaasif79,

your problem is that the fl variable you defined is an empty string.

if your goal is simply to build a VCF file, I found a different vcf file in the extdata directory of the package that you could read in using readVcf.

Try:

fl <- system.file('extdata', 'gl_chr1.vcf', package="VariantAnnotation")
vcf <- readVcf(fl, "human_hg19.fa")

Kind Regards,

Federico

0
4.2 years ago by
United States
James W. MacDonald48k wrote:

The call to

fl <- system.file('UVAcl1_haplotype_snps_only.vcf', package="VariantAnnotation")

assumes that there is a file called UVAcl1_haplotype_snps_only.vcf in the package directory of VariantAnnotation, which I seriously doubt is true. Perhaps the file is in your working directory? In which case you just want to do

vcf <- readVcf("UVAcl1_haplotype_snps_only.vcf", "hg19")

Note that you cannot follow the vignettes verbatim, as they often rely on data that are supplied with the package to give examples. You must also consult the help pages for the function you are using to see what the expected inputs should be.