Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
CovarianceMatrix function error (Read 1495 times)
Danny
YaBB Newbies
*
Offline



Posts: 31
CovarianceMatrix function error
Apr 13th, 2009 at 12:18pm
 
From ErrorEstimates.bf with the command

CovarianceMatrix (cmx, theLnLik);

I'm getting the error:

Error:LUDecompose doesn't work on singular matrices (row 33)
Call stack
1 : Calculate the Covariance Matrix storing into, cmx, the following likelihood function:theLnLik ;

I get the same error with the MPI, MP and SP versions of the svn downloaded around March 3, 2009.

I reduced the number of taxa from around 80 to around 25 and I didn't get the error.  I'm guessing there is a branch that has an omega at a high or low extreme with zero variance and covariance?
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: CovarianceMatrix function error
Reply #1 - Apr 13th, 2009 at 1:07pm
 
Dear Danny,

Yeah, sounds like the numerical Hessian is singular.
Try using the profile likelihood approach.

Cheers,
Sergei
Back to top
 

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego
WWW WWW  
IP Logged