Entering edit mode
In the DiffBind package, how can I extract the t-statistic from the dba.report?
In the DiffBind package, how can I extract the t-statistic from the dba.report?
DiffBind does not use t-statistics to determine the confidence in assessing differential binding. T-tests assume that the read counts can be accurately modelled using a normal distribution, which is an invalid assumption for these data. Instead we use an FDR-corrected p-value computed using exact tests from the DESeq2 and/or edgeR packages.
These confidence statistics can be retrieved as the $FDR values in the report returned by dba.report():
data(tamoxifen_analysis)
report <- dba.report(tamoxifen, th=.01)
report$FDR
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.