[R] Ancestral state reconstruction under a split-tree for BiSSE

Hello! 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 nodes.i<-match(nodes,phy$node.label)+length(phy$tip.label) pars.b<-c(428.597, 90.777, 421.878, 81.815, 0.201, 2.900) #p

Re: [R] Discrete trait Ornstein–Uhlenbeck in R?

Thanks for the feedback. By discrete, I meant a categorical trait (diel activity). I've looked at the discrete models in geiger, but I had to ask about Ornstein-Uhlenbeck because one of my advisors suggested it. -- View this message in context: http://r.789695.n4.nabble.com/Discrete-trait-Ornst

[R] Discrete trait Ornstein–Uhlenbeck in R?

Is there a package that will allow me to fit Brownian motion and Ornstein–Uhlenbeck models of evolution for discrete traits? I know that geiger and ouch have commands for fitting these models for continuous traits, but these aren't suitable for discrete trait evolution, correct? -- View this