powered by
MAF file is more recommended. In this function, we will mimic the MAF object from the key c(1, 2, 4, 5, 7) columns of VCF file.
c(1, 2, 4, 5, 7)
read_vcf( vcfs, samples = NULL, genome_build = c("hg19", "hg38", "T2T", "mm10", "mm9", "ce11"), keep_only_pass = FALSE, verbose = TRUE )
a MAF.
VCF file paths.
sample names for VCF files.
genome build version like "hg19".
if TRUE, keep only 'PASS' mutation for analysis.
TRUE
if TRUE, print extra info.
read_maf, read_copynumber
vcfs <- list.files(system.file("extdata", package = "sigminer"), "*.vcf", full.names = TRUE) # \donttest{ maf <- read_vcf(vcfs) maf <- read_vcf(vcfs, keep_only_pass = TRUE) # }
Run the code above in your browser using DataLab