Skip to content

phylogenetic correction and MCMC model

2 messages · David Costantini, Bert Gunter

#
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 on your website 
Thanks in advance 
David
#
Such specialized questions are usually better posted on appropriate R-sigs,
in this case, https://stat.ethz.ch/mailman/listinfo/r-sig-phylo I presume.

Bert Gunter

"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )


On Mon, May 17, 2021 at 7:33 AM David Costantini <david.costantini at mnhn.fr>
wrote: