Bayesian phylogeography finds its roots.
Directory of Open Access Journals (Sweden)
Philippe Lemey
2009-09-01
Full Text Available As a key factor in endemic and epidemic dynamics, the geographical distribution of viruses has been frequently interpreted in the light of their genetic histories. Unfortunately, inference of historical dispersal or migration patterns of viruses has mainly been restricted to model-free heuristic approaches that provide little insight into the temporal setting of the spatial dynamics. The introduction of probabilistic models of evolution, however, offers unique opportunities to engage in this statistical endeavor. Here we introduce a Bayesian framework for inference, visualization and hypothesis testing of phylogeographic history. By implementing character mapping in a Bayesian software that samples time-scaled phylogenies, we enable the reconstruction of timed viral dispersal patterns while accommodating phylogenetic uncertainty. Standard Markov model inference is extended with a stochastic search variable selection procedure that identifies the parsimonious descriptions of the diffusion process. In addition, we propose priors that can incorporate geographical sampling distributions or characterize alternative hypotheses about the spatial dynamics. To visualize the spatial and temporal information, we summarize inferences using virtual globe software. We describe how Bayesian phylogeography compares with previous parsimony analysis in the investigation of the influenza A H5N1 origin and H5N1 epidemiological linkage among sampling localities. Analysis of rabies in West African dog populations reveals how virus diffusion may enable endemic maintenance through continuous epidemic cycles. From these analyses, we conclude that our phylogeographic framework will make an important asset in molecular epidemiology that can be easily generalized to infer biogeogeography from genetic data for many organisms.
DEFF Research Database (Denmark)
Marske, Katharine Ann; Rahbek, Carsten; Nogues, David Bravo
2013-01-01
Synthesis of ecological and evolutionary concepts and tools has led to improved understanding of how diversification, dispersal, community assembly, long-term coexistence and extinction shape patterns of biological diversity. Phylogeography, with its focus on Quaternary interactions within and...... range reduction and its effects on the strength and temporal flexibility of networks of species interactions. We conclude with an outlook on the data-gathering protocols necessary for this collaborative, interdisciplinary research agenda....
Bayesian Hypothesis Testing for Planet Finding
Braems, I.; Kasdin, N. J.
2003-12-01
One of the most important performance metrics of any space planet finding system is integration time. The time needed to make a positive detection of an extrasolar planet determines the number of systems we can observe for the life of the mission and the stability requirements of the spacecraft and optical control systems. Most astronomical detection approaches rely on fairly simple signal-to-noise calculations and a threshold determined by the ability of the human eye to extract the planet image from the background (usually a signal-to-noise ratio of five). In this paper we present an alternative approach to detection using Bayesian hypothesis testing. This optimal approach provides a quantitative measure of the probability of detection under various conditions and integration times (such as known or unknown background levels) and under different prior assumptions. We also show how the technique allows for a much higher probability of detection for shorter integration times than the previous photometric approaches. We gratefully acknowledge the support of the Jet Propulsion Laboratory of the National Aeronautics and Space Administration for this work and Institut National de Recherche en Informatique et Automatique (INRIA) for its support of Ms. Braems.
Bayesian methods for finding sparse representations
Wipf, David Paul
2006-01-01
Finding the sparsest or minimum L0-norm representation of a signal given a (possibly) overcomplete dictionary of basis vectors is an important problem in many application domains, including neuroelectromagnetic source localization, compressed sensing, sparse component analysis, feature selection, image restoration/compression, and neural coding. Unfortunately, the required optimization is typically NP-hard, and so approximate procedures that succeed with high probability are sought. Nearly al...
Finding Clocks in Genes: A Bayesian Approach to Estimate Periodicity
Directory of Open Access Journals (Sweden)
Yan Ren
2016-01-01
Full Text Available Identification of rhythmic gene expression from metabolic cycles to circadian rhythms is crucial for understanding the gene regulatory networks and functions of these biological processes. Recently, two algorithms, JTK_CYCLE and ARSER, have been developed to estimate periodicity of rhythmic gene expression. JTK_CYCLE performs well for long or less noisy time series, while ARSER performs well for detecting a single rhythmic category. However, observing gene expression at high temporal resolution is not always feasible, and many scientists are interested in exploring both ultradian and circadian rhythmic categories simultaneously. In this paper, a new algorithm, named autoregressive Bayesian spectral regression (ABSR, is proposed. It estimates the period of time-course experimental data and classifies gene expression profiles into multiple rhythmic categories simultaneously. Through the simulation studies, it is shown that ABSR substantially improves the accuracy of periodicity estimation and clustering of rhythmic categories as compared to JTK_CYCLE and ARSER for the data with low temporal resolution. Moreover, ABSR is insensitive to rhythmic patterns. This new scheme is applied to existing time-course mouse liver data to estimate period of rhythms and classify the genes into ultradian, circadian, and arrhythmic categories. It is observed that 49.2% of the circadian profiles detected by JTK_CYCLE with 1-hour resolution are also detected by ABSR with only 4-hour resolution.
Finding the Most Distant Quasars Using Bayesian Selection Methods
Mortlock, Daniel
2014-01-01
Quasars, the brightly glowing disks of material that can form around the super-massive black holes at the centres of large galaxies, are amongst the most luminous astronomical objects known and so can be seen at great distances. The most distant known quasars are seen as they were when the Universe was less than a billion years old (i.e., $\\sim\\!7%$ of its current age). Such distant quasars are, however, very rare, and so are difficult to distinguish from the billions of other comparably-bright sources in the night sky. In searching for the most distant quasars in a recent astronomical sky survey (the UKIRT Infrared Deep Sky Survey, UKIDSS), there were $\\sim\\!10^3$ apparently plausible candidates for each expected quasar, far too many to reobserve with other telescopes. The solution to this problem was to apply Bayesian model comparison, making models of the quasar population and the dominant contaminating population (Galactic stars) to utilise the information content in the survey measurements. The result wa...
A Bayesian Method For Finding Galaxies That Cause Quasar Absorption Lines
Shoemaker, Emileigh Suzanne; Laubner, David Andrew; Scott, Jennifer E.
2016-01-01
We present a study of candidate absorber-galaxy pairs for 39 low redshift quasar sightlines (0.06 Digital Sky Survey (SDSS). We downloaded the COS linelists for these quasar spectra from MAST and queried the SDSS DR12 database for photometric data on all galaxies within 1 Mpc of each of these quasar lines of sight. We calculated photometric redshifts for all the SDSS galaxies using the Bayesian Photometric Redshift code. We used all these absorber and galaxy data as input into an absorber-galaxy matching code which also employs a Bayesian scheme, along with known statistics of the intergalactic medium and circumgalactic media of galaxies, for finding the most probable galaxy match for each absorber. We compare our candidate absorber-galaxy matches to existing studies in the literature and explore trends in the absorber and galaxy properties among the matched and non-matched populations. This method of matching absorbers and galaxies can be used to find targets for follow up spectroscopic studies.
Breathnach, Michelle; Moore, Elizabeth
2013-06-10
The Bayesian Approach allows forensic scientists to evaluate the significance of scientific evidence in light of two conflicting hypothesis. This aids the investigator to calculate a numerical value of the probability that the scientific findings support one hypothesis over conflicting opinions. In the case where oral intercourse is alleged, α-amylase, an indicator of saliva, is detected on penile swabs. The value of this finding is unknown as it may indicate the presence of saliva resulting from oral intercourse however it may also represent the presence of saliva due to innocent means such as background levels of salivary-α-amylase in the male population due to secondary transfer. Therefore, it is difficult to attach significance to this finding without background information and knowledge. A population study of the background levels of salivary-α-amylase was performed by analysing items of underwear worn under normal circumstances by 69 male volunteers. The Phadebas press test was used to screen the garments for amylase-containing stains and the positive areas were subjected to further confirmation of saliva by the RSID-Saliva kit. 44% of underwear screened had stains containing amylase. This study determined the background level of salivary-α-amylase and DNA on the inside front of male underwear which has potential implications on the interpretation of evidence in alleged oral intercourse. PMID:23683908
Herschtal, A.; Foroudi, F.; Greer, P. B.; Eade, T. N.; Hindson, B. R.; Kron, T.
2012-05-01
Early approaches to characterizing errors in target displacement during a fractionated course of radiotherapy assumed that the underlying fraction-to-fraction variability in target displacement, known as the ‘treatment error’ or ‘random error’, could be regarded as constant across patients. More recent approaches have modelled target displacement allowing for differences in random error between patients. However, until recently it has not been feasible to compare the goodness of fit of alternate models of random error rigorously. This is because the large volumes of real patient data necessary to distinguish between alternative models have only very recently become available. This work uses real-world displacement data collected from 365 patients undergoing radical radiotherapy for prostate cancer to compare five candidate models for target displacement. The simplest model assumes constant random errors across patients, while other models allow for random errors that vary according to one of several candidate distributions. Bayesian statistics and Markov Chain Monte Carlo simulation of the model parameters are used to compare model goodness of fit. We conclude that modelling the random error as inverse gamma distributed provides a clearly superior fit over all alternatives considered. This finding can facilitate more accurate margin recipes and correction strategies.
Gasparini, Mauro; Eisele, J
2003-01-01
Compositional random vectors are fundamental tools in the Bayesian analysis of categorical data. Many of the issues that are discussed with reference to the statistical analysis of compositional data have a natural counterpart in the construction of a Bayesian statistical model for categorical data. This note builds on the idea of cross-fertilization of the two areas recommended by Aitchison (1986) in his seminal book on compositional data. Particular emphasis is put on the pro...
Crandell, Jamie L.; Voils, Corrine I.; Chang, YunKyung; Sandelowski, Margarete
2011-01-01
The possible utility of Bayesian methods for the synthesis of qualitative and quantitative research has been repeatedly suggested but insufficiently investigated. In this project, we developed and used a Bayesian method for synthesis, with the goal of identifying factors that influence adherence to HIV medication regimens. We investigated the effect of 10 factors on adherence. Recognizing that not all factors were examined in all studies, we considered standard methods for dealing with missin...
Limited, episodic diversification and contrasting phylogeography in a New Zealand cicada radiation
DEFF Research Database (Denmark)
Marshall, David; Hill, Kathy; Marske, Katharine;
2012-01-01
four extant species (Amphipsalta - 3 spp. + Notopsalta - 1 sp.) and has been little studied. We examined mitochondrial and nuclear-gene phylogenies and phylogeography, Bayesian relaxed-clock divergence timing (incorporating literature-based uncertainty of molecular clock estimates) and ecological niche...... models of the species from the smaller radiation....
Schöniger, Anneli; Illman, Walter A.; Wöhling, Thomas; Nowak, Wolfgang
2015-12-01
Groundwater modelers face the challenge of how to assign representative parameter values to the studied aquifer. Several approaches are available to parameterize spatial heterogeneity in aquifer parameters. They differ in their conceptualization and complexity, ranging from homogeneous models to heterogeneous random fields. While it is common practice to invest more effort into data collection for models with a finer resolution of heterogeneities, there is a lack of advice which amount of data is required to justify a certain level of model complexity. In this study, we propose to use concepts related to Bayesian model selection to identify this balance. We demonstrate our approach on the characterization of a heterogeneous aquifer via hydraulic tomography in a sandbox experiment (Illman et al., 2010). We consider four increasingly complex parameterizations of hydraulic conductivity: (1) Effective homogeneous medium, (2) geology-based zonation, (3) interpolation by pilot points, and (4) geostatistical random fields. First, we investigate the shift in justified complexity with increasing amount of available data by constructing a model confusion matrix. This matrix indicates the maximum level of complexity that can be justified given a specific experimental setup. Second, we determine which parameterization is most adequate given the observed drawdown data. Third, we test how the different parameterizations perform in a validation setup. The results of our test case indicate that aquifer characterization via hydraulic tomography does not necessarily require (or justify) a geostatistical description. Instead, a zonation-based model might be a more robust choice, but only if the zonation is geologically adequate.
DEFF Research Database (Denmark)
Marske, Katharine Ann; Leschen, Richard; Buckley, Thomas
2011-01-01
stochastic search variable selection incorporated in BEAST to identify historical dispersal patterns via ancestral state reconstruction. Ecological niche models (ENMs) were incorporated to reconstruct the potential geographic distribution of each species during the Last Glacial Maximum (LGM). Coalescent......Mitochondrial DNA (cox1) sequence data and recently developed coalescent phylogeography models were used to construct geo-spatial histories for the New Zealand fungus beetles Epistranus lawsoni and Pristoderus bakewelli (Zopheridae). These methods utilize continuous-time Markov chains and Bayesian...
A Bayesian Dose-finding Design for Oncology Clinical Trials of Combinational Biological Agents.
Cai, Chunyan; Yuan, Ying; Ji, Yuan
2014-01-01
Treating patients with novel biological agents is becoming a leading trend in oncology. Unlike cytotoxic agents, for which efficacy and toxicity monotonically increase with dose, biological agents may exhibit non-monotonic patterns in their dose-response relationships. Using a trial with two biological agents as an example, we propose a dose-finding design to identify the biologically optimal dose combination (BODC), which is defined as the dose combination of the two agents with the highest efficacy and tolerable toxicity. A change-point model is used to reflect the fact that the dose-toxicity surface of the combinational agents may plateau at higher dose levels, and a flexible logistic model is proposed to accommodate the possible non-monotonic pattern for the dose-efficacy relationship. During the trial, we continuously update the posterior estimates of toxicity and efficacy and assign patients to the most appropriate dose combination. We propose a novel dose-finding algorithm to encourage sufficient exploration of untried dose combinations in the two-dimensional space. Extensive simulation studies show that the proposed design has desirable operating characteristics in identifying the BODC under various patterns of dose-toxicity and dose-efficacy relationships. PMID:24511160
Human phylogeography and diversity.
Harcourt, Alexander H
2016-07-19
Homo sapiens phylogeography begins with the species' origin nearly 200 kya in Africa. First signs of the species outside Africa (in Arabia) are from 125 kya. Earliest dates elsewhere are now 100 kya in China, 45 kya in Australia and southern Europe (maybe even 60 kya in Australia), 32 kya in northeast Siberia, and maybe 20 kya in the Americas. Humans reached arctic regions and oceanic islands last-arctic North America about 5 kya, mid- and eastern Pacific islands about 2-1 kya, and New Zealand about 700 y ago. Initial routes along coasts seem the most likely given abundant and easily harvested shellfish there as indicated by huge ancient oyster shell middens on all continents. Nevertheless, the effect of geographic barriers-mountains and oceans-is clear. The phylogeographic pattern of diasporas from several single origins-northeast Africa to Eurasia, southeast Eurasia to Australia, and northeast Siberia to the Americas-allows the equivalent of a repeat experiment on the relation between geography and phylogenetic and cultural diversity. On all continents, cultural diversity is high in productive low latitudes, presumably because such regions can support populations of sustainable size in a small area, therefore allowing a high density of cultures. Of course, other factors operate. South America has an unusually low density of cultures in its tropical latitudes. A likely factor is the phylogeographic movement of peoples from the Old World bringing novel and hence, lethal diseases to the New World, a foretaste, perhaps, of present day global transport of tropical diseases. PMID:27432967
Graph hierarchies for phylogeography.
Cybis, Gabriela B; Sinsheimer, Janet S; Lemey, Philippe; Suchard, Marc A
2013-03-19
Bayesian phylogeographic methods simultaneously integrate geographical and evolutionary modelling, and have demonstrated value in assessing spatial spread patterns of measurably evolving organisms. We improve on existing phylogeographic methods by combining information from multiple phylogeographic datasets in a hierarchical setting. Consider N exchangeable datasets or strata consisting of viral sequences and locations, each evolving along its own phylogenetic tree and according to a conditionally independent geographical process. At the hierarchical level, a random graph summarizes the overall dispersion process by informing which migration rates between sampling locations are likely to be relevant in the strata. This approach provides an efficient and improved framework for analysing inherently hierarchical datasets. We first examine the evolutionary history of multiple serotypes of dengue virus in the Americas to showcase our method. Additionally, we explore an application to intrahost HIV evolution across multiple patients. PMID:23382428
Nolan, L A; Kaban, Ata; Raychaudhuri, S
2006-01-01
We present the results of a novel application of Bayesian modelling techniques, which, although purely data driven, have a physically interpretable result, and will be useful as an efficient data mining tool. We base our studies on the UV-to-optical spectra (observed and synthetic) of early-type galaxies. A probabilistic latent variable architecture is formulated, and a rigorous Bayesian methodology is employed for solving the inverse modelling problem from the available data. A powerful aspect of our formalism is that it allows us to recover a limited fraction of missing data due to incomplete spectral coverage, as well as to handle observational errors in a principled way. We apply this method to a sample of 21 well-studied early-type spectra, with known star-formation histories. We find that our data-driven Bayesian modelling allows us to identify those early-types which contain a significant stellar population <~ 1 Gyr old. This method would therefore be a very useful tool for automatically discovering...
Comparative phylogeography of the ocean planet.
Bowen, Brian W; Gaither, Michelle R; DiBattista, Joseph D; Iacchei, Matthew; Andrews, Kimberly R; Grant, W Stewart; Toonen, Robert J; Briggs, John C
2016-07-19
Understanding how geography, oceanography, and climate have ultimately shaped marine biodiversity requires aligning the distributions of genetic diversity across multiple taxa. Here, we examine phylogeographic partitions in the sea against a backdrop of biogeographic provinces defined by taxonomy, endemism, and species composition. The taxonomic identities used to define biogeographic provinces are routinely accompanied by diagnostic genetic differences between sister species, indicating interspecific concordance between biogeography and phylogeography. In cases where individual species are distributed across two or more biogeographic provinces, shifts in genotype frequencies often align with biogeographic boundaries, providing intraspecific concordance between biogeography and phylogeography. Here, we provide examples of comparative phylogeography from (i) tropical seas that host the highest marine biodiversity, (ii) temperate seas with high productivity but volatile coastlines, (iii) migratory marine fauna, and (iv) plankton that are the most abundant eukaryotes on earth. Tropical and temperate zones both show impacts of glacial cycles, the former primarily through changing sea levels, and the latter through coastal habitat disruption. The general concordance between biogeography and phylogeography indicates that the population-level genetic divergences observed between provinces are a starting point for macroevolutionary divergences between species. However, isolation between provinces does not account for all marine biodiversity; the remainder arises through alternative pathways, such as ecological speciation and parapatric (semiisolated) divergences within provinces and biodiversity hotspots. PMID:27432963
Directory of Open Access Journals (Sweden)
Shao-Yu Chen
Full Text Available BACKGROUND: The Siberian salamander (Ranodon sibiricus, distributed in geographically isolated areas of Central Asia, is an ideal alpine species for studies of conservation and phylogeography. However, there are few data regarding the genetic diversity in R. sibiricus populations. METHODOLOGY/PRINCIPAL FINDINGS: We used two genetic markers (mtDNA and microsatellites to survey all six populations of R. sibiricus in China. Both of the markers revealed extreme genetic uniformity among these populations. There were only three haplotypes in the mtDNA, and the overall nucleotide diversity in the mtDNA was 0.00064, ranging from 0.00000 to 0.00091 for the six populations. Although we recovered 70 sequences containing microsatellite repeats, there were only two loci that displayed polymorphism. We used the approximate Bayesian computation (ABC method to study the demographic history of the populations. This analysis suggested that the extant populations diverged from the ancestral population approximately 120 years ago and that the historical population size was much larger than the present population size; i.e., R. sibiricus has experienced dramatic population declines. CONCLUSION/SIGNIFICANCE: Our findings suggest that the genetic diversity in the R. sibiricus populations is the lowest among all investigated amphibians. We conclude that the isolation of R. sibiricus populations occurred recently and was a result of recent human activity and/or climatic changes. The Pleistocene glaciation oscillations may have facilitated intraspecies genetic homogeneity rather than enhanced divergence. A low genomic evolutionary rate and elevated inbreeding frequency may have also contributed to the low genetic variation observed in this species. Our findings indicate the urgency of implementing a protection plan for this endangered species.
2011-01-01
Predicting the reliability of software systems based on a component-based approach is inherently difficult, in particular due to failure dependencies between software components. One possible way to assess and include dependency aspects in software reliability models is to find upper bounds for probabilities that software components fail simultaneously and then include these into the reliability models. In earlier research, it has been shown that including partial dependency information may g...
Bayesian optimization for materials design
Frazier, Peter I.; Wang, Jialei
2015-01-01
We introduce Bayesian optimization, a technique developed for optimizing time-consuming engineering simulations and for fitting machine learning models on large datasets. Bayesian optimization guides the choice of experiments during materials design and discovery to find good material designs in as few experiments as possible. We focus on the case when materials designs are parameterized by a low-dimensional vector. Bayesian optimization is built on a statistical technique called Gaussian pro...
Phylogeography of genus Squalius in Albania
Radek Šanda; Miroslav Švátora
2015-01-01
This study is focused on the issue of the Squalius genus phylogeography in Albania in the Balkan region. Phylogenetic analyses of sequence variation at mitochondrial DNA (cytochrome b gene) were used to examine these issues for the freshwater fish of the genus Squalius from various river systems in the Adriatic Sea region. There were identified three genetic lineages of unclear taxonomic position, where the genetic variation between identified clades range from 1.6 to 2.1 %. The first lineage...
... Brain George Hightower searches for genetic mutations that affect HIV's ability to infect the brain. Read Issue All Issues Explore Findings by Topic Cell Biology Cellular Structures, Functions, Processes, Imaging, Stress Response Chemistry and Biochemistry Enzymes, Molecular Probes, Metabolic ...
Lesaffre, Emmanuel
2012-01-01
The growth of biostatistics has been phenomenal in recent years and has been marked by considerable technical innovation in both methodology and computational practicality. One area that has experienced significant growth is Bayesian methods. The growing use of Bayesian methodology has taken place partly due to an increasing number of practitioners valuing the Bayesian paradigm as matching that of scientific discovery. In addition, computational advances have allowed for more complex models to be fitted routinely to realistic data sets. Through examples, exercises and a combination of introd
Connectivity in the deep: Phylogeography of the velvet belly lanternshark
Gubili, Chrysoula; Macleod, Kirsty; Perry, William; Hanel, Pia; Batzakas, Ioannis; Farrell, Edward D.; Lynghammar, Arve; Mancusi, Cecilia; Mariani, Stefano; Menezes, Gui M.; Neat, Francis; Scarcella, Giuseppe; Griffiths, Andrew M.
2016-09-01
The velvet belly lanternshark, Etmopterus spinax, is a deep-sea bioluminescent squaloid shark, found predominantly in the Northeast Atlantic and Mediterranean Sea. It has been exposed to relatively high levels of mortality associated with by-catch in some regions. Its late maturity and low fecundity potentially renders it vulnerable to over-exploitation, although little remains known about processes of connectivity between key habitats/regions. This study utilised DNA sequencing of partial regions of the mitochondrial control region and nuclear ribosomal internal transcribed spacer 2 to investigate population structure and phylogeography of this species across the Northeast Atlantic and Mediterranean Basin. Despite the inclusion of samples from the range edges or remote locations, no evidence of significant population structure was detected. An important exception was identified using the control region sequence, with much greater (and statistically significant) levels of genetic differentiation between the Mediterranean and Atlantic. This suggests that the Strait of Gibraltar may represent an important bathymetric barrier, separating regions with very low levels of female dispersal. Bayesian estimation of divergence time also places the separation between the Mediterranean and Atlantic lineages within the last 100,000 years, presumably connected with perturbations during the last Glacial Period. These results demonstrate population subdivision at a much smaller geographic distance than has generally been identified in previous work on deep-sea sharks. This highlights a very significant role for shallow bathymetry in promoting genetic differentiation in deepwater taxa. It acts as an important exception to a general paradigm of marine species being connected by high levels of gene-flow, representing single stocks over large scales. It may also have significant implications for the fisheries management of this species.
Draper, D.
2001-01-01
© 2012 Springer Science+Business Media, LLC. All rights reserved. Article Outline: Glossary Definition of the Subject and Introduction The Bayesian Statistical Paradigm Three Examples Comparison with the Frequentist Statistical Paradigm Future Directions Bibliography
Comparative Phylogeography in Rainforest Trees from Lower Guinea, Africa
Myriam Heuertz; Jérôme Duminil; Gilles Dauby; Vincent Savolainen; Olivier J Hardy
2014-01-01
Comparative phylogeography is an effective approach to assess the evolutionary history of biological communities. We used comparative phylogeography in fourteen tree taxa from Lower Guinea (Atlantic Equatorial Africa) to test for congruence with two simple evolutionary scenarios based on physio-climatic features 1) the W-E environmental gradient and 2) the N-S seasonal inversion, which determine climatic and seasonality differences in the region. We sequenced the trnC-ycf6 plastid DNA region ...
Epidemiological history and phylogeography of West Nile virus lineage 2.
Ciccozzi, Massimo; Peletto, Simone; Cella, Eleonora; Giovanetti, Marta; Lai, Alessia; Gabanelli, Elena; Acutis, Pier Luigi; Modesto, Paola; Rezza, Giovanni; Platonov, Alexander E; Lo Presti, Alessandra; Zehender, Gianguglielmo
2013-07-01
West Nile virus (WNV) was first isolated in Uganda. In Europe WNV was sporadically detected until 1996, since then the virus has been regularly isolated from birds and mosquitoes and caused several outbreaks in horses and humans. Phylogenetic analysis showed two main different WNV lineages. The lineage 1 is widespread and segregates into different subclades (1a-c). WNV-1a includes numerous strains from Africa, America, and Eurasia. The spatio-temporal history of WNV-1a in Europe was recently described, identifying two main routes of dispersion, one in Eastern and the second in Western Europe. The West Nile lineage 2 (WNV-2) is mainly present in sub-Saharan Africa but has been recently emerged in Eastern and Western European countries. In this study we reconstruct the phylogeny of WNV-2 on a spatio-temporal scale in order to estimate the time of origin and patterns of geographical dispersal of the different isolates, particularly in Europe. Phylogeography findings obtained from E and NS5 gene analyses suggest that there were at least two separate introductions of WNV-2 from the African continent dated back approximately to the year 1999 (Central Europe) and 2000 (Russia), respectively. The epidemiological implications and clinical consequences of lineage 1 and 2 cocirculation deserve further investigations. PMID:23542457
Kirstein, Roland
2005-01-01
This paper presents a modification of the inspection game: The ?Bayesian Monitoring? model rests on the assumption that judges are interested in enforcing compliant behavior and making correct decisions. They may base their judgements on an informative but imperfect signal which can be generated costlessly. In the original inspection game, monitoring is costly and generates a perfectly informative signal. While the inspection game has only one mixed strategy equilibrium, three Perfect Bayesia...
Comparative phylogeography: concepts, methods and general patterns in neotropical birds
International Nuclear Information System (INIS)
Understanding the patterns and processes involved in intraspecific lineages diversification in time and space is the aim of phylogeography. The comparison of those phylogeographic patterns among co-distributed species shows insights of a community history. Here I review the concepts and methodologies of comparative phylogeography, an active research field that has heterogeneous analytical methods. In order to present a framework for phylogeography in the neotropics, I comment the general phylogeographic patterns of the birds from this region. this review is based on more than 100 studies conducted during the last 25 years and indicate that despite different co-distributed species seem to share some points in their phylogeographic pattern they have idiosyncratic aspects, indicating an unique history for each one.
Bessiere, Pierre; Ahuactzin, Juan Manuel; Mekhnacha, Kamel
2013-01-01
Probability as an Alternative to Boolean LogicWhile logic is the mathematical foundation of rational reasoning and the fundamental principle of computing, it is restricted to problems where information is both complete and certain. However, many real-world problems, from financial investments to email filtering, are incomplete or uncertain in nature. Probability theory and Bayesian computing together provide an alternative framework to deal with incomplete and uncertain data. Decision-Making Tools and Methods for Incomplete and Uncertain DataEmphasizing probability as an alternative to Boolean
Liley, James; Wallace, Chris
2015-01-01
Genome-wide association studies (GWAS) have been successful in identifying single nucleotide polymorphisms (SNPs) associated with many traits and diseases. However, at existing sample sizes, these variants explain only part of the estimated heritability. Leverage of GWAS results from related phenotypes may improve detection without the need for larger datasets. The Bayesian conditional false discovery rate (cFDR) constitutes an upper bound on the expected false discovery rate (FDR) across a set of SNPs whose p values for two diseases are both less than two disease-specific thresholds. Calculation of the cFDR requires only summary statistics and have several advantages over traditional GWAS analysis. However, existing methods require distinct control samples between studies. Here, we extend the technique to allow for some or all controls to be shared, increasing applicability. Several different SNP sets can be defined with the same cFDR value, and we show that the expected FDR across the union of these sets may exceed expected FDR in any single set. We describe a procedure to establish an upper bound for the expected FDR among the union of such sets of SNPs. We apply our technique to pairwise analysis of p values from ten autoimmune diseases with variable sharing of controls, enabling discovery of 59 SNP-disease associations which do not reach GWAS significance after genomic control in individual datasets. Most of the SNPs we highlight have previously been confirmed using replication studies or larger GWAS, a useful validation of our technique; we report eight SNP-disease associations across five diseases not previously declared. Our technique extends and strengthens the previous algorithm, and establishes robust limits on the expected FDR. This approach can improve SNP detection in GWAS, and give insight into shared aetiology between phenotypically related conditions. PMID:25658688
Local phylogeography colonisation of the British Isles by small mammals
Czech Academy of Sciences Publication Activity Database
Searle, J.; Barnett, R.; Dobney, K. M.; Herman, J.S.; Hoelzel, A. R.; Jones, E. P.; Kotlík, Petr; McDevitt, A.D.; Martínková, Natália; Rambau, R.V.
Mendoza : Biodiversity Research Group, 2009, s. 9-10. [International Mammalogical Congress /10./. Mendoza (AR), 09.08.2009-14.08.2009] R&D Projects: GA AV ČR IAA600450701 Institutional research plan: CEZ:AV0Z50450515; CEZ:AV0Z60930519 Keywords : phylogeography * colonisation * small mammals Subject RIV: EH - Ecology, Behaviour
Phylogeography of Pipistrellus pipistrellus species group: radiation and reticulation
Czech Academy of Sciences Publication Activity Database
Hulva, P.; Fornůsková, Alena; Chudárková, A.; Allegrini, B.; Evin, A.; Benda, P.; Horáček, I.; Bartonička, T.; Bryja, Josef
Kostelec nad Černými lesy: Lesnická práce, 2010 - (Horáček, I.; Benda, P.). s. 177-178 ISBN 978-80-87154-46-5. [International Bat Research Conference /15./. 23.08.2010-27.08.2010, Praha] Institutional research plan: CEZ:AV0Z60930519 Keywords : phylogeography * bats Subject RIV: EG - Zoology
Phylogeography of genus Squalius in Albania
Directory of Open Access Journals (Sweden)
Radek Šanda
2015-11-01
Full Text Available This study is focused on the issue of the Squalius genus phylogeography in Albania in the Balkan region. Phylogenetic analyses of sequence variation at mitochondrial DNA (cytochrome b gene were used to examine these issues for the freshwater fish of the genus Squalius from various river systems in the Adriatic Sea region. There were identified three genetic lineages of unclear taxonomic position, where the genetic variation between identified clades range from 1.6 to 2.1 %. The first lineage is distributed in the Neretva River drainage in Bosnia and Hercegovina, i.e. north of Albania and in the whole remaining Periadriatic regions, whereas the second lineage is especially spread in the northern part of Albania and the third lineage occurs especially in the zone of the European ancient lake systems on the Balkan Peninsula (lakes Ohrid and Prespa, from where expands to the southern part of Albania. Both lineages sympatrically coincide in the hydrological river-lake system of Ohrid-Drin-Skadar. The phylogenetic and taxonomic position of the Squalius genus in the region seems to be interesting topic for subsequent and more detailed study.
Evolutionary lessons from California plant phylogeography.
Sork, Victoria L; Gugger, Paul F; Chen, Jin-Ming; Werth, Silke
2016-07-19
Phylogeography documents the spatial distribution of genetic lineages that result from demographic processes, such as population expansion, population contraction, and gene movement, shaped by climate fluctuations and the physical landscape. Because most phylogeographic studies have used neutral markers, the role of selection may have been undervalued. In this paper, we contend that plants provide a useful evolutionary lesson about the impact of selection on spatial patterns of neutral genetic variation, when the environment affects which individuals can colonize new sites, and on adaptive genetic variation, when environmental heterogeneity creates divergence at specific loci underlying local adaptation. Specifically, we discuss five characteristics found in plants that intensify the impact of selection: sessile growth form, high reproductive output, leptokurtic dispersal, isolation by environment, and the potential to evolve longevity. Collectively, these traits exacerbate the impact of environment on movement between populations and local selection pressures-both of which influence phylogeographic structure. We illustrate how these unique traits shape these processes with case studies of the California endemic oak, Quercus lobata, and the western North American lichen, Ramalina menziesii Obviously, the lessons we learn from plant traits are not unique to plants, but they highlight the need for future animal, plant, and microbe studies to incorporate its impact. Modern tools that generate genome-wide sequence data are now allowing us to decipher how evolutionary processes affect the spatial distribution of different kinds of genes and also to better model future spatial distribution of species in response to climate change. PMID:27432984
Chloroplast DNA Phylogeography of Holy Basil (Ocimum tenuiflorum in Indian Subcontinent
Directory of Open Access Journals (Sweden)
Felix Bast
2014-01-01
Full Text Available Ocimum tenuiflorum L., holy basil “Tulsi”, is an important medicinal plant that is being grown and traditionally revered throughout Indian Subcontinent for thousands of years; however, DNA sequence-based genetic diversity of this aromatic herb is not yet known. In this report, we present our studies on the phylogeography of this species using trnL-trnF intergenic spacer of plastid genome as the DNA barcode for isolates from Indian subcontinent. Our pairwise distance analyses indicated that genetic heterogeneity of isolates remained quite low, with overall mean nucleotide p-distance of 5×10-4. However, our sensitive phylogenetic analysis using maximum likelihood framework was able to reveal subtle intraspecific molecular evolution of this species within the subcontinent. All isolates except that from North-Central India formed a distinct phylogenetic clade, notwithstanding low bootstrap support and collapse of the clade in Bayesian Inference. North-Central isolates occupied more basal position compared to other isolates, which is suggestive of its evolutionarily primitive status. Indian isolates formed a monophyletic and well-supported clade within O. tenuiflorum clade, which indicates a distinct haplotype. Given the vast geographical area of more than 3 million km2 encompassing many exclusive biogeographical and ecological zones, relatively low rate of evolution of this herb at this locus in India is particularly interesting.
Chloroplast DNA phylogeography of holy basil (Ocimum tenuiflorum) in Indian subcontinent.
Bast, Felix; Rani, Pooja; Meena, Devendra
2014-01-01
Ocimum tenuiflorum L., holy basil "Tulsi", is an important medicinal plant that is being grown and traditionally revered throughout Indian Subcontinent for thousands of years; however, DNA sequence-based genetic diversity of this aromatic herb is not yet known. In this report, we present our studies on the phylogeography of this species using trnL-trnF intergenic spacer of plastid genome as the DNA barcode for isolates from Indian subcontinent. Our pairwise distance analyses indicated that genetic heterogeneity of isolates remained quite low, with overall mean nucleotide p-distance of 5 × 10(-4). However, our sensitive phylogenetic analysis using maximum likelihood framework was able to reveal subtle intraspecific molecular evolution of this species within the subcontinent. All isolates except that from North-Central India formed a distinct phylogenetic clade, notwithstanding low bootstrap support and collapse of the clade in Bayesian Inference. North-Central isolates occupied more basal position compared to other isolates, which is suggestive of its evolutionarily primitive status. Indian isolates formed a monophyletic and well-supported clade within O. tenuiflorum clade, which indicates a distinct haplotype. Given the vast geographical area of more than 3 million km(2) encompassing many exclusive biogeographical and ecological zones, relatively low rate of evolution of this herb at this locus in India is particularly interesting. PMID:24523650
Gray, Rebecca R; Salemi, Marco
2012-12-01
The rate of new emerging infectious diseases entering the human population has increased over the past century, with pathogens originating from animals or from products of animal origin accounting for the vast majority. Primary risk factors for the emergence and spread of emerging zoonoses include expansion and intensification of animal agriculture and long-distance live animal transport, live animal markets, bushmeat consumption and habitat destruction. Developing effective control strategies is contingent upon the ability to test causative hypotheses of disease transmission within a statistical framework. Broadly speaking, molecular phylogeography offers a framework in which specific hypotheses regarding pathogen gene flow and dispersal within an ecological context can be compared. A number of different methods has been developed for this application. Here, our intent is firstly to discuss the application of a wide variety of statistically based methods (including Bayesian reconstruction, network parsimony analysis and regression) to specific viruses (influenza, salmon anaemia virus, foot and mouth disease and Rift Valley Fever) that have been associated with animal farming/movements; and secondly to place them in the larger framework of the threat of potential zoonotic events as well as the economic and biosecurity implications of pathogen outbreaks among our animal food sources. PMID:22931895
Mitochondrial DNA under siege in avian phylogeography.
Zink, Robert M; Barrowclough, George F
2008-05-01
Mitochondrial DNA (mtDNA) has been the workhorse of research in phylogeography for almost two decades. However, concerns with basing evolutionary interpretations on mtDNA results alone have been voiced since the inception of such studies. Recently, some authors have suggested that the potential problems with mtDNA are so great that inferences about population structure and species limits are unwarranted unless corroborated by other evidence, usually in the form of nuclear gene data. Here we review the relative merits of mitochondrial and nuclear phylogeographical studies, using birds as an exemplar class of organisms. A review of population demographic and genetic theory indicates that mitochondrial and nuclear phylogeographical results ought to concur for both geographically unstructured populations and for populations that have long histories of isolation. However, a relatively common occurrence will be shallow, but geographically structured mtDNA trees--without nuclear gene corroboration--for populations with relatively shorter periods of isolation. This is expected because of the longer coalescence times of nuclear genes (approximately four times that of mtDNA); such cases do not contradict the mtDNA inference of recent isolation and evolutionary divergence. Rather, the nuclear markers are more lagging indicators of changes in population structure. A review of the recent literature on birds reveals the existence of relatively few cases in which nuclear markers contradict mitochondrial markers in a fashion not consistent with coalescent theory. Preliminary information from nuclear genes suggests that mtDNA patterns will prove to be robust indicators of patterns of population history and species limits. At equilibrium, mitochondrial loci are generally a more sensitive indicator of population structure than are nuclear loci, and mitochondrial estimates of F(ST)-like statistics are generally expected to exceed nuclear ones. Hence, invoking behavioural or ecological
Bayesian artificial intelligence
Korb, Kevin B
2010-01-01
Updated and expanded, Bayesian Artificial Intelligence, Second Edition provides a practical and accessible introduction to the main concepts, foundation, and applications of Bayesian networks. It focuses on both the causal discovery of networks and Bayesian inference procedures. Adopting a causal interpretation of Bayesian networks, the authors discuss the use of Bayesian networks for causal modeling. They also draw on their own applied research to illustrate various applications of the technology.New to the Second EditionNew chapter on Bayesian network classifiersNew section on object-oriente
DEFF Research Database (Denmark)
Jensen, Finn Verner; Nielsen, Thomas Dyhre
2016-01-01
Mathematically, a Bayesian graphical model is a compact representation of the joint probability distribution for a set of variables. The most frequently used type of Bayesian graphical models are Bayesian networks. The structural part of a Bayesian graphical model is a graph consisting of nodes and...... largely due to the availability of efficient inference algorithms for answering probabilistic queries about the states of the variables in the network. Furthermore, to support the construction of Bayesian network models, learning algorithms are also available. We give an overview of the Bayesian network...
Directory of Open Access Journals (Sweden)
Chiari Ylenia
2007-01-01
Full Text Available Abstract Background The genus Mantella, endemic poison frogs of Madagascar with 16 described species, are known in the field of international pet trade and entered under the CITES control for the last four years. The phylogeny and phylogeography of this genus have been recently subject of study for conservation purposes. Here we report on the studies of the phylogeography of the Mantella cowani group using a fragment of 453 bp of the mitochondrial cytochrome b gene from 195 individuals from 21 localities. This group is represented by five forms: M. cowani, a critically endangered species, a vulnerable species, M. haraldmeieri, and the non-threatened M. baroni, M. aff. baroni, and M. nigricans. Results The Bayesian phylogenetic and haplotype network analyses revealed the presence of three separated haplotype clades: (1 M. baroni, M. aff. baroni, M. nigricans, and putative hybrids of M. cowani and M. baroni, (2 M. cowani and putative hybrids of M. cowani and M. baroni, and (3 M. haraldmeieri. The putative hybrids were collected from sites where M. cowani and M. baroni live in sympatry. Conclusion These results suggest (a a probable hybridization between M. cowani and M. baroni, (b a lack of genetic differentiation between M. baroni/M. aff. baroni and M. nigricans, (c evidence of recent gene-flow between the northern (M. nigricans, eastern (M. baroni, and south-eastern (M. aff. baroni forms of distinct coloration, and (d the existence of at least three units for conservation in the Mantella cowani group.
Gelman, Andrew; Stern, Hal S; Dunson, David B; Vehtari, Aki; Rubin, Donald B
2013-01-01
FUNDAMENTALS OF BAYESIAN INFERENCEProbability and InferenceSingle-Parameter Models Introduction to Multiparameter Models Asymptotics and Connections to Non-Bayesian ApproachesHierarchical ModelsFUNDAMENTALS OF BAYESIAN DATA ANALYSISModel Checking Evaluating, Comparing, and Expanding ModelsModeling Accounting for Data Collection Decision AnalysisADVANCED COMPUTATION Introduction to Bayesian Computation Basics of Markov Chain Simulation Computationally Efficient Markov Chain Simulation Modal and Distributional ApproximationsREGRESSION MODELS Introduction to Regression Models Hierarchical Linear
Yuan, Ying; MacKinnon, David P.
2009-01-01
This article proposes Bayesian analysis of mediation effects. Compared to conventional frequentist mediation analysis, the Bayesian approach has several advantages. First, it allows researchers to incorporate prior information into the mediation analysis, thus potentially improving the efficiency of estimates. Second, under the Bayesian mediation analysis, inference is straightforward and exact, which makes it appealing for studies with small samples. Third, the Bayesian approach is conceptua...
Directory of Open Access Journals (Sweden)
Joaquín Muñoz
2013-11-01
Full Text Available Since Darwin’s time, waterbirds have been considered an important vector for the dispersal of continental aquatic invertebrates. Bird movements have facilitated the worldwide invasion of the American brine shrimp Artemia franciscana, transporting cysts (diapausing eggs, and favouring rapid range expansions from introduction sites. Here we address the impact of bird migratory flyways on the population genetic structure and phylogeography of A. franciscana in its native range in the Americas. We examined sequence variation for two mitochondrial gene fragments (COI and 16S for a subset of the data in a large set of population samples representing the entire native range of A. franciscana. Furthermore, we performed Mantel tests and redundancy analyses (RDA to test the role of flyways, geography and human introductions on the phylogeography and population genetic structure at a continental scale. A. franciscana mitochondrial DNA was very diverse, with two main clades, largely corresponding to Pacific and Atlantic populations, mirroring American bird flyways. There was a high degree of regional endemism, with populations subdivided into at least 12 divergent, geographically restricted and largely allopatric mitochondrial lineages, and high levels of population structure (ΦST of 0.92, indicating low ongoing gene flow. We found evidence of human-mediated introductions in nine out of 39 populations analysed. Once these populations were removed, Mantel tests revealed a strong association between genetic variation and geographic distance (i.e., isolation-by-distance pattern. RDA showed that shared bird flyways explained around 20% of the variance in genetic distance between populations and this was highly significant, once geographic distance was controlled for. The variance explained increased to 30% when the factor human introduction was included in the model. Our findings suggest that bird-mediated transport of brine shrimp propagules does not result
Phylogeography and Molecular Epidemiology of Yersinia pestis in Madagascar
Vogler, A.; F. Chan; D. Wagner; Roumagnac, P.; J. Lee; Nera, R.; Eppinger, M; Ravel, J; Rahalison, L.; Rasoamanana, B; Beckstrom-Sternberg, S; Achtman, M; Chanteau, S.; Keim, P
2011-01-01
BACKGROUND: Plague was introduced to Madagascar in 1898 and continues to be a significant human health problem. It exists mainly in the central highlands, but in the 1990s was reintroduced to the port city of Mahajanga, where it caused extensive human outbreaks. Despite its prevalence, the phylogeography and molecular epidemiology of Y. pestis in Madagascar has been difficult to study due to the great genetic similarity among isolates. We examine island-wide geographic-genetic patterns based ...
Bayesian Games with Intentions
Bjorndahl, Adam; Halpern, Joseph Y.; Pass, Rafael
2016-01-01
We show that standard Bayesian games cannot represent the full spectrum of belief-dependent preferences. However, by introducing a fundamental distinction between intended and actual strategies, we remove this limitation. We define Bayesian games with intentions, generalizing both Bayesian games and psychological games, and prove that Nash equilibria in psychological games correspond to a special class of equilibria as defined in our setting.
A Bayesian Framework for Active Artificial Perception
Ferreira, Joao; Lobo, Jorge; Bessiere, Pierre; Castelo-Branco, M; Dias, Jorge
2012-01-01
In this text, we present a Bayesian framework for active multimodal perception of 3D structure and motion. The design of this framework finds its inspiration in the role of the dorsal perceptual pathway of the human brain. Its composing models build upon a common egocentric spatial configuration that is naturally fitting for the integration of readings from multiple sensors using a Bayesian approach. In the process, we will contribute with efficient and robust probabilistic solutions for cycl...
Li, Ning; Chen, Xiao; Sun, Dianrong; Song, Na; Lin, Qin; Gao, Tianxiang
2015-08-01
The red stingray Dasyatis akajei is distributed in both marine and freshwater, but little is known about its phylogeography and population structure. We sampled 107 individuals from one freshwater region and 6 coastal localities within the distribution range of D. akajei. Analyses of the first hypervariable region of mitochondrial DNA control region of 474 bp revealed only 17 polymorphism sites that defined 28 haplotypes, with no unique haplotype for the freshwater population. A high level of haplotype diversity and low nucleotide diversity were observed in both marine (h = 0.9393 ± 0.0104, π = 0.0069 ± 0.0040) and freshwater populations (h = 0.8333 ± 0.2224, π = 0.0084 ± 0.0063). Significant level of genetic structure was detected between four marine populations (TZ, WZ, ND and ZZ) via both hierarchical molecular variance analysis (AMOVA) and pairwise FST (with two exceptions), which is unusual for elasmobranchs detected previously over such short geographical distance. However, limited sampling suggested that the freshwater population was not particularly distinct (p > 0.05), but additional samples would be needed to confirm it. Demersal and slow-moving characters likely have contributed to the genetically heterogeneous population structure. The demographic history of D. akajei examined by mismatch distribution analyses, neutrality tests and Bayesian skyline analyses suggested a sudden population expansion dating to upper Pleistocene. The information on genetic diversity and genetic structure will have implications for the management of fisheries and conservation efforts. PMID:24409898
Bayesian Seismology of the Sun
Gruberbauer, Michael
2013-01-01
We perform a Bayesian grid-based analysis of the solar l=0,1,2 and 3 p modes obtained via BiSON in order to deliver the first Bayesian asteroseismic analysis of the solar composition problem. We do not find decisive evidence to prefer either of the contending chemical compositions, although the revised solar abundances (AGSS09) are more probable in general. We do find indications for systematic problems in standard stellar evolution models, unrelated to the consequences of inadequate modelling of the outer layers on the higher-order modes. The seismic observables are best fit by solar models that are several hundred million years older than the meteoritic age of the Sun. Similarly, meteoritic age calibrated models do not adequately reproduce the observed seismic observables. Our results suggest that these problems will affect any asteroseismic inference that relies on a calibration to the Sun.
Rubin, Donald B.
1981-01-01
The Bayesian bootstrap is the Bayesian analogue of the bootstrap. Instead of simulating the sampling distribution of a statistic estimating a parameter, the Bayesian bootstrap simulates the posterior distribution of the parameter; operationally and inferentially the methods are quite similar. Because both methods of drawing inferences are based on somewhat peculiar model assumptions and the resulting inferences are generally sensitive to these assumptions, neither method should be applied wit...
Bayesian analysis of rare events
Straub, Daniel; Papaioannou, Iason; Betz, Wolfgang
2016-06-01
In many areas of engineering and science there is an interest in predicting the probability of rare events, in particular in applications related to safety and security. Increasingly, such predictions are made through computer models of physical systems in an uncertainty quantification framework. Additionally, with advances in IT, monitoring and sensor technology, an increasing amount of data on the performance of the systems is collected. This data can be used to reduce uncertainty, improve the probability estimates and consequently enhance the management of rare events and associated risks. Bayesian analysis is the ideal method to include the data into the probabilistic model. It ensures a consistent probabilistic treatment of uncertainty, which is central in the prediction of rare events, where extrapolation from the domain of observation is common. We present a framework for performing Bayesian updating of rare event probabilities, termed BUS. It is based on a reinterpretation of the classical rejection-sampling approach to Bayesian analysis, which enables the use of established methods for estimating probabilities of rare events. By drawing upon these methods, the framework makes use of their computational efficiency. These methods include the First-Order Reliability Method (FORM), tailored importance sampling (IS) methods and Subset Simulation (SuS). In this contribution, we briefly review these methods in the context of the BUS framework and investigate their applicability to Bayesian analysis of rare events in different settings. We find that, for some applications, FORM can be highly efficient and is surprisingly accurate, enabling Bayesian analysis of rare events with just a few model evaluations. In a general setting, BUS implemented through IS and SuS is more robust and flexible.
Bayesian statistics an introduction
Lee, Peter M
2012-01-01
Bayesian Statistics is the school of thought that combines prior beliefs with the likelihood of a hypothesis to arrive at posterior beliefs. The first edition of Peter Lee’s book appeared in 1989, but the subject has moved ever onwards, with increasing emphasis on Monte Carlo based techniques. This new fourth edition looks at recent techniques such as variational methods, Bayesian importance sampling, approximate Bayesian computation and Reversible Jump Markov Chain Monte Carlo (RJMCMC), providing a concise account of the way in which the Bayesian approach to statistics develops as wel
Understanding Computational Bayesian Statistics
Bolstad, William M
2011-01-01
A hands-on introduction to computational statistics from a Bayesian point of view Providing a solid grounding in statistics while uniquely covering the topics from a Bayesian perspective, Understanding Computational Bayesian Statistics successfully guides readers through this new, cutting-edge approach. With its hands-on treatment of the topic, the book shows how samples can be drawn from the posterior distribution when the formula giving its shape is all that is known, and how Bayesian inferences can be based on these samples from the posterior. These ideas are illustrated on common statistic
Pan-African phylogeography of a model organism, the African clawed frog "Xenopus laevis"
Czech Academy of Sciences Publication Activity Database
Furman, B. L. S.; Bewick, A. J.; Harrison, T. L.; Greenbaum, E.; Gvoždík, Václav; Kusamba, C.; Evans, B. J.
2015-01-01
Roč. 24, č. 4 (2015), s. 909-925. ISSN 0962-1083 Institutional support: RVO:68081766 Keywords : gene flow * phylogeography * population genetics * species limits Subject RIV: EG - Zoology Impact factor: 6.494, year: 2014
An Intuitive Dashboard for Bayesian Network Inference
Reddy, Vikas; Charisse Farr, Anna; Wu, Paul; Mengersen, Kerrie; Yarlagadda, Prasad K. D. V.
2014-03-01
Current Bayesian network software packages provide good graphical interface for users who design and develop Bayesian networks for various applications. However, the intended end-users of these networks may not necessarily find such an interface appealing and at times it could be overwhelming, particularly when the number of nodes in the network is large. To circumvent this problem, this paper presents an intuitive dashboard, which provides an additional layer of abstraction, enabling the end-users to easily perform inferences over the Bayesian networks. Unlike most software packages, which display the nodes and arcs of the network, the developed tool organises the nodes based on the cause-and-effect relationship, making the user-interaction more intuitive and friendly. In addition to performing various types of inferences, the users can conveniently use the tool to verify the behaviour of the developed Bayesian network. The tool has been developed using QT and SMILE libraries in C++.
An Intuitive Dashboard for Bayesian Network Inference
International Nuclear Information System (INIS)
Current Bayesian network software packages provide good graphical interface for users who design and develop Bayesian networks for various applications. However, the intended end-users of these networks may not necessarily find such an interface appealing and at times it could be overwhelming, particularly when the number of nodes in the network is large. To circumvent this problem, this paper presents an intuitive dashboard, which provides an additional layer of abstraction, enabling the end-users to easily perform inferences over the Bayesian networks. Unlike most software packages, which display the nodes and arcs of the network, the developed tool organises the nodes based on the cause-and-effect relationship, making the user-interaction more intuitive and friendly. In addition to performing various types of inferences, the users can conveniently use the tool to verify the behaviour of the developed Bayesian network. The tool has been developed using QT and SMILE libraries in C++
Frühwirth-Schnatter, Sylvia
1990-01-01
In the paper at hand we apply it to Bayesian statistics to obtain "Fuzzy Bayesian Inference". In the subsequent sections we will discuss a fuzzy valued likelihood function, Bayes' theorem for both fuzzy data and fuzzy priors, a fuzzy Bayes' estimator, fuzzy predictive densities and distributions, and fuzzy H.P.D .-Regions. (author's abstract)
Yuan, Ying; MacKinnon, David P.
2009-01-01
In this article, we propose Bayesian analysis of mediation effects. Compared with conventional frequentist mediation analysis, the Bayesian approach has several advantages. First, it allows researchers to incorporate prior information into the mediation analysis, thus potentially improving the efficiency of estimates. Second, under the Bayesian…
A Bayesian Method for Estimating Evolutionary History
Kim, Joungyoun; Anthony, Nicola M; Larget, Bret R.
2012-01-01
Phylogeography is the study of evolutionary history among populations in a species associated with geographic genetic variation. This paper examines the phylogeography of three African gorilla subspecies based on two types of DNA sequence data. One type is HV1, the first hyper-variable region in the control region of the mitochondrial genome. The other type is nuclear mitochondrial DNA (Numt DNA), which results from the introgression of a copy of HV1 from the mitochondrial genome into the nuc...
Kadane, Joseph B.
2009-01-01
This paper reviews the maxims used by three early modern fictional detectives: Monsieur Lecoq, C. Auguste Dupin and Sherlock Holmes. It find similarities between these maxims and Bayesian thought. Poe's Dupin uses ideas very similar to Bayesian game theory. Sherlock Holmes' statements also show thought patterns justifiable in Bayesian terms.
Kadane, Joseph B
2010-01-01
This paper reviews the maxims used by three early modern fictional detectives: Monsieur Lecoq, C. Auguste Dupin and Sherlock Holmes. It find similarities between these maxims and Bayesian thought. Poe's Dupin uses ideas very similar to Bayesian game theory. Sherlock Holmes' statements also show thought patterns justifiable in Bayesian terms.
Directory of Open Access Journals (Sweden)
Jasna Puizina
2013-08-01
Full Text Available Eobania vermiculata (O.F. Müller, 1774, is a typical Mediterranean species of large land snails. Nonindigenous populations of this species, however, are already established in the USA, Australia and elsewhere in the world, where this species is considered to represent a potentially serious threat as a pest and invasive species. The aims of this study were: 1 to determine the pattern of genetic variation within the Croatian E. vermiculata populations based on analyses of sequence diversity of two mitochondrial genes, 16S rDNA and the cytochrome oxidase I (COI, and 2 to shed more light upon the phylogeography of E. vermiculata in this area. Seventy-seven specimens of land snail Eobania vermiculata were sampled at 19 sampling sites along Croatian coastal region and islands. The partial 16S rRNA gene sequences (379 bp grouped into 14 haplotypes, whereas the partial COI gene sequences (523 bp grouped into 13 haplotypes. The overall population is characterized by relatively high haplotype (gene diversity (0.719±0.042 for 16S rDNA and 0.869±0.020 for COI. Demographic Fu F’s tests and Tajima's D value indicated no significant change in the population size, thus suggesting long historical presence of E. vermiculata in this region. Maximum likelihood phylogenetic analysis, Bayesian inference and median joining haplotype network showed a genetic splitting of Croatian 16S rRNA and COI sequences, with a clear distinction between south-Adriatic and north-Adriatic haplotypes. A possible explanation for the observed phylogeography of E. vermiculata, could be related to the climate change, glaciations and the Adriatic Sea level oscillations during the Quaternary
Granade, Christopher; Cory, D G
2015-01-01
In recent years, Bayesian methods have been proposed as a solution to a wide range of issues in quantum state and process tomography. State-of- the-art Bayesian tomography solutions suffer from three problems: numerical intractability, a lack of informative prior distributions, and an inability to track time-dependent processes. Here, we solve all three problems. First, we use modern statistical methods, as pioneered by Husz\\'ar and Houlsby and by Ferrie, to make Bayesian tomography numerically tractable. Our approach allows for practical computation of Bayesian point and region estimators for quantum states and channels. Second, we propose the first informative priors on quantum states and channels. Finally, we develop a method that allows online tracking of time-dependent states and estimates the drift and diffusion processes affecting a state. We provide source code and animated visual examples for our methods.
Bayesian exploratory factor analysis
Gabriella Conti; Sylvia Frühwirth-Schnatter; James Heckman; Rémi Piatek
2014-01-01
This paper develops and applies a Bayesian approach to Exploratory Factor Analysis that improves on ad hoc classical approaches. Our framework relies on dedicated factor models and simultaneously determines the number of factors, the allocation of each measurement to a unique factor, and the corresponding factor loadings. Classical identifi cation criteria are applied and integrated into our Bayesian procedure to generate models that are stable and clearly interpretable. A Monte Carlo study c...
Bayesian Exploratory Factor Analysis
Conti, Gabriella; Frühwirth-Schnatter, Sylvia; Heckman, James J.; Piatek, Rémi
2014-01-01
This paper develops and applies a Bayesian approach to Exploratory Factor Analysis that improves on ad hoc classical approaches. Our framework relies on dedicated factor models and simultaneously determines the number of factors, the allocation of each measurement to a unique factor, and the corresponding factor loadings. Classical identification criteria are applied and integrated into our Bayesian procedure to generate models that are stable and clearly interpretable. A Monte Carlo study co...
Bayesian Exploratory Factor Analysis
Gabriella Conti; Sylvia Fruehwirth-Schnatter; Heckman, James J.; Remi Piatek
2014-01-01
This paper develops and applies a Bayesian approach to Exploratory Factor Analysis that improves on \\emph{ad hoc} classical approaches. Our framework relies on dedicated factor models and simultaneously determines the number of factors, the allocation of each measurement to a unique factor, and the corresponding factor loadings. Classical identification criteria are applied and integrated into our Bayesian procedure to generate models that are stable and clearly interpretable. A Monte Carlo s...
Bayesian exploratory factor analysis
Conti, Gabriella; Frühwirth-Schnatter, Sylvia; Heckman, James J.; Piatek, Rémi
2014-01-01
This paper develops and applies a Bayesian approach to Exploratory Factor Analysis that improves on ad hoc classical approaches. Our framework relies on dedicated factor models and simultaneously determines the number of factors, the allocation of each measurement to a unique factor, and the corresponding factor loadings. Classical identification criteria are applied and integrated into our Bayesian procedure to generate models that are stable and clearly interpretable. A Monte Carlo st...
Bayesian exploratory factor analysis
Conti, Gabriella; Frühwirth-Schnatter, Sylvia; Heckman, James; Piatek, Rémi
2014-01-01
This paper develops and applies a Bayesian approach to Exploratory Factor Analysis that improves on ad hoc classical approaches. Our framework relies on dedicated factor models and simultaneously determines the number of factors, the allocation of each measurement to a unique factor, and the corresponding factor loadings. Classical identification criteria are applied and integrated into our Bayesian procedure to generate models that are stable and clearly interpretable. A Monte Carlo study co...
Carbonetto, Peter; Kisynski, Jacek; De Freitas, Nando; Poole, David L
2012-01-01
The Bayesian Logic (BLOG) language was recently developed for defining first-order probability models over worlds with unknown numbers of objects. It handles important problems in AI, including data association and population estimation. This paper extends BLOG by adopting generative processes over function spaces - known as nonparametrics in the Bayesian literature. We introduce syntax for reasoning about arbitrary collections of objects, and their properties, in an intuitive manner. By expl...
Bayesian default probability models
Andrlíková, Petra
2014-01-01
This paper proposes a methodology for default probability estimation for low default portfolios, where the statistical inference may become troublesome. The author suggests using logistic regression models with the Bayesian estimation of parameters. The piecewise logistic regression model and Box-Cox transformation of credit risk score is used to derive the estimates of probability of default, which extends the work by Neagu et al. (2009). The paper shows that the Bayesian models are more acc...
Numeracy, frequency, and Bayesian reasoning
Directory of Open Access Journals (Sweden)
Gretchen B. Chapman
2009-02-01
Full Text Available Previous research has demonstrated that Bayesian reasoning performance is improved if uncertainty information is presented as natural frequencies rather than single-event probabilities. A questionnaire study of 342 college students replicated this effect but also found that the performance-boosting benefits of the natural frequency presentation occurred primarily for participants who scored high in numeracy. This finding suggests that even comprehension and manipulation of natural frequencies requires a certain threshold of numeracy abilities, and that the beneficial effects of natural frequency presentation may not be as general as previously believed.
Parkinson, C L; Zamudio, K R; Greene, H W
2000-04-01
We used mitochondrial DNA sequences from three gene regions and two tRNAs (ND4, tRNA-HIS-SER, 12S, and 16S rDNA) to investigate the historical ecology of the New World pitviper clade Agkistrodon, with emphasis on the disjunct subspecies of the cantil, A. bilineatus. We found strong evidence that the copperhead (A. contortrix) is basal to its congeners, and that the cottonmouth (A. piscivorus) is basal to cantils. Phylogeography and natural history of the living terminal taxa imply that Agkistrodon primitively occupied relatively temperate habitats, with subsequent evolution of tropicality in ancestral A. bilineatus. Our best supported phylogeny rejects three gulf arc scenarios for the biogeography of A. bilineatus. We find significant statistical support for an initial divergence between populations on the east and west coasts of México and subsequent occupancy of the Yucatán Peninsula, by way of subhumid corridors in northern Central America. Based on phylogenetic relationships, morphological and molecular divergence, and allopatry we elevate A. b. taylori of northeastern México to species status. Taylor's cantil is likely threatened by habitat destruction and small geographical range, and we offer recommendations for its conservation and management. PMID:10736044
Smith, Andrea L; Monteiro, Luis; Hasegawa, Osamu; Friesen, Vicki L
2007-06-01
Factors shaping population differentiation in low latitude seabirds are not well-understood. In this study, we examined global patterns of DNA sequence variation in the mitochondrial control region of the band-rumped storm-petrel (Oceanodroma castro), a highly pelagic seabird distributed across the sub-tropical and tropical Atlantic and Pacific Oceans. Despite previous classification as a single, monotypic species, fixed haplotype differences occurred between Atlantic and Pacific populations, and among all Pacific populations. In addition, Cape Verde and Galapagos birds formed distinct clades, estimated to have diverged from all other populations at least 150,000years ago. Azores hot season populations were also genetically distinct, lending support to previous phenotypic evidence that they be recognized as a separate species. Seasonal populations in Madeira probably represent separate genetic management units. The phylogeography of the band-rumped storm-petrel appears to have been shaped by both nonphysical barriers to gene flow and Pleistocene oceanographic conditions. Ancestral populations likely expanded through contiguous range expansion and infrequent long-distance colonization into their current breeding range. These findings suggest several possible revisions to the taxonomy of the band-rumped storm-petrel. PMID:17408975
Comparative phylogeography of African savannah ungulates
DEFF Research Database (Denmark)
Lorenzen, Eline; Heller, Rasmus; Siegismund, Hans Redlef
2012-01-01
The savannah biome of sub-Saharan Africa harbours the highest diversity of ungulates (hoofed mammals) on Earth. In this review, we compile population genetic data from 19 codistributed ungulate taxa of the savannah biome and find striking concordance in the phylogeographic structuring of species....
A review of marine phylogeography in southern Africa
Directory of Open Access Journals (Sweden)
Nigel P. Barker
2011-05-01
Full Text Available The southern African marine realm is located at the transition zone between the Atlantic and Indo-Pacific biomes. Its biodiversity is particularly rich and comprises faunal and floral elements from the two major oceanic regions, as well as a large number of endemics. Within this realm, strikingly different biota occur in close geographic proximity to each other, and many of the species with distributions spanning two or more of the region’s marine biogeographic provinces are divided into evolutionary units that can often only be distinguished on the basis of genetic data. In this review, we describe the state of marine phylogeography in southern Africa, that is, the study of evolutionary relationships at the species level, or amongst closely related species, in relation to the region’s marine environment. We focus particularly on coastal phylogeography, where much progress has recently been made in identifying phylogeographic breaks and explaining how they originated and are maintained. We also highlight numerous shortcomings that should be addressed in the near future. These include: the limited data available for commercially important organisms, particularly offshore species; the paucity of oceanographic data for nearshore areas; a dearth of studies based on multilocus data; and the fact that studying the role of diversifying selection in speciation has been limited to physiological approaches to the exclusion of genetics. It is becoming apparent that the southern African marine realm is one of the world’s most interesting environments in which to study the evolutionary processes that shape not only regional, but also global patterns of marine biodiversity.
Comparative phylogeography in rainforest trees from Lower Guinea, Africa.
Heuertz, Myriam; Duminil, Jérôme; Dauby, Gilles; Savolainen, Vincent; Hardy, Olivier J
2014-01-01
Comparative phylogeography is an effective approach to assess the evolutionary history of biological communities. We used comparative phylogeography in fourteen tree taxa from Lower Guinea (Atlantic Equatorial Africa) to test for congruence with two simple evolutionary scenarios based on physio-climatic features 1) the W-E environmental gradient and 2) the N-S seasonal inversion, which determine climatic and seasonality differences in the region. We sequenced the trnC-ycf6 plastid DNA region using a dual sampling strategy: fourteen taxa with small sample sizes (dataset 1, mean n = 16/taxon), to assess whether a strong general pattern of allele endemism and genetic differentiation emerged; and four taxonomically well-studied species with larger sample sizes (dataset 2, mean n = 109/species) to detect the presence of particular shared phylogeographic patterns. When grouping the samples into two alternative sets of two populations, W and E, vs. N and S, neither dataset exhibited a strong pattern of allelic endemism, suggesting that none of the considered regions consistently harboured older populations. Differentiation in dataset 1 was similarly strong between W and E as between N and S, with 3-5 significant F ST tests out of 14 tests in each scenario. Coalescent simulations indicated that, given the power of the data, this result probably reflects idiosyncratic histories of the taxa, or a weak common differentiation pattern (possibly with population substructure) undetectable across taxa in dataset 1. Dataset 2 identified a common genetic break separating the northern and southern populations of Greenwayodendron suaveolens subsp. suaveolens var. suaveolens, Milicia excelsa, Symphonia globulifera and Trichoscypha acuminata in Lower Guinea, in agreement with differentiation across the N-S seasonal inversion. Our work suggests that currently recognized tree taxa or suspected species complexes can contain strongly differentiated genetic lineages, which could lead
Comparative phylogeography in rainforest trees from Lower Guinea, Africa.
Directory of Open Access Journals (Sweden)
Myriam Heuertz
Full Text Available Comparative phylogeography is an effective approach to assess the evolutionary history of biological communities. We used comparative phylogeography in fourteen tree taxa from Lower Guinea (Atlantic Equatorial Africa to test for congruence with two simple evolutionary scenarios based on physio-climatic features 1 the W-E environmental gradient and 2 the N-S seasonal inversion, which determine climatic and seasonality differences in the region. We sequenced the trnC-ycf6 plastid DNA region using a dual sampling strategy: fourteen taxa with small sample sizes (dataset 1, mean n = 16/taxon, to assess whether a strong general pattern of allele endemism and genetic differentiation emerged; and four taxonomically well-studied species with larger sample sizes (dataset 2, mean n = 109/species to detect the presence of particular shared phylogeographic patterns. When grouping the samples into two alternative sets of two populations, W and E, vs. N and S, neither dataset exhibited a strong pattern of allelic endemism, suggesting that none of the considered regions consistently harboured older populations. Differentiation in dataset 1 was similarly strong between W and E as between N and S, with 3-5 significant F ST tests out of 14 tests in each scenario. Coalescent simulations indicated that, given the power of the data, this result probably reflects idiosyncratic histories of the taxa, or a weak common differentiation pattern (possibly with population substructure undetectable across taxa in dataset 1. Dataset 2 identified a common genetic break separating the northern and southern populations of Greenwayodendron suaveolens subsp. suaveolens var. suaveolens, Milicia excelsa, Symphonia globulifera and Trichoscypha acuminata in Lower Guinea, in agreement with differentiation across the N-S seasonal inversion. Our work suggests that currently recognized tree taxa or suspected species complexes can contain strongly differentiated genetic lineages
Dimensionality reduction in Bayesian estimation algorithms
Directory of Open Access Journals (Sweden)
G. W. Petty
2013-03-01
Full Text Available An idealized synthetic database loosely resembling 3-channel passive microwave observations of precipitation against a variable background is employed to examine the performance of a conventional Bayesian retrieval algorithm. For this dataset, algorithm performance is found to be poor owing to an irreconcilable conflict between the need to find matches in the dependent database versus the need to exclude inappropriate matches. It is argued that the likelihood of such conflicts increases sharply with the dimensionality of the observation space of real satellite sensors, which may utilize 9 to 13 channels to retrieve precipitation, for example. An objective method is described for distilling the relevant information content from N real channels into a much smaller number (M of pseudochannels while also regularizing the background (geophysical plus instrument noise component. The pseudochannels are linear combinations of the original N channels obtained via a two-stage principal component analysis of the dependent dataset. Bayesian retrievals based on a single pseudochannel applied to the independent dataset yield striking improvements in overall performance. The differences between the conventional Bayesian retrieval and reduced-dimensional Bayesian retrieval suggest that a major potential problem with conventional multichannel retrievals – whether Bayesian or not – lies in the common but often inappropriate assumption of diagonal error covariance. The dimensional reduction technique described herein avoids this problem by, in effect, recasting the retrieval problem in a coordinate system in which the desired covariance is lower-dimensional, diagonal, and unit magnitude.
Tactile length contraction as Bayesian inference.
Tong, Jonathan; Ngo, Vy; Goldreich, Daniel
2016-08-01
To perceive, the brain must interpret stimulus-evoked neural activity. This is challenging: The stochastic nature of the neural response renders its interpretation inherently uncertain. Perception would be optimized if the brain used Bayesian inference to interpret inputs in light of expectations derived from experience. Bayesian inference would improve perception on average but cause illusions when stimuli violate expectation. Intriguingly, tactile, auditory, and visual perception are all prone to length contraction illusions, characterized by the dramatic underestimation of the distance between punctate stimuli delivered in rapid succession; the origin of these illusions has been mysterious. We previously proposed that length contraction illusions occur because the brain interprets punctate stimulus sequences using Bayesian inference with a low-velocity expectation. A novel prediction of our Bayesian observer model is that length contraction should intensify if stimuli are made more difficult to localize. Here we report a tactile psychophysical study that tested this prediction. Twenty humans compared two distances on the forearm: a fixed reference distance defined by two taps with 1-s temporal separation and an adjustable comparison distance defined by two taps with temporal separation t ≤ 1 s. We observed significant length contraction: As t was decreased, participants perceived the two distances as equal only when the comparison distance was made progressively greater than the reference distance. Furthermore, the use of weaker taps significantly enhanced participants' length contraction. These findings confirm the model's predictions, supporting the view that the spatiotemporal percept is a best estimate resulting from a Bayesian inference process. PMID:27121574
Bayesian Inference Methods for Sparse Channel Estimation
DEFF Research Database (Denmark)
Pedersen, Niels Lovmand
2013-01-01
inference algorithms based on the proposed prior representation for sparse channel estimation in orthogonal frequency-division multiplexing receivers. The inference algorithms, which are mainly obtained from variational Bayesian methods, exploit the underlying sparse structure of wireless channel responses......This thesis deals with sparse Bayesian learning (SBL) with application to radio channel estimation. As opposed to the classical approach for sparse signal representation, we focus on the problem of inferring complex signals. Our investigations within SBL constitute the basis for the development of...... Bayesian inference algorithms for sparse channel estimation. Sparse inference methods aim at finding the sparse representation of a signal given in some overcomplete dictionary of basis vectors. Within this context, one of our main contributions to the field of SBL is a hierarchical representation of...
Bayesian least squares deconvolution
Asensio Ramos, A.; Petit, P.
2015-11-01
Aims: We develop a fully Bayesian least squares deconvolution (LSD) that can be applied to the reliable detection of magnetic signals in noise-limited stellar spectropolarimetric observations using multiline techniques. Methods: We consider LSD under the Bayesian framework and we introduce a flexible Gaussian process (GP) prior for the LSD profile. This prior allows the result to automatically adapt to the presence of signal. We exploit several linear algebra identities to accelerate the calculations. The final algorithm can deal with thousands of spectral lines in a few seconds. Results: We demonstrate the reliability of the method with synthetic experiments and we apply it to real spectropolarimetric observations of magnetic stars. We are able to recover the magnetic signals using a small number of spectral lines, together with the uncertainty at each velocity bin. This allows the user to consider if the detected signal is reliable. The code to compute the Bayesian LSD profile is freely available.
Bayesian least squares deconvolution
Ramos, A Asensio
2015-01-01
Aims. To develop a fully Bayesian least squares deconvolution (LSD) that can be applied to the reliable detection of magnetic signals in noise-limited stellar spectropolarimetric observations using multiline techniques. Methods. We consider LSD under the Bayesian framework and we introduce a flexible Gaussian Process (GP) prior for the LSD profile. This prior allows the result to automatically adapt to the presence of signal. We exploit several linear algebra identities to accelerate the calculations. The final algorithm can deal with thousands of spectral lines in a few seconds. Results. We demonstrate the reliability of the method with synthetic experiments and we apply it to real spectropolarimetric observations of magnetic stars. We are able to recover the magnetic signals using a small number of spectral lines, together with the uncertainty at each velocity bin. This allows the user to consider if the detected signal is reliable. The code to compute the Bayesian LSD profile is freely available.
Loredo, T J
2004-01-01
I describe a framework for adaptive scientific exploration based on iterating an Observation--Inference--Design cycle that allows adjustment of hypotheses and observing protocols in response to the results of observation on-the-fly, as data are gathered. The framework uses a unified Bayesian methodology for the inference and design stages: Bayesian inference to quantify what we have learned from the available data and predict future data, and Bayesian decision theory to identify which new observations would teach us the most. When the goal of the experiment is simply to make inferences, the framework identifies a computationally efficient iterative ``maximum entropy sampling'' strategy as the optimal strategy in settings where the noise statistics are independent of signal properties. Results of applying the method to two ``toy'' problems with simulated data--measuring the orbit of an extrasolar planet, and locating a hidden one-dimensional object--show the approach can significantly improve observational eff...
Bayesian and frequentist inequality tests
David M. Kaplan; Zhuo, Longhao
2016-01-01
Bayesian and frequentist criteria are fundamentally different, but often posterior and sampling distributions are asymptotically equivalent (and normal). We compare Bayesian and frequentist hypothesis tests of inequality restrictions in such cases. For finite-dimensional parameters, if the null hypothesis is that the parameter vector lies in a certain half-space, then the Bayesian test has (frequentist) size $\\alpha$; if the null hypothesis is any other convex subspace, then the Bayesian test...
Bayesian multiple target tracking
Streit, Roy L
2013-01-01
This second edition has undergone substantial revision from the 1999 first edition, recognizing that a lot has changed in the multiple target tracking field. One of the most dramatic changes is in the widespread use of particle filters to implement nonlinear, non-Gaussian Bayesian trackers. This book views multiple target tracking as a Bayesian inference problem. Within this framework it develops the theory of single target tracking, multiple target tracking, and likelihood ratio detection and tracking. In addition to providing a detailed description of a basic particle filter that implements
Bayesian Exploratory Factor Analysis
DEFF Research Database (Denmark)
Conti, Gabriella; Frühwirth-Schnatter, Sylvia; Heckman, James J.;
2014-01-01
This paper develops and applies a Bayesian approach to Exploratory Factor Analysis that improves on ad hoc classical approaches. Our framework relies on dedicated factor models and simultaneously determines the number of factors, the allocation of each measurement to a unique factor, and the...... corresponding factor loadings. Classical identification criteria are applied and integrated into our Bayesian procedure to generate models that are stable and clearly interpretable. A Monte Carlo study confirms the validity of the approach. The method is used to produce interpretable low dimensional aggregates...
Bayesian Geostatistical Design
DEFF Research Database (Denmark)
Diggle, Peter; Lophaven, Søren Nymand
2006-01-01
locations to, or deletion of locations from, an existing design, and prospective design, which consists of choosing positions for a new set of sampling locations. We propose a Bayesian design criterion which focuses on the goal of efficient spatial prediction whilst allowing for the fact that model...
Czech Academy of Sciences Publication Activity Database
Krejsa, Jiří; Věchet, S.
Bratislava: Slovak University of Technology in Bratislava, 2010, s. 217-222. ISBN 978-80-227-3353-3. [Robotics in Education . Bratislava (SK), 16.09.2010-17.09.2010] Institutional research plan: CEZ:AV0Z20760514 Keywords : mobile robot localization * bearing only beacons * Bayesian filters Subject RIV: JD - Computer Applications, Robotics
DEFF Research Database (Denmark)
Antoniou, Constantinos; Harrison, Glenn W.; Lau, Morten I.;
2015-01-01
A large literature suggests that many individuals do not apply Bayes’ Rule when making decisions that depend on them correctly pooling prior information and sample data. We replicate and extend a classic experimental study of Bayesian updating from psychology, employing the methods of experimenta...
Bayesian Independent Component Analysis
DEFF Research Database (Denmark)
Winther, Ole; Petersen, Kaare Brandt
2007-01-01
In this paper we present an empirical Bayesian framework for independent component analysis. The framework provides estimates of the sources, the mixing matrix and the noise parameters, and is flexible with respect to choice of source prior and the number of sources and sensors. Inside the engine...
Noncausal Bayesian Vector Autoregression
DEFF Research Database (Denmark)
Lanne, Markku; Luoto, Jani
We propose a Bayesian inferential procedure for the noncausal vector autoregressive (VAR) model that is capable of capturing nonlinearities and incorporating effects of missing variables. In particular, we devise a fast and reliable posterior simulator that yields the predictive distribution as a...
Loredo, Thomas J.
2004-04-01
I describe a framework for adaptive scientific exploration based on iterating an Observation-Inference-Design cycle that allows adjustment of hypotheses and observing protocols in response to the results of observation on-the-fly, as data are gathered. The framework uses a unified Bayesian methodology for the inference and design stages: Bayesian inference to quantify what we have learned from the available data and predict future data, and Bayesian decision theory to identify which new observations would teach us the most. When the goal of the experiment is simply to make inferences, the framework identifies a computationally efficient iterative ``maximum entropy sampling'' strategy as the optimal strategy in settings where the noise statistics are independent of signal properties. Results of applying the method to two ``toy'' problems with simulated data-measuring the orbit of an extrasolar planet, and locating a hidden one-dimensional object-show the approach can significantly improve observational efficiency in settings that have well-defined nonlinear models. I conclude with a list of open issues that must be addressed to make Bayesian adaptive exploration a practical and reliable tool for optimizing scientific exploration.
Bayesian logistic regression analysis
Van Erp, H.R.N.; Van Gelder, P.H.A.J.M.
2012-01-01
In this paper we present a Bayesian logistic regression analysis. It is found that if one wishes to derive the posterior distribution of the probability of some event, then, together with the traditional Bayes Theorem and the integrating out of nuissance parameters, the Jacobian transformation is an
Lee, Kyung Min; Yang, Eun Chan; Coyer, James A.; Zuccarello, Giuseppe C.; Wang, Wei-Lui; Choi, Chang Geun; Boo, Sung Min
2012-01-01
Although benthic marine algae are essential components of marine coastal systems that have been influenced profoundly by past and present climate change, our knowledge of seaweed phylogeography is limited. The brown alga Ishige okamurae Yendo occurs in the northwest Pacific, where it occupies a char
Czech Academy of Sciences Publication Activity Database
Kindler, C.; Böhme, W.; Corti, C.; Gvoždík, Václav; Jablonski, D.; Jandzik, D.; Metallinou, M.; Široký, P.; Fritz, U.
2013-01-01
Roč. 42, č. 5 (2013), s. 458-472. ISSN 0300-3256 Institutional support: RVO:67985904 Keywords : TURTLES EMYS-ORBICULARIS * NUCLEAR-DNA SEQUENCES * MOLECULAR PHYLOGEOGRAPHY Subject RIV: EG - Zoology Impact factor: 2.922, year: 2013
Antunes, Jorge T.; Pedro N. Leão; Vasconcelos, Vítor M.
2015-01-01
Cylindrospermopsis raciborskii is a cyanobacterial species extensively studied for its toxicity, bloom formation and invasiveness potential, which have consequences to public and environmental health. Its current geographical distribution, spanning different climates, suggests that C. raciborskii has acquired the status of a cosmopolitan species. From phylogeography studies, a tropical origin for this species seems convincing, with different conjectural routes of expansion toward temperate cl...
Garcia-Rodriguez, A. I.; Bowen, B.W.; Domning, D.; Mignucci-Giannoni, A. A.; Marmontel, M.; Montoya-Ospina, R. A.; Morales-Vela, B.; Rudin, M.; Bonde, R.K.; McGuire, P.M.
1998-01-01
To resolve the population genetic structure and phylogeography of the West Indian manatee (Trichechus manatus), mitochondrial (mt) DNA control region sequences were compared among eight locations across the western Atlantic region. Fifteen haplotypes were identified among 86 individuals from Florida, Puerto Rico, the Dominican Republic, Mexico, Colombia, Venezuela, Guyana and Brazil. Despite the manatee's ability to move thousands of kilometres along continental margins, strong population separations between most locations were demonstrated with significant haplotype frequency shifts. These findings are consistent with tagging studies which indicate that stretches of open water and unsuitable coastal habitats constitute substantial barriers to gene flow and colonization. Low levels of genetic diversity within Florida and Brazilian samples might be explained by recent colonization into high latitudes or bottleneck effects. Three distinctive mtDNA lineages were observed in an intraspecific phylogeny of T. manatus, corresponding approximately to: (i) Florida and the West Indies; (ii) the Gulf of Mexico to the Caribbean rivers of South America; and (iii) the northeast Atlantic coast of South America. These lineages, which are not concordant with previous subspecies designations, are separated by sequence divergence estimates of d = 0.04-0.07, approximately the same level of divergence observed between T. manatus and the Amazonian manatee (T. inunguis, n = 16). Three individuals from Guyana, identified as T. manatus, had mtDNA haplotypes which are affiliated with the endemic Amazon form T. inunguis. The three primary T. manatus lineages and the T. inunguis lineage may represent relatively deep phylogeographic partitions which have been bridged recently due to changes in habitat availability (after the Wisconsin glacial period, 10 000 BP), natural colonization, and human-mediated transplantation.
Directory of Open Access Journals (Sweden)
Bin Wang
Full Text Available BACKGROUND: The influence of Pleistocene climatic fluctuations on intraspecific diversification in the Qinling-Daba Mountains of East Asia remains poorly investigated. We tested hypotheses concerning refugia during the last glacial maximum (LGM in this region by examining the phylogeography of the swelled vent frog (Feirana quadranus; Dicroglossidae, Anura, Amphibia. METHODOLOGY/PRINCIPAL FINDINGS: We obtained complete mitochondrial ND2 gene sequences of 224 individuals from 34 populations of Feirana quadranus for phylogeographic analyses. Additionally, we obtained nuclear tyrosinase gene sequences of 68 F. quadranus, one F. kangxianensis and three F. taihangnica samples to test for mitochondrial introgression among them. Phylogenetic analyses based on all genes revealed no introgression among them. Phylogenetic analyses based on ND2 datasets revealed that F. quadranus was comprised of six lineages which were separated by deep valleys; the sole exception is that the Main Qinling and Micang-Western Qinling lineages overlap in distribution. Analyses of population structure indicated restricted gene flow among lineages. Coalescent simulations and divergence dating indicated that the basal diversification within F. quadranus may be associated with the dramatic uplifts of the Tibetan Plateau during the Pliocene. Coalescent simulations indicated that Wuling, Daba, and Western Qinling-Micang-Longmen Mountains were refugia for F. quadranus during the LGM. Demographic analyses indicated that the Daba lineage experienced population size increase prior to the LGM but the Main Qinling and the Micang-Western Qinling lineages expanded in population size and range after the LGM, and the other lineages almost have stable population size or slight slow population size decline. CONCLUSIONS/SIGNIFICANCE: The Qinling-Daba Mountains hosted three refugia for F. quadranus during the LGM. Populations that originated in the Daba Mountains colonized the Main Qinling
Quantum Bayesianism at the Perimeter
Fuchs, Christopher A
2010-01-01
The author summarizes the Quantum Bayesian viewpoint of quantum mechanics, developed originally by C. M. Caves, R. Schack, and himself. It is a view crucially dependent upon the tools of quantum information theory. Work at the Perimeter Institute for Theoretical Physics continues the development and is focused on the hard technical problem of a finding a good representation of quantum mechanics purely in terms of probabilities, without amplitudes or Hilbert-space operators. The best candidate representation involves a mysterious entity called a symmetric informationally complete quantum measurement. Contemplation of it gives a way of thinking of the Born Rule as an addition to the rules of probability theory, applicable when one gambles on the consequences of interactions with physical systems. The article ends by outlining some directions for future work.
State Information in Bayesian Games
Cuff, Paul
2009-01-01
Two-player zero-sum repeated games are well understood. Computing the value of such a game is straightforward. Additionally, if the payoffs are dependent on a random state of the game known to one, both, or neither of the players, the resulting value of the game has been analyzed under the framework of Bayesian games. This investigation considers the optimal performance in a game when a helper is transmitting state information to one of the players. Encoding information for an adversarial setting (game) requires a different result than rate-distortion theory provides. Game theory has accentuated the importance of randomization (mixed strategy), which does not find a significant role in most communication modems and source coding codecs. Higher rates of communication, used in the right way, allow the message to include the necessary random component useful in games.
Probability and Bayesian statistics
1987-01-01
This book contains selected and refereed contributions to the "Inter national Symposium on Probability and Bayesian Statistics" which was orga nized to celebrate the 80th birthday of Professor Bruno de Finetti at his birthplace Innsbruck in Austria. Since Professor de Finetti died in 1985 the symposium was dedicated to the memory of Bruno de Finetti and took place at Igls near Innsbruck from 23 to 26 September 1986. Some of the pa pers are published especially by the relationship to Bruno de Finetti's scientific work. The evolution of stochastics shows growing importance of probability as coherent assessment of numerical values as degrees of believe in certain events. This is the basis for Bayesian inference in the sense of modern statistics. The contributions in this volume cover a broad spectrum ranging from foundations of probability across psychological aspects of formulating sub jective probability statements, abstract measure theoretical considerations, contributions to theoretical statistics an...
Bayesian Magic in Asteroseismology
Kallinger, T.
2015-09-01
Only a few years ago asteroseismic observations were so rare that scientists had plenty of time to work on individual data sets. They could tune their algorithms in any possible way to squeeze out the last bit of information. Nowadays this is impossible. With missions like MOST, CoRoT, and Kepler we basically drown in new data every day. To handle this in a sufficient way statistical methods become more and more important. This is why Bayesian techniques started their triumph march across asteroseismology. I will go with you on a journey through Bayesian Magic Land, that brings us to the sea of granulation background, the forest of peakbagging, and the stony alley of model comparison.
Approximate Bayesian Computation: a nonparametric perspective
Blum, Michael
2010-01-01
Approximate Bayesian Computation is a family of likelihood-free inference techniques that are well-suited to models defined in terms of a stochastic generating mechanism. In a nutshell, Approximate Bayesian Computation proceeds by computing summary statistics s_obs from the data and simulating summary statistics for different values of the parameter theta. The posterior distribution is then approximated by an estimator of the conditional density g(theta|s_obs). In this paper, we derive the asymptotic bias and variance of the standard estimators of the posterior distribution which are based on rejection sampling and linear adjustment. Additionally, we introduce an original estimator of the posterior distribution based on quadratic adjustment and we show that its bias contains a fewer number of terms than the estimator with linear adjustment. Although we find that the estimators with adjustment are not universally superior to the estimator based on rejection sampling, we find that they can achieve better perfor...
Bayesian Nonparametric Graph Clustering
Banerjee, Sayantan; Akbani, Rehan; Baladandayuthapani, Veerabhadran
2015-01-01
We present clustering methods for multivariate data exploiting the underlying geometry of the graphical structure between variables. As opposed to standard approaches that assume known graph structures, we first estimate the edge structure of the unknown graph using Bayesian neighborhood selection approaches, wherein we account for the uncertainty of graphical structure learning through model-averaged estimates of the suitable parameters. Subsequently, we develop a nonparametric graph cluster...
Approximate Bayesian recursive estimation
Czech Academy of Sciences Publication Activity Database
Kárný, Miroslav
2014-01-01
Roč. 285, č. 1 (2014), s. 100-111. ISSN 0020-0255 R&D Projects: GA ČR GA13-13502S Institutional support: RVO:67985556 Keywords : Approximate parameter estimation * Bayesian recursive estimation * Kullback–Leibler divergence * Forgetting Subject RIV: BB - Applied Statistics, Operational Research Impact factor: 4.038, year: 2014 http://library.utia.cas.cz/separaty/2014/AS/karny-0425539.pdf
Bayesian Benchmark Dose Analysis
Fang, Qijun; Piegorsch, Walter W.; Barnes, Katherine Y.
2014-01-01
An important objective in environmental risk assessment is estimation of minimum exposure levels, called Benchmark Doses (BMDs) that induce a pre-specified Benchmark Response (BMR) in a target population. Established inferential approaches for BMD analysis typically involve one-sided, frequentist confidence limits, leading in practice to what are called Benchmark Dose Lower Limits (BMDLs). Appeal to Bayesian modeling and credible limits for building BMDLs is far less developed, however. Indee...
Bayesian Generalized Rating Curves
Helgi Sigurðarson 1985
2014-01-01
A rating curve is a curve or a model that describes the relationship between water elevation, or stage, and discharge in an observation site in a river. The rating curve is fit from paired observations of stage and discharge. The rating curve then predicts discharge given observations of stage and this methodology is applied as stage is substantially easier to directly observe than discharge. In this thesis a statistical rating curve model is proposed working within the framework of Bayesian...
Heteroscedastic Treed Bayesian Optimisation
Assael, John-Alexander M.; Wang, Ziyu; Shahriari, Bobak; De Freitas, Nando
2014-01-01
Optimising black-box functions is important in many disciplines, such as tuning machine learning models, robotics, finance and mining exploration. Bayesian optimisation is a state-of-the-art technique for the global optimisation of black-box functions which are expensive to evaluate. At the core of this approach is a Gaussian process prior that captures our belief about the distribution over functions. However, in many cases a single Gaussian process is not flexible enough to capture non-stat...
Efficient Bayesian Phase Estimation
Wiebe, Nathan; Granade, Chris
2016-07-01
We introduce a new method called rejection filtering that we use to perform adaptive Bayesian phase estimation. Our approach has several advantages: it is classically efficient, easy to implement, achieves Heisenberg limited scaling, resists depolarizing noise, tracks time-dependent eigenstates, recovers from failures, and can be run on a field programmable gate array. It also outperforms existing iterative phase estimation algorithms such as Kitaev's method.
Brody, Samuel; Lapata, Mirella
2009-01-01
Sense induction seeks to automatically identify word senses directly from a corpus. A key assumption underlying previous work is that the context surrounding an ambiguous word is indicative of its meaning. Sense induction is thus typically viewed as an unsupervised clustering problem where the aim is to partition a word’s contexts into different classes, each representing a word sense. Our work places sense induction in a Bayesian context by modeling the contexts of the ambiguous word as samp...
Bayesian Neural Word Embedding
Barkan, Oren
2016-01-01
Recently, several works in the domain of natural language processing presented successful methods for word embedding. Among them, the Skip-gram (SG) with negative sampling, known also as Word2Vec, advanced the state-of-the-art of various linguistics tasks. In this paper, we propose a scalable Bayesian neural word embedding algorithm that can be beneficial to general item similarity tasks as well. The algorithm relies on a Variational Bayes solution for the SG objective and a detailed step by ...
Wiegerinck, Wim; Schoenaker, Christiaan; Duane, Gregory
2016-04-01
Recently, methods for model fusion by dynamically combining model components in an interactive ensemble have been proposed. In these proposals, fusion parameters have to be learned from data. One can view these systems as parametrized dynamical systems. We address the question of learnability of dynamical systems with respect to both short term (vector field) and long term (attractor) behavior. In particular we are interested in learning in the imperfect model class setting, in which the ground truth has a higher complexity than the models, e.g. due to unresolved scales. We take a Bayesian point of view and we define a joint log-likelihood that consists of two terms, one is the vector field error and the other is the attractor error, for which we take the L1 distance between the stationary distributions of the model and the assumed ground truth. In the context of linear models (like so-called weighted supermodels), and assuming a Gaussian error model in the vector fields, vector field learning leads to a tractable Gaussian solution. This solution can then be used as a prior for the next step, Bayesian attractor learning, in which the attractor error is used as a log-likelihood term. Bayesian attractor learning is implemented by elliptical slice sampling, a sampling method for systems with a Gaussian prior and a non Gaussian likelihood. Simulations with a partially observed driven Lorenz 63 system illustrate the approach.
Bayesian theory and applications
Dellaportas, Petros; Polson, Nicholas G; Stephens, David A
2013-01-01
The development of hierarchical models and Markov chain Monte Carlo (MCMC) techniques forms one of the most profound advances in Bayesian analysis since the 1970s and provides the basis for advances in virtually all areas of applied and theoretical Bayesian statistics. This volume guides the reader along a statistical journey that begins with the basic structure of Bayesian theory, and then provides details on most of the past and present advances in this field. The book has a unique format. There is an explanatory chapter devoted to each conceptual advance followed by journal-style chapters that provide applications or further advances on the concept. Thus, the volume is both a textbook and a compendium of papers covering a vast range of topics. It is appropriate for a well-informed novice interested in understanding the basic approach, methods and recent applications. Because of its advanced chapters and recent work, it is also appropriate for a more mature reader interested in recent applications and devel...
Opatova, Vera; Arnedo, Miquel A.
2014-01-01
Studies conducted on volcanic islands have greatly contributed to our current understanding of how organisms diversify. The Canary Islands archipelago, located northwest of the coast of northern Africa, harbours a large number of endemic taxa. Because of their low vagility, mygalomorph spiders are usually absent from oceanic islands. The spider Titanidiops canariensis, which inhabits the easternmost islands of the archipelago, constitutes an exception to this rule. Here, we use a multi-locus approach that combines three mitochondrial and four nuclear genes to investigate the origins and phylogeography of this remarkable trap-door spider. We provide a timeframe for the colonisation of the Canary Islands using two alternative approaches: concatenation and species tree inference in a Bayesian relaxed clock framework. Additionally, we investigate the existence of cryptic species on the islands by means of a Bayesian multi-locus species delimitation method. Our results indicate that T. canariensis colonised the Canary Islands once, most likely during the Miocene, although discrepancies between the timeframes from different approaches make the exact timing uncertain. A complex evolutionary history for the species in the archipelago is revealed, which involves two independent colonisations of Fuerteventura from the ancestral range of T. canariensis in northern Lanzarote and a possible back colonisation of southern Lanzarote. The data further corroborate a previously proposed volcanic refugium, highlighting the impact of the dynamic volcanic history of the island on the phylogeographic patterns of the endemic taxa. T. canariensis includes at least two different species, one inhabiting the Jandia peninsula and central Fuerteventura and one spanning from central Fuerteventura to Lanzarote. Our data suggest that the extant northern African Titanidiops lineages may have expanded to the region after the islands were colonised and, hence, are not the source of colonisation. In
A Bayesian Modelling of Wildfires in Portugal
Silva, Giovani L.; Soares, Paulo; Marques, Susete; Dias, Inês M.; Oliveira, Manuela M.; Borges, Guilherme J.
2015-01-01
In the last decade wildfires became a serious problem in Portugal due to different issues such as climatic characteristics and nature of Portuguese forest. In order to analyse wildfire data, we employ beta regression for modelling the proportion of burned forest area, under a Bayesian perspective. Our main goal is to find out fire risk factors that influence the proportion of area burned and what may make a forest type susceptible or resistant to fire. Then, we analyse wildfire...
Summary Statistics in Approximate Bayesian Computation
Prangle, Dennis
2015-01-01
This document is due to appear as a chapter of the forthcoming Handbook of Approximate Bayesian Computation (ABC) edited by S. Sisson, Y. Fan, and M. Beaumont. Since the earliest work on ABC, it has been recognised that using summary statistics is essential to produce useful inference results. This is because ABC suffers from a curse of dimensionality effect, whereby using high dimensional inputs causes large approximation errors in the output. It is therefore crucial to find low dimensional ...
Improving Environmental Scanning Systems Using Bayesian Networks
Simon Welter; Jörg H. Mayer; Reiner Quick
2013-01-01
As companies’ environment is becoming increasingly volatile, scanning systems gain in importance. We propose a hybrid process model for such systems' information gathering and interpretation tasks that combines quantitative information derived from regression analyses and qualitative knowledge from expert interviews. For the latter, we apply Bayesian networks. We derive the need for such a hybrid process model from a literature review. We lay out our model to find a suitable set of business e...
Unbounded Bayesian Optimization via Regularization
Shahriari, Bobak; Bouchard-Côté, Alexandre; De Freitas, Nando
2015-01-01
Bayesian optimization has recently emerged as a popular and efficient tool for global optimization and hyperparameter tuning. Currently, the established Bayesian optimization practice requires a user-defined bounding box which is assumed to contain the optimizer. However, when little is known about the probed objective function, it can be difficult to prescribe such bounds. In this work we modify the standard Bayesian optimization framework in a principled way to allow automatic resizing of t...
Bayesian ensemble refinement by replica simulations and reweighting
Hummer, Gerhard
2015-01-01
We describe different Bayesian ensemble refinement methods, examine their interrelation, and discuss their practical application. With ensemble refinement, the properties of dynamic and partially disordered (bio)molecular structures can be characterized by integrating a wide range of experimental data, including measurements of ensemble-averaged observables. We start from a Bayesian formulation in which the posterior is a functional that ranks different configuration space distributions. By maximizing this posterior, we derive an optimal Bayesian ensemble distribution. For discrete configurations, this optimal distribution is identical to that obtained by the maximum entropy "ensemble refinement of SAXS" (EROS) formulation. Bayesian replica ensemble refinement enhances the sampling of relevant configurations by imposing restraints on averages of observables in coupled replica molecular dynamics simulations. We find that the strength of the restraint scales with the number of replicas and we show that this sca...
Bayesian classification and regression trees for predicting incidence of cryptosporidiosis.
Directory of Open Access Journals (Sweden)
Wenbiao Hu
Full Text Available BACKGROUND: Classification and regression tree (CART models are tree-based exploratory data analysis methods which have been shown to be very useful in identifying and estimating complex hierarchical relationships in ecological and medical contexts. In this paper, a Bayesian CART model is described and applied to the problem of modelling the cryptosporidiosis infection in Queensland, Australia. METHODOLOGY/PRINCIPAL FINDINGS: We compared the results of a Bayesian CART model with those obtained using a Bayesian spatial conditional autoregressive (CAR model. Overall, the analyses indicated that the nature and magnitude of the effect estimates were similar for the two methods in this study, but the CART model more easily accommodated higher order interaction effects. CONCLUSIONS/SIGNIFICANCE: A Bayesian CART model for identification and estimation of the spatial distribution of disease risk is useful in monitoring and assessment of infectious diseases prevention and control.
Exploiting correlation and budget constraints in Bayesian multi-armed bandit optimization
Hoffman, Matthew W.; Shahriari, Bobak; De Freitas, Nando
2013-01-01
We address the problem of finding the maximizer of a nonlinear smooth function, that can only be evaluated point-wise, subject to constraints on the number of permitted function evaluations. This problem is also known as fixed-budget best arm identification in the multi-armed bandit literature. We introduce a Bayesian approach for this problem and show that it empirically outperforms both the existing frequentist counterpart and other Bayesian optimization methods. The Bayesian approach place...
Preliminary investigation of a Bayesian network for mammographic diagnosis of breast cancer.
Kahn, C. E.; Roberts, L. M.; K. Wang; Jenks, D.; Haddawy, P.
1995-01-01
Bayesian networks use the techniques of probability theory to reason under conditions of uncertainty. We investigated the use of Bayesian networks for radiological decision support. A Bayesian network for the interpretation of mammograms (MammoNet) was developed based on five patient-history features, two physical findings, and 15 mammographic features extracted by experienced radiologists. Conditional-probability data, such as sensitivity and specificity, were derived from peer-reviewed jour...
Bayesian network learning for natural hazard assessments
Vogel, Kristin
2016-04-01
Even though quite different in occurrence and consequences, from a modelling perspective many natural hazards share similar properties and challenges. Their complex nature as well as lacking knowledge about their driving forces and potential effects make their analysis demanding. On top of the uncertainty about the modelling framework, inaccurate or incomplete event observations and the intrinsic randomness of the natural phenomenon add up to different interacting layers of uncertainty, which require a careful handling. Thus, for reliable natural hazard assessments it is crucial not only to capture and quantify involved uncertainties, but also to express and communicate uncertainties in an intuitive way. Decision-makers, who often find it difficult to deal with uncertainties, might otherwise return to familiar (mostly deterministic) proceedings. In the scope of the DFG research training group „NatRiskChange" we apply the probabilistic framework of Bayesian networks for diverse natural hazard and vulnerability studies. The great potential of Bayesian networks was already shown in previous natural hazard assessments. Treating each model component as random variable, Bayesian networks aim at capturing the joint distribution of all considered variables. Hence, each conditional distribution of interest (e.g. the effect of precautionary measures on damage reduction) can be inferred. The (in-)dependencies between the considered variables can be learned purely data driven or be given by experts. Even a combination of both is possible. By translating the (in-)dependences into a graph structure, Bayesian networks provide direct insights into the workings of the system and allow to learn about the underlying processes. Besides numerous studies on the topic, learning Bayesian networks from real-world data remains challenging. In previous studies, e.g. on earthquake induced ground motion and flood damage assessments, we tackled the problems arising with continuous variables
Bayesian nonparametric data analysis
Müller, Peter; Jara, Alejandro; Hanson, Tim
2015-01-01
This book reviews nonparametric Bayesian methods and models that have proven useful in the context of data analysis. Rather than providing an encyclopedic review of probability models, the book’s structure follows a data analysis perspective. As such, the chapters are organized by traditional data analysis problems. In selecting specific nonparametric models, simpler and more traditional models are favored over specialized ones. The discussed methods are illustrated with a wealth of examples, including applications ranging from stylized examples to case studies from recent literature. The book also includes an extensive discussion of computational methods and details on their implementation. R code for many examples is included in on-line software pages.
Decentralized Distributed Bayesian Estimation
Czech Academy of Sciences Publication Activity Database
Dedecius, Kamil; Sečkárová, Vladimíra
Praha: ÚTIA AVČR, v.v.i, 2011 - (Janžura, M.; Ivánek, J.). s. 16-16 [7th International Workshop on Data–Algorithms–Decision Making. 27.11.2011-29.11.2011, Mariánská] R&D Projects: GA ČR 102/08/0567; GA ČR GA102/08/0567 Institutional research plan: CEZ:AV0Z10750506 Keywords : estimation * distributed estimation * model Subject RIV: BB - Applied Statistics, Operational Research http://library.utia.cas.cz/separaty/2011/AS/dedecius-decentralized distributed bayesian estimation.pdf
Congdon, Peter
2014-01-01
This book provides an accessible approach to Bayesian computing and data analysis, with an emphasis on the interpretation of real data sets. Following in the tradition of the successful first edition, this book aims to make a wide range of statistical modeling applications accessible using tested code that can be readily adapted to the reader's own applications. The second edition has been thoroughly reworked and updated to take account of advances in the field. A new set of worked examples is included. The novel aspect of the first edition was the coverage of statistical modeling using WinBU
Computationally efficient Bayesian tracking
Aughenbaugh, Jason; La Cour, Brian
2012-06-01
In this paper, we describe the progress we have achieved in developing a computationally efficient, grid-based Bayesian fusion tracking system. In our approach, the probability surface is represented by a collection of multidimensional polynomials, each computed adaptively on a grid of cells representing state space. Time evolution is performed using a hybrid particle/grid approach and knowledge of the grid structure, while sensor updates use a measurement-based sampling method with a Delaunay triangulation. We present an application of this system to the problem of tracking a submarine target using a field of active and passive sonar buoys.
Improved iterative Bayesian unfolding
D'Agostini, G
2010-01-01
This paper reviews the basic ideas behind a Bayesian unfolding published some years ago and improves their implementation. In particular, uncertainties are now treated at all levels by probability density functions and their propagation is performed by Monte Carlo integration. Thus, small numbers are better handled and the final uncertainty does not rely on the assumption of normality. Theoretical and practical issues concerning the iterative use of the algorithm are also discussed. The new program, implemented in the R language, is freely available, together with sample scripts to play with toy models.
The NIFTY way of Bayesian signal inference
International Nuclear Information System (INIS)
We introduce NIFTY, 'Numerical Information Field Theory', a software package for the development of Bayesian signal inference algorithms that operate independently from any underlying spatial grid and its resolution. A large number of Bayesian and Maximum Entropy methods for 1D signal reconstruction, 2D imaging, as well as 3D tomography, appear formally similar, but one often finds individualized implementations that are neither flexible nor easily transferable. Signal inference in the framework of NIFTY can be done in an abstract way, such that algorithms, prototyped in 1D, can be applied to real world problems in higher-dimensional settings. NIFTY as a versatile library is applicable and already has been applied in 1D, 2D, 3D and spherical settings. A recent application is the D3PO algorithm targeting the non-trivial task of denoising, deconvolving, and decomposing photon observations in high energy astronomy
QBism, the Perimeter of Quantum Bayesianism
Fuchs, Christopher A
2010-01-01
This article summarizes the Quantum Bayesian point of view of quantum mechanics, with special emphasis on the view's outer edges---dubbed QBism. QBism has its roots in personalist Bayesian probability theory, is crucially dependent upon the tools of quantum information theory, and most recently, has set out to investigate whether the physical world might be of a type sketched by some false-started philosophies of 100 years ago (pragmatism, pluralism, nonreductionism, and meliorism). Beyond conceptual issues, work at Perimeter Institute is focused on the hard technical problem of finding a good representation of quantum mechanics purely in terms of probabilities, without amplitudes or Hilbert-space operators. The best candidate representation involves a mysterious entity called a symmetric informationally complete quantum measurement. Contemplation of it gives a way of thinking of the Born Rule as an addition to the rules of probability theory, applicable when an agent considers gambling on the consequences of...
Bayesian Magnetohydrodynamic Seismology of Coronal Loops
Arregui, Inigo
2011-01-01
We perform a Bayesian parameter inference in the context of resonantly damped transverse coronal loop oscillations. The forward problem is solved in terms of parametric results for kink waves in one-dimensional flux tubes in the thin tube and thin boundary approximations. For the inverse problem, we adopt a Bayesian approach to infer the most probable values of the relevant parameters, for given observed periods and damping times, and to extract their confidence levels. The posterior probability distribution functions are obtained by means of Markov Chain Monte Carlo simulations, incorporating observed uncertainties in a consistent manner. We find well localized solutions in the posterior probability distribution functions for two of the three parameters of interest, namely the Alfven travel time and the transverse inhomogeneity length-scale. The obtained estimates for the Alfven travel time are consistent with previous inversion results, but the method enables us to additionally constrain the transverse inho...
Distributed Detection via Bayesian Updates and Consensus
Liu, Qipeng; Wang, Xiaofan
2014-01-01
In this paper, we discuss a class of distributed detection algorithms which can be viewed as implementations of Bayes' law in distributed settings. Some of the algorithms are proposed in the literature most recently, and others are first developed in this paper. The common feature of these algorithms is that they all combine (i) certain kinds of consensus protocols with (ii) Bayesian updates. They are different mainly in the aspect of the type of consensus protocol and the order of the two operations. After discussing their similarities and differences, we compare these distributed algorithms by numerical examples. We focus on the rate at which these algorithms detect the underlying true state of an object. We find that (a) The algorithms with consensus via geometric average is more efficient than that via arithmetic average; (b) The order of consensus aggregation and Bayesian update does not apparently influence the performance of the algorithms; (c) The existence of communication delay dramatically slows do...
Bayesian Inference on Gravitational Waves
Directory of Open Access Journals (Sweden)
Asad Ali
2015-12-01
Full Text Available The Bayesian approach is increasingly becoming popular among the astrophysics data analysis communities. However, the Pakistan statistics communities are unaware of this fertile interaction between the two disciplines. Bayesian methods have been in use to address astronomical problems since the very birth of the Bayes probability in eighteenth century. Today the Bayesian methods for the detection and parameter estimation of gravitational waves have solid theoretical grounds with a strong promise for the realistic applications. This article aims to introduce the Pakistan statistics communities to the applications of Bayesian Monte Carlo methods in the analysis of gravitational wave data with an overview of the Bayesian signal detection and estimation methods and demonstration by a couple of simplified examples.
Adaptive Dynamic Bayesian Networks
Energy Technology Data Exchange (ETDEWEB)
Ng, B M
2007-10-26
A discrete-time Markov process can be compactly modeled as a dynamic Bayesian network (DBN)--a graphical model with nodes representing random variables and directed edges indicating causality between variables. Each node has a probability distribution, conditional on the variables represented by the parent nodes. A DBN's graphical structure encodes fixed conditional dependencies between variables. But in real-world systems, conditional dependencies between variables may be unknown a priori or may vary over time. Model errors can result if the DBN fails to capture all possible interactions between variables. Thus, we explore the representational framework of adaptive DBNs, whose structure and parameters can change from one time step to the next: a distribution's parameters and its set of conditional variables are dynamic. This work builds on recent work in nonparametric Bayesian modeling, such as hierarchical Dirichlet processes, infinite-state hidden Markov networks and structured priors for Bayes net learning. In this paper, we will explain the motivation for our interest in adaptive DBNs, show how popular nonparametric methods are combined to formulate the foundations for adaptive DBNs, and present preliminary results.
Bayesian analysis toolkit - BAT
International Nuclear Information System (INIS)
Statistical treatment of data is an essential part of any data analysis and interpretation. Different statistical methods and approaches can be used, however the implementation of these approaches is complicated and at times inefficient. The Bayesian analysis toolkit (BAT) is a software package developed in C++ framework that facilitates the statistical analysis of the data using Bayesian theorem. The tool evaluates the posterior probability distributions for models and their parameters using Markov Chain Monte Carlo which in turn provide straightforward parameter estimation, limit setting and uncertainty propagation. Additional algorithms, such as simulated annealing, allow extraction of the global mode of the posterior. BAT sets a well-tested environment for flexible model definition and also includes a set of predefined models for standard statistical problems. The package is interfaced to other software packages commonly used in high energy physics, such as ROOT, Minuit, RooStats and CUBA. We present a general overview of BAT and its algorithms. A few physics examples are shown to introduce the spectrum of its applications. In addition, new developments and features are summarized.
Exon-primed intron-crossing (EPIC markers as a tool for ant phylogeography
Directory of Open Access Journals (Sweden)
Patrícia R. Ströher
2013-12-01
Full Text Available Exon-primed intron-crossing (EPIC markers as a tool for ant phylogeography. Due to their local abundance, diversity of adaptations and worldwide distribution, ants are a classic example of adaptive radiation. Despite this evolutionary and ecological importance, phylogeographical studies on ants have relied largely on mitochondrial markers. In this study we design and test exon-primed intron-crossing (EPIC markers, which can be widely used to uncover ant intraspecific variation. Candidate markers were obtained through screening the available ant genomes for unlinked conserved exonic regions interspersed with introns. A subset of 15 markers was tested in vitro and showed successful amplification in several phylogenetically distant ant species. These markers represent an important step forward in ant phylogeography and population genetics, allowing for more extensive characterization of variation in ant nuclear DNA without the need to develop species-specific markers.
Czech Academy of Sciences Publication Activity Database
Mazoch, V.; Bryja, Josef; Patzenhauerová, Hana; Šumbera, R.
Brno: Ústav biologie obratlovců AV ČR, 2012 - (Bryja, J.; Albrechtová, J.; Tkadlec, E.). s. 131-132 ISBN 978-80-87189-11-5. [Zoologické dny. 09.02.2012-10.02.2012, Olomouc] R&D Projects: GA ČR GAP506/10/0983 Institutional research plan: CEZ:AV0Z60930519 Keywords : phylogeography * rodents * Africa Subject RIV: EG - Zoology
Czech Academy of Sciences Publication Activity Database
Kodandaramaiah, U.; Konvička, Martin; Tammaru, T.; Wahlberg, N.; Gotthard, K.
2012-01-01
Roč. 16, č. 2 (2012), s. 305-313. ISSN 1366-638X R&D Projects: GA MŠk LC06073; GA ČR GAP505/10/2167 Institutional support: RVO:60077344 Keywords : Lopinga achine * phylogeography * conservation Subject RIV: EH - Ecology, Behaviour Impact factor: 1.801, year: 2012 http://link.springer.com/article/10.1007/s10841-012-9465-4?null
Directory of Open Access Journals (Sweden)
Hung Kuo-Hsiang
2011-01-01
Full Text Available Abstract Background Tetraena mongolica (Zygophyllaceae, an endangered endemic species in western Inner Mongolia, China. For endemic species with a limited geographical range and declining populations, historical patterns of demography and hierarchical genetic structure are important for determining population structure, and also provide information for developing effective and sustainable management plans. In this study, we assess genetic variation, population structure, and phylogeography of T. mongolica from eight populations. Furthermore, we evaluate the conservation and management units to provide the information for conservation. Results Sequence variation and spatial apportionment of the atpB-rbcL noncoding spacer region of the chloroplast DNA were used to reconstruct the phylogeography of T. mongolica. A total of 880 bp was sequenced from eight extant populations throughout the whole range of its distribution. At the cpDNA locus, high levels of genetic differentiation among populations and low levels of genetic variation within populations were detected, indicating that most seed dispersal was restricted within populations. Conclusions Demographic fluctuations, which led to random losses of genetic polymorphisms from populations, due to frequent flooding of the Yellow River and human disturbance were indicated by the analysis of BEAST skyline plot. Nested clade analysis revealed that restricted gene flow with isolation by distance plus occasional long distance dispersal is the main evolutionary factor affecting the phylogeography and population structure of T. mongolica. For setting a conservation management plan, each population of T. mongolica should be recognized as a conservation unit.
Phylogeography of the inshore fish, Bostrychus sinensis, along the Pacific coastline of China.
Qiu, Fan; Li, Hai; Lin, Hungdu; Ding, Shaoxiong; Miyamoto, Michael M
2016-03-01
This study assesses the phylogeography of the Chinese four-eyed sleeper (Bostrychus sinensis) with one mitochondrial and one nuclear genes and two morphological characters. Phylogenetic and population genetic analyses of the sequences reveals two phylogeographic lineages from the East and South China Seas, which are corroborated by the morphological data. The vicariance of the two lineages is attributed to the Pleistocene Ice Age exposure of the Taiwan Strait and consequent connection of Taiwan to the mainland, which thereby introduced an ecological barrier to gene flow between populations in the East and South China Seas. The distributions of the two lineages now overlap in the East China Sea and this secondary contact is attributed to biased northward migration along the two main currents of the Taiwan Strait following its interglacial re-flooding. In conclusion, this study reinforces the importance of "vicariance, then secondary contact" due to Late Pliocene and Pleistocene sea-level changes to the phylogeography of marine species. Specifically, it corroborates the importance of Pleistocene sea-level changes in the Taiwan Strait to the phylogeography of Chinese inshore species. PMID:26732489
Book review: Bayesian analysis for population ecology
Link, William A.
2011-01-01
Brian Dennis described the field of ecology as “fertile, uncolonized ground for Bayesian ideas.” He continued: “The Bayesian propagule has arrived at the shore. Ecologists need to think long and hard about the consequences of a Bayesian ecology. The Bayesian outlook is a successful competitor, but is it a weed? I think so.” (Dennis 2004)
DEFF Research Database (Denmark)
Hartelius, Karsten; Carstensen, Jens Michael
2003-01-01
A method for locating distorted grid structures in images is presented. The method is based on the theories of template matching and Bayesian image restoration. The grid is modeled as a deformable template. Prior knowledge of the grid is described through a Markov random field (MRF) model which...... represents the spatial coordinates of the grid nodes. Knowledge of how grid nodes are depicted in the observed image is described through the observation model. The prior consists of a node prior and an arc (edge) prior, both modeled as Gaussian MRFs. The node prior models variations in the positions of grid...... nodes and the arc prior models variations in row and column spacing across the grid. Grid matching is done by placing an initial rough grid over the image and applying an ensemble annealing scheme to maximize the posterior distribution of the grid. The method can be applied to noisy images with missing...
Phylogeography of Francisella tularensis subspecies holarctica from the country of Georgia
Directory of Open Access Journals (Sweden)
Busch Joseph D
2011-06-01
lineages previously described from the B.Br.013 group. This finding suggests that additional phylogenetic studies of F. tularensis subsp. holarctica populations in Eastern Europe and Asia have the potential to yield important new insights into the evolutionary history and phylogeography of this broadly dispersed F. tularensis subspecies.
Current trends in Bayesian methodology with applications
Upadhyay, Satyanshu K; Dey, Dipak K; Loganathan, Appaia
2015-01-01
Collecting Bayesian material scattered throughout the literature, Current Trends in Bayesian Methodology with Applications examines the latest methodological and applied aspects of Bayesian statistics. The book covers biostatistics, econometrics, reliability and risk analysis, spatial statistics, image analysis, shape analysis, Bayesian computation, clustering, uncertainty assessment, high-energy astrophysics, neural networking, fuzzy information, objective Bayesian methodologies, empirical Bayes methods, small area estimation, and many more topics.Each chapter is self-contained and focuses on
Bayesian analysis of cosmic structures
Kitaura, Francisco-Shu
2011-01-01
We revise the Bayesian inference steps required to analyse the cosmological large-scale structure. Here we make special emphasis in the complications which arise due to the non-Gaussian character of the galaxy and matter distribution. In particular we investigate the advantages and limitations of the Poisson-lognormal model and discuss how to extend this work. With the lognormal prior using the Hamiltonian sampling technique and on scales of about 4 h^{-1} Mpc we find that the over-dense regions are excellent reconstructed, however, under-dense regions (void statistics) are quantitatively poorly recovered. Contrary to the maximum a posteriori (MAP) solution which was shown to over-estimate the density in the under-dense regions we obtain lower densities than in N-body simulations. This is due to the fact that the MAP solution is conservative whereas the full posterior yields samples which are consistent with the prior statistics. The lognormal prior is not able to capture the full non-linear regime at scales ...
Applications of Bayesian spectrum representation in acoustics
Botts, Jonathan M.
framework. The application to reflection data is useful for representing frequency-dependent impedance boundaries in finite difference acoustic simulations. Furthermore, since the filter transfer function is a parametric model, it can be modified to incorporate arbitrary frequency weighting and account for the band-limited nature of measured reflection spectra. Finally, the model is modified to compensate for dispersive error in the finite difference simulation, from the filter design process. Stemming from the filter boundary problem, the implementation of pressure sources in finite difference simulation is addressed in order to assure that schemes properly converge. A class of parameterized source functions is proposed and shown to offer straightforward control of residual error in the simulation. Guided by the notion that the solution to be approximated affects the approximation error, sources are designed which reduce residual dispersive error to the size of round-off errors. The early part of a room impulse response can be characterized by a series of isolated plane waves. Measured with an array of microphones, plane waves map to a directional response of the array or spatial intensity map. Probabilistic inversion of this response results in estimates of the number and directions of image source arrivals. The model-based inversion is shown to avoid ambiguities associated with peak-finding or inspection of the spatial intensity map. For this problem, determining the number of arrivals in a given frame is critical for properly inferring the state of the sound field. This analysis is effectively compression of the spatial room response, which is useful for analysis or encoding of the spatial sound field. Parametric, model-based formulations of these problems enhance the solution in all cases, and a Bayesian interpretation provides a principled approach to model comparison and parameter estimation. v
Phylogeography of foot-and-mouth disease virus serotype O in Ecuador.
de Carvalho, Luiz Max Fagundes; Santos, Leonardo Bacelar Lima; Faria, Nuno Rodrigues; de Castro Silveira, Waldemir
2013-01-01
Foot-and-mouth disease virus (FMDV) is the causative agent of the most important disease of domestic cattle, foot-and-mouth disease. In Ecuador, FMDV is maintained at an endemic state, with sporadic outbreaks. To unravel the tempo and mode of FMDV spread within the country we conducted a Bayesian phylogeographic analysis using a continuous time Markov chain (CTMC) to model the diffusion of FMDV between Ecuadorian provinces. We implement this framework through Markov chain Monte Carlo available in the BEAST package to study 90 FMDV serotype O isolates from 17 provinces in the period 2002-2010. The Bayesian approach also allowed us to test hypotheses on the mechanisms of viral spread by incorporating environmental and epidemiological data in our prior distributions and perform Bayesian model selection. Our analyses suggest an intense flow of viral strains throughout the country that is possibly coupled to animal movements and ecological factors, since most of inferred viral spread routes were in Coast and Highland regions. Moreover, our results suggest that both short- and long-range spread occur within Ecuador. The province of Esmeraldas, in the border with Colombia and where most animal commerce is done, was found to be the most probable origin of the circulating strains, pointing to a transboundary behavior of FMDV in South America. These findings suggest that uncontrolled animal movements can create a favorable environment for FMDV maintenance and pose a serious threat to control programmes. Also, we show that phylogeographic modeling can be a powerful tool in unraveling the spatial dynamics of viruses and potentially in controlling the spread of these pathogens. PMID:22985683
Bayesian Methods and Universal Darwinism
Campbell, John
2010-01-01
Bayesian methods since the time of Laplace have been understood by their practitioners as closely aligned to the scientific method. Indeed a recent champion of Bayesian methods, E. T. Jaynes, titled his textbook on the subject Probability Theory: the Logic of Science. Many philosophers of science including Karl Popper and Donald Campbell have interpreted the evolution of Science as a Darwinian process consisting of a 'copy with selective retention' algorithm abstracted from Darwin's theory of...
Portfolio Allocation for Bayesian Optimization
Brochu, Eric; Hoffman, Matthew W.; De Freitas, Nando
2010-01-01
Bayesian optimization with Gaussian processes has become an increasingly popular tool in the machine learning community. It is efficient and can be used when very little is known about the objective function, making it popular in expensive black-box optimization scenarios. It uses Bayesian methods to sample the objective efficiently using an acquisition function which incorporates the model's estimate of the objective and the uncertainty at any given point. However, there are several differen...
Neuronanatomy, neurology and Bayesian networks
Bielza Lozoya, Maria Concepcion
2014-01-01
Bayesian networks are data mining models with clear semantics and a sound theoretical foundation. In this keynote talk we will pinpoint a number of neuroscience problems that can be addressed using Bayesian networks. In neuroanatomy, we will show computer simulation models of dendritic trees and classification of neuron types, both based on morphological features. In neurology, we will present the search for genetic biomarkers in Alzheimer's disease and the prediction of health-related qualit...
Bayesian Networks and Influence Diagrams
DEFF Research Database (Denmark)
Kjærulff, Uffe Bro; Madsen, Anders Læsø
Probabilistic networks, also known as Bayesian networks and influence diagrams, have become one of the most promising technologies in the area of applied artificial intelligence, offering intuitive, efficient, and reliable methods for diagnosis, prediction, decision making, classification......, troubleshooting, and data mining under uncertainty. Bayesian Networks and Influence Diagrams: A Guide to Construction and Analysis provides a comprehensive guide for practitioners who wish to understand, construct, and analyze intelligent systems for decision support based on probabilistic networks. Intended...
Bayesian logistic betting strategy against probability forecasting
Kumon, Masayuki; Takemura, Akimichi; Takeuchi, Kei
2012-01-01
We propose a betting strategy based on Bayesian logistic regression modeling for the probability forecasting game in the framework of game-theoretic probability by Shafer and Vovk (2001). We prove some results concerning the strong law of large numbers in the probability forecasting game with side information based on our strategy. We also apply our strategy for assessing the quality of probability forecasting by the Japan Meteorological Agency. We find that our strategy beats the agency by exploiting its tendency of avoiding clear-cut forecasts.
Recovery of shapes: hypermodels and Bayesian learning
International Nuclear Information System (INIS)
We discuss the problem of recovering an image from its blurred and noisy copy with the additional information that the image consists of simple shapes with sharp edges. An iterative algorithm is given, based on the idea of updating the Tikhonov type smoothness penalty on the basis of the previous estimate. This algorithm is discussed in the framework of Bayesian hypermodels and it is shown that the approach can be justified as a sequential iterative scheme for finding the mode of the posterior density. An effective numerical algorithm based on preconditioned Krylov subspace iterations is suggested and demonstrated with a computed example
Dale Poirier
2008-01-01
This paper provides Bayesian rationalizations for White’s heteroskedastic consistent (HC) covariance estimator and various modifications of it. An informed Bayesian bootstrap provides the statistical framework.
Nonparametric Bayesian Classification
Coram, M A
2002-01-01
A Bayesian approach to the classification problem is proposed in which random partitions play a central role. It is argued that the partitioning approach has the capacity to take advantage of a variety of large-scale spatial structures, if they are present in the unknown regression function $f_0$. An idealized one-dimensional problem is considered in detail. The proposed nonparametric prior uses random split points to partition the unit interval into a random number of pieces. This prior is found to provide a consistent estimate of the regression function in the $\\L^p$ topology, for any $1 \\leq p < \\infty$, and for arbitrary measurable $f_0:[0,1] \\rightarrow [0,1]$. A Markov chain Monte Carlo (MCMC) implementation is outlined and analyzed. Simulation experiments are conducted to show that the proposed estimate compares favorably with a variety of conventional estimators. A striking resemblance between the posterior mean estimate and the bagged CART estimate is noted and discussed. For higher dimensions, a ...
BAT - Bayesian Analysis Toolkit
International Nuclear Information System (INIS)
One of the most vital steps in any data analysis is the statistical analysis and comparison with the prediction of a theoretical model. The many uncertainties associated with the theoretical model and the observed data require a robust statistical analysis tool. The Bayesian Analysis Toolkit (BAT) is a powerful statistical analysis software package based on Bayes' Theorem, developed to evaluate the posterior probability distribution for models and their parameters. It implements Markov Chain Monte Carlo to get the full posterior probability distribution that in turn provides a straightforward parameter estimation, limit setting and uncertainty propagation. Additional algorithms, such as Simulated Annealing, allow to evaluate the global mode of the posterior. BAT is developed in C++ and allows for a flexible definition of models. A set of predefined models covering standard statistical cases are also included in BAT. It has been interfaced to other commonly used software packages such as ROOT, Minuit, RooStats and CUBA. An overview of the software and its algorithms is provided along with several physics examples to cover a range of applications of this statistical tool. Future plans, new features and recent developments are briefly discussed.
Exploiting sensitivity analysis in Bayesian networks for consumer satisfaction study
Jaronski, W.; Bloemer, J.M.M.; Vanhoof, K.; Wets, G.
2004-01-01
The paper presents an application of Bayesian network technology in a empirical customer satisfaction study. The findings of the study should provide insight as to the importance of product/service dimensions in terms of the strength of their influence on overall satisfaction. To this end we apply a
The bootstrap and Bayesian bootstrap method in assessing bioequivalence
International Nuclear Information System (INIS)
Parametric method for assessing individual bioequivalence (IBE) may concentrate on the hypothesis that the PK responses are normal. Nonparametric method for evaluating IBE would be bootstrap method. In 2001, the United States Food and Drug Administration (FDA) proposed a draft guidance. The purpose of this article is to evaluate the IBE between test drug and reference drug by bootstrap and Bayesian bootstrap method. We study the power of bootstrap test procedures and the parametric test procedures in FDA (2001). We find that the Bayesian bootstrap method is the most excellent.
Bayesian seismic AVO inversion
Energy Technology Data Exchange (ETDEWEB)
Buland, Arild
2002-07-01
A new linearized AVO inversion technique is developed in a Bayesian framework. The objective is to obtain posterior distributions for P-wave velocity, S-wave velocity and density. Distributions for other elastic parameters can also be assessed, for example acoustic impedance, shear impedance and P-wave to S-wave velocity ratio. The inversion algorithm is based on the convolutional model and a linearized weak contrast approximation of the Zoeppritz equation. The solution is represented by a Gaussian posterior distribution with explicit expressions for the posterior expectation and covariance, hence exact prediction intervals for the inverted parameters can be computed under the specified model. The explicit analytical form of the posterior distribution provides a computationally fast inversion method. Tests on synthetic data show that all inverted parameters were almost perfectly retrieved when the noise approached zero. With realistic noise levels, acoustic impedance was the best determined parameter, while the inversion provided practically no information about the density. The inversion algorithm has also been tested on a real 3-D dataset from the Sleipner Field. The results show good agreement with well logs but the uncertainty is high. The stochastic model includes uncertainties of both the elastic parameters, the wavelet and the seismic and well log data. The posterior distribution is explored by Markov chain Monte Carlo simulation using the Gibbs sampler algorithm. The inversion algorithm has been tested on a seismic line from the Heidrun Field with two wells located on the line. The uncertainty of the estimated wavelet is low. In the Heidrun examples the effect of including uncertainty of the wavelet and the noise level was marginal with respect to the AVO inversion results. We have developed a 3-D linearized AVO inversion method with spatially coupled model parameters where the objective is to obtain posterior distributions for P-wave velocity, S
Brochu, Eric; Cora, Vlad M.; De Freitas, Nando
2010-01-01
We present a tutorial on Bayesian optimization, a method of finding the maximum of expensive cost functions. Bayesian optimization employs the Bayesian technique of setting a prior over the objective function and combining it with evidence to get a posterior function. This permits a utility-based selection of the next observation to make on the objective function, which must take into account both exploration (sampling from areas of high uncertainty) and exploitation (sampling areas likely to...
Czech Academy of Sciences Publication Activity Database
Hulva, P.; Fornůsková, Alena; Chudárková, A.; Evin, A.; Allegrini, B.; Benda, P.; Bryja, Josef
2010-01-01
Roč. 19, č. 24 (2010), s. 5417-5431. ISSN 0962-1083 R&D Projects: GA MŠk LC06073 Institutional research plan: CEZ:AV0Z60930519 Keywords : hybrid speciation * introgression * Mediterranean * microsatellites * mitochondrial DNA * phylogeography * bats * radiation Subject RIV: EG - Zoology Impact factor: 6.457, year: 2010
Czech Academy of Sciences Publication Activity Database
Hahn, M.W.; Koll, U.; Jezberová, Jitka; Camacho, A.
2015-01-01
Roč. 17, č. 3 (2015), s. 829-840. ISSN 1462-2912 R&D Projects: GA ČR(CZ) GEEEF/10/E011 Institutional support: RVO:60077344 Keywords : polynucleobacter * phylogeography * microbiology * bacteria Subject RIV: EE - Microbiology, Virology Impact factor: 6.201, year: 2014
Lemmon, Alan R; Lemmon, Emily Moriarty
2012-10-01
One of the major challenges for researchers studying phylogeography and shallow-scale phylogenetics is the identification of highly variable and informative nuclear loci for the question of interest. Previous approaches to locus identification have generally required extensive testing of anonymous nuclear loci developed from genomic libraries of the target taxon, testing of loci of unknown utility from other systems, or identification of loci from the nearest model organism with genomic resources. Here, we present a fast and economical approach to generating thousands of variable, single-copy nuclear loci for any system using next-generation sequencing. We performed Illumina paired-end sequencing of three reduced-representation libraries (RRLs) in chorus frogs (Pseudacris) to identify orthologous, single-copy loci across libraries and to estimate sequence divergence at multiple taxonomic levels. We also conducted PCR testing of these loci across the genus Pseudacris and outgroups to determine whether loci developed for phylogeography can be extended to deeper phylogenetic levels. Prior to sequencing, we conducted in silico digestion of the most closely related reference genome (Xenopus tropicalis) to generate expectations for the number of loci and degree of coverage for a particular experimental design. Using the RRL approach, we: (i) identified more than 100,000 single-copy nuclear loci, 6339 of which were obtained for divergent conspecifics and 904 of which were obtained for heterospecifics; (ii) estimated average nuclear sequence divergence at 0.1% between alleles within an individual, 1.1% between conspecific individuals that represent two different clades, and 1.8% between species; and (iii) determined from PCR testing that 53% of the loci successfully amplify within-species and also many amplify to the genus-level and deeper in the phylogeny (16%). Our study effectively identified nuclear loci present in the genome that have levels of sequence divergence on
DEFF Research Database (Denmark)
Muñoz, Joaquin; Gómez, Africa; Green, Andy J.;
2008-01-01
There has been a recent appreciation of the ecological impacts of zooplanktonic species invasions. The North American brine shrimp Artemia franciscana is one such alien invader in hyper-saline water ecosystems at a global scale. It has been shown to outcompete native Artemia species, leading to...... their local extinction. We used partial sequences of the mitochondrial Cytochrome c Oxidase Subunit 1 (COI or cox1) gene to investigate the genetic diversity and phylogeography of A. salina, an extreme halophilic sexual brine shrimp, over its known distribution range (Mediterranean Basin and South...
Bayesian modeling using WinBUGS
Ntzoufras, Ioannis
2009-01-01
A hands-on introduction to the principles of Bayesian modeling using WinBUGS Bayesian Modeling Using WinBUGS provides an easily accessible introduction to the use of WinBUGS programming techniques in a variety of Bayesian modeling settings. The author provides an accessible treatment of the topic, offering readers a smooth introduction to the principles of Bayesian modeling with detailed guidance on the practical implementation of key principles. The book begins with a basic introduction to Bayesian inference and the WinBUGS software and goes on to cover key topics, including: Markov Chain Monte Carlo algorithms in Bayesian inference Generalized linear models Bayesian hierarchical models Predictive distribution and model checking Bayesian model and variable evaluation Computational notes and screen captures illustrate the use of both WinBUGS as well as R software to apply the discussed techniques. Exercises at the end of each chapter allow readers to test their understanding of the presented concepts and all ...
Risks Analysis of Logistics Financial Business Based on Evidential Bayesian Network
Bin Suo; Ying Yan
2013-01-01
Risks in logistics financial business are identified and classified. Making the failure of the business as the root node, a Bayesian network is constructed to measure the risk levels in the business. Three importance indexes are calculated to find the most important risks in the business. And more, considering the epistemic uncertainties in the risks, evidence theory associate with Bayesian network is used as an evidential network in the risk analysis of logistics finance. To find how much un...
Probability biases as Bayesian inference
Directory of Open Access Journals (Sweden)
Andre; C. R. Martins
2006-11-01
Full Text Available In this article, I will show how several observed biases in human probabilistic reasoning can be partially explained as good heuristics for making inferences in an environment where probabilities have uncertainties associated to them. Previous results show that the weight functions and the observed violations of coalescing and stochastic dominance can be understood from a Bayesian point of view. We will review those results and see that Bayesian methods should also be used as part of the explanation behind other known biases. That means that, although the observed errors are still errors under the be understood as adaptations to the solution of real life problems. Heuristics that allow fast evaluations and mimic a Bayesian inference would be an evolutionary advantage, since they would give us an efficient way of making decisions. %XX In that sense, it should be no surprise that humans reason with % probability as it has been observed.
Bayesian Methods and Universal Darwinism
Campbell, John
2010-01-01
Bayesian methods since the time of Laplace have been understood by their practitioners as closely aligned to the scientific method. Indeed a recent champion of Bayesian methods, E. T. Jaynes, titled his textbook on the subject Probability Theory: the Logic of Science. Many philosophers of science including Karl Popper and Donald Campbell have interpreted the evolution of Science as a Darwinian process consisting of a 'copy with selective retention' algorithm abstracted from Darwin's theory of Natural Selection. Arguments are presented for an isomorphism between Bayesian Methods and Darwinian processes. Universal Darwinism, as the term has been developed by Richard Dawkins, Daniel Dennett and Susan Blackmore, is the collection of scientific theories which explain the creation and evolution of their subject matter as due to the operation of Darwinian processes. These subject matters span the fields of atomic physics, chemistry, biology and the social sciences. The principle of Maximum Entropy states that system...
Phylogeography, hybridization and Pleistocene refugia of the kob antelope (Kobus kob)
DEFF Research Database (Denmark)
Lorenzen, Eline Deidre; De Neergaard, Rikke; Arctander, Peter; Siegismund, Hans R
2007-01-01
Mitochondrial DNA control region sequences and seven microsatellites were used to estimate the genetic structuring, evolutionary history and historic migration patterns of the kob antelope (Kobus kob). Ten populations were analysed, representing the three recognized K. kob subspecies: K. k. kob in...... distinct K. k. leucotis populations in Sudan and Ethiopia. This was regardless of marker type. Pairwise comparisons and genetic distances between populations grouped Murchison with K. k. leucotis, as did the Bayesian analysis, which failed to find any genetic structuring within the group. We propose that...... the divergent phenotype and life-history adaptations of K. k. leucotis reflect the isolation of kob populations in refugia in west and east Africa during the Pleistocene. Subsequent dispersal has led to secondary contact and hybridization in northern Uganda between lineages, which was supported by...
Bayesian methods for proteomic biomarker development
Directory of Open Access Journals (Sweden)
Belinda Hernández
2015-12-01
In this review we provide an introduction to Bayesian inference and demonstrate some of the advantages of using a Bayesian framework. We summarize how Bayesian methods have been used previously in proteomics and other areas of bioinformatics. Finally, we describe some popular and emerging Bayesian models from the statistical literature and provide a worked tutorial including code snippets to show how these methods may be applied for the evaluation of proteomic biomarkers.
Bayesian test and Kuhn's paradigm
Institute of Scientific and Technical Information of China (English)
Chen Xiaoping
2006-01-01
Kuhn's theory of paradigm reveals a pattern of scientific progress,in which normal science alternates with scientific revolution.But Kuhn underrated too much the function of scientific test in his pattern,because he focuses all his attention on the hypothetico-deductive schema instead of Bayesian schema.This paper employs Bayesian schema to re-examine Kuhn's theory of paradigm,to uncover its logical and rational components,and to illustrate the tensional structure of logic and belief,rationality and irrationality,in the process of scientific revolution.
3D Bayesian contextual classifiers
DEFF Research Database (Denmark)
Larsen, Rasmus
2000-01-01
We extend a series of multivariate Bayesian 2-D contextual classifiers to 3-D by specifying a simultaneous Gaussian distribution for the feature vectors as well as a prior distribution of the class variables of a pixel and its 6 nearest 3-D neighbours.......We extend a series of multivariate Bayesian 2-D contextual classifiers to 3-D by specifying a simultaneous Gaussian distribution for the feature vectors as well as a prior distribution of the class variables of a pixel and its 6 nearest 3-D neighbours....
Phylogeography of the current rabies viruses in Indonesia.
Dibia, I Nyoman; Sumiarto, Bambang; Susetya, Heru; Putra, Anak Agung Gde; Scott-Orr, Helen; Mahardika, Gusti Ngurah
2015-01-01
Rabies is a major fatal zoonotic disease in Indonesia. This study was conducted to determine the recent dynamics of rabies virus (RABV) in various areas and animal species throughout Indonesia. A total of 27 brain samples collected from rabid animals of various species in Bali, Sumatra, Kalimantan, Sulawesi, Java, and Flores in 2008 to 2010 were investigated. The cDNA of the nucleoprotein gene from each sample was generated and amplified by one-step reverse transcription-PCR, after which the products were sequenced and analyzed. The symmetric substitution model of a Bayesian stochastic search variable selection extension of the discrete phylogeographic model of the social network was applied in BEAST ver. 1.7.5 software. The spatial dispersal was visualized in Cartographica using Spatial Phylogenetic Reconstruction of Evolutionary Dynamics. We demonstrated inter-island introduction and reintroduction, and dog was found to be the only source of infection of other animals. Ancestors of Indonesian RABVs originated in Java and its descendants were transmitted to Kalimantan, then further to Sumatra, Flores, and Bali. The Flores descendent was subsequently transmitted to Sulawesi and back to Kalimantan. The viruses found in various animal species were transmitted by the dog. PMID:25643792
Bayesian Model Averaging for Propensity Score Analysis
Kaplan, David; Chen, Jianshen
2013-01-01
The purpose of this study is to explore Bayesian model averaging in the propensity score context. Previous research on Bayesian propensity score analysis does not take into account model uncertainty. In this regard, an internally consistent Bayesian framework for model building and estimation must also account for model uncertainty. The…
Bayesian networks and food security - An introduction
Stein, A.
2004-01-01
This paper gives an introduction to Bayesian networks. Networks are defined and put into a Bayesian context. Directed acyclical graphs play a crucial role here. Two simple examples from food security are addressed. Possible uses of Bayesian networks for implementation and further use in decision sup
Bayesian variable order Markov models: Towards Bayesian predictive state representations
C. Dimitrakakis
2009-01-01
We present a Bayesian variable order Markov model that shares many similarities with predictive state representations. The resulting models are compact and much easier to specify and learn than classical predictive state representations. Moreover, we show that they significantly outperform a more st
REMITTANCES, DUTCH DISEASE, AND COMPETITIVENESS: A BAYESIAN ANALYSIS
FARID MAKHLOUF; MAZHAR MUGHAL
2013-01-01
The paper studies symptoms of Dutch disease in the Pakistani economy arising from international remittances. An IV Bayesian analysis is carried out to take care of the endogeneity and uncertainty due to the managed float of Pakistani Rupee. We find evidence for both spending and resource movement effects in both the short and the long-run. These impacts are stronger and different from those the Official Development Assistance and the FDI exert. We find that while aggregate remittances and the...
Bayesian network learning with cutting planes
Cussens, James
2012-01-01
The problem of learning the structure of Bayesian networks from complete discrete data with a limit on parent set size is considered. Learning is cast explicitly as an optimisation problem where the goal is to find a BN structure which maximises log marginal likelihood (BDe score). Integer programming, specifically the SCIP framework, is used to solve this optimisation problem. Acyclicity constraints are added to the integer program (IP) during solving in the form of cutting planes. Finding good cutting planes is the key to the success of the approach -the search for such cutting planes is effected using a sub-IP. Results show that this is a particularly fast method for exact BN learning.
A Bayesian approach to earthquake source studies
Minson, Sarah
Bayesian sampling has several advantages over conventional optimization approaches to solving inverse problems. It produces the distribution of all possible models sampled proportionally to how much each model is consistent with the data and the specified prior information, and thus images the entire solution space, revealing the uncertainties and trade-offs in the model. Bayesian sampling is applicable to both linear and non-linear modeling, and the values of the model parameters being sampled can be constrained based on the physics of the process being studied and do not have to be regularized. However, these methods are computationally challenging for high-dimensional problems. Until now the computational expense of Bayesian sampling has been too great for it to be practicable for most geophysical problems. I present a new parallel sampling algorithm called CATMIP for Cascading Adaptive Tempered Metropolis In Parallel. This technique, based on Transitional Markov chain Monte Carlo, makes it possible to sample distributions in many hundreds of dimensions, if the forward model is fast, or to sample computationally expensive forward models in smaller numbers of dimensions. The design of the algorithm is independent of the model being sampled, so CATMIP can be applied to many areas of research. I use CATMIP to produce a finite fault source model for the 2007 Mw 7.7 Tocopilla, Chile earthquake. Surface displacements from the earthquake were recorded by six interferograms and twelve local high-rate GPS stations. Because of the wealth of near-fault data, the source process is well-constrained. I find that the near-field high-rate GPS data have significant resolving power above and beyond the slip distribution determined from static displacements. The location and magnitude of the maximum displacement are resolved. The rupture almost certainly propagated at sub-shear velocities. The full posterior distribution can be used not only to calculate source parameters but also
Theory-independent limits on correlations from generalized Bayesian networks
International Nuclear Information System (INIS)
Bayesian networks provide a powerful tool for reasoning about probabilistic causation, used in many areas of science. They are, however, intrinsically classical. In particular, Bayesian networks naturally yield the Bell inequalities. Inspired by this connection, we generalize the formalism of classical Bayesian networks in order to investigate non-classical correlations in arbitrary causal structures. Our framework of ‘generalized Bayesian networks’ replaces latent variables with the resources of any generalized probabilistic theory, most importantly quantum theory, but also, for example, Popescu–Rohrlich boxes. We obtain three main sets of results. Firstly, we prove that all of the observable conditional independences required by the classical theory also hold in our generalization; to obtain this, we extend the classical d-separation theorem to our setting. Secondly, we find that the theory-independent constraints on probabilities can go beyond these conditional independences. For example we find that no probabilistic theory predicts perfect correlation between three parties using only bipartite common causes. Finally, we begin a classification of those causal structures, such as the Bell scenario, that may yield a separation between classical, quantum and general-probabilistic correlations. (paper)
Bayesian Analysis of Experimental Data
Directory of Open Access Journals (Sweden)
Lalmohan Bhar
2013-10-01
Full Text Available Analysis of experimental data from Bayesian point of view has been considered. Appropriate methodology has been developed for application into designed experiments. Normal-Gamma distribution has been considered for prior distribution. Developed methodology has been applied to real experimental data taken from long term fertilizer experiments.
Bayesian image restoration, using configurations
DEFF Research Database (Denmark)
Thorarinsdottir, Thordis Linda
2006-01-01
configurations are expressed in terms of the mean normal measure of the random set. These probabilities are used as prior probabilities in a Bayesian image restoration approach. Estimation of the remaining parameters in the model is outlined for the salt and pepper noise. The inference in the model is discussed...
Bayesian image restoration, using configurations
DEFF Research Database (Denmark)
Thorarinsdottir, Thordis
configurations are expressed in terms of the mean normal measure of the random set. These probabilities are used as prior probabilities in a Bayesian image restoration approach. Estimation of the remaining parameters in the model is outlined for salt and pepper noise. The inference in the model is discussed in...
ANALYSIS OF BAYESIAN CLASSIFIER ACCURACY
Directory of Open Access Journals (Sweden)
Felipe Schneider Costa
2013-01-01
Full Text Available The naÃ¯ve Bayes classifier is considered one of the most effective classification algorithms today, competing with more modern and sophisticated classifiers. Despite being based on unrealistic (naÃ¯ve assumption that all variables are independent, given the output class, the classifier provides proper results. However, depending on the scenario utilized (network structure, number of samples or training cases, number of variables, the network may not provide appropriate results. This study uses a process variable selection, using the chi-squared test to verify the existence of dependence between variables in the data model in order to identify the reasons which prevent a Bayesian network to provide good performance. A detailed analysis of the data is also proposed, unlike other existing work, as well as adjustments in case of limit values between two adjacent classes. Furthermore, variable weights are used in the calculation of a posteriori probabilities, calculated with mutual information function. Tests were applied in both a naÃ¯ve Bayesian network and a hierarchical Bayesian network. After testing, a significant reduction in error rate has been observed. The naÃ¯ve Bayesian network presented a drop in error rates from twenty five percent to five percent, considering the initial results of the classification process. In the hierarchical network, there was not only a drop in fifteen percent error rate, but also the final result came to zero.
Bayesian Agglomerative Clustering with Coalescents
Teh, Yee Whye; Daumé III, Hal; Roy, Daniel
2009-01-01
We introduce a new Bayesian model for hierarchical clustering based on a prior over trees called Kingman's coalescent. We develop novel greedy and sequential Monte Carlo inferences which operate in a bottom-up agglomerative fashion. We show experimentally the superiority of our algorithms over others, and demonstrate our approach in document clustering and phylolinguistics.
Bayesian Networks and Influence Diagrams
DEFF Research Database (Denmark)
Kjærulff, Uffe Bro; Madsen, Anders Læsø
Bayesian Networks and Influence Diagrams: A Guide to Construction and Analysis, Second Edition, provides a comprehensive guide for practitioners who wish to understand, construct, and analyze intelligent systems for decision support based on probabilistic networks. This new edition contains six new...
Antunes, Jorge T; Leão, Pedro N; Vasconcelos, Vítor M
2015-01-01
Cylindrospermopsis raciborskii is a cyanobacterial species extensively studied for its toxicity, bloom formation and invasiveness potential, which have consequences to public and environmental health. Its current geographical distribution, spanning different climates, suggests that C. raciborskii has acquired the status of a cosmopolitan species. From phylogeography studies, a tropical origin for this species seems convincing, with different conjectural routes of expansion toward temperate climates. This expansion may be a result of the species physiological plasticity, or of the existence of different ecotypes with distinct environmental requirements. In particular, C. raciborskii is known to tolerate wide temperature and light regimes and presents diverse nutritional strategies. This cyanobacterium is also thought to have benefited from climate change conditions, regarding its invasiveness into temperate climates. Other factors, recently put forward, such as allelopathy, may also be important to its expansion. The effect of C. raciborskii in the invaded communities is still mostly unknown but may strongly disturb species diversity at different trophic levels. In this review we present an up-to-date account of the distribution, phylogeography, ecophysiology, as well some preliminary reports of the impact of C. raciborskii in different organisms. PMID:26042108
Topics in Bayesian statistics and maximum entropy
International Nuclear Information System (INIS)
Notions of Bayesian decision theory and maximum entropy methods are reviewed with particular emphasis on probabilistic inference and Bayesian modeling. The axiomatic approach is considered as the best justification of Bayesian analysis and maximum entropy principle applied in natural sciences. Particular emphasis is put on solving the inverse problem in digital image restoration and Bayesian modeling of neural networks. Further topics addressed briefly include language modeling, neutron scattering, multiuser detection and channel equalization in digital communications, genetic information, and Bayesian court decision-making. (author)
Fully Bayesian Experimental Design for Pharmacokinetic Studies
Directory of Open Access Journals (Sweden)
Elizabeth G. Ryan
2015-03-01
Full Text Available Utility functions in Bayesian experimental design are usually based on the posterior distribution. When the posterior is found by simulation, it must be sampled from for each future dataset drawn from the prior predictive distribution. Many thousands of posterior distributions are often required. A popular technique in the Bayesian experimental design literature, which rapidly obtains samples from the posterior, is importance sampling, using the prior as the importance distribution. However, importance sampling from the prior will tend to break down if there is a reasonable number of experimental observations. In this paper, we explore the use of Laplace approximations in the design setting to overcome this drawback. Furthermore, we consider using the Laplace approximation to form the importance distribution to obtain a more efficient importance distribution than the prior. The methodology is motivated by a pharmacokinetic study, which investigates the effect of extracorporeal membrane oxygenation on the pharmacokinetics of antibiotics in sheep. The design problem is to find 10 near optimal plasma sampling times that produce precise estimates of pharmacokinetic model parameters/measures of interest. We consider several different utility functions of interest in these studies, which involve the posterior distribution of parameter functions.
A Bayesian framework for active artificial perception.
Ferreira, João Filipe; Lobo, Jorge; Bessière, Pierre; Castelo-Branco, Miguel; Dias, Jorge
2013-04-01
In this paper, we present a Bayesian framework for the active multimodal perception of 3-D structure and motion. The design of this framework finds its inspiration in the role of the dorsal perceptual pathway of the human brain. Its composing models build upon a common egocentric spatial configuration that is naturally fitting for the integration of readings from multiple sensors using a Bayesian approach. In the process, we will contribute with efficient and robust probabilistic solutions for cyclopean geometry-based stereovision and auditory perception based only on binaural cues, modeled using a consistent formalization that allows their hierarchical use as building blocks for the multimodal sensor fusion framework. We will explicitly or implicitly address the most important challenges of sensor fusion using this framework, for vision, audition, and vestibular sensing. Moreover, interaction and navigation require maximal awareness of spatial surroundings, which, in turn, is obtained through active attentional and behavioral exploration of the environment. The computational models described in this paper will support the construction of a simultaneously flexible and powerful robotic implementation of multimodal active perception to be used in real-world applications, such as human-machine interaction or mobile robot navigation. PMID:23014760
Phylogeography of Japanese encephalitis virus: genotype is associated with climate.
Directory of Open Access Journals (Sweden)
Amy J Schuh
Full Text Available The circulation of vector-borne zoonotic viruses is largely determined by the overlap in the geographical distributions of virus-competent vectors and reservoir hosts. What is less clear are the factors influencing the distribution of virus-specific lineages. Japanese encephalitis virus (JEV is the most important etiologic agent of epidemic encephalitis worldwide, and is primarily maintained between vertebrate reservoir hosts (avian and swine and culicine mosquitoes. There are five genotypes of JEV: GI-V. In recent years, GI has displaced GIII as the dominant JEV genotype and GV has re-emerged after almost 60 years of undetected virus circulation. JEV is found throughout most of Asia, extending from maritime Siberia in the north to Australia in the south, and as far as Pakistan to the west and Saipan to the east. Transmission of JEV in temperate zones is epidemic with the majority of cases occurring in summer months, while transmission in tropical zones is endemic and occurs year-round at lower rates. To test the hypothesis that viruses circulating in these two geographical zones are genetically distinct, we applied Bayesian phylogeographic, categorical data analysis and phylogeny-trait association test techniques to the largest JEV dataset compiled to date, representing the envelope (E gene of 487 isolates collected from 12 countries over 75 years. We demonstrated that GIII and the recently emerged GI-b are temperate genotypes likely maintained year-round in northern latitudes, while GI-a and GII are tropical genotypes likely maintained primarily through mosquito-avian and mosquito-swine transmission cycles. This study represents a new paradigm directly linking viral molecular evolution and climate.
Phylogeography of mtDNA haplogroup R7 in the Indian peninsula
Directory of Open Access Journals (Sweden)
Shukla Parul
2008-08-01
Full Text Available Abstract Background Human genetic diversity observed in Indian subcontinent is second only to that of Africa. This implies an early settlement and demographic growth soon after the first 'Out-of-Africa' dispersal of anatomically modern humans in Late Pleistocene. In contrast to this perspective, linguistic diversity in India has been thought to derive from more recent population movements and episodes of contact. With the exception of Dravidian, which origin and relatedness to other language phyla is obscure, all the language families in India can be linked to language families spoken in different regions of Eurasia. Mitochondrial DNA and Y chromosome evidence has supported largely local evolution of the genetic lineages of the majority of Dravidian and Indo-European speaking populations, but there is no consensus yet on the question of whether the Munda (Austro-Asiatic speaking populations originated in India or derive from a relatively recent migration from further East. Results Here, we report the analysis of 35 novel complete mtDNA sequences from India which refine the structure of Indian-specific varieties of haplogroup R. Detailed analysis of haplogroup R7, coupled with a survey of ~12,000 mtDNAs from caste and tribal groups over the entire Indian subcontinent, reveals that one of its more recently derived branches (R7a1, is particularly frequent among Munda-speaking tribal groups. This branch is nested within diverse R7 lineages found among Dravidian and Indo-European speakers of India. We have inferred from this that a subset of Munda-speaking groups have acquired R7 relatively recently. Furthermore, we find that the distribution of R7a1 within the Munda-speakers is largely restricted to one of the sub-branches (Kherwari of northern Munda languages. This evidence does not support the hypothesis that the Austro-Asiatic speakers are the primary source of the R7 variation. Statistical analyses suggest a significant correlation between
Bayesian methods for measures of agreement
Broemeling, Lyle D
2009-01-01
Using WinBUGS to implement Bayesian inferences of estimation and testing hypotheses, Bayesian Methods for Measures of Agreement presents useful methods for the design and analysis of agreement studies. It focuses on agreement among the various players in the diagnostic process.The author employs a Bayesian approach to provide statistical inferences based on various models of intra- and interrater agreement. He presents many examples that illustrate the Bayesian mode of reasoning and explains elements of a Bayesian application, including prior information, experimental information, the likelihood function, posterior distribution, and predictive distribution. The appendices provide the necessary theoretical foundation to understand Bayesian methods as well as introduce the fundamentals of programming and executing the WinBUGS software.Taking a Bayesian approach to inference, this hands-on book explores numerous measures of agreement, including the Kappa coefficient, the G coefficient, and intraclass correlation...
Plug & Play object oriented Bayesian networks
DEFF Research Database (Denmark)
Bangsø, Olav; Flores, J.; Jensen, Finn Verner
2003-01-01
Object oriented Bayesian networks have proven themselves useful in recent years. The idea of applying an object oriented approach to Bayesian networks has extended their scope to larger domains that can be divided into autonomous but interrelated entities. Object oriented Bayesian networks have...... been shown to be quite suitable for dynamic domains as well. However, processing object oriented Bayesian networks in practice does not take advantage of their modular structure. Normally the object oriented Bayesian network is transformed into a Bayesian network and, inference is performed...... by constructing a junction tree from this network. In this paper we propose a method for translating directly from object oriented Bayesian networks to junction trees, avoiding the intermediate translation. We pursue two main purposes: firstly, to maintain the original structure organized in an instance tree...
Bayesian Networks as a Decision Tool for O&M of Offshore Wind Turbines
DEFF Research Database (Denmark)
Nielsen, Jannie Jessen; Sørensen, John Dalsgaard
2010-01-01
Costs to operation and maintenance (O&M) of offshore wind turbines are large. This paper presents how influence diagrams can be used to assist in rational decision making for O&M. An influence diagram is a graphical representation of a decision tree based on Bayesian Networks. Bayesian Networks...... offer efficient Bayesian updating of a damage model when imperfect information from inspections/monitoring is available. The extension to an influence diagram offers the calculation of expected utilities for decision alternatives, and can be used to find the optimal strategy among different alternatives...
Bayesian analysis for the Burr type XII distribution based on record values
Directory of Open Access Journals (Sweden)
Mustafa Nadar
2013-05-01
Full Text Available In this paper we reviewed and extended some results that have been derived on record values from the two parameters Burr Type XII distribution. The two parameters were assumed to be random variables and Bayes estimates were derived on the basis of a linear exponential (LINEX loss function. Estimates for future record values were derived using non Bayesian and Bayesian approaches. In the Bayesian approach we reviewed the estimators obtained by Ahmedi and Doostparast (2006 using the well known squared error loss (SEL function and we derived estimate for the future record value under LINEX loss function. A numerical example with tables and figures illustrated the findings.
Flexible Bayesian Nonparametric Priors and Bayesian Computational Methods
Zhu, Weixuan
2016-01-01
The definition of vectors of dependent random probability measures is a topic of interest in Bayesian nonparametrics. They represent dependent nonparametric prior distributions that are useful for modelling observables for which specific covariate values are known. Our first contribution is the introduction of novel multivariate vectors of two-parameter Poisson-Dirichlet process. The dependence is induced by applying a L´evy copula to the marginal L´evy intensities. Our attenti...
Theoretical Analysis of Bayesian Optimisation with Unknown Gaussian Process Hyper-Parameters
Wang, Ziyu; De Freitas, Nando
2014-01-01
Bayesian optimisation has gained great popularity as a tool for optimising the parameters of machine learning algorithms and models. Somewhat ironically, setting up the hyper-parameters of Bayesian optimisation methods is notoriously hard. While reasonable practical solutions have been advanced, they can often fail to find the best optima. Surprisingly, there is little theoretical analysis of this crucial problem in the literature. To address this, we derive a cumulative regret bound for Baye...
Modeling Social Annotation: a Bayesian Approach
Plangprasopchok, Anon
2008-01-01
Collaborative tagging systems, such as del.icio.us, CiteULike, and others, allow users to annotate objects, e.g., Web pages or scientific papers, with descriptive labels called tags. The social annotations, contributed by thousands of users, can potentially be used to infer categorical knowledge, classify documents or recommend new relevant information. Traditional text inference methods do not make best use of socially-generated data, since they do not take into account variations in individual users' perspectives and vocabulary. In a previous work, we introduced a simple probabilistic model that takes interests of individual annotators into account in order to find hidden topics of annotated objects. Unfortunately, our proposed approach had a number of shortcomings, including overfitting, local maxima and the requirement to specify values for some parameters. In this paper we address these shortcomings in two ways. First, we extend the model to a fully Bayesian framework. Second, we describe an infinite ver...
Learning Bayesian network structure with immune algorithm
Institute of Scientific and Technical Information of China (English)
Zhiqiang Cai; Shubin Si; Shudong Sun; Hongyan Dui
2015-01-01
Finding out reasonable structures from bulky data is one of the difficulties in modeling of Bayesian network (BN), which is also necessary in promoting the application of BN. This pa-per proposes an immune algorithm based method (BN-IA) for the learning of the BN structure with the idea of vaccination. Further-more, the methods on how to extract the effective vaccines from local optimal structure and root nodes are also described in details. Final y, the simulation studies are implemented with the helicopter convertor BN model and the car start BN model. The comparison results show that the proposed vaccines and the BN-IA can learn the BN structure effectively and efficiently.
Bayesian variable selection with spherically symmetric priors
De Kock, M B
2014-01-01
We propose that Bayesian variable selection for linear parametrisations with Gaussian iid likelihoods be based on the spherical symmetry of the diagonalised parameter space. This reduces the multidimensional parameter space problem to one dimension without the need for conjugate priors. Combining this likelihood with what we call the r-prior results in a framework in which we can derive closed forms for the evidence, posterior and characteristic function for four different r-priors, including the hyper-g prior and the Zellner-Siow prior, which are shown to be special cases of our r-prior. Two scenarios of a single variable dispersion parameter and of fixed dispersion are studied separately, and asymptotic forms comparable to the traditional information criteria are derived. In a simple simulation exercise, we find that model comparison based on our uniform r-prior appears to fare better than the current model comparison schemes.
Bayesian Overlapping Community Detection in Dynamic Networks
Ghorbani, Mahsa; Khodadadi, Ali
2016-01-01
Detecting community structures in social networks has gained considerable attention in recent years. However, lack of prior knowledge about the number of communities, and their overlapping nature have made community detection a challenging problem. Moreover, many of the existing methods only consider static networks, while most of real world networks are dynamic and evolve over time. Hence, finding consistent overlapping communities in dynamic networks without any prior knowledge about the number of communities is still an interesting open research problem. In this paper, we present an overlapping community detection method for dynamic networks called Dynamic Bayesian Overlapping Community Detector (DBOCD). DBOCD assumes that in every snapshot of network, overlapping parts of communities are dense areas and utilizes link communities instead of common node communities. Using Recurrent Chinese Restaurant Process and community structure of the network in the last snapshot, DBOCD simultaneously extracts the numbe...
Bayesian optimization for tuning chaotic systems
Directory of Open Access Journals (Sweden)
M. Abbas
2014-08-01
Full Text Available In this work, we consider the Bayesian optimization (BO approach for tuning parameters of complex chaotic systems. Such problems arise, for instance, in tuning the sub-grid scale parameterizations in weather and climate models. For such problems, the tuning procedure is generally based on a performance metric which measures how well the tuned model fits the data. This tuning is often a computationally expensive task. We show that BO, as a tool for finding the extrema of computationally expensive objective functions, is suitable for such tuning tasks. In the experiments, we consider tuning parameters of two systems: a simplified atmospheric model and a low-dimensional chaotic system. We show that BO is able to tune parameters of both the systems with a low number of objective function evaluations and without the need of any gradient information.
Safety culture in Bayesian and legal contexts
International Nuclear Information System (INIS)
While contemplating the similarities between the law of torts and concepts of safety, the author realized that there was a close correspondence between the law of negligence and the way safety ought to be generally defined. This definition of safety is provided herein. A safety culture must have an adequate definition of safety in order to function most effectively. This paper provides a practical definition of safety that answers the question 'How safe is safe enough? The development rests on two bases: the subjectivistic-Bayesian definition of probability and certain legal definitions primarily from the tort law of negligence. The development also leads to the conclusion that one cannot generally expect greater specificity in determining how safe is safe enough than one finds in the legal definition of liability under the tort of negligence. It then follows that some of the public's aversion to complex technical undertakings is rooted in its typically intuitive and vague notions concerning safety
Bayesian approach to rough set
Marwala, Tshilidzi
2007-01-01
This paper proposes an approach to training rough set models using Bayesian framework trained using Markov Chain Monte Carlo (MCMC) method. The prior probabilities are constructed from the prior knowledge that good rough set models have fewer rules. Markov Chain Monte Carlo sampling is conducted through sampling in the rough set granule space and Metropolis algorithm is used as an acceptance criteria. The proposed method is tested to estimate the risk of HIV given demographic data. The results obtained shows that the proposed approach is able to achieve an average accuracy of 58% with the accuracy varying up to 66%. In addition the Bayesian rough set give the probabilities of the estimated HIV status as well as the linguistic rules describing how the demographic parameters drive the risk of HIV.
Attention in a bayesian framework
DEFF Research Database (Denmark)
Whiteley, Louise Emma; Sahani, Maneesh
2012-01-01
include both selective phenomena, where attention is invoked by cues that point to particular stimuli, and integrative phenomena, where attention is invoked dynamically by endogenous processing. However, most previous Bayesian accounts of attention have focused on describing relatively simple experimental...... settings, where cues shape expectations about a small number of upcoming stimuli and thus convey "prior" information about clearly defined objects. While operationally consistent with the experiments it seeks to describe, this view of attention as prior seems to miss many essential elements of both its......The behavioral phenomena of sensory attention are thought to reflect the allocation of a limited processing resource, but there is little consensus on the nature of the resource or why it should be limited. Here we argue that a fundamental bottleneck emerges naturally within Bayesian models of...
Bayesian Sampling using Condition Indicators
DEFF Research Database (Denmark)
Faber, Michael H.; Sørensen, John Dalsgaard
2002-01-01
allows for a Bayesian formulation of the indicators whereby the experience and expertise of the inspection personnel may be fully utilized and consistently updated as frequentistic information is collected. The approach is illustrated on an example considering a concrete structure subject to corrosion......The problem of control quality of components is considered for the special case where the acceptable failure rate is low, the test costs are high and where it may be difficult or impossible to test the condition of interest directly. Based on the classical control theory and the concept of...... condition indicators introduced by Benjamin and Cornell (1970) a Bayesian approach to quality control is formulated. The formulation is then extended to the case where the quality control is based on sampling of indirect information about the condition of the components, i.e. condition indicators. This...
BAYESIAN IMAGE RESTORATION, USING CONFIGURATIONS
Directory of Open Access Journals (Sweden)
Thordis Linda Thorarinsdottir
2011-05-01
Full Text Available In this paper, we develop a Bayesian procedure for removing noise from images that can be viewed as noisy realisations of random sets in the plane. The procedure utilises recent advances in configuration theory for noise free random sets, where the probabilities of observing the different boundary configurations are expressed in terms of the mean normal measure of the random set. These probabilities are used as prior probabilities in a Bayesian image restoration approach. Estimation of the remaining parameters in the model is outlined for salt and pepper noise. The inference in the model is discussed in detail for 3 X 3 and 5 X 5 configurations and examples of the performance of the procedure are given.
Bayesian priors for transiting planets
Kipping, David M
2016-01-01
As astronomers push towards discovering ever-smaller transiting planets, it is increasingly common to deal with low signal-to-noise ratio (SNR) events, where the choice of priors plays an influential role in Bayesian inference. In the analysis of exoplanet data, the selection of priors is often treated as a nuisance, with observers typically defaulting to uninformative distributions. Such treatments miss a key strength of the Bayesian framework, especially in the low SNR regime, where even weak a priori information is valuable. When estimating the parameters of a low-SNR transit, two key pieces of information are known: (i) the planet has the correct geometric alignment to transit and (ii) the transit event exhibits sufficient signal-to-noise to have been detected. These represent two forms of observational bias. Accordingly, when fitting transits, the model parameter priors should not follow the intrinsic distributions of said terms, but rather those of both the intrinsic distributions and the observational ...
Bayesian Inference for Radio Observations
Lochner, Michelle; Zwart, Jonathan T L; Smirnov, Oleg; Bassett, Bruce A; Oozeer, Nadeem; Kunz, Martin
2015-01-01
(Abridged) New telescopes like the Square Kilometre Array (SKA) will push into a new sensitivity regime and expose systematics, such as direction-dependent effects, that could previously be ignored. Current methods for handling such systematics rely on alternating best estimates of instrumental calibration and models of the underlying sky, which can lead to inaccurate uncertainty estimates and biased results because such methods ignore any correlations between parameters. These deconvolution algorithms produce a single image that is assumed to be a true representation of the sky, when in fact it is just one realisation of an infinite ensemble of images compatible with the noise in the data. In contrast, here we report a Bayesian formalism that simultaneously infers both systematics and science. Our technique, Bayesian Inference for Radio Observations (BIRO), determines all parameters directly from the raw data, bypassing image-making entirely, by sampling from the joint posterior probability distribution. Thi...
Bayesian inference on proportional elections.
Brunello, Gabriel Hideki Vatanabe; Nakano, Eduardo Yoshio
2015-01-01
Polls for majoritarian voting systems usually show estimates of the percentage of votes for each candidate. However, proportional vote systems do not necessarily guarantee the candidate with the most percentage of votes will be elected. Thus, traditional methods used in majoritarian elections cannot be applied on proportional elections. In this context, the purpose of this paper was to perform a Bayesian inference on proportional elections considering the Brazilian system of seats distribution. More specifically, a methodology to answer the probability that a given party will have representation on the chamber of deputies was developed. Inferences were made on a Bayesian scenario using the Monte Carlo simulation technique, and the developed methodology was applied on data from the Brazilian elections for Members of the Legislative Assembly and Federal Chamber of Deputies in 2010. A performance rate was also presented to evaluate the efficiency of the methodology. Calculations and simulations were carried out using the free R statistical software. PMID:25786259
A Bayesian Nonparametric IRT Model
Karabatsos, George
2015-01-01
This paper introduces a flexible Bayesian nonparametric Item Response Theory (IRT) model, which applies to dichotomous or polytomous item responses, and which can apply to either unidimensional or multidimensional scaling. This is an infinite-mixture IRT model, with person ability and item difficulty parameters, and with a random intercept parameter that is assigned a mixing distribution, with mixing weights a probit function of other person and item parameters. As a result of its flexibility...
Bayesian segmentation of hyperspectral images
Mohammadpour, Adel; Mohammad-Djafari, Ali
2007-01-01
In this paper we consider the problem of joint segmentation of hyperspectral images in the Bayesian framework. The proposed approach is based on a Hidden Markov Modeling (HMM) of the images with common segmentation, or equivalently with common hidden classification label variables which is modeled by a Potts Markov Random Field. We introduce an appropriate Markov Chain Monte Carlo (MCMC) algorithm to implement the method and show some simulation results.
Bayesian segmentation of hyperspectral images
Mohammadpour, Adel; Féron, Olivier; Mohammad-Djafari, Ali
2004-11-01
In this paper we consider the problem of joint segmentation of hyperspectral images in the Bayesian framework. The proposed approach is based on a Hidden Markov Modeling (HMM) of the images with common segmentation, or equivalently with common hidden classification label variables which is modeled by a Potts Markov Random Field. We introduce an appropriate Markov Chain Monte Carlo (MCMC) algorithm to implement the method and show some simulation results.
Bayesian Stable Isotope Mixing Models
Parnell, Andrew C.; Phillips, Donald L.; Bearhop, Stuart; Semmens, Brice X.; Ward, Eric J.; Moore, Jonathan W.; Andrew L Jackson; Inger, Richard
2012-01-01
In this paper we review recent advances in Stable Isotope Mixing Models (SIMMs) and place them into an over-arching Bayesian statistical framework which allows for several useful extensions. SIMMs are used to quantify the proportional contributions of various sources to a mixture. The most widely used application is quantifying the diet of organisms based on the food sources they have been observed to consume. At the centre of the multivariate statistical model we propose is a compositional m...
Bayesian Network--Response Regression
WANG, LU; Durante, Daniele; Dunson, David B.
2016-01-01
There is an increasing interest in learning how human brain networks vary with continuous traits (e.g., personality, cognitive abilities, neurological disorders), but flexible procedures to accomplish this goal are limited. We develop a Bayesian semiparametric model, which combines low-rank factorizations and Gaussian process priors to allow flexible shifts of the conditional expectation for a network-valued random variable across the feature space, while including subject-specific random eff...
Bayesian estimation of turbulent motion
Héas, P.; Herzet, C.; Mémin, E.; Heitz, D.; P. D. Mininni
2013-01-01
International audience Based on physical laws describing the multi-scale structure of turbulent flows, this article proposes a regularizer for fluid motion estimation from an image sequence. Regularization is achieved by imposing some scale invariance property between histograms of motion increments computed at different scales. By reformulating this problem from a Bayesian perspective, an algorithm is proposed to jointly estimate motion, regularization hyper-parameters, and to select the ...
Elements of Bayesian experimental design
Energy Technology Data Exchange (ETDEWEB)
Sivia, D.S. [Rutherford Appleton Lab., Oxon (United Kingdom)
1997-09-01
We consider some elements of the Bayesian approach that are important for optimal experimental design. While the underlying principles used are very general, and are explained in detail in a recent tutorial text, they are applied here to the specific case of characterising the inferential value of different resolution peakshapes. This particular issue was considered earlier by Silver, Sivia and Pynn (1989, 1990a, 1990b), and the following presentation confirms and extends the conclusions of their analysis.
Skill Rating by Bayesian Inference
Di Fatta, Giuseppe; Haworth, Guy McCrossan; Regan, Kenneth W.
2009-01-01
Systems Engineering often involves computer modelling the behaviour of proposed systems and their components. Where a component is human, fallibility must be modelled by a stochastic agent. The identification of a model of decision-making over quantifiable options is investigated using the game-domain of Chess. Bayesian methods are used to infer the distribution of players’ skill levels from the moves they play rather than from their competitive results. The approach is used on large sets of ...
Topics in Nonparametric Bayesian Statistics
2003-01-01
The intersection set of Bayesian and nonparametric statistics was almost empty until about 1973, but now seems to be growing at a healthy rate. This chapter gives an overview of various theoretical and applied research themes inside this field, partly complementing and extending recent reviews of Dey, Müller and Sinha (1998) and Walker, Damien, Laud and Smith (1999). The intention is not to be complete or exhaustive, but rather to touch on research areas of interest, partly by example.
Cover Tree Bayesian Reinforcement Learning
Tziortziotis, Nikolaos; Dimitrakakis, Christos; Blekas, Konstantinos
2013-01-01
This paper proposes an online tree-based Bayesian approach for reinforcement learning. For inference, we employ a generalised context tree model. This defines a distribution on multivariate Gaussian piecewise-linear models, which can be updated in closed form. The tree structure itself is constructed using the cover tree method, which remains efficient in high dimensional spaces. We combine the model with Thompson sampling and approximate dynamic programming to obtain effective exploration po...
Bayesian kinematic earthquake source models
Minson, S. E.; Simons, M.; Beck, J. L.; Genrich, J. F.; Galetzka, J. E.; Chowdhury, F.; Owen, S. E.; Webb, F.; Comte, D.; Glass, B.; Leiva, C.; Ortega, F. H.
2009-12-01
Most coseismic, postseismic, and interseismic slip models are based on highly regularized optimizations which yield one solution which satisfies the data given a particular set of regularizing constraints. This regularization hampers our ability to answer basic questions such as whether seismic and aseismic slip overlap or instead rupture separate portions of the fault zone. We present a Bayesian methodology for generating kinematic earthquake source models with a focus on large subduction zone earthquakes. Unlike classical optimization approaches, Bayesian techniques sample the ensemble of all acceptable models presented as an a posteriori probability density function (PDF), and thus we can explore the entire solution space to determine, for example, which model parameters are well determined and which are not, or what is the likelihood that two slip distributions overlap in space. Bayesian sampling also has the advantage that all a priori knowledge of the source process can be used to mold the a posteriori ensemble of models. Although very powerful, Bayesian methods have up to now been of limited use in geophysical modeling because they are only computationally feasible for problems with a small number of free parameters due to what is called the "curse of dimensionality." However, our methodology can successfully sample solution spaces of many hundreds of parameters, which is sufficient to produce finite fault kinematic earthquake models. Our algorithm is a modification of the tempered Markov chain Monte Carlo (tempered MCMC or TMCMC) method. In our algorithm, we sample a "tempered" a posteriori PDF using many MCMC simulations running in parallel and evolutionary computation in which models which fit the data poorly are preferentially eliminated in favor of models which better predict the data. We present results for both synthetic test problems as well as for the 2007 Mw 7.8 Tocopilla, Chile earthquake, the latter of which is constrained by InSAR, local high
Bayesian Kernel Mixtures for Counts
Canale, Antonio; David B Dunson
2011-01-01
Although Bayesian nonparametric mixture models for continuous data are well developed, there is a limited literature on related approaches for count data. A common strategy is to use a mixture of Poissons, which unfortunately is quite restrictive in not accounting for distributions having variance less than the mean. Other approaches include mixing multinomials, which requires finite support, and using a Dirichlet process prior with a Poisson base measure, which does not allow smooth deviatio...
Bayesian Optimization for Adaptive MCMC
Mahendran, Nimalan; Wang, Ziyu; Hamze, Firas; De Freitas, Nando
2011-01-01
This paper proposes a new randomized strategy for adaptive MCMC using Bayesian optimization. This approach applies to non-differentiable objective functions and trades off exploration and exploitation to reduce the number of potentially costly objective function evaluations. We demonstrate the strategy in the complex setting of sampling from constrained, discrete and densely connected probabilistic graphical models where, for each variation of the problem, one needs to adjust the parameters o...
Inference in hybrid Bayesian networks
DEFF Research Database (Denmark)
Lanseth, Helge; Nielsen, Thomas Dyhre; Rumí, Rafael;
2009-01-01
and reliability block diagrams). However, limitations in the BNs' calculation engine have prevented BNs from becoming equally popular for domains containing mixtures of both discrete and continuous variables (so-called hybrid domains). In this paper we focus on these difficulties, and summarize some of the last...... decade's research on inference in hybrid Bayesian networks. The discussions are linked to an example model for estimating human reliability....
Quantile pyramids for Bayesian nonparametrics
2009-01-01
P\\'{o}lya trees fix partitions and use random probabilities in order to construct random probability measures. With quantile pyramids we instead fix probabilities and use random partitions. For nonparametric Bayesian inference we use a prior which supports piecewise linear quantile functions, based on the need to work with a finite set of partitions, yet we show that the limiting version of the prior exists. We also discuss and investigate an alternative model based on the so-called substitut...
Space Shuttle RTOS Bayesian Network
Morris, A. Terry; Beling, Peter A.
2001-01-01
With shrinking budgets and the requirements to increase reliability and operational life of the existing orbiter fleet, NASA has proposed various upgrades for the Space Shuttle that are consistent with national space policy. The cockpit avionics upgrade (CAU), a high priority item, has been selected as the next major upgrade. The primary functions of cockpit avionics include flight control, guidance and navigation, communication, and orbiter landing support. Secondary functions include the provision of operational services for non-avionics systems such as data handling for the payloads and caution and warning alerts to the crew. Recently, a process to selection the optimal commercial-off-the-shelf (COTS) real-time operating system (RTOS) for the CAU was conducted by United Space Alliance (USA) Corporation, which is a joint venture between Boeing and Lockheed Martin, the prime contractor for space shuttle operations. In order to independently assess the RTOS selection, NASA has used the Bayesian network-based scoring methodology described in this paper. Our two-stage methodology addresses the issue of RTOS acceptability by incorporating functional, performance and non-functional software measures related to reliability, interoperability, certifiability, efficiency, correctness, business, legal, product history, cost and life cycle. The first stage of the methodology involves obtaining scores for the various measures using a Bayesian network. The Bayesian network incorporates the causal relationships between the various and often competing measures of interest while also assisting the inherently complex decision analysis process with its ability to reason under uncertainty. The structure and selection of prior probabilities for the network is extracted from experts in the field of real-time operating systems. Scores for the various measures are computed using Bayesian probability. In the second stage, multi-criteria trade-off analyses are performed between the scores
Bayesian analysis of contingency tables
Gómez Villegas, Miguel A.; González Pérez, Beatriz
2005-01-01
The display of the data by means of contingency tables is used in different approaches to statistical inference, for example, to broach the test of homogeneity of independent multinomial distributions. We develop a Bayesian procedure to test simple null hypotheses versus bilateral alternatives in contingency tables. Given independent samples of two binomial distributions and taking a mixed prior distribution, we calculate the posterior probability that the proportion of successes in the first...
Bayesian Credit Ratings (new version)
Paola Cerchiello; Paolo Giudici
2013-01-01
In this contribution we aim at improving ordinal variable selection in the context of causal models. In this regard, we propose an approach that provides a formal inferential tool to compare the explanatory power of each covariate, and, therefore, to select an effective model for classification purposes. Our proposed model is Bayesian nonparametric, and, thus, keeps the amount of model specification to a minimum. We consider the case in which information from the covariates is at the ordinal ...
Bayesian second law of thermodynamics
Bartolotta, Anthony; Carroll, Sean M.; Leichenauer, Stefan; Pollack, Jason
2016-08-01
We derive a generalization of the second law of thermodynamics that uses Bayesian updates to explicitly incorporate the effects of a measurement of a system at some point in its evolution. By allowing an experimenter's knowledge to be updated by the measurement process, this formulation resolves a tension between the fact that the entropy of a statistical system can sometimes fluctuate downward and the information-theoretic idea that knowledge of a stochastically evolving system degrades over time. The Bayesian second law can be written as Δ H (ρm,ρ ) + F |m≥0 , where Δ H (ρm,ρ ) is the change in the cross entropy between the original phase-space probability distribution ρ and the measurement-updated distribution ρm and F |m is the expectation value of a generalized heat flow out of the system. We also derive refined versions of the second law that bound the entropy increase from below by a non-negative number, as well as Bayesian versions of integral fluctuation theorems. We demonstrate the formalism using simple analytical and numerical examples.
Quantum Inference on Bayesian Networks
Yoder, Theodore; Low, Guang Hao; Chuang, Isaac
2014-03-01
Because quantum physics is naturally probabilistic, it seems reasonable to expect physical systems to describe probabilities and their evolution in a natural fashion. Here, we use quantum computation to speedup sampling from a graphical probability model, the Bayesian network. A specialization of this sampling problem is approximate Bayesian inference, where the distribution on query variables is sampled given the values e of evidence variables. Inference is a key part of modern machine learning and artificial intelligence tasks, but is known to be NP-hard. Classically, a single unbiased sample is obtained from a Bayesian network on n variables with at most m parents per node in time (nmP(e) - 1 / 2) , depending critically on P(e) , the probability the evidence might occur in the first place. However, by implementing a quantum version of rejection sampling, we obtain a square-root speedup, taking (n2m P(e) -1/2) time per sample. The speedup is the result of amplitude amplification, which is proving to be broadly applicable in sampling and machine learning tasks. In particular, we provide an explicit and efficient circuit construction that implements the algorithm without the need for oracle access.
Levy, Hila; Clucas, Gemma V; Rogers, Alex D; Leaché, Adam D; Ciborowski, Kate L; Polito, Michael J; Lynch, Heather J; Dunn, Michael J; Hart, Tom
2016-03-01
Climate change, fisheries' pressure on penguin prey, and direct human disturbance of wildlife have all been implicated in causing large shifts in the abundance and distribution of penguins in the Southern Ocean. Without mark-recapture studies, understanding how colonies form and, by extension, how ranges shift is challenging. Genetic studies, particularly focused on newly established colonies, provide a snapshot of colonization and can reveal the extent to which shifts in abundance and occupancy result from changes in demographic rates (e.g., reproduction and survival) or migration among suitable patches of habitat. Here, we describe the population structure of a colonial seabird breeding across a large latitudinal range in the Southern Ocean. Using multilocus microsatellite genotype data from 510 Gentoo penguin (Pygoscelis papua) individuals from 14 colonies along the Scotia Arc and Antarctic Peninsula, together with mitochondrial DNA data, we find strong genetic differentiation between colonies north and south of the Polar Front, that coincides geographically with the taxonomic boundary separating the subspecies P. p. papua and P. p. ellsworthii. Using a discrete Bayesian phylogeographic approach, we show that southern Gentoos expanded from a possible glacial refuge in the center of their current range, colonizing regions to the north and south through rare, long-distance dispersal. Our findings show that this dispersal is important for new colony foundation and range expansion in a seabird species that ordinarily exhibits high levels of natal philopatry, though persistent oceanographic features serve as barriers to movement. PMID:26933489
12th Brazilian Meeting on Bayesian Statistics
Louzada, Francisco; Rifo, Laura; Stern, Julio; Lauretto, Marcelo
2015-01-01
Through refereed papers, this volume focuses on the foundations of the Bayesian paradigm; their comparison to objectivistic or frequentist Statistics counterparts; and the appropriate application of Bayesian foundations. This research in Bayesian Statistics is applicable to data analysis in biostatistics, clinical trials, law, engineering, and the social sciences. EBEB, the Brazilian Meeting on Bayesian Statistics, is held every two years by the ISBrA, the International Society for Bayesian Analysis, one of the most active chapters of the ISBA. The 12th meeting took place March 10-14, 2014 in Atibaia. Interest in foundations of inductive Statistics has grown recently in accordance with the increasing availability of Bayesian methodological alternatives. Scientists need to deal with the ever more difficult choice of the optimal method to apply to their problem. This volume shows how Bayes can be the answer. The examination and discussion on the foundations work towards the goal of proper application of Bayesia...
A Decomposition Algorithm for Learning Bayesian Network Structures from Data
DEFF Research Database (Denmark)
Zeng, Yifeng; Cordero Hernandez, Jorge
2008-01-01
the complete network. The new learning algorithm firstly finds local components from the data, and then recover the complete network by joining the learned components. We show the empirical performance of the decomposition algorithm in several benchmark networks.......It is a challenging task of learning a large Bayesian network from a small data set. Most conventional structural learning approaches run into the computational as well as the statistical problems. We propose a decomposition algorithm for the structure construction without having to learn...
Directory of Open Access Journals (Sweden)
Tania Anaid Gutiérrez-García
2013-12-01
Full Text Available Central America is an ideal region for comparative phylogeographic studies because of its intricate geologic and biogeographic history, diversity of habitats and dynamic climatic and tectonic history. The aim of this work was to assess the phylogeography of two rodents codistributed throughout Central America, in order to identify if they show concordant genetic and phylogeographic patterns. The synopsis includes four parts: (1 an overview of the field of comparative phylogeography; (2 a detailed review that describes how genetic and geologic studies can be combined to elucidate general patterns of the biogeographic and evolutionary history of Central America; and a phylogeographic analysis of two species at both the (3 intraspecific and (4 comparative phylogeographic levels. The last incorporates specific ecological features and evaluates their influence on the species’ genetic patterns. Results showed a concordant genetic structure influenced by geographic distance for both rodents, but dissimilar dispersal patterns due to ecological features and life history.
Bayesian Posterior Distributions Without Markov Chains
Cole, Stephen R.; Chu, Haitao; Greenland, Sander; Hamra, Ghassan; Richardson, David B.
2012-01-01
Bayesian posterior parameter distributions are often simulated using Markov chain Monte Carlo (MCMC) methods. However, MCMC methods are not always necessary and do not help the uninitiated understand Bayesian inference. As a bridge to understanding Bayesian inference, the authors illustrate a transparent rejection sampling method. In example 1, they illustrate rejection sampling using 36 cases and 198 controls from a case-control study (1976–1983) assessing the relation between residential ex...
Joaquín Muñoz; Francisco Amat; Green, Andy J.; Jordi Figuerola; Africa Gómez
2013-01-01
Since Darwin’s time, waterbirds have been considered an important vector for the dispersal of continental aquatic invertebrates. Bird movements have facilitated the worldwide invasion of the American brine shrimp Artemia franciscana, transporting cysts (diapausing eggs), and favouring rapid range expansions from introduction sites. Here we address the impact of bird migratory flyways on the population genetic structure and phylogeography of A. franciscana in its native range in the Americas. ...
Meyer Axel; Barluenga Marta
2010-01-01
Abstract Background Elucidation of the mechanisms driving speciation requires detailed knowledge about the phylogenetic relationships and phylogeography of the incipient species within their entire ranges as well as their colonization history. The Midas cichlid species complex Amphilophus spp. has been proven to be a powerful model system for the study of ecological specialization, sexual selection and the mechanisms of sympatric speciation. Here we present a comprehensive and integrative phy...
Arteaga, Alejandro; Pyron, R. Alexander; Peñafiel, Nicolás; Romero-Barreto, Paulina; Culebras, Jaime; Bustamante, Lucas; Yánez-Muñoz, Mario H.; Guayasamin, Juan M.
2016-01-01
Comparative phylogeography allow us to understand how shared historical circumstances have shaped the formation of lineages, by examining a broad spectrum of co-distributed populations of different taxa. However, these types of studies are scarce in the Neotropics, a region that is characterized by high diversity, complex geology, and poorly understood biogeography. Here, we investigate the diversification patterns of five lineages of amphibians and reptiles, co-distributed across the Choco a...
Honorio Coronado, Eurídice N.; Dexter, Kyle G.; Poelchau, Monica F; Hollingsworth, Peter M; Phillips, Oliver L.; Pennington, R Toby
2014-01-01
Aim: To examine the phylogeography of Ficus insipida subsp. insipida in order to investigate patterns of spatial genetic structure across the Neotropics and within Amazonia. Location: Neotropics. Methods: Plastid DNA (trnH-psbA; 410 individuals from 54 populations) and nuclear ribosomal internal transcribed spacer (ITS; 85 individuals from 27 populations) sequences were sampled from Mexico to Bolivia, representing the full extent of the taxon's distribution. Divergence of plastid lineages was...
Sezonlin, Michel; Dupas, Stéphane; Le Ru, Bruno; Le Gall, Philippe; Moyal, Pascal; Calatayud, Paul-André; Giffard, I; Faure, N; Silvain, Jean-François
2006-01-01
The population genetics and phylogeography of African phytophagous insects have received little attention. Some, such as the maize stalk borer Busseola fusca, display significant geographic differences in ecological preferences that may be congruent with patterns of molecular variation. To test this, we collected 307 individuals of this species from maize and cultivated sorghum at 52 localities in West, Central and East Africa during the growing season. For all collected individuals, we seque...
Tahsin, Tasnia; Beard, Rachel; Rivera, Robert; Lauder, Rob; Wallstrom, Garrick; Scotch, Matthew; Gonzalez, Graciela
2014-01-01
Zoonotic viruses represent emerging or re-emerging pathogens that pose significant public health threats throughout the world. It is therefore crucial to advance current surveillance mechanisms for these viruses through outlets such as phylogeography. Despite the abundance of zoonotic viral sequence data in publicly available databases such as GenBank, phylogeographic analysis of these viruses is often limited by the lack of adequate geographic metadata. However, many GenBank records include references to articles with more detailed information and automated systems may help extract this information efficiently and effectively. In this paper, we describe our efforts to determine the proportion of GenBank records with "insufficient" geographic metadata for seven well-studied viruses. We also evaluate the performance of four different Named Entity Recognition (NER) systems for automatically extracting related entities using a manually created gold-standard. PMID:25717409
Bayesian networks with applications in reliability analysis
Langseth, Helge
2002-01-01
A common goal of the papers in this thesis is to propose, formalize and exemplify the use of Bayesian networks as a modelling tool in reliability analysis. The papers span work in which Bayesian networks are merely used as a modelling tool (Paper I), work where models are specially designed to utilize the inference algorithms of Bayesian networks (Paper II and Paper III), and work where the focus has been on extending the applicability of Bayesian networks to very large domains (Paper IV and ...
Mokrousov, Igor; Vyazovaya, Anna; Iwamoto, Tomotada; Skiba, Yuriy; Pole, Ilva; Zhdanova, Svetlana; Arikawa, Kentaro; Sinkov, Viacheslav; Umpeleva, Tatiana; Valcheva, Violeta; Alvarez Figueroa, Maria; Ranka, Renate; Jansone, Inta; Ogarkov, Oleg; Zhuravlev, Viacheslav; Narvskaya, Olga
2016-06-01
Currently, Mycobacterium tuberculosis isolates of Latin-American Mediterranean (LAM) family may be detected far beyond the geographic areas that coined its name 15years ago. Here, we established the framework phylogeny of this geographically intriguing and pathobiologically important mycobacterial lineage and hypothesized how human demographics and migration influenced its phylogeography. Phylogenetic analysis of LAM isolates from all continents based on 24 variable number of tandem repeats (VNTR) loci and other markers identified three global sublineages with certain geographic affinities and defined by large deletions RD115, RD174, and by spoligotype SIT33. One minor sublineage (spoligotype SIT388) appears endemic in Japan. One-locus VNTR signatures were established for sublineages and served for their search in published literature and geographic mapping. We suggest that the LAM family originated in the Western Mediterranean region. The most widespread RD115 sublineage seems the most ancient and encompasses genetically and geographically distant branches, including extremely drug resistant KZN in South Africa and LAM-RUS recently widespread across Northern Eurasia. The RD174 sublineage likely started its active spread in Brazil; its earlier branch is relatively dominated by isolates from South America and the derived one is dominated by Portuguese and South/Southeastern African isolates. The relatively most recent SIT33-sublineage is marked with enigmatic gaps and peaks across the Americas and includes South African clade F11/RD761, which likely emerged within the SIT33 subpopulation after its arrival to Africa. In addition to SIT388-sublineage, other deeply rooted, endemic LAM sublineages may exist that remain to be discovered. As a general conclusion, human mass migration appears to be the major factor that shaped the M. tuberculosis phylogeography over large time-spans. PMID:27001605
Bayesian Methods and Universal Darwinism
Campbell, John
2009-12-01
Bayesian methods since the time of Laplace have been understood by their practitioners as closely aligned to the scientific method. Indeed a recent Champion of Bayesian methods, E. T. Jaynes, titled his textbook on the subject Probability Theory: the Logic of Science. Many philosophers of science including Karl Popper and Donald Campbell have interpreted the evolution of Science as a Darwinian process consisting of a `copy with selective retention' algorithm abstracted from Darwin's theory of Natural Selection. Arguments are presented for an isomorphism between Bayesian Methods and Darwinian processes. Universal Darwinism, as the term has been developed by Richard Dawkins, Daniel Dennett and Susan Blackmore, is the collection of scientific theories which explain the creation and evolution of their subject matter as due to the Operation of Darwinian processes. These subject matters span the fields of atomic physics, chemistry, biology and the social sciences. The principle of Maximum Entropy states that Systems will evolve to states of highest entropy subject to the constraints of scientific law. This principle may be inverted to provide illumination as to the nature of scientific law. Our best cosmological theories suggest the universe contained much less complexity during the period shortly after the Big Bang than it does at present. The scientific subject matter of atomic physics, chemistry, biology and the social sciences has been created since that time. An explanation is proposed for the existence of this subject matter as due to the evolution of constraints in the form of adaptations imposed on Maximum Entropy. It is argued these adaptations were discovered and instantiated through the Operations of a succession of Darwinian processes.
Range Wide Phylogeography of Dactylopius coccus (Hemiptera: Dactylopiidae)
DEFF Research Database (Denmark)
Van Dam, Alex; Portillo Martinez, Liberato; Jeri Chavez, Antonio;
2015-01-01
studies have been inconclusive. Here, we fill in the remaining gaps in the ecological record and look for corroboration from mtDNA markers as to the origin of this species. We use three mtDNA genes (CO1, tRNA-Leucine, and CO2) spanning 1294 bp, along with climate niche modeling of Holocene and Pleistocene...... cochineal distributions. We find the center of origin of D. coccus to be Oaxaca Mexico based on mtDNA data and climate niche modeling. Further meta-genomic data are needed to rule out selective sweeps from past and present endosymbionts for these results to be definitive....
Brochu, Eric; de Freitas, Nando
2010-01-01
We present a tutorial on Bayesian optimization, a method of finding the maximum of expensive cost functions. Bayesian optimization employs the Bayesian technique of setting a prior over the objective function and combining it with evidence to get a posterior function. This permits a utility-based selection of the next observation to make on the objective function, which must take into account both exploration (sampling from areas of high uncertainty) and exploitation (sampling areas likely to offer improvement over the current best observation). We also present two detailed extensions of Bayesian optimization, with experiments---active user modelling with preferences, and hierarchical reinforcement learning---and a discussion of the pros and cons of Bayesian optimization based on our experiences.