[#6719] Propagation of NAs sometimes differing from base R

Date:
2021-04-13 04:19 Priority:
3

State:
Open

Submitted by:
David Cortes (david.cortes) Assigned to:
Martin Maechler (mmaechler)

Hardware:
None Product:
None

Operating System:
None Component:
None

Version:
None Severity:
None

Resolution:
None

URL:

Summary:
Propagation of NAs sometimes differing from base R

Detailed description
Not sure if this is worth fixing, but there are some mismatches in the propagation of NAs (particularly NA0) between base R and Matrix when multiplying a "dgCMatrix" by a numeric vector with number of elements matching to the number of rows of the matrix. I'm attaching serialized files as I haven't been able to pinpoint what characteristics does the matrix need to have in order to trigger this bug. library(Matrix) X <- readRDS("X.rds") v <- readRDS("v.rds") res_Matrix <- (X v)[1,1:5] res_R <- (as.matrix(X) * v)[1,1:5] Matrix: [1] 0 0 0 0 0 Base R: [1] NA NA NA NA NA

Comments:

Message
Date: 2022-06-16 15:38 Sender: Martin Maechler Moved from Bugs to Feature Requests
Date: 2021-04-17 08:21 Sender: Martin Maechler Thank you David. This is actually quite on purpose: A somewhat conscious decision to keep sparse matrices sparse. It needs to be documented however, and if users / collaborators are interested I'd open a "wishlist" item to ask for a global option (still turned off by default) that would make Matrix pkg behavior less sparse but compatible with base R, even in this boundary case.

Attachments:

Size	Name	Date	By	Download
1 KiB	X.rds	2021-04-13 04:19	david.cortes	X.rds
504 bytes	v.rds	2021-04-13 04:19	david.cortes	v.rds

Field	Old Value	Date	By
type	Bugs	2022-06-16 15:38	mmaechler
Operating System	None	2021-04-17 08:21	mmaechler
Severity	None	2021-04-17 08:21	mmaechler
assigned_to	none	2021-04-17 08:21	mmaechler
File Added	5205: v.rds	2021-04-13 04:19	david.cortes
File Added	5204: X.rds	2021-04-13 04:19	david.cortes