[R] question on suppressing error messages with Rmath library
Luke Tierney
luke at stat.uiowa.edu
Wed Mar 21 16:04:05 CET 2007
You might get less noise in the replies if you were explicit about
using Rmath stand-alone and asked on r-devel.
As far as I can see you would need to compile a version of the
stand-alone library that defines the macros for handling of warning
messages differently -- the current one just calls printf in the
stand-alone library. (You might be able to trick the linker into using
a version of printf for calls from within Rmath that does nothing, but
I suspect recompiling the sourses is easier.) We will probably be
rethinking this soon in conjunction with some other changes to
vectorized math in R.
Best,
luke
On Wed, 21 Mar 2007, Ranjan Maitra wrote:
> Dear list,
>
> I have been using the Rmath library for quite a while: in the current instance, I am calling dnt (non-central t density function) repeatedly for several million. When the argument is small, I get the warning message:
>
> full precision was not achieved in 'pnt'
>
> which is nothing unexpected. (The density calls pnt, if you look at the function dnt.) However, to have this happen a huge number of times, when the optimizer is churning through the dataset is bothersome, but more importantly, a bottleneck in terms of speed. Is it possible to switch this off? Is there an setting somewhere that I am missing?
>
> Many thanks and best wishes,
> Ranjan
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
--
Luke Tierney
Chair, Statistics and Actuarial Science
Ralph E. Wareham Professor of Mathematical Sciences
University of Iowa Phone: 319-335-3386
Department of Statistics and Fax: 319-335-3017
Actuarial Science
241 Schaeffer Hall email: luke at stat.uiowa.edu
Iowa City, IA 52242 WWW: http://www.stat.uiowa.edu
More information about the R-help
mailing list