[R-sig-phylo] using ace() with fixed internal node - singularity error

Dear R-sig-phylo members, I am trying to use ace() with a fixed internal node as described here: https://stat.ethz.ch/pipermail/r-sig-phylo/2011-June/001465.html My aim is to provide some confidence in the ancestral state estimation of the root node by setting the root node to have specific

[R-sig-phylo] Error estimating RMA in phytools

Hi I am working on the correlation among traits of 36 plants. I am using the function phyl.RMA to estimate relationship between pairs of traits but I keep having a error message CNvrSRLphylo-phyl.RMA(Chem2mean[,1],Morpho2mean[,1],treeARP,method=lambda) Error in solve.default(kronecker(R, C)) :

Re: [R-sig-phylo] Error estimating RMA in phytools

Hi Oscar. My educated guess after some investigating is this is because the objects Chem2mean and Morpho2mean are data frames not matrices; and thus Chem2mean[,1] does not contain names. Try first doing: Chem2mean-as.matrix(Chem2mean), and the same with Morpho2mean, and then repeat the