An Error In The Estimation Of Ancestral State With R Package(Ape)

Question

[Ontheway]

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?

[David W]

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.

[zx8754]

http://stackoverflow.com/questions/4442518/general-suggestions-for-debugging-r