The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.

fpCompare

Alex M. Chubaty

August 13, 2022

Reliable comparison of floating point numbers

Comparisons of floating point numbers are problematic due to errors associated with the binary representation of decimal numbers. Computer scientists and programmers are aware of these problems (.e.g., Goldberg 1991) and yet people still use numerical methods which fail to account for floating point errors (this pitfall is the first to be highlighted in the First Circle of “The R Inferno” (Burns 2012)).

To avoid these and other numerical rounding issues, R’s help file for relational operators (e.g., ?'>') suggests using identical and all.equal when making numerical comparisons:

x1 <- 0.5 - 0.3
x2 <- 0.3 - 0.1
x1 == x2                           # FALSE on most machines
identical(all.equal(x1, x2), TRUE) # TRUE everywhere

Inspired by R FAQ 7.31 and this Stack Overflow answer, this package provides new relational operators useful for performing floating point number comparisons with a set tolerance:

fpCompare1 base
%>=% >=
%>>% >
%<=% <=
%<<% <
%==% ==
%!=% !=

These functions use the base relational operators to make comparisons, but incorporate a tolerance value (fpCompare.tolerance) similar to all.equal. The default fpCompare.tolerance value is .Machine$double.eps^0.5, set via options. This is the same default used in all.equal for numeric comparisons.

# set telorance value
tol = .Machine$double.eps^0.5       # default value
options(fpCompare.tolerance = tol)

# perform comparisons
x1 <- 0.5 - 0.3
x2 <- 0.3 - 0.1
x1 == x2         # FALSE on most machines
x1 %==% x2       # TRUE everywhere

Installation

From CRAN

install.packages("fpCompare")

From GitHub

library(devtools)
install_github("PredictiveEcology/fpCompare")

Bug Reports

https://github.com/PredictiveEcology/fpCompare/issues

References

Burns, Patrick. 2012. The R Inferno. 2nd ed. https://www.burns-stat.com/pages/Tutor/R_inferno.pdf.
Goldberg, David. 1991. What every computer scientist should know about floating-point arithmetic.” ACM Computing Surveys, 171–264. https://dl.acm.org/doi/10.1145/103162.103163.

  1. The %<<% and %>>% symbols are used instead of %<% and %>% to avoid a conflict with magrittr’s pipe operator (%>%).↩︎

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.