From dpleydel at univ-fcomte.fr  Mon Feb 20 14:54:53 2006
From: dpleydel at univ-fcomte.fr (dpleydel@univ-fcomte.fr)
Date: Mon, 20 Feb 2006 14:54:53 +0100 (CET)
Subject: [Rd] Matrix / SparseM conflict (PR#8618)
Message-ID: <20060220135453.F08793F201@slim.kubism.ku.dk>

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.

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
>



[Rd] Matrix / SparseM conflict (PR#8618)
dpleydel@univ-fcomte.fr