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.