Priors and Posteriors in Bayesian Timing of Divergence Analyses: The Age of Butterflies Revisited.
Ontology highlight
ABSTRACT:
SUBMITTER: Chazot N
PROVIDER: S-EPMC6893297 | biostudies-literature | 2019 Sep
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Molecular divergence time analyses often rely on the age of fossil lineages to calibrate node age estimates. Most divergence time analyses are now performed in a Bayesian framework, where fossil calibrations are incorporated as parametric prior probabilities on node ages. It is widely accepted that an ideal parameterization of such node age prior probabilities should be based on a comprehensive analysis of the fossil record of the clade of interest, but there is currently no generally applicable approach for calculating such informative priors. We provide here a simple and easily implemented method that employs fossil data to estimate the likely amount of missing history prior to the oldest fossil occurrence of a clade, which can be used to fit an informative parametric prior probability distribution on a node age. Specifically, our method uses the extant diversity and the stratigraphic distribution of fossil lineages confidently assigned to a clade to fit a branching model of lineage diversification. Conditioning this on a simple model of fossil preservation, we estimate the likely amount of missing history prior to the oldest fossil occurrence of a clade. The likelihood surface of missing history can then be translated into a parametric prior probability distribution on the age of the clade of interest. We show that the method performs well with simulated fossil distribution data, but that the likelihood surface of missing history can at times be too complex for the distribution-fitting algorithm employed by our software tool. An empirical example of the application of our method is performed to estimate echinoid node ages. A simulation-based sensitivity analysis using the echinoid data set shows that node age prior distributions estimated under poor preservation rates are significantly less informative than those estimated under high preservation rates.
| S-EPMC3673923 | biostudies-literature
Project description:Sparse-data problems are common, and approaches are needed to evaluate the sensitivity of parameter estimates based on sparse data. We propose a Bayesian approach that uses weakly informative priors to quantify sensitivity of parameters to sparse data. The weakly informative prior is based on accumulated evidence regarding the expected magnitude of relationships using relative measures of disease association. We illustrate the use of weakly informative priors with an example of the association of lifetime alcohol consumption and head and neck cancer. When data are sparse and the observed information is weak, a weakly informative prior will shrink parameter estimates toward the prior mean. Additionally, the example shows that when data are not sparse and the observed information is not weak, a weakly informative prior is not influential. Advancements in implementation of Markov Chain Monte Carlo simulation make this sensitivity analysis easily accessible to the practicing epidemiologist.
| S-EPMC3607322 | biostudies-literature
Project description:Bayesian methods are an important set of tools for performing meta-analyses. They avoid some potentially unrealistic assumptions that are required by conventional frequentist methods. More importantly, meta-analysts can incorporate prior information from many sources, including experts' opinions and prior meta-analyses. Nevertheless, Bayesian methods are used less frequently than conventional frequentist methods, primarily because of the need for nontrivial statistical coding, while frequentist approaches can be implemented via many user-friendly software packages. This article aims at providing a practical review of implementations for Bayesian meta-analyses with various prior distributions. We present Bayesian methods for meta-analyses with the focus on odds ratio for binary outcomes. We summarize various commonly used prior distribution choices for the between-studies heterogeneity variance, a critical parameter in meta-analyses. They include the inverse-gamma, uniform, and half-normal distributions, as well as evidence-based informative log-normal priors. Five real-world examples are presented to illustrate their performance. We provide all of the statistical code for future use by practitioners. Under certain circumstances, Bayesian methods can produce markedly different results from those by frequentist methods, including a change in decision on statistical significance. When data information is limited, the choice of priors may have a large impact on meta-analytic results, in which case sensitivity analyses are recommended. Moreover, the algorithm for implementing Bayesian analyses may not converge for extremely sparse data; caution is needed in interpreting respective results. As such, convergence should be routinely examined. When select statistical assumptions that are made by conventional frequentist methods are violated, Bayesian methods provide a reliable alternative to perform a meta-analysis.
| S-EPMC8036799 | biostudies-literature
Project description:Bayesian analyses with the arm-based (AB) network meta-analysis (NMA) model require researchers to specify a prior distribution for the covariance matrix of the treatment-specific event rates in a transformed scale, for example, the treatment-specific log-odds when a logit transformation is used. The commonly used conjugate prior for the covariance matrix, the inverse-Wishart (IW) distribution, has several limitations. For example, although the IW distribution is often described as noninformative or weakly informative, it may in fact provide strong information when some variance components are small (eg, when the standard deviation of study-specific log-odds of a treatment is smaller than 1/2), as is common in NMAs with binary outcomes. In addition, the IW prior generally leads to underestimation of correlations between treatment-specific log-odds, which are critical for borrowing strength across treatment arms to estimate treatment effects efficiently and to reduce potential bias. Alternatively, several separation strategies (ie, separate priors on variances and correlations) can be considered. To study the IW prior's impact on NMA results and compare it with separation strategies, we did simulation studies under different missing-treatment mechanisms. A separation strategy with appropriate priors for the correlation matrix and variances performs better than the IW prior, and should be recommended as the default vague prior in the AB NMA approach. Finally, we reanalyzed three case studies and illustrated the importance, when performing AB-NMA, of sensitivity analyses with different prior specifications on variances.
| S-EPMC7486995 | biostudies-literature
Project description:In this paper, we investigate whether timing in monolingual acquisition interacts with age of onset and input effects in child bilingualism. Six different morpho-syntactic and semantic phenomena acquired early, late or very late are considered, with their timing in L1 acquisition varying between age 3 (subject-verb agreement) and after age 6 (case marking). Data from simultaneous bilingual children (2L1) whose mean age of onset to German was 3 months are compared with data from early second language learners of German (eL2) whose mean age of onset to German was 35 months as well as with data from monolingual children. To explore change over time, children were tested twice at the ages of 4;4 and 5;8 years. The main findings were that 2L1 children had an advantage over their eL2 peers in early acquired phenomena, which disappeared with time, whereas in late acquired phenomena 2L1 and eL2 children did not differ. Moreover, 2L1 children performed like monolingual children in early acquired phenomena but had a disadvantage in the late acquired phenomena with the amount of delay decreasing with time. We conclude that age of onset effects are modulated by effects of timing in monolingual acquisition. Contrary to expectation, input in terms of language dominance, measured as the dominant language used at home, did not affect simultaneous bilingual children's performance in any of the phenomena. We discuss the implications of our findings for the hypothesis that acquisition of late phenomena is determined by input alone and suggest an alternative concept: the learner's internal need for time to master a phenomenon, which is determined by its complexity and cross-linguistic robustness.
| S-EPMC6339886 | biostudies-literature
Project description:Simple heuristics are often regarded as tractable decision strategies because they ignore a great deal of information in the input data. One puzzle is why heuristics can outperform full-information models, such as linear regression, which make full use of the available information. These "less-is-more" effects, in which a relatively simpler model outperforms a more complex model, are prevalent throughout cognitive science, and are frequently argued to demonstrate an inherent advantage of simplifying computation or ignoring information. In contrast, we show at the computational level (where algorithmic restrictions are set aside) that it is never optimal to discard information. Through a formal Bayesian analysis, we prove that popular heuristics, such as tallying and take-the-best, are formally equivalent to Bayesian inference under the limit of infinitely strong priors. Varying the strength of the prior yields a continuum of Bayesian models with the heuristics at one end and ordinary regression at the other. Critically, intermediate models perform better across all our simulations, suggesting that down-weighting information with the appropriate prior is preferable to entirely ignoring it. Rather than because of their simplicity, our analyses suggest heuristics perform well because they implement strong priors that approximate the actual structure of the environment. We end by considering how new heuristics could be derived by infinitely strengthening the priors of other Bayesian models. These formal results have implications for work in psychology, machine learning and economics.
| S-EPMC5886040 | biostudies-literature
Project description:Magnetoencephalography (MEG) records brain activity with excellent temporal and good spatial resolution, while functional magnetic resonance imaging (fMRI) offers good temporal and excellent spatial resolution. The aim of this study is to implement a Bayesian framework to use fMRI data as spatial priors for MEG inverse solutions. We used simulated MEG data with both evoked and induced activity and experimental MEG data from sixteen participants to examine the effectiveness of using fMRI spatial priors in MEG source reconstruction. For simulated MEG data, incorporating the prior information from fMRI increased the spatial resolution of MEG source reconstruction by 3 mm on average. For experimental MEG data, fMRI spatial information reduced the spurious clusters for evoked activity and showed more left-lateralized activation pattern for induced activity. The use of fMRI spatial priors greatly reduced location error for induced source in MEG data. Our results provide empirical evidence that the use of fMRI spatial priors improves the accuracy of MEG source reconstruction. The combined MEG and fMRI approach can provide neuroimaging data with better spatial and temporal resolutions to add another perspective to our understanding of the neurobiology of language. The potential clinical applications include pre-surgical evaluation of language function for epilepsy patients and evaluation of language network for children with language disorders.
| S-EPMC9007727 | biostudies-literature
Project description:MotivationThe flexibility of a Bayesian framework is promising for GWAS, but current approaches can benefit from more informative prior models. We introduce a novel Bayesian approach to GWAS, called Structured and Non-Local Priors (SNLPs) GWAS, that improves over existing methods in two important ways. First, we describe a model that allows for a marker's gene-parent membership and other characteristics to influence its probability of association with an outcome. Second, we describe a non-local alternative model for differential minor allele rates at each marker, in which the null and alternative hypotheses have no common support.ResultsWe employ a non-parametric model that allows for clustering of the genes in tandem with a regression model for marker-level covariates, and demonstrate how incorporating these additional characteristics can improve power. We further demonstrate that our non-local alternative model gives symmetric rates of convergence for the null and alternative hypotheses, whereas commonly used local alternative models have asymptotic rates that favor the alternative hypothesis over the null. We demonstrate the robustness and flexibility of our structured and non-local model for different data generating scenarios and signal-to-noise ratios. We apply our Bayesian GWAS method to single nucleotide polymorphisms data collected from a pool of Alzheimer's disease and cognitively normal patients from the Alzheimer's Database Neuroimaging Initiative.Availability and implementationR code to perform the SNLPs method is available at https://github.com/lockEF/BayesianScreening.
| S-EPMC6956774 | biostudies-literature
Project description:Understanding patterns of species occurrence and the processes underlying these patterns is fundamental to the study of ecology. One of the more commonly used approaches to investigate species occurrence patterns is occupancy modeling, which can account for imperfect detection of a species during surveys. In recent years, there has been a proliferation of Bayesian modeling in ecology, which includes fitting Bayesian occupancy models. The Bayesian framework is appealing to ecologists for many reasons, including the ability to incorporate prior information through the specification of prior distributions on parameters. While ecologists almost exclusively intend to choose priors so that they are "uninformative" or "vague", such priors can easily be unintentionally highly informative. Here we report on how the specification of a "vague" normally distributed (i.e., Gaussian) prior on coefficients in Bayesian occupancy models can unintentionally influence parameter estimation. Using both simulated data and empirical examples, we illustrate how this issue likely compromises inference about species-habitat relationships. While the extent to which these informative priors influence inference depends on the data set, researchers fitting Bayesian occupancy models should conduct sensitivity analyses to ensure intended inference, or employ less commonly used priors that are less informative (e.g., logistic or t prior distributions). We provide suggestions for addressing this issue in occupancy studies, and an online tool for exploring this issue under different contexts.
| S-EPMC5826699 | biostudies-literature
Project description:In archaeological research, changes in material culture and the evolution of styles are taken as major indicators for socio-cultural transformation. They form the basis for typo-chronological classification and the establishment of phases and periods. Central European Bronze Age material culture from burials reveals changes during the Bronze Age and represents a perfect case study for analyzing phenomena of cultural change and the adoption of innovation in the societies of prehistoric Europe. Our study focuses on the large-scale change in material culture which took place in the second millennium BC and the emergence at the same period of new burial rites: the shift from inhumation burials in flat graves to complex mounds and simple cremation burials. Paul Reinecke was the first to divide the European Bronze Age (EBA) into two phases, Bz A1 and A2. The shift from the first to the second phase has so far been ascribed to technical advances. Our study adopted an innovative approach to quantifying this phenomenon. Through regressive reciprocal averaging and Bayesian analysis of radiocarbon-dated grave contexts located in Switzerland and southern Germany, we modelled chronological changes in the material culture and changes in burial rites in these regions in a probabilistic way. We used kernel density models to summarize radiocarbon dates, with the aim of visualizing cultural changes in the third and second millennium BC. In 2015, Stockhammer et al. cast doubt on the chronological sequence of the Reinecke phases of the EBA on the basis of newly collected radiocarbon dates from southern Germany. Our intervention is a direct response to the results of that study. We fully agree with Stockhammer's et al. dating of the start of EBA, but propose a markedly different dating of the EBA/MBA transition. Our modelling of radiocarbon data demonstrates a statistically significant typological sequence of phases Bz A1, Bz A2 and Bz B and disproves their postulated chronological overlap. The linking of the archaeological relative-chronological system with absolute dates is of major importance to understanding the temporal dimension of the EBA phases.
| S-EPMC7769557 | biostudies-literature