Hi Kimberly.
What version of phytools are you using? Neither make.simmap nor
describe.simmap call MatrixExp.
All the best, Liam
Liam J. Revell
Associate Professor, University of Massachusetts Boston
Profesor Asistente, Universidad Católica de la Ssma Concepción
web: http://faculty.umb.edu/liam.revell/, http://www.phytools.org
On 9/15/2018 10:38 AM, Kimberly Lynn Foster wrote:
Hello, I am looking at a family of fish and habitat transitions using SIMMAP in
phytools.
I keep running into the error below when using a subsample of 500 trees. I used
compare.tree and check.names to make sure all my tips match my data and visa
versa. Any help would be greatly appreciated.
simtrees<-list(500)
YY<-list(500)
for (i in 1:500){
+ simtrees[[i]]<-make.simmap(CrenuchidTreesLadderized[[i]], habitattrait,
model="ER", nsim=1000)
+ YY[[i]]<- describe.simmap(simtree[[i]], plot=F)
+ }
Error in MatrixExp(Y) : "mat" must be a square matrix
Kim
[[alternative HTML version deleted]]
___
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
___
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/