From henrik@bengt@@on @ending from gm@il@com  Tue Oct  2 18:57:14 2018
From: henrik@bengt@@on @ending from gm@il@com (Henrik Bengtsson)
Date: Tue, 2 Oct 2018 09:57:14 -0700
Subject: [Rd] maximum matrix size
In-Reply-To: <bb6fa2$a7n2cg@ironport10.mayo.edu>
References: <bb6fa2$a7n2cg@ironport10.mayo.edu>
Message-ID: <CAFDcVCQgXAzupyHB2ufvM=3oWi_h0dpjxvf42VV55CmOyUp9HQ@mail.gmail.com>

On Tue, Oct 2, 2018 at 9:43 AM Therneau, Terry M., Ph.D. via R-devel
<r-devel at r-project.org> wrote:
>
> I am now getting the occasional complaint about survival routines that are not able to
> handle big data.   I looked in the manuals to try and update my understanding of max
> vector size, max matrix, max data set, etc; but it is either not there or I missed it (the
> latter more likely).   Is it still .Machine$integer.max for everything?   Will that
> change?   Found where?

FWIW, this is the reference I've decided to follow for matrixStats:

"* For now, keep 2^31-1 limit on matrix rows and columns."

from Slide 5 in Luke Tierney's 'Some new developments for the R
engine', June 24, 2012
(http://homepage.stat.uiowa.edu/~luke/talks/purdue12.pdf).

/Henrik

>
> I am going to need to go through the survival package and put specific checks in front
> some or all of my .Call() statements, in order to give a sensible message whenever a
> bounday is struck.  A well meaning person just posted a suggested "bug fix" to the github
> source of one routine where my .C call allocates a scratch vector, suggesting  "resid =
> double( as.double(n) *nvar)" to prevent a "NA produced by integer overflow" message,  in
> the code below.   A fix is obvously not quite that easy :-)
>
>          resid <- .C(Ccoxscore, as.integer(n),
>                  as.integer(nvar),
>                  as.double(y),
>                  x=as.double(x),
>                  as.integer(newstrat),
>                  as.double(score),
>                  as.double(weights[ord]),
>                  as.integer(method=='efron'),
>                  resid= double(n*nvar),
>                  double(2*nvar))$resid
>
> Terry T.
>
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel



[Rd] maximum matrix size
Henrik Bengtsson