-
Notifications
You must be signed in to change notification settings - Fork 9
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
Prob of shift proportional to branch length maybe broke on dev version? #11
Comments
Hello, I am running into the same error. I am using version 2.1. If I just keep re-running it, it will manage to run for a period of time and then throw an error again. Would really like to use branch lengths in the prior but the errors worrying me. See below for a very long example of the output. Thanks for your help! `Checking inputs for errors:...........
|
Resolved this issue in version 2.1.3 (hopefully) |
Runs that previously worked on the CRAN version of bayou are no longer working for me on the dev version. These runs use the probability of a shift proportional branch length prior outlined in some of the documentation. The runs work for a few generations, but ultimately crash. For instance,
prior2 <- make.prior(prunedMBC, dists=list(dalpha="dhalfcauchy", dsig2="dhalfcauchy",dsb="dsb", dk="cdpois", dtheta="dnorm"), param=list(dalpha=list(scale=1), dsig2=list(scale=1), dk=list(lambda=15, kmax=200), dsb=list(bmax=Inf,prob=prunedMBC$edge.length), dtheta=list(mean=mean(avMBC), sd=2)))
mcmc2 <- bayou.makeMCMC(prunedMBC, avMBC, SE=cvMBC, model="OU", prior=prior2, new.dir=TRUE, samp=10^2, chunk=10^3, ticker.freq=10^4, plot.freq=NULL)
system.time(mcmc2$run(ngen=10^7))
returns
Perhaps something needs to get fixed with this model on the dev version? Or, perhaps I've misspecified something in the code (entirely possible).
The text was updated successfully, but these errors were encountered: