dba.analyze {DiffBind} | R Documentation |
Performs differential binding affinity analysis
dba.analyze(DBA, method=DBA$config$AnalysisMethod, bSubControl=TRUE, bFullLibrarySize=TRUE, bTagwise=TRUE, filter=0, filterFun=max, bCorPlot=DBA$config$bCorPlot, bReduceObjects=TRUE, bParallel=DBA$config$RunParallel)
DBA |
DBA object. If no contrasts are specified (DBA$contrast is NULL), default contrasts will be added via a call to |
method |
method, or vector of methods, by which to analyze differential binding affinity. Supported methods: also, for backward compatibility: Additionally, if this value is set to |
bSubControl |
logical indicating whether Control read counts are subtracted for each site in each sample before performing analysis. |
bFullLibrarySize |
logical indicating if the full library size (total number of reads in BAM/SAM/BED file) for each sample is used for scaling normalization. If FALSE, the total number of reads present in the peaks for each sample is used (generally preferable if overall biding levels are expected to be similar between samples). |
bTagwise |
logical indicating if dispersion should be calculated on a tagwise (or per-condition) basis. If there are only a very few members of each group in a contrast (e.g. no replicates), this should be set to FALSE. |
filter |
value to use for filtering intervals with low read counts.
Each contrast will be filtered separately.
The |
filterFun |
function that will be invoked for each interval with a vector of scores for each sample. Returns a score that will be evaluated against the |
bCorPlot |
logical indicating whether to plot a correlation heatmap for the analyzed data (first contrast only). If no sites are significantly differentially bound using the default thresholds, no heatmap will be plotted. |
bReduceObjects |
logical indicating whether strip the analysis objects of unnecessary fields to save memory. If it is desired to use the |
bParallel |
logical indicating that the analyses is to be done in parallel using multicore (one process for each contrast for each method, plus an additional process per method). |
See the DBA User Guide for more details on how the edgeR
and DESeq2
analyses are carried out.
DBA object with results of analysis added to DBA$contrasts.
If there is a blocking factor for the contrast(s) specified using a previous call to dba.contrast
, a multi-factor analysis will automatically be carried out in addition to a single factor analysis.
Rory Stark
data(tamoxifen_counts) tamoxifen <- dba.analyze(tamoxifen) tamoxifen data(tamoxifen_counts) tamoxifen <- dba.contrast(tamoxifen,categories=DBA_CONDITION,block=tamoxifen$masks$MCF7) tamoxifen <- dba.analyze(tamoxifen,method=DBA_ALL_METHODS) tamoxifen