Ancestral state reconstruction under a split-tree for BiSSE
I am trying to do ancestral reconstruction under a split BiSSE model.
#phy is my tree
nodes<-c(755,620,602,448,6,340) #vector of nodes at which to split the phylogeny
pars.b<-c(428.597, 90.777, 421.878, 81.815, 0.201, 2.900) #pars were taken from the non-split model
All of that performs as expected, but when I try to perform ancestral reconstruction on my split model, I get this error:
Why does it expect 6 parameters for my split model instead of the 42? Am I using the wrong
function for the ancestral reconstruction, or is the code incorrect? The updates for version 0.9-1
of diversitree state that reconstruction under split models is availabe.
It was also suggested to me that I could average the parameters from fit.split and do asr with an unsplit tree.
Is this viable?