Find research datasets worth reusing

opencc-by-4.0Jul 2021View details →

Data and software associated with the paper "Bayesian Inference of Joint Coalescence Times of Sampled Sequences"

1. Data files and run logs produced for the paper "Bayesian Inference of Joint Coalescence Times of Sampled Sequences". 2. Software script versions used in the above.

opencc-by-4.0May 2016View details →

Fig. 6. Bayesian inference tree for 5,179 in A new species of Tritetrabdella (Hirudinida: Hirudiniformes: Haemadipsidae) from northern Indochina

Fig. 6. Bayesian inference tree for 5,179 bp alignment positions of nuclear 18S rRNA and 28S rRNA and mitochondrial cytochrome c oxidase subunit I markers. Numbers on nodes indicate bootstrap values for maximum likelihood and Bayesian posterior probabilities.

opencc-by-4.0Apr 2023View details →

Data of "Bayesian inference of high-dimensional finite-strain visco-elastic-visco-plastic model parameters for additive manufactured polymers and neural network based material parameters generator"

General Data of <a href="http://doi.org/10.1016/j.ijsolstr.2023.112470">https://doi.org/10.1016/j.ijsolstr.2023.112470</a> related to MOAMMM project. Data related to the publication (we would be grateful if you could cite the paper in the case in which you are using the data): title = "Bayesian inference of high-dimensional finite-strain visco-elastic-visco-plastic model parameters for additive manufactured polymers and neural network based material parameters generator.", journal = "International Journal of Solids and Structures", year = "2023", volume = "283", pages = "112470", doi = "10.1016/j.ijsolstr.2023.112470", author = "Ling Wu, Cyrielle Anglade, Lucia Cobian, Miguel Monclus, Javier Segurado, Fatma Karayagiz, Ubiratan Freitas, and Ludovic Noels" This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 862015. Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union. Neither the European Union nor the granting authority can be held responsible for them. Description BI code and results of the inference of a pressure-dependent visco-elastic visco-plastic model developed in [NGU16] with a umat implementation in <a href="https://gitlab.uliege.be/moammm/moammmPublic/code/-/tree/main/MaterialModels/FiniteStrain/Finite_VEVP">https://gitlab.uliege.be/moammm/moammmPublic/code/-/tree/main/MaterialModels/FiniteStrain/Finite_VEVP</a>. The BI is described in [WU23] .The experimental results used in the BI are reported in [COB22,COB22b]. To run the BI you need the open source code <a href="https://gitlab.onelab.info/cm3/cm3Libraries">https://gitlab.onelab.info/cm3/cm3Libraries</a> If you use these data or model, we would be grateful if you could cite the related papers. Bibliography <ul> <li>[WU23] L. Wu, C. Anglade, L. Cobian, M. Monclus, J. Segurado, F. Karayagiz, U. Santos Freitas, L. Noels, Bayesian inference of high-dimensional finite-strain visco-elastic-visco-plastic model parameters for additive manufactured polymers and neural network based material parameters generator, International Journal of Solids and Structures (2023) 112470: https://doi.org/10.1016/j.ijsolstr.2023.112470</li> <li>[COB22] L. Cobian, M. Rueda-Ruiz, J.P. Fernandez-Blazquez, V. Martinez, F. Galvez, F. Karayagiz, T. Lück, J. Segurado, M.A. Monclus, Micromechanical characterization of the material response in a PA12-SLS fabricated lattice structure and its correlation with bulk behaviour, Polymer Testing 110 (2022) 107556: https://doi.org/10.1016/j.polymertesting.2022.107556 (in Open access)</li> <li>[COB22b] Data of “. Cobian, M. Rueda-Ruiz, J.P. Fernandez-Blazquez, V. Martinez, F. Galvez, F. Karayagiz, T. Lück, J. Segurado, M.A. Monclus, Micromechanical characterization of the material response in a PA12-SLS fabricated lattice structure and its correlation with bulk behaviour, Polymer Testing 110 (2022) 107556” http://dx.doi.org/10.5281/zenodo.6136935 (in Open access)</li> <li>[NGU16] V. D. Nguyen, F. Lani, T. Pardoen, X. Morelle, L. Noels, A large strain hyperelastic viscoelastic-viscoplastic-damage constitutive model based on a multi-mechanism non-local damage continuum for amorphous glassy polymers. International Journal of Solids and Structures 96 (2016): 192-216; https://dx.doi.org/10.1016/j.ijsolstr.2016.06.008, Open access: https://orbi.uliege.be/handle/2268/197898</li> </ul> Directories All the codes and experimental results are in five directories: <ol> <li>experimentalTests: experimental data, see the README.txt in each subdirectory for details</li> <li>BayesianVE: BI of the visco-elastic parameters <ol> <li>PlotExperimentalCurves: to vizualize the experimental curves and prepare the observations for the BI in the VE range <ol> <li>Loadcase_H.py and Loadcase_V.py read experimental results and create Load_ExpVE_H.dat and Load_ExpVE_V.dat, which keep the experimental observations and loading conditions to perform the BI.</li> <li>PrintDir_H & PrintDir_V subdirectories with the functions called by Loadcase_H.py and Loadcase_V.py</li> <li>Load_ExpVE_H.dat and Load_ExpVE_V.dat created files with the observations and loading conditions to perform the BI</li> </ol> </li> <li>VE_V2Step and VE_H: BI for viscoelastic properties of "V" specimen (VE_V2Step) and "H" specimen (VE_H) <ol> <li>BI_allpos_sequence.py runs the BI using Predict_VETest.py and creates the MCMC_VE_....dat</li> <li>WarmStart = True is used to restart an inference</li> <li>MCMC_VE_....dat in the VE_V2Step and VE_H directories are the BI results</li> <li>When proceeding in two steps in VE_V2Step, a first step generates MCMC_VE_VN8_1st.dat whose posterior is used as prior in the second step to generate MCMC_VE_VN8_2nd.dat</li> </ol> </li> <li>CheckBayRes: to visualize predictions of a BI sample and experimental curves <ol> <li>MCMCRes.py is used to check the numerical predictions of a BI parameter sample (read last sample by default, V or H direction can be selected at line</li> <li>ResKGEmu.py plots the evolution of elastic properties with time</li> <li>uses as input VE_V2Step/MCMC_VE_....dat or VE_H/MCMC_VE_....dat</li> <li>uses local ViscoElasticTest.py, line.geo, line. msh as interface with https://gitlab.onelab.info/cm3/cm3Libraries code</li> <li>uses local functions plotExpLoad_Unload.py, plotExp.py</li> </ol> </li> <li>ViscoElasticTest.py, line.geo, line.msh: interface with <a href="https://gitlab.onelab.info/cm3/cm3Libraries">https://gitlab.onelab.info/cm3/cm3Libraries</a> code used by VE_V2Step and VE_H to call the VEVP model</li> </ol> </li> <li>BayesianVEVP: BI of the visco-elastic and visco-plastic parameters <ol> <li>PlotExperimentalCurves: to vizualize the experimental curves and prepare the observations for the BI in the VE-VP ranges <ol> <li>Loadcase_H.py and Loadcase_V.py read experimental results and create Load_ExpVEVP_H.dat and Load_ExpVEVP_V.dat, which keep the experimental observations and loading conditions to perform BI at the viscoplastic stage.</li> <li>PrintDir_H & PrintDir_V subdirectories with the functions called by Loadcase_H.py and Loadcase_V.py</li> <li>Load_ExpVEVP_H.dat and Load_ExpVEVP_V.dat created files with the observations and loading conditions to perform the BI</li> </ol> </li> <li>VP_V2step and VP_H2step: BI for viscoelastic-viscoplastic properties of "V" specimen (VP_V2Step) and "H" specimen (VP_H2Step) <ol> <li>BI_allpos_sequence.py runs the BI using Predict_VETest.py and creates the MCMC_VP_....dat</li> <li>WarmStart = True is used to restart an inference</li> <li>It starts from the VE prosterior as prior, see point 2, and generates a MCMC_VP_?_1of2Steps.dat (? being H or V)</li> <li>Then using MCMC_VP_?_1of2Steps.dat posterior to get a new prior, it generates MCMC_VP_?_2of2Steps.dat (? being H or V)</li> </ol> </li> <li>CheckBayRes: visualize predictions of a BI sample and experimental curves <ol> <li>MCMCRes.py is used to check the numerical predictions with 3 BI parameter samples ([28000, 45000,70000] by default, V or H direction can be selected at line 12) using the samples of BayesianVEVP/VP_?2step/MCMC_VP_?_2of2Steps.dat (? being H or V)</li> <li>plot_hist.py is used to plot histograms of all the inferred parameters using the samples of BayesianVEVP/VP_?2step/MCMC_VP_?_2of2Steps.dat (? being H or V)</li> <li>Plot_Prop.py plots joints histograms of the inferred parameters using the samples of BayesianVEVP/VP_?2step/MCMC_VP_?_2of2Steps.dat (? being H or V)</li> <li>ResKGEmu.py plots the evolution of elastic properties with time</li> </ol> </li> <li>VEVPTest.py: interface with https://gitlab.onelab.info/cm3/cm3Libraries code used by VP_V2Step and VP_H2Step to call the VEVP model</li> </ol> </li> <li>RandomParametersGenerator: used to generate the parameters from the BI samples, with the same statistical content <ol> <li>Generator <ol> <li>DataProcess.py: creates normalized data for training from final inferred parameters in ../MCMC_ResData and creates ?_dirNormData (? being H or V)</li> <li>KmeanDataProcess.py: performs clustering for the data of H_dirNormDat and creates H_dirNormData_2cluster (no need for V direction because not bimodal)</li> <li>Gan_V.py and Gan_H.py are used to train the random material parameter generators and create the VDir_Gan or HDir_Gan200_0/HDir_Gan200_1</li> <li>GenerateParameters.py generates random parameters using the Gan files VDir_Gan or HDir_Gan200_0/HDir_Gan200_1 and checks the joint histograms of generated parameters, generated parameters are in V_GenData and H_GenData</li> <li>Ganlib.py is used by the generator</li> </ol> </li> <li>CheckRes <ol> <li>GenDataRes.py is used to check the numerical predictions with the generated parameter samples, see point 4) (using V_GenData and H_GenData).</li> <li>Plot_PropGen.py plots joints histograms of the generated parameters using the samples of V_GenData or H_GenData</li> </ol> </li> </ol> </li> <li>MCMC_ResData:All final data used in the paper (they can substitute the ones used here above) <ol> <li>H_direction and V_direction keep the MCMC random walk results of BI.</li> <li>RandomParameterGenerator keeps results of the generator Paper</li> </ol> </li> </ol> Figures of [WU23] <ul> <li>Fig. 5: From directory BayesianVE/PlotExperimentalCurves, run python3 ./PrintDir_V/plotExp_T.py or ./PrintDir_V/plotExp_C.py or ./PrintDir_V/plotExp_R.py</li> <li>Fig. 7: BayesianVEVP/CheckBayRes/Plot_Prop.py with direct = "V" and then with direct = "H" and with Var = [0,1,20,24,28,29,30,31]</li> <li>Fig. 8: BayesianVEVP/CheckBayRes/MCMCRes.py with direct = "V" (requires <a href="https://gitlab.onelab.info/cm3/cm3Libraries">https://gitlab.onelab.info/cm3/cm3Libraries</a> code)</li> <li>Fig. 9: BayesianVEVP/CheckBayRes/MCMCRes.py with direct = "H" (requires<a href="https://gitlab.onelab.info/cm3/cm3Libraries"> https://gitlab.onelab.info/cm3/cm3Libraries</a> code)</li> <li>Fig. 11: RandomParametersGenerator/CheckRes/Plot_PropGen.py with direct = "V" and then with direct = "H" and with Var = [0,1,20,24,28,29,30,31]</li> <li>Fig. 12: RandomParametersGenerator/CheckRes/GenDataRes.py with direct = "V" (requires <a href="https://gitlab.onelab.info/cm3/cm3Libraries">https://gitlab.onelab.info/cm3/cm3Libraries</a> code)</li> <li>Fig. 13: RandomParametersGenerator/CheckRes/GenDataRes.py with direct = "H" (requires <a href="https://gitlab.onelab.info/cm3/cm3Libraries">https://gitlab.onelab.info/cm3/cm3Libraries</a> code)</li> <li>Fig. 14A: From directory BayesianVE/PlotExperimentalCurves, run python3 ./PrintDir_H/plotExp_T.py or ./PrintDir_H/plotExp_C.py or ./PrintDir_H/plotExp_R.py</li> <li>Fig. 15C: BayesianVEVP/CheckBayRes/Plot_Prop.py with direct = "V", Var = [2,3,8,9,14,15,18,19] and [20,21,22,23,24,25,26,27]</li> <li>Fig. 16C: BayesainVEVP/CheckBayRes/plot_hist.py with direct = "V"</li> <li>Fig. 17C: BayesianVEVP/CheckBayRes/plot_hist.py with direct = "V"</li> <li>Fig. 18C: BayesianVEVP/CheckBayRes/Plot_Prop.py with direct = "H", Var = [2,3,8,9,14,15,18,19] and [20,21,22,23,24,25,26,27]</li> <li>Fig. 19C: BayesianVEVP/CheckBayRes/plot_hist.py with direct = "H"</li> <li>Fig. 20C: BayesianVEVP/CheckBayRes/plot_hist.py with direct = "H"</li> <li>Fig. 21D: RandomParametersGenerator/CheckRes/Plot_PropGen.py with direct = "V" , Var = [2,3,8,9,14,15,18,19] and [20,21,22,23,24,25,26,27]</li> <li>Fig. 22D: RandomParametersGenerator/CheckRes/Plot_PropGen.py with direct = "H" , Var = [2,3,8,9,14,15,18,19] and [20,21,22,23,24,25,26,27]</li> </ul>

opencc-by-4.0Apr 2023View details →

The equation of state for neutron star matter has been obtained through Bayesian inference utilizing a relativistic mean field model with a non-linear mesonic interaction.

The equation of state for matter in neutron stars has been obtained through Bayesian inference utilizing a relativistic mean field model with a non-linear mesonic interaction. ---------------------------- Dr. Tuhin Malik Department of Physics, University of Coimbra tm@uc.pt Date: 22 Apr 2023 ----------------------------- The high density behavior of nuclear matter is analyzed within a relativistic mean field description with non-linear meson interactions. To assess  the model parameters and their output, a Bayesian inference technique is used. The Bayesian setup is limited only by a few nuclear saturation properties, the neutron star maximum mass larger than 2 M$_\odot$, and the low-density pure neutron matter equation of state (EOS) produced by an accurate N$^3$LO calculation in chiral effective field theory. Depending on the strength of the non-linear  scalar vector  field contribution, we have found three distinct classes of EOSs, each one correlated to different star properties distributions. If  the non-linear vector  field contribution is absent, the gravitational maximum mass and the sound velocity at high densities are the greatest. However, it also gives the smallest speed of sound at  densities below three times saturation density. On the other hand,  models with the strongest  non-linear vector  field contribution,  predict the largest radii and tidal deformabilities for 1.4 M$_\odot$ stars, together with  the smallest mass for the onset of the nucleonic direct Urca processes and the smallest central baryonic densities for the maximum mass configuration.  {These models have the largest speed of sound below three times saturation density, but the smallest at high densities, in particular, above four times saturation density the speed of sound decreases approaching approximately $\sqrt{0.4}c$ at the center of the maximum mass star. On the contrary, a weak non-linear vector contribution gives a monotonically increasing speed of sound.} {A 2.75 M$_\odot$ NS maximum mass was obtained in the tail of the posterior with a weak non-linear vector field interaction. This indicates that the secondary object in GW190814 could also be an NS. {The possible onset of hyperons  and the compatibility of the different sets of models with pQCD are discussed. It is shown that pQCD favors models with a large contribution from the non-linear vector  field term or which include hyperons.}} The article e-Print:  <a href="https://arxiv.org/abs/2301.08169">2301.08169</a> We release model parameters, its nuclear saturation properties, equation of state,  and TOV solutions derived from Bayesian Inference with Prior Set 0, 1, 2, and 3. We also share Set 0 with Hyperon.  For every Set, our data release packet contains four CSV files, namely "set{X}_prop.csv", "set{X}_eos.csv", "set{X}_tov.csv", and "set{X}_cs2.csv", where X in [0,1,2,3 and 0_hyp]. set{X}_prop.csv: The file contains the parameters for the RMF model, as well as a few NS properties and nuclear saturation properties. It has the following columns: model name,gs,gv,gr,B,C,xi,lam,rho0,e0,k0,q0,z0,jsym0,lsym0, ksym0,qsym0,zsym0,m_max,r_max,r14,lam14,cs2_max, ec,rhoc,rho_durca. It is to be noted that the parameter B and C are the 10^3*b and 10^3 c (see article for details).  set{X}_eos.csv: For those models in set{X}_prop.csv, it is the NS matter EOS file. It has the following columns: model name, baryon number density, energy density and pressure. The units for baryon number density is fm-3 and MeV/fm3 is for both energy density and pressure. The EOS is for the core only. The crust is not added.  set{X}_tov.csv: For those models in set{X}_prop.csv, it is the TOV solution. It has the following columns: model name, ns radius (km), ns mass (msun),  and dimensionless tidal deformability lambda.  set{X}_cs2.csv: For those models in set{X}_prop.csv, it is the square of the speed of sound over density. It has the following columns: model name, number density fm-3, and square of the speed of sound c2.  -------------------------------------------------------------------------

opencc-by-4.0Nov 2020View details →

Figure 4. Bayesian Inference phylogenetic tree inferred from 1039 in Genetic Relationships of Long-nosed Potoroos Potorous tridactylus (Kerr, 1792) from the Bass Strait Islands, with Notes on the Subspecies Potorous tridactylus benormi Courtney, 1963

Figure 4. Bayesian Inference phylogenetic tree inferred from 1039 bp of concatenated CO1 and ND2 mitochondrial DNA sequence data. Posterior probabilities for major lineages are shown. A similar tree topology was also inferred from Maximum Likelihood.

dryad40/100

Inferring human immunodeficiency virus 1 proviral integration dates with Bayesian Inference

Human immunodeficiency virus 1 (HIV) proviruses archived in the persistent reservoir currently pose the greatest obstacle to HIV cure due to their evasion of combined antiretroviral therapy and ability to reseed HIV infection. Understanding the dynamics of the HIV persistent reservoir is imperative for discovering a durable HIV cure. Here, we explore Bayesian methods using the software BEAST2 to estimate HIV proviral integration dates. We started with within-host longitudinal HIV sequences collected prior to therapy, along with sequences collected from the persistent reservoir during suppressive therapy. We built a BEAST2 model to estimate integration dates of proviral sequences collected during suppressive therapy, implementing a tip date random walker to adjust the sequence tip dates and a latency-specific prior to inform the dates. To validate our method, we implemented it on both simulated and empirical data sets. Consistent with previous studies, we found that proviral integration dates were spread throughout active infection. Path sampling to select an alternative prior for date estimation in place of the latency-specific prior produced unrealistic results in one empirical data set, while on another data set the latency-specific prior was selected as best-fitting. Our Bayesian method outperforms current date estimation techniques with a root mean squared error of 0.89 years on simulated data relative to 1.23–1.89 years with previously developed methods. Bayesian methods offer an adaptable framework for inferring proviral integration dates.

opencc-zeroAug 2023View details →

opencc-by-4.0Aug 2023View details →

Data for AJ paper: Mass ratio of single-line spectroscopic binaries with visual orbits using Bayesian inference and suitable priors

Data and plots for paper "Mass ratio of single-line spectroscopic binaries with visual orbits using Bayesian inference and suitable priors" accepted for publication in The Astronomical Journal.

dryad40/100

Inferring human immunodeficiency virus 1 proviral integration dates with Bayesian Inference

Open the record for dataset details and reuse information.

publicAug 2023View details →

opencc-by-4.0Apr 2020View details →

Data of Bayesian inference of non-linear multiscale model parameters accelerated by a Deep Neural Network

<pre>Data from title = "Bayesian inference of non-linear multiscale model parameters accelerated by a Deep Neural Network", journal = "Computer Methods in Applied Mechanics and Engineering", pages = "112693", year = "2020", issn = "0045-7825", doi = "https://doi.org/10.1016/j.cma.2019.112693", author = "Wu, Ling and Zulueta, Kepa and Major, Zoltan and Arriaga, Aitor and Noels, Ludovic" </pre>

opencc-by-4.0Nov 2020View details →

Data for "Bayesian inference for biophysical neuron models enables stimulus optimization for retinal neuroprosthetics"

Experimental and precomputed data for the paper "Bayesian inference for biophysical neuron models enables stimulus optimization for retinal neuroprosthetics" by Oesterle et al. 2020 (DOI: <a href="https://doi.org/10.7554/eLife.54997">10.7554/eLife.54997</a>). The cone bipolar cell data has been described and published in the paper "Inhibition decorrelates visual feature representations in the inner retina" by Franke et al. 2017 (DOI: <a href="https://doi.org/10.1038/nature21394">10.1038/nature21394</a>).  This data is both a supplement to the Oesterle et al. paper and the code for this paper. The code is available in this <a href="http://github.com/berenslab/CBC_inference">GitHub repository</a>. We recommend downloading the GitHub repository and to follow the instructions there.

opencc-zeroOct 2019View details →

Data from: A Bayesian approach for inferring the impact of a discrete character on rates of continuous-character evolution in the presence of background-rate variation

Understanding how and why rates of character evolution vary across the Tree of Life is central to many evolutionary questions; e.g., does the trophic apparatus (a set of continuous characters) evolve at a higher rate in fish lineages that dwell in reef versus non-reef habitats (a discrete character)? Existing approaches for inferring the relationship between a discrete character and rates of continuous-character evolution rely on comparing a null model (in which rates of continuous-character evolution are constant across lineages) to an alternative model (in which rates of continuous-character evolution depend on the state of the discrete character under consideration). However, these approaches are susceptible to a "straw-man" effect: the influence of the discrete character is inflated because the null model is extremely unrealistic. Here, we describe MuSSCRat, a Bayesian approach for inferring the impact of a discrete trait on rates of continuous-character evolution in the presence of alternative sources of rate variation ("background-rate variation"). We demonstrate by simulation that our method is able to reliably infer the degree of state-dependent rate variation, and show that ignoring background-rate variation leads to biased inferences regarding the degree of state-dependent rate variation in grunts (the fish group Haemulidae).

opencc-zeroApr 2020View details →

Bayesian inference of ancestral host-parasite interactions under a phylogenetic model of host repertoire evolution

Intimate ecological interactions, such as those between parasites and their hosts, may persist over long time spans, coupling the evolutionary histories of the lineages involved. Most methods that reconstruct the coevolutionary history of such interactions make the simplifying assumption that parasites have a single host. Many methods also focus on congruence between host and parasite phylogenies, using cospeciation as the null model. However, there is an increasing body of evidence suggesting that the host ranges of parasites are more complex: that host ranges often include more than one host and evolve via gains and losses of hosts rather than through cospeciation alone. Here, we develop a Bayesian approach for inferring coevolutionary history based on a model accommodating these complexities. Specifically, a parasite is assumed to have a host repertoire, which includes both potential hosts and one or more actual hosts. Over time, potential hosts can be added or lost, and potential hosts can develop into actual hosts or vice versa. Thus, host colonization is modeled as a two-step process that may potentially be influenced by host relatedness. We first explore the statistical behavior of our model by simulating evolution of host-parasite interactions under a range of parameter values. We then use our approach, implemented in the program RevBayes, to infer the coevolutionary history between 34 Nymphalini butterfly species and 25 angiosperm families. Our analysis suggests that host relatedness among angiosperm families influences how easily Nymphalini lineages gain new hosts.

opencc-by-4.0Feb 2017View details →

Figure 2. - Phylogenetic relationships among Dicronocephalus species reconstructed with Bayesian inference using COI sequences. Numbers above branches indicate ML bootstrap values and Bayesian posterior probabilities. Numbers below branches are bootstrap, symmetric resampling, and jacknife support from parsimony searches, respectively. Scale bar represents 10% nucleotide mutation rate.

Figure 2. - Phylogenetic relationships among Dicronocephalus species reconstructed with Bayesian inference using COI sequences. Numbers above branches indicate ML bootstrap values and Bayesian posterior probabilities. Numbers below branches are bootstrap, symmetric resampling, and jacknife support from parsimony searches, respectively. Scale bar represents 10% nucleotide mutation rate.

opencc-by-4.0Feb 2017View details →

Figure 4. - Phylogenetic relationships among Dicronocephalus species reconstructed with Bayesian inference using COI and 16S rRNA sequences. Numbers above branches indicate ML bootstrap values and Bayesian posterior probabilities. Numbers below branches are bootstrap, symmetric resampling, and jacknife support from parsimony searches, respectively. Scale bar represents 10% nucleotide mutation rate.

Figure 4. - Phylogenetic relationships among Dicronocephalus species reconstructed with Bayesian inference using COI and 16S rRNA sequences. Numbers above branches indicate ML bootstrap values and Bayesian posterior probabilities. Numbers below branches are bootstrap, symmetric resampling, and jacknife support from parsimony searches, respectively. Scale bar represents 10% nucleotide mutation rate.

opencc-by-4.0Feb 2017View details →

Figure 3. - Phylogenetic relationships among Dicronocephalus species reconstructed with Bayesian inference using 16S rRNA sequences. Numbers above branches indicate ML bootstrap values and Bayesian posterior probabilities. Numbers below branches are bootstrap, symmetric resampling, and jacknife support from parsimony searches, respectively. Scale bar represents 10% nucleotide mutation rate.

Figure 3. - Phylogenetic relationships among Dicronocephalus species reconstructed with Bayesian inference using 16S rRNA sequences. Numbers above branches indicate ML bootstrap values and Bayesian posterior probabilities. Numbers below branches are bootstrap, symmetric resampling, and jacknife support from parsimony searches, respectively. Scale bar represents 10% nucleotide mutation rate.

opencc-zeroJan 2024View details →

Fast Bayesian inference of phylogenies from multiple continuous characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remains challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

opencc-zeroJun 2022View details →

Dodonaphy - a Software using Hyperbolic Space for Bayesian Phylogenetic Inference

Bayesian inference for phylogenetics is a gold standard for computing distributions of phylogenies. It faces the challenging problem of moving throughout the high-dimensional space of trees. However, hyperbolic space offers a low dimensional representation of tree-like data. In this paper, we embed genomic sequences into hyperbolic space and perform hyperbolic Markov Chain Monte Carlo for Bayesian inference. The posterior probability is computed by decoding a neighbour joining tree from proposed embedding locations. We empirically demonstrate the fidelity of this method on eight data sets. The sampled posterior distribution recovers the splits and branch lengths to a high degree. We investigated the effects of curvature and embedding dimension on the Markov Chain's performance. Finally, we discuss the prospects for adapting this method to navigate tree space with gradients. This software embeds phylogenetic taxa in hyerbolic space to perform Bayesian inference. Version 1.0.0 includes a Markov Chain Monte Carlo (MCMC) that we compared to the state-of-art on eight datasets (previously published elsewhere). This package is implemented in Python3.9 with a simle command line interface provided.

opencc-zeroJun 2024View details →

Data and supplementary information from: Sequential bayesian phylogenetic inference

The ideal approach to Bayesian phylogenetic inference is to estimate all parameters of interest jointly in a single hierarchical model. However, this is often not feasible in practice due to the high computational cost. Instead, phylogenetic pipelines generally consist of sequential analyses, whereby a single point estimate from a given analysis is used as input for the next analysis (e.g., a single multiple sequence alignment is used to estimate a gene tree). In this framework, uncertainty is not propagated from step to step, which can lead to inaccurate or spuriously certain results. Here, we formally develop and test the sequential approach for Bayesian phylogenetic inference, which uses importance sampling to generate observations for the next step of an analysis pipeline from the posterior produced in the previous step. The sequential approach presented here not only accounts for uncertainty between analysis steps, but also allows for greater flexibility in software choice (and hence model availability) and can be more efficient computationally than the traditional joint approach when multiple models are being tested. We show that the sequential approach is identical in practice to the joint approach only if sufficient information in the data is present (a narrow posterior distribution) and/or sufficiently many importance samples are used. Conversely, we show that the common practice of using a single point estimate can be biased, e.g., a single phylogeny estimate to transform an unrooted phylogeny into time-calibrate phylogeny. We demonstrate the theory of sequential Bayesian inference using both a toy example and an empirical case study of insect divergence times estimation using a relaxed clock model from transcriptome data. In the empirical example, we estimate three posterior distributions of branch lengths from the same data (DNA character matrix with a GTR+Gamma+I substitution model, an amino acid data matrix with empirical substitution models, and an amino acid data matrix with the PhyloBayes CAT-GTR model). Finally, we apply three different-node calibration strategies and show that both, the data source and underlying substitution process to estimate branch lengths as well as the node-calibration strategies, impact divergence time estimates. Thus, our new sequential Bayesian phylogenetic inference provides the opportunity to efficiently test different approach for divergence time estimation, including branch lengths estimation from other software.

opencc-by-4.0Feb 2019View details →

Bayesian Inference of Atomistic Structure in Functional Materials

<pre>Configurational 5D DFT dataset of C60/TiO2 material accompanying the manuscript "Bayesian Inference of Atomistic Structure in Functional Materials" by M. Todorovic, M. U. Gutmann, J. Corander and P. Rinke.</pre>