[R] Using Rmpfr to round with precision in R

Martin Maechler maechler at stat.math.ethz.ch
Thu Nov 27 10:10:39 CET 2014

```>>>>> Artur Augusto <arturaugusto at gmail.com>
>>>>>     on Wed, 19 Nov 2014 15:26:12 -0200 writes:

> I'm trying to use the Rmpfr library with the round() function to apply
> the round half to even rule and achieve correct results, without
> errors due the finite precision of float point values.

> Example of problem:

> round(1.225,2)
> #[1] 1.23

> So far, this is what I've achieved:

> library(Rmpfr)
> x <- c(1.225, 1.2225, 1.22225, 1.222225)
> n <- c(2, 3, 4, 5)
> round2 <- function(x, n){
> sprintf(paste("%#.", n, "f", sep=""), round(mpfr(as.character(x), 200), n))
> }
> mapply(round2, x, n)
> #[1] "1.22"    "1.222"   "1.2222"  "1.22222"

The above  round2() function is really pretty strange [what do you really want?]
but below,  you have a very good point :

> But in some cases I don't get the desired results:

> round2(1.152, 2)# Should be 1.15
> #[1] "1.16"

> Reading the Rmpfr docs, at the roundMpfr() function, it says:

Actually,  roundMpfr() is *not* used  by  round(<mpfr>, *)
It has somewhat different semantics. That's why I used a new
function name, roundMpfr(),  whereas the round() method for
"mpfr" objects is what you get from

selectMethod(round, "mpfr")

> The mpfr class group method Math2 implements a method for round(x,
> digits) which rounds to decimal digits.

Indeed... and that is different from the roundMpfr() as
mentioned above.

> But I can't figure how to use it.

> How can I achieve desired rounded results?
> Artur

Wait a day or two [depending on how quickly the new Rmpfr version 0.5-7
gets published on CRAN], and update your Rmpfr package.
The problem is fixed in Rmpfr 0.5-7.

If you had CC'ed your e-mail to
maintainer("Rmpfr")
I would have heard of the problem earlier and almost surely
have fixed it earlier.

Best regards,
Martin Maechler  (maintainer of 'Rmpfr').

```

More information about the R-help mailing list