Diffbind - "dba.report" showing values with less significant figures?
1
0
Entering edit mode
wesley.cai • 0
@wesleycai-11862
Last seen 6.8 years ago

Hi,

I have been using diffbind for a while now, and I use this code to generate a results file:

d_all.DB <- dba.report(d_all, contrast=1, method=DBA_DESEQ2, th=1, DataType=DBA_DATA_FRAME)
write.table(d_all.DB, "results.txt", sep="\t", eol="\n", quote=FALSE, row.names=FALSE)

I used to get 8+ significant figures in the resulting values, but now it only shows 3. Is there a way to specify the sigfigs?

Thanks!

Wes

diffbind dba.report • 763 views
ADD COMMENT
0
Entering edit mode
Rory Stark ★ 5.2k
@rory-stark-5741
Last seen 14 days ago
Cambridge, UK

Hi Wes-

Yes indeed, a change was made to limit the number of significant digits for the main stats in the report. This was done for two reasons:

  1. Reporting 8+ significant digits is inherently misleading. While the computer is perfectly happy to calculate these digits, it gives a false sense of precision in the calculations. When we are dealing with whole read counts, often <100, calculating fold changes (especially log fold changes) out that far is not really meaningful. For confidence statistics, while we retain the magnitude in the form of "e-07" etc, the specific digits after that can't really be taken seriously. (Indeed, calculations of p-values beyond about the e-04 and lower rate should be taken with a grain of salt). Basically, too many digits encourages an illusion of precision that is not real.
  2. When each number has 8 digits to the right of the decimal place, the vast majority of information contained in the report is actually insignificant, making it more difficult to focus on the meaningful aspect of the values. It is certainly a lot easier on the eye to look at tables of numbers with fewer digits!

I would consider adding a feature in future allowing you to specify the number of significant digits as this would be an easy feature to add.

Regards-

Rory

ADD COMMENT

Login before adding your answer.

Traffic: 1037 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6