Hi Graeme, I am not sure what you want to do with this but...it is straightforward to disallow reverse transitions in diversitree
## Simulate a dollo character library(diversitree) set.seed(1) t <- tree.bd(c(1,0), max.taxa=100) d <- sim.character(t, c(0.1, 0), model="mk2") ## make a likelihood function lik <- make.mk2(t, d) ## constrain the likelihood so that it only allows forward transitions lik_con <- constrain(lik, q10~0) ## estimate the MLE fit <- find.mle(lik_con, 0.1) ## use parameters to do ASR using ML parameters. Other operations also possible. asr <- asr.joint(lik_con, fit$par) Hope this helps, Matt On Mon, Sep 14, 2015 at 5:00 PM, Graeme Lloyd <[email protected]> wrote: > Hi Everyone, > > I was wondering if anyone knows if it is possible to set up a custom rate > model to represent a binary discrete Dollo character with ace or > make.simmap? > > Thanks, > > > Graeme > _______________________________________________ > R-sig-phylo mailing list - [email protected] > https://stat.ethz.ch/mailman/listinfo/r-sig-phylo > Searchable archive at > http://www.mail-archive.com/[email protected]/ > [[alternative HTML version deleted]] _______________________________________________ R-sig-phylo mailing list - [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-phylo Searchable archive at http://www.mail-archive.com/[email protected]/
