Re: [R] Matrix::solve() with 1-d arrays -- treating "array" as "numeric"?

2021-04-19 Thread Martin Maechler
> Deepayan Sarkar > on Mon, 19 Apr 2021 09:56:58 +0530 writes: > On Sat, Apr 17, 2021 at 9:08 PM Martin Maechler > wrote: >> >> > Deepayan Sarkar > on Fri, 16 Apr 2021 11:34:20 >> +0530 writes: >> >> > I get what I initially thought was

Re: [R] Matrix::solve() with 1-d arrays -- treating "array" as "numeric"?

2021-04-18 Thread Deepayan Sarkar
On Sat, Apr 17, 2021 at 9:08 PM Martin Maechler wrote: > > > Deepayan Sarkar > > on Fri, 16 Apr 2021 11:34:20 +0530 writes: > > > I get what I initially thought was unexpected behaviour from: > > > x <- tapply(runif(100), sample(5, 100, TRUE), mean) > > solve(Diagonal(5),

Re: [R] Matrix::solve() with 1-d arrays -- treating "array" as "numeric"?

2021-04-17 Thread Martin Maechler
> Deepayan Sarkar > on Fri, 16 Apr 2021 11:34:20 +0530 writes: > I get what I initially thought was unexpected behaviour from: > x <- tapply(runif(100), sample(5, 100, TRUE), mean) > solve(Diagonal(5), x) > # Error: not-yet-implemented method for solve(, ). > #

[R] Matrix::solve() with 1-d arrays

2021-04-16 Thread Deepayan Sarkar
I get what I initially thought was unexpected behaviour from: x <- tapply(runif(100), sample(5, 100, TRUE), mean) solve(Diagonal(5), x) # Error: not-yet-implemented method for solve(, ). # ->> Ask the package authors to implement the missing feature. This is because x is a 1-D array, so the