Hi David, It looks like phylo_ultra might be a list? Is phylo_ultra[[1]] a tree?
Also, don't use nodes="TIPS"; this is just to demonstrate how poor the algorithm is when you don't use the expanded inverse. I see people using nodes="TIPS" a lot - where does this code come from? Cheers, Jarrod On 17/05/2021 16:06, David Costantini wrote:
This email was sent to you by someone outside the University. You should only click on links or attachments if you are certain that the email is genuine and the content is safe. Dear All I am trying to apply a phylogenetic correction to an MCMC model, but I have problems in making the inverse matrix. I can visualise the treeplot very well, but when I use the script: inv.phylo<-inverseA(phylo_ultra,nodes="TIPS",scale=TRUE) R tells me that there is an error: Error in pedigree[, 2] : incorrect number of dimensions In addition: Warning message: In if (attr(pedigree, "class") == "phylo") { : Do you have any experience with this? I couldn't find a solution so far. Thanks in advance David [[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/
The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. Is e buidheann carthannais a th’ ann an Oilthigh Dhùn Èideann, clàraichte an Alba, àireamh clàraidh SC005336. _______________________________________________ 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/