Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Bug in info of musse #27

Open
KYPpa opened this issue Oct 29, 2020 · 0 comments
Open

Bug in info of musse #27

KYPpa opened this issue Oct 29, 2020 · 0 comments

Comments

@KYPpa
Copy link

KYPpa commented Oct 29, 2020

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.)
name.pretty="MuSSE", # pretty name for printing
## Parameters:
np=as.integer(k * (k + 2)), # number of ODE parameters

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.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant