search for: brobdingnag

Hi. I am having difficulty following section 1.6.6 of the R-extensions manual. I am trying to update the Brobdingnag package to include a NAMESPACE file (the untb package requires the Brobdingnag package). Without the NAMESPACE file, the package passes R CMD check cleanly. However, if I include a NAMESPACE file, even an empty one, R CMD check gives the following error in 00install.out: wt118:~/packages% cat...

Hello List. please find uploaded to CRAN a new package, Brobdingnag. This package does two things: (1) allows computation of very large numbers using a logarithmic representation. (2) provides a "Hello, World" example of S4 methods in use: there are two classes of object (brob and glub) and one virtual class (swift). The package includes a vi...

...[1] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN <snip> > is.numeric(x^B) [1] TRUE > is.real(x^B) [1] TRUE > is.infinite(x^B) [1] FALSE FALSE FALSE FALSE FALSE I am lost, I checked the R mailing help, but could not find anything directly. I loaded package Brobdingnag and tried: as.brob(x^B) [1] NAexp(NaN) NAexp(NaN) NAexp(NaN) NAexp(NaN) NAexp(NaN) > as.brob(x)^B [1] NAexp(187.67) NAexp(187.65) NAexp(187.63) NAexp(187.61) I guess I must be misunderstanding something fundamental. Any clues would be really appreciated. Willem

I think Robin knows about FAQ 7.31/floating point (author of 'Brobdingnag', among other numerical packages). I agree that this is surprising (to me). To reframe this question: is there way to get an *exact* ASCII representation of a numeric value (i.e., guaranteeing the restored value is identical() to the original) ? .deparseOpts has ?"digits17"?: R...

HI R Gurus, > This is my first foray into using c-code with R, so ... > I had a look at the archives and did not find anything on this, so > hopefully I am not doubling up. > I have previously used the following approach where I needed some very small numbers/large (using Brobdingnag): surfacewithdiff <- function(t, y, p) { const=p["const"] kay =p["kay"] psii=p["psii"] m1=p["m1"] m2=p["m2"] M=p["M"] N=p["N"] T=p["T"] bexp<-as.brob((const*(psii-m1*y[1]-m...

Hi R used with C-code experts, I had a look at the archives and did not find anything on this, so hopefully I am not doubling up. I have previously used the following approach where I needed some very small/large numbers (using Brobdingnag): surfacewithdiff <- function(t, y, p) { const=p["const"] kay =p["kay"] psii=p["psii"] m1=p["m1"] m2=p["m2"] M=p["M"] N=p["N"] T=p["T"] bexp<-as.brob((const*(psii-m1*y[1]-m...

...far as I know I can't use "apply") -- this produces a list of brobs. My problem is, how do I act on, plot this list, or do vector calculations with the elements? From the other emails I have read there does not seem to be a means of indexing a range of elements from a list. library(Brobdingnag) t<-c(as.brob(5), as.brob(6), as.brob(7), as.brob(8), as.brob(10)) > t [[1]] [1] +exp(1.6094) [[2]] [1] +exp(1.7918) [[3]] [1] +exp(1.9459) [[4]] [1] +exp(2.0794) [[5]] [1] +exp(2.3026) > t[2:4] [[1]] [1] +exp(1.7918) [[2]] [1] +exp(1.9459) [[3]] [1] +exp(2.0794) So if you try p...

...2.6.0 (it is major.minor.patchlevel). This is most easily done by > update.packages(checkBuilt=TRUE, ask=FALSE) However, if you don't want to do that yet, be aware that - Certain S4-using packages must be reinstalled, and using old versions can make R malfunction. I believe these are Brobdingnag, Matrix, NADA and kappalab. - Any S4 package that makes an existing function S4-generic will grab that function as it existed in the version of R under which they were installed. If it has changed, there is potential trouble. - Packages that create or change character strings at C level will...

...od and I will experiment with this. Best wishes Robin On Sun, Mar 1, 2020 at 6:22 AM Duncan Murdoch <murdoch.duncan at gmail.com> wrote: > > On 29/02/2020 4:19 a.m., Ben Bolker wrote: > > > > I think Robin knows about FAQ 7.31/floating point (author of > > 'Brobdingnag', among other numerical packages). I agree that this is > > surprising (to me). > > > > To reframe this question: is there way to get an *exact* ASCII > > representation of a numeric value (i.e., guaranteeing the restored value > > is identical() to the origina...

Dear R, I have two questions: 1, Why both R and Matlab give 0*Inf==NaN? To my knowledge, it should be zero mathematically. Am I right? 2, I need to calculate e.g. exp(a)/(exp(b)+c), where both a and b are very large numbers (>>1000, e.g a=1000, b=1007, and c=5). R gives me NaN when I use the following command: > exp(1000)/(exp(1007)+5) [1] NaN I am pretty sure this should be close to

...6:22 AM Duncan Murdoch > > <murdoch.duncan at gmail.com> wrote: > >> > >> On 29/02/2020 4:19 a.m., Ben Bolker wrote: > >> > > >> > I think Robin knows about FAQ 7.31/floating point > >> (author of > 'Brobdingnag', among other numerical > >> packages). I agree that this is > surprising (to me). > >> > > >> > To reframe this question: is there way to get an > >> *exact* ASCII > representation of a numeric value (i.e., > >>...

On 29/02/2020 4:19 a.m., Ben Bolker wrote: > > I think Robin knows about FAQ 7.31/floating point (author of > 'Brobdingnag', among other numerical packages). I agree that this is > surprising (to me). > > To reframe this question: is there way to get an *exact* ASCII > representation of a numeric value (i.e., guaranteeing the restored value > is identical() to the original) ? > > .depar...

I'm looking for some tutorial on S4. I only find the following one, which is not in English. Can somebody let me know if there is any introductory material? I'm very familiar with OO and C++. If there is some material that suits my background, it will be great. https://stat.ethz.ch/pipermail/r-help/2009-January/184108.html

...mewhere). More generally, there are methods for computing sums on the log scale, e.g. https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.misc.logsumexp.html I don't know where this has been implemented in the R ecosystem, but this sort of computation is the basis of the "Brobdingnag" package for operating on very large ("Brobdingnagian") and very small ("Lilliputian") numbers. On 2019-06-21 6:58 p.m., jing hua zhao wrote: > Hi Peter, Rui, Chrstophe and Gabriel, > > Thanks for your inputs -- the use of qnorm(., log=TRUE) is a good point in...

Hi everyone I have a problem about "infinite". If I type 10^308, R shows "1e+308" When I type 10^309, R shows "Inf" So, we know if a value is large than 1.XXXe+308, R will show "Inf" How can i do let the value, like "10^400" ,typed in R to show the word "1e+400" not "Inf" -- View this message in context:

...n > On Sun, Mar 1, 2020 at 6:22 AM Duncan Murdoch > <murdoch.duncan at gmail.com> wrote: >> >> On 29/02/2020 4:19 a.m., Ben Bolker wrote: >> > >> > I think Robin knows about FAQ 7.31/floating point >> (author of > 'Brobdingnag', among other numerical >> packages). I agree that this is > surprising (to me). >> > >> > To reframe this question: is there way to get an >> *exact* ASCII > representation of a numeric value (i.e., >> guaranteeing the restored val...

Hi the list. I try to learn the S4 programming. I find the wiki and several doc. But I still have few questions... 1. To define 'representation', we can use two syntax : - representation=list(temps = 'numeric',traj = 'matrix') - representation(temps = 'numeric',traj = 'matrix') Is there any difference ? 2. 'validityMethod' check the

...doch >> > <murdoch.duncan at gmail.com> wrote: >> >> >> >> On 29/02/2020 4:19 a.m., Ben Bolker wrote: >> >> > >> >> > I think Robin knows about FAQ 7.31/floating point >> >> (author of > 'Brobdingnag', among other numerical >> >> packages). I agree that this is > surprising (to me). >> >> > >> >> > To reframe this question: is there way to get an >> >> *exact* ASCII > representation of a numeric value (i.e., &g...

Hi, I'm trying to calculate p-value to findout definitely expressed genes compare A to B situation. I got this data(this is a part of data) from whole organism , and each number means each expression values (that means, we could think 'a' gene is 13 in A situation, and it turns 30 in B situation) To findout probability, I'm going to use Audic - Claverie Method. ( The significance

Hello, I have a text file with one column containing long number but stored as string. I download the file with read.table (and colClass) and the first row of this column is : "095842087016731010" As I need to make some calculations with these numbers, I tried to convert them using as.numeric. But then I get as.numeric("095842087016731010") 95842087016731008 I understand