vcf.venn: Create a venn diagram of SNV overlap in a group of VCF files.

Description Usage Arguments Value Examples

Description

There is a maximum number of samples (nine) that can be included in a venn diagram (this limit is imposed by the R package producing the figures).

Usage

1
vcf.venn(vcf.files, genome, sample.names = NULL, ...)

Arguments

vcf.files

A character vector of VCF file to process.

genome

A character string identifying the genome

sample.names

A character vector of samples names for each of the VCF files. Must be the same length as vcf.files. By default the VCF files names are used as sample names.

...

Additionally parameters to pass to read.vcf.

Value

A list containing:

venn

An object of class Venn, see documentation of package Vennerable for further details and methods.

data

The data list that was used to create the Venn instance.

Examples

1
2
3
4
5
6
7
8
9
## Not run: 
f1 <- system.file('extdata', 'LC1_TUMOUR_A.vcf', package = 'GenomicVis')
f2 <- system.file('extdata', 'LC1_TUMOUR_B.vcf', package = 'GenomicVis')
vcf.files <- c(f1, f2)
sample.names <- c('LC1_A', 'LC1_B')
v <- vcf.venn(vcf.files, 'GRCh37', sample.names)
plot(v$venn)

## End(Not run)

jjellis/GenomicVis documentation built on May 19, 2019, 11:39 a.m.