An Error In The Estimation Of Ancestral State With R Package(Ape)
Hi, I am trying to estimate the ancestral characters in phylogenetic analysis with the R package(ape). it did well in the estimation of continuous characters, only with a warning as follows:
> #*In sqrt(1/out$hessian) : 产生了NaNs*
but returned an error when estimating a discrete characters. the error is showed as follows:
#*> ace(data_input[,2],tree_input,'d')
#错误于E(Q * phy$edge.length[i]) , this means "error in E(Q * phy$edge.length[i])"
# BLAS/LAPACK 例行程序'DGEEV '产生了错误-13信息*, this means there is an error message(-13) in the program DGEEV of BLAS/LAPACK.
the phylogeny inference tree was selected randomly in our previous work and the trait data was created rndomly too.
can some give some advice to solve this problem?
• 2,367 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Hard to tell much from this error message, other than the fact it's a LAPACK library throwing the error. It might become clearer if you describe/show some of your data. You might also ask at the R special-interest group mailing list for PCMs.
http://stackoverflow.com/questions/4442518/general-suggestions-for-debugging-r