Re: [Rd] na.omit inconsistent with is.na on list

2021-08-13 Thread Iñaki Ucar
On Thu, 12 Aug 2021 at 22:20, Gabriel Becker wrote: > > Hi Toby, > > This definitely appears intentional, the first expression of > stats:::na.omit.default is > >if (!is.atomic(object)) > > return(object) I don't follow your point. This only means that the *default* method is not int

Re: [Rd] na.omit inconsistent with is.na on list

2021-08-13 Thread Hugh Parsonage
The data.frame method deliberately skips non-atomic columns before invoking is.na(x) so I think it is fair to assume this behaviour is intentional and assumed. Not so clear to me that there is a sensible answer for list columns. (List columns seem to collide with the expectation that in each varia

[Rd] svd For Large Matrix

2021-08-13 Thread Dario Strbenac via R-devel
Good day, I have a real scenario involving 45 million biological cells (samples) and 60 proteins (variables) which leads to a segmentation fault for svd. I thought this might be a good example of why it might benefit from a long vector upgrade. test <- matrix(rnorm(4500*60), ncol = 60) test

Re: [Rd] [External] svd For Large Matrix

2021-08-13 Thread luke-tierney
[copying the list] svd() does support matrices with long vector data. Your example works fine for me on a machine with enough memory with either the reference BLAS/LAPACK or the BLAS/LAPACK used on Fedora 33 (flexiblas backed, I believe, by a version of openBLAS). Take a look at sessionInfo() to

Re: [Rd] [External] svd For Large Matrix

2021-08-13 Thread Prof Brian Ripley
On 13/08/2021 15:58, luke-tier...@uiowa.edu wrote: [copying the list] svd() does support matrices with long vector data. Your example works fine for me on a machine with enough memory with either the reference BLAS/LAPACK or the BLAS/LAPACK used on Fedora 33 (flexiblas backed, I believe, by a ve

[Rd] , Problem in random number generation for Marsaglia-Multicarry +, Kinderman-Ramage

2021-08-13 Thread André GILLIBERT
Thank you. It looks like the M-M generator is deeply flawed, although R implementation seems to be faithful. I can understand that R keeps this generator for historical reasons. However, R should at least warn that this generator is buggy and should not be used! That could be as simple as

Re: [Rd] [External] svd For Large Matrix

2021-08-13 Thread Dario Strbenac via R-devel
Good day, Ah, I was confident it wouldn't be environment-specific but it is. My environment is R version 4.1.0 (2021-05-18) Platform: x86_64-pc-linux-gnu (64-bit) Running under: Debian GNU/Linux 10 (buster) Matrix products: default BLAS: /usr/lib/x86_64-linux-gnu/atlas/libblas.so.3.10.3 LAPAC

Re: [Rd] [External] svd For Large Matrix

2021-08-13 Thread Dirk Eddelbuettel
Dario, On 14 August 2021 at 00:00, Dario Strbenac via R-devel wrote: | Good day, | | Ah, I was confident it wouldn't be environment-specific but it is. My environment is | | R version 4.1.0 (2021-05-18) | Platform: x86_64-pc-linux-gnu (64-bit) | Running under: Debian GNU/Linux 10 (buster) |