View source: R/HoneyBADGER_gexp.R
HMM to identify CNV boundaries using normalized gene expression data
1 2 3 | calcGexpCnvBoundaries(gexp.norm, genes, m = 0.15, chrs = paste0("chr",
c(1:22)), min.traverse = 3, t = 1e-06, min.num.genes = 3, trim = 0.1,
verbose = FALSE)
|
gexp.norm |
Normalized gene expression matrix |
genes |
GRanges annotation of gene names and coordinates |
m |
Expression magnitude deviation needed to distinguish CNV from neutral |
chrs |
List of chromosome names. Genes not mapping to these chromosomes will be excluded. Default autosomes only: paste0('chr', c(1:22)) |
min.traverse |
Depth traversal to look for subclonal CNVs. Higher depth, potentially smaller subclones detectable. (default: 2) |
t |
HMM transition parameter. Higher number, more transitions. (default: 1e-6) |
min.num.genes |
Minimum number of genes within a CNV. (default: 3) |
trim |
Trim boundary SNPs |
verbose |
Verbosity (default: FALSE) |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.