Your matrix is not symmetric, positive definite. If you don't know what this means, you shouldn't be using chol()
This may be because it isn't to begin with, or due to numerical error, it doesn't behave as one in the decomposition. My relative ignorance of numeric methods for linear algebra prevents me from saying more than that. -- Bert On Thu, Jun 14, 2012 at 4:23 AM, <nata...@orchidpharma.com> wrote: > Dear friends, > > When I do Cholesky decomposition for a 15x15 matrix using the function > chol(), I get the following error for which I do not understand the meaning > of the error > > Error in chol.default(M_cov) : > the leading minor of order 10 is not positive definite > > When I searched online for similar error reported earlier I could get few > hits but not of much help to resolve my error and one post suggested to use > different function called sechol() from accuracy package but that did not > work and it leads to different errors. So I want to stick to function chol() > itself. > > Could you please help me to find where things are going wrong in my matrix? > > > Thanks and regards, > B.Natarj > > ______________________________________________ > 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. -- Bert Gunter Genentech Nonclinical Biostatistics Internal Contact Info: Phone: 467-7374 Website: http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm ______________________________________________ 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.
