try Cholesky() in package Matrix. albyn
On Tue, Mar 10, 2009 at 02:33:01PM -0700, Manli Yan wrote: > Hi everyone: > I try to use r to do the Cholesky Decomposition,which is A=LDL',so far I > only found how to decomposite A in to LL' by using chol(A),the function > Cholesky(A) doesnt work,any one know other command to decomposte A in to > LDL' > > My r code is: > library(Matrix) > A=matrix(c(1,1,1,1,5,5,1,5,14),nrow=3) > > > chol(A) > [,1] [,2] [,3] > [1,] 1 1 1 > [2,] 0 2 2 > [3,] 0 0 3 > > > Cholesky(A) > Error in function (classes, fdef, mtable) : > unable to find an inherited method for function "Cholesky", for signature > "matrix" > > whatz wrong??? > thanks~ > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.