rational solution in R [Software]

posted by mittyri – Russia, 2020-12-21 14:49 (777 d 14:38 ago) – Posting: # 22152
Views: 6,293

Hi Helmut,

❝ At least consistently giving not zero. ;-)

library(Rmpfr)

❝ prec <- 64

❝ x <- mpfr(0.5, prec) - mpfr(0.4, prec) - mpfr(0.1, prec)

❝ x

1 'mpfr' number of precision  64   bits

❝ [1] -2.77555756156289135106e-17

❝ y <- mpfr(-0.4, prec) - mpfr(0.1, prec) + mpfr(0.5, prec)

❝ y

1 'mpfr' number of precision  64   bits

❝ [1] -2.77555756156289135106e-17



I think the right solution is
install.packages("remotes")
remotes::install_github("bertcarnell/rational")
library(rational)
(rational(5L, 10L) - rational(4L, 10L) - rational(1L, 10L)) == 0
[1] TRUE


Interestingly but I did not find good solution using CRAN

Kind regards,
Mittyri

Complete thread:

UA Flag
Activity
 Admin contact
22,481 posts in 4,710 threads, 1,603 registered users;
19 visitors (0 registered, 19 guests [including 8 identified bots]).
Forum time: 05:28 CET (Europe/Vienna)

The difference between a surrogate and a true endpoint
is like the difference between a cheque and cash.
You can get the cheque earlier but then,
of course, it might bounce.    Stephen Senn

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5