Re: [R-sig-phylo] BiSSE models on multiphylo objects

2015-12-06 Thread Daniel Caetano
Hi Kevin, Although diversitree was written in a modular and efficient way, it is very challenging to make modifications to the code. Functions in the package are usually very general and make it hard to understand their role in the MCMC algorithm. My opinion is that it is better to use the lik

Re: [R-sig-phylo] Stuck with Phylomatic tree!

2015-03-05 Thread Daniel Caetano
)alsophila)cyatheaceae)))cyatheales)))monilophyte)euphyllophyte;" Then you can use the function 'read.newick' in phytools followed by the function 'collapse.singles'. Please also note that the newick tree need to be in a single line (same string). Hope it helps, *Daniel Caetano

[R-sig-phylo] Simulate multivariate OU process

2014-11-12 Thread Daniel Caetano
Theoretical Biology) processes and I also searched for something in 'ouch'. Someone know a function that simulate such a process? At the moment I am considering implementing the mvOU simulation using the already implemented mvBM simulation in 'mvSLOUCH'. Thanks for the help,

Re: [R-sig-phylo] Error in diversitree with MuSSE models (Arbuckle, Kevin)

2014-04-03 Thread Daniel Caetano
f elements equal to the number of parameters in your full model. To find the MLE of your constrained likelihood function you need to point to the free parameters in 'start'. Try this, check if start has only "lambda1","mu1","q12": fit.allcon <- find.mle(lik.allc