ADWGS_test | R Documentation |
ADWGS_test executes the statistical test for ActiveDriverWGS
ADWGS_test( id, gr_element_coords, gr_site_coords, gr_maf, win_size, this_genome, detect_depleted_mutations = FALSE )
id |
A string used to identify the element of interest. |
gr_element_coords |
A GenomicRanges object that describes the elements of interest containing the chromosome, start and end coordinates, and an mcols column corresponding to id |
gr_site_coords |
A GenomicRanges object that describes the sites of interest which reside in the elements of interest containing the chromosome, start and end coordinates, and an mcols column corresponding to id. Examples of sites include transcription factor binding sites in promoter regions or phosphosites in exons of protein coding genes. An empty GenomicRanges object nullifies the requirement for sites to exist. |
gr_maf |
A GenomicRanges object that describes the mutations in the dataset containing the chromosome, start and end coordinates, patient id, and trinucleotide context |
win_size |
An integer indicating the size of the background window in base pairs that is used to establish the expected mutation rate and respective null model. The default is 50000bps |
this_genome |
The reference genome object of BSgenome, for example BSgenome.Hsapiens.UCSC.hg19::Hsapiens |
detect_depleted_mutations |
if TRUE, detect elements with significantly fewer than expected mutations. FALSE by default |
A data frame containing the following columns
A string identifying the element of interest
The p-value of the element
The number of patients with a mutation in the element
The expected number of patients with a mutation in the element with respect to background
A boolean indicating whether the element is enriched in mutations
The p-value of the site
The number of patients with a mutation in the site
The expected number of patients with a mutation in the site with respect to element
A boolean indicating whether the site is enriched in mutations
A numeric indicator denoting the order in which the results were calculated
The FDR corrected p-value of the element
The FDR corrected p-value of the site
A V indicates the presence of site mutations
library(GenomicRanges) # Regions data(cancer_genes) gr_element_coords = GRanges(seqnames = cancer_genes$chr, IRanges(start = cancer_genes$start, end = cancer_genes$end), mcols = cancer_genes$id) # Sites (NULL) gr_site_coords = GRanges(c(seqnames=NULL,ranges=NULL,strand=NULL)) # Reference genome this_genome = BSgenome.Hsapiens.UCSC.hg19::Hsapiens # Mutations data(cll_mutations) cll_mutations = format_muts(cll_mutations, this_genome = this_genome) gr_maf = GRanges(cll_mutations$chr, IRanges(cll_mutations$pos1, cll_mutations$pos2), mcols=cll_mutations[,c("patient", "tag")]) # ADWGS_test id = "ATM" result = ADWGS_test(id, gr_element_coords, gr_site_coords, gr_maf, win_size = 50000, this_genome = this_genome)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.