Question: how to read in vcf.gz and tbi file for filtervcf using VariantAnnotation in R
0
gravatar for ansakim20
14 days ago by
ansakim200
ansakim200 wrote:

Hi I am new to this forum and bioinformatics, I am trying to understand the functionality of the VariantAnnotation package in Bioconductor R, in the filterVcf vignette the files are read in from the package as follows

###################################################
### code chunk number 5: createFilteredFile
###################################################
file.gz     <- system.file("extdata", "chr7-sub.vcf.gz", 
                           package="VariantAnnotation")
file.gz.tbi <- system.file("extdata", "chr7-sub.vcf.gz.tbi", 
                           package="VariantAnnotation")

can someone please explain, how do I read in my own vcf and tbi file?

Following is the link to the R code

https://bioconductor.org/packages/release/bioc/vignettes/VariantAnnotation/inst/doc/VariantAnnotation.R

snp R • 74 views
ADD COMMENTlink modified 13 days ago by genomax91k • written 14 days ago by ansakim200
0
gravatar for ansakim20
13 days ago by
ansakim200
ansakim200 wrote:

I figured it out since the file was in my machine, I just had to assign the path of that file

myfile <- path to my file
ADD COMMENTlink modified 13 days ago by genomax91k • written 13 days ago by ansakim200
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: 1712 users visited in the last hour