diversitree
diversitree copied to clipboard
Bug in info of musse
Edit of the argnames of musse objects fails like below.
argnames(lik)[1]<-c("l") Error in
argnames<-.dtlik
(*tmp*
, value = c("l", "lambda2", "lambda3", : Wrong argnames length: expected 15, got 12
Then I checked the source and found in "model-musse.R"
make.info.musse <- function(k, phy) {
list(name="musse", # canonical name (as in make.
The parameters should be composed of k lambda, k mu and k*(k-1) transition rates. Then np=as.integer(k * (k + 1)) seems to be correct.