Re: [Rd] Matrix / SparseM conflict (PR#8618)

From: Duncan Murdoch <murdoch_at_stats.uwo.ca>
Date: Mon 20 Feb 2006 - 14:05:44 GMT

On 2/20/2006 8:54 AM, dpleydel@univ-fcomte.fr wrote:
> Full_Name: David Pleydell
> Version: 2.2.1
> OS: Debian Etch
> Submission from: (NULL) (193.55.70.206)
>
>
> There appears to be a conflict between the chol functions from the Matrix and
> the SparseM packages. chol() can only be applied to a matrix of class dspMatrix
> if SparseM is not in the path.

This isn't a bug, it's simply that both Matrix and SparseM define a generic for chol (though Matrix actually gets it by converting the function in base).

If you want access to the version in Matrix after attaching SparseM ahead of it in the search path, use Matrix::chol.

 > library(Matrix)
 > sm <- as(as(Matrix(diag(5) + 1), "dsyMatrix"), "dspMatrix")
 > library(SparseM)

[1] "SparseM library loaded"
 > chol(sm)
Error in chol(sm) : no applicable method for "chol"  > Matrix::chol(sm)
5 x 5 Matrix of class "pCholesky"

      [,1] [,2] [,3] [,4] [,5]

[1,] 1.4142136 0.7071068 0.7071068 0.7071068 0.7071068
[2,]         . 1.2247449 0.4082483 0.4082483 0.4082483
[3,]         .         . 1.1547005 0.2886751 0.2886751
[4,]         .         .         . 1.1180340 0.2236068
[5,]         .         .         .         . 1.0954451

You might want to discuss this with the SparseM package maintainers; perhaps they don't really need to define their own generic.

Duncan Murdoch

>
> with gratitude
> David
>
>

>> library(Matrix)
>> sm <- as(as(Matrix(diag(5) + 1), "dsyMatrix"), "dspMatrix")
>> chol(sm)

> 5 x 5 Matrix of class "pCholesky"
> [,1] [,2] [,3] [,4] [,5]
> [1,] 1.4142136 0.7071068 0.7071068 0.7071068 0.7071068
> [2,] . 1.2247449 0.4082483 0.4082483 0.4082483
> [3,] . . 1.1547005 0.2886751 0.2886751
> [4,] . . . 1.1180340 0.2236068
> [5,] . . . . 1.0954451
>> >> library(SparseM)

> [1] "SparseM library loaded"
>> chol(sm)

> Error in chol(sm) : no applicable method for "chol"
>> 
>> detach("package:SparseM")
>> chol(sm)

> 5 x 5 Matrix of class "pCholesky"
> [,1] [,2] [,3] [,4] [,5]
> [1,] 1.4142136 0.7071068 0.7071068 0.7071068 0.7071068
> [2,] . 1.2247449 0.4082483 0.4082483 0.4082483
> [3,] . . 1.1547005 0.2886751 0.2886751
> [4,] . . . 1.1180340 0.2236068
> [5,] . . . . 1.0954451
>>

>
> ______________________________________________
> R-devel@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel


R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Tue Feb 21 01:08:16 2006

This archive was generated by hypermail 2.1.8 : Tue 21 Feb 2006 - 03:44:41 GMT