Version 3.4.0 uses a separate calculation for 1 - E, e.g. one minus the local extinction rates to avoid numerical abberations. These are only used in the CLA calculations.
3.3.0Version 3.3.0 normalizes evaluation of the loglikelihood at every integration step. This strongly reduces any numerical abberations that could previously occur. These mainly occured along long branches, where values could become so close to zero that they became too close to numerical precision. If absolutely desired, this normalization can be disabled by setting âuse_normalizationâ to FALSE.
3.2.0Added the function âplot_idparslistâ, which provides code to visualize the rates in the idparslist setup, using ggplot.
3.1.0Version 3.1.0 fixes a bug in the simulation code that caused trait changes during speciation to not be tracked appropriately. This could, for instance, interfere with conditioning and this bug especially impacted ClaSSE-type simulations.
Minor changesVersion 3.0.1 patches some inaccuracies in simulation functions, and deprecates expand_q_matrix, as this was making some incorrect assumptions.
Breaking changespool_init_states
Version 3.0.0 extends the C++ code base used for the standard likelihood to the âcla_â likelihood, harnessing the same computation improvement.
Breaking changescreate_lambda_matrices()
is now called create_lambda_list()
create_transition_matrix()
is now called create_q_matrix()
create_mus()
is now called create_mu_vector()
create_default_q_list()
is now called create_default_shift_matrix()
create_default_lambda_list ()
is now called create_default_lambda_transition_matrix()
create_default_q_list()
is now called create_default_shift_matrix()
phylo_Vign
is now called phylo_vignette
traitinfo
is now called traits
phy
is now called example_phy_GeoSSE
plot_state_exact()
argument steps
renamed to num_steps
and argument focal_tree
renamed to phy
for consistency with other functions.NEWS.md
file to track changes to the package.default_params_doc()
.vignette("complete_tree", package = "secsse")
vignette("starting_secsse", package = "secsse")
.secsse_sim()
argument conditioning
now defaults to "obs_states"
from "none"
.secsse_sim()
fix bug causing error when simulating trees with extinct species.fill_in()
, create_default_q_list()
, create_default_transition_list()
, create_mus()
Version 2.5.0 appeared in 2021 on GitHub and was published in May 2023 on CRAN. Version 2.5.0 marks the first version using C++ to perform the integration, and it used tbb (from the RcppParallel package) to perform multithreading. This marks a ten fold increase in speed over previous versions. Secondly, 2.5.0 introduces the function secsse_sim()
to simulate a diversification process using the (cla) secsse framework. Lastly, in version 2.5.0 functions were added to allow visualisation of inferred rates of speciation across the tree (e.g. plot_state_exact()
and secsse_loglik_eval()
).
Version 2.0.0 appeared in June of 2019 on CRAN and extended the package with the cla framework, e.g. including state shifts during speciation / asymmetric inheritance during speciation.
1.0.0The first version of secsse appeared in January of 2019 on CRAN. It used the package deSolve to solve all integrations, and could switch between either using a fully R based evaluation, or use FORTRAN to speed up calculations. Furthermore, using the foreach package, within-R parallelization was implemented. However, parallelization only situationally improved computation times, and generally, computation was relatively slow.
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4