Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Bayesian inference”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23Linked to original sources

Bayesian class discovery in microarray datasets.

A novel approach to class discovery in gene expression datasets is presented. In the context of clinical diagnosis, the central goal of class discovery algorithms is to simultaneously find putative (sub-)types of diseases and to identify informative subsets of genes with disease-type specific expression profile. Contrary to many other approaches in the literature, the method presented implements a wrapper strategy for feature selection, in the sense that the features are directly selected by optimizing the discriminative power of the used partitioning algorithm. The usual combinatorial problems associated with wrapper approaches are overcome by a Bayesian inference mechanism. On the technical side, we present an efficient optimization algorithm with guaranteed local convergence property. The only free parameter of the optimization method is selected by a resampling-based stability analysis. Experiments with Leukemia and Lymphoma datasets demonstrate that our method is able to correctly infer partitions and corresponding subsets of genes which both are relevant in a biological sense. Moreover, the frequently observed problem of ambiguities caused by different but equally high-scoring partitions is successfully overcome by the model selection method proposed.

Algorithms↗

Control of neuromuscular blockade in the presence of sensor faults.

The problem of embedding sensor fault tolerance in feedback control of neuromuscular blockade is considered. For tackling interruptions of feedback measurements, a structure based upon Bayesian inference as well as a predictive filter is proposed. This algorithm is general and can be applied to different situations. Here, it is incorporated in an adaptive automatic system for feedback control of neuromuscular blockade using continuous infusion of muscle relaxants. A significant contribution consists in the experimental clinical testing of the algorithm in patients undergoing surgery.

Algorithms↗

Variational learning for Gaussian mixture models.

This paper proposes a joint maximum likelihood and Bayesian methodology for estimating Gaussian mixture models. In Bayesian inference, the distributions of parameters are modeled, characterized by hyperparameters. In the case of Gaussian mixtures, the distributions of parameters are considered as Gaussian for the mean, Wishart for the covariance, and Dirichlet for the mixing probability. The learning task consists of estimating the hyperparameters characterizing these distributions. The integration in the parameter space is decoupled using an unsupervised variational methodology entitled variational expectation-maximization (VEM). This paper introduces a hyperparameter initialization procedure for the training algorithm. In the first stage, distributions of parameters resulting from successive runs of the expectation-maximization algorithm are formed. Afterward, maximum-likelihood estimators are applied to find appropriate initial values for the hyperparameters. The proposed initialization provides faster convergence, more accurate hyperparameter estimates, and better generalization for the VEM training algorithm. The proposed methodology is applied in blind signal detection and in color image segmentation.

Algorithms↗

Mapping Sub-National Respiratory Virus Circulation in Cambodia Using Metatranscriptomic Sequencing: A Multi-Center Hospital-Based Surveillance Study.

BACKGROUND: Genomic surveillance can guide early detection of and response to emerging epidemics. Metatranscriptomic sequencing was used to investigate sub-national respiratory virus circulation in Cambodia from 2020 to 2023. METHODS: Nasopharyngeal swabs were collected from individuals aged 2 months to 65 years with influenza-like illness in four Cambodian hospitals. Metatranscriptomic data were generated by short-read RNA sequencing. Bernoulli space-time scan statistics were used to identify temporal virus clusters. Bayesian inference of phylogenetic trees was used to compute divergence times for temporally clustered, highly represented viruses (influenza A/H3N2 and B, Betacoronavirus 1, respiratory syncytial virus [RSV] A and B), and publicly available global influenza virus genomes. RESULTS: Of 1093 individuals, 499 (45.7%) had detectable respiratory viruses belonging to 68 distinct species. Moderate (N > 20) discrete time-clusters were noted of RSV-A (37 cases), Betacoronavirus 1 (21 cases), RSV-B (22 cases), and A/H3N2 (30 cases). The posterior median of time to most recent common ancestor ranged from 0.71 years (95% HPD 0.38-1.10) for Betacoronavirus 1 and 1.31 years (95% HPD 0.60-3.20) for A/H3N2, to 2.75 years (1.82-4.26) for RSV-A and 4.79 years (2.39-7.74) for RSV-B. A/H3N2 and influenza B virus genomes mapped to clades 3C.2a1b.2a.2a and Victoria 1A.3a.2, respectively, and inter-mixed with concurrent global strains. CONCLUSIONS: Multiple respiratory viruses circulated at a sub-national level in Cambodia from 2020 to 2023 despite pandemic disruptions. Influenza virus population diversity decreased during the height of lockdown but recovered in mid-2022. Re-emerging influenza strains were distinct from historically circulating strains and clustered with contemporaneous global variants, suggesting multiple external introductions.

Humans↗

Stochastic algorithms for Markov models estimation with intermittent missing data.

Multistate Markov models are frequently used to characterize disease processes, but their estimation from longitudinal data is often hampered by complex patterns of incompleteness. Two algorithms for estimating Markov chain models in the case of intermittent missing data in longitudinal studies, a stochastic EM algorithm and the Gibbs sampler, are described. The first can be viewed as a random perturbation of the EM algorithm and is appropriate when the M step is straightforward but the E step is computationally burdensome. It leads to a good approximation of the maximum likelihood estimates. The Gibbs sampler is used for a full Bayesian inference. The performances of the two algorithms are illustrated on two simulated data sets. A motivating example concerned with the modelling of the evolution of parasitemia by Plasmodium falciparum (malaria) in a cohort of 105 young children in Cameroon is described and briefly analyzed.

Algorithms↗

Bayesian semiparametric models for survival data with a cure fraction.

We propose methods for Bayesian inference for a new class of semiparametric survival models with a cure fraction. Specifically, we propose a semiparametric cure rate model with a smoothing parameter that controls the degree of parametricity in the right tail of the survival distribution. We show that such a parameter is crucial for these kinds of models and can have an impact on the posterior estimates. Several novel properties of the proposed model are derived. In addition, we propose a class of improper noninformative priors based on this model and examine the properties of the implied posterior. Also, a class of informative priors based on historical data is proposed and its theoretical properties are investigated. A case study involving a melanoma clinical trial is discussed in detail to demonstrate the proposed methodology.

Bayes Theorem↗

Dynamic conditionally linear mixed models for longitudinal data.

We develop a new class of models, dynamic conditionally linear mixed models, for longitudinal data by decomposing the within-subject covariance matrix using a special Cholesky decomposition. Here 'dynamic' means using past responses as covariates and 'conditional linearity' means that parameters entering the model linearly may be random, but nonlinear parameters are nonrandom. This setup offers several advantages and is surprisingly similar to models obtained from the first-order linearization method applied to nonlinear mixed models. First, it allows for flexible and computationally tractable models that include a wide array of covariance structures; these structures may depend on covariates and hence may differ across subjects. This class of models includes, e.g., all standard linear mixed models, antedependence models, and Vonesh-Carter models. Second, it guarantees the fitted marginal covariance matrix of the data is positive definite. We develop methods for Bayesian inference and motivate the usefulness of these models using a series of longitudinal depression studies for which the features of these new models are well suited.

Antidepressive Agents↗

Phylogeographical patterns in the American Pacific Northwest: lessons from the arionid slug Prophysaon coeruleum.

Some of the earliest studies of phylogenetic concordance involve native plants from the Pacific Northwest where many taxa showed clear genetic breaks between southern and northern populations. To test whether similar breaks also occur in invertebrate species with low dispersal capacities, genetic data from two mitochondrial genes are assembled for individuals of the arionid slug Prophysaon coeruleum throughout the species' range. Bayesian inference revealed three major clades and a total of eight subclades. It is argued that the demographic and genealogical past of P. coeruleum has resulted in a deep and shallow phylogeographical structure. The deep structure is at least 2.6-5.9 million years old and therefore clearly predates the Pleistocene period. Superimposed on this structure is a shallow structure that is most likely less than 2 million years old and probably the result of Pleistocene perturbations. Molecular analyses revealed that the three known colour traits in P. coeruleum do not represent monophyletic groups and that they do not match the patterns of genetic structure found. It is argued that the colour traits are perhaps a response to different levels of UV-radiation. The study adds to the increasing evidence that the phylogeographical structure of some taxa is more complex than previously thought. Moreover, it shows that genealogical concordance should not be deduced from phylogeographical patterns alone but should be based on an understanding of timing and causes of historical processes that lead to those patterns.

Analysis of Variance↗

Parallels in the evolution of the two largest New and Old World seed-beetle genera (Coleoptera, Bruchidae).

This study provides the first phylogenetic analysis of a large sample of the two largest genera of seed-beetles, Acanthoscelides Schilsky and Bruchidius Schilsky, which mostly feed on legumes (Fabaceae). The goal of this study was to investigate evolutionary patterns in relation to biogeography and host-plant associations. We used three mitochondrial molecular markers and parsimony and Bayesian inference methods to reconstruct the phylogeny of 76 species. In addition, we critically reviewed host-plant records in the literature for these two bruchid genera. Our results demonstrated the existence of two major clades, one New World and one largely Old World, which generally correspond to the two genera. Yet, current classification of several species is erroneous, so that both genera as currently defined are paraphyletic. We highlighted a strong trend toward specialization (with high taxonomic conservatism in host-plant use) exhibited by the two studied genera. However, we showed the existence of several host shifts during the evolution of this group of bruchids. Our phylogenetic hypotheses and our evaluation of host-plant associations both suggest that the two genera have undergone parallel evolution, as they have independently colonized similar host plants in their respective areas of distribution. Our estimation of divergence times indicated a more ancient origin for bruchids than that suggested by the fossil records. Interestingly, the suggested timing of diversification is consistent with the hypothesis of a radiation that could have occurred contemporaneously with the diversification of their legume hosts.

Animals↗

On avoiding statistical bias in linkage-based counselling.

Using the Succession Rule of Laplace (1795) and related reasoning, this paper shows how to give unbiased counselling to patients when predictions are to be based on small samples. The recombination fraction can be regarded as a probability parameter, theta, which itself has a probability distribution between the limits of 0 and 1/2. The probability of a recombinant, P(Rec), is not numerically equal to the maximum likelihood estimate of theta, nor is it numerically equal to the maximum posterior probability estimate in Bayesian inference. Rather it is equal to the infinite sum of all possible theta values, each weighted according to its probability density p(theta) which denotes the relative probability that that theta value is the true one. The various published proposals for obtaining an unbiased estimate of theta are shown to be equivalent one to another, except for the simplifying approximations used.

Bias↗

Phylogenetic relationships among A-genome species of the genus Oryza revealed by intron sequences of four nuclear genes.

The A-genome group in Oryza consists of eight diploid species and is distributed world-wide. Here we reconstructed the phylogeny among the A-genome species based on sequences of nuclear genes and MITE (miniature inverted-repeat transposable elements) insertions. Thirty-seven accessions representing two cultivated and six wild species from the A-genome group were sampled. Introns of four nuclear single-copy genes on different chromosomes were sequenced and analysed by both maximum parsimony (MP) and Bayesian inference methods. All the species except for Oryza rufipogon and Oryza nivara formed a monophyletic group and the Australian endemic Oryza meridionalis was the earliest divergent lineage. Two subspecies of Oryza sativa (ssp. indica and ssp. japonica) formed two separate monophyletic groups, suggestive of their polyphyletic origin. Based on molecular clock approach, we estimated that the divergence of the A-genome group occurred c. 2.0 million years ago (mya) while the two subspecies (indica and japonica) separated c. 0.4 mya. Intron sequences of nuclear genes provide sufficient resolution and are informative for phylogenetic inference at lower taxonomic levels.

Base Sequence↗

Phylogenetic analysis of the myostatin gene sub-family and the differential expression of a novel member in zebrafish.

The myostatin (MSTN)-null phenotype in mammals is characterized by extreme gains in skeletal muscle mass or "double muscling" as the cytokine negatively regulates skeletal muscle growth. Recent attempts, however, to reproduce a comparable phenotype in zebrafish have failed. Several aspects of MSTN biology in the fishes differ significantly from those in mammals and at least two distinct paralogs have been identified in some species, which possibly suggests functional divergence between the different vertebrate classes or between fish paralogs. We therefore conducted a phylogenetic analysis of the entire MSTN gene sub-family. Maximum likelihood, Bayesian inference, and bootstrap analyses indicated a monophyletic distribution of all MSTN genes with two distinct fish clades: MSTN-1 and -2. These analyses further indicated that all Salmonid genes described are actually MSTN-1 orthologs and that additional MSTN-2 paralogs may be present in most, if not all, teleosts. An additional zebrafish homolog was identified by BLAST searches of the zebrafish Hierarchical Tets Generation System database and was subsequently cloned. Comparative sequence analysis of both genes (zebrafish MSTN (zfMSTN)-1 and -2) revealed many differences, primarily within the latency-associated peptide regions, but also within the bioactive domains. The 2-kb promoter region of zfMSTN-2 contained many putative cis regulatory elements that are active during myogenesis, but are lacking in the zfMSTN-1 promoter. In fact, zfMSTN-2 expression was limited to the early stages of somitogenesis, whereas zfMSTN-1 was expressed throughout embryogenesis. These data suggest that zfMSTN-2 may be more closely associated with skeletal muscle growth and development. They also resolve the previous ambiguity in classification of fish MSTN genes.

Animals↗

A molecular phylogenetic investigation of zoothamnium (ciliophora, peritrichia, sessilida).

The gene coding for 18S small subunit ribosomal RNA (ssu rRNA) was sequenced in seven free-living, marine species of the sessiline peritrich genus Zoothamnium. These were Zoothamnium niveum, Zoothamnium alternans, Zoothamnium pelagicum, and four unidentified species. The ssu rRNA gene also was sequenced in Vorticella convallaria, Vorticella microstoma, and in an unidentified, freshwater species of Vorticella. Phylogenetic trees were constructed using these new sequences to test a previously published phylogenetic association between Zoothamnium arbuscula, currently in the family Zoothamniidae, and peritrichs in the family Vorticellidae. Trees constructed by means of neighbor-joining, maximum parsimony, maximum likelihood, and Bayesian inference methods all had similar topologies. The seven new sequences of Zoothamnium species grouped into three well-supported clades, each of which contained a diversity of morphological types. The three clades formed a poorly supported, larger clade that was deeply divergent from Z. arbuscula, which remained more closely associated with vorticellid peritrichs. It is apparent that Zoothamnium is a richly diverse genus and that a much more intensive investigation, involving both morphological and molecular data and a wider selection of species, will be necessary to resolve its phylogeny. A greater amount of molecular diversity than is predicted by morphological data exists within all major clades of sessiline peritrichs that have been included in molecular phylogenies, indicating that characteristics of stalk and peristomial structure traditionally used to differentiate taxa at the generic level and above may not be uniformly reliable.

Animals↗

Using the acetabulum to estimate age at death of adult males.

The acetabular region is often present and adequately preserved in adult human skeletal remains. Close morphological examination of the 242 left male os coxae from the identified collection of Coimbra (Portugal) has enabled the recognition of seven variables that can be used to estimate age at death. This paper describes these variables and argues their appropriateness by analyzing the correlation between these criteria and the age, the intra- and interobserver consistence, and the accuracy in age prediction using Bayesian inference to estimate age of identified specimens. Results show significant close correlation between the acetabular criteria and age, nonsignificant differences in intra- and interobserver test, and 89% accuracy in Bayes prediction. Obtained estimated age of the specimens had similar accuracy in all ages. These results indicate that these seven variables, based on the acetabular area, are potentially useful to estimate age at death for adult specimens.

Acetabulum↗

Phytophthora polonica, a new species isolated from declining Alnus glutinosa stands in Poland.

In a survey of Phytophthora associated with alder decline in Poland, several isolates of a homothallic Phytophthora sp., which could not be assigned to other taxa including Phytophthora alni subspecies, were consistently recovered from rhizosphere soil samples. Their morphology and pathogenicity, as well as sequence data for three nuclear regions (internal transcribed spacer rDNA, elongation factor-1alpha and beta-tubulin) and a coding mitochondrial DNA region (nadh1), were examined. The new Phytophthora species is characterized by the moderate to slow growth rate of its colony in carrot agar at 20 degrees C, high optimal (c. 30 degrees C) and maximum (c. 38 degrees C) growth temperatures, formation of catenulate, often lateral, hyphal swellings, large chlamydospores in agar media and in soil extract, persistent, ovoid to ellipsoid nonpapillate sporangia and large oogonia with paragynous and sometimes amphigynous antheridia. Phytophthora polonica was slightly pathogenic to alder twigs and not pathogenic to trunks of several tree species. In a phylogenetic analysis using either Bayesian inference or maximum likelihood methods, P. polonica falls in clade 8 'sensu Kroon et al. (2004)' of Phytophthora.

Alnus↗

Contrasting Patterns of Connectivity Between Populations of Euphotic and Mesophotic Hydroids in Reunion Island Support the Deep Reef Refuge Hypothesis.

In the context of coral reef decline, mesophotic coral ecosystems (MCEs, 30-150 m) offer hope for the recovery of degraded euphotic reefs. The Deep Reef Refuge Hypothesis (DRRH) postulates the potential of mesophotic reefs to reseed euphotic reefs. This hypothesis needs to be further tested by estimating connectivity along the depth gradient. Mesophotic data are lacking worldwide, particularly in the southwestern Indian Ocean (SWIO). Here, using a total of 2218 samples collected at depths ranging from 10 to 103 m, we estimated the connectivity of 7 hydroid species sampled at euphotic, upper, and lower mesophotic depths around Reunion Island using a multi-species comparative framework. Population genetic analyses using 8-17 microsatellite markers per species (80 markers in total) as well as Bayesian inference were performed to estimate population structure and contemporary migration rates to highlight connectivity patterns and directionality of gene flow between depths. The results revealed three main genetic patterns depending on the species: a horizontal stepping stone pattern between areas around the island, a vertical stepping stone pattern between adjacent depths, and a quasi-panmictic pattern. Each species showed some specificity within these patterns, but overall, at least 4 of the 7 species support the assumption of vertical connectivity from the Deep Reef Refuge Hypothesis, highlighting the importance of studying multiple species. The existence of vertical connectivity between euphotic and mesophotic depths in the southwestern Indian Ocean confirms the importance of mesophotic coral ecosystems for conservation efforts and our global understanding of coral reef ecosystem dynamics.

Animals↗

Hyperspectral techniques in analysis of oral dosage forms.

Pharmaceutical oral dosage forms are used in this paper to test the sensitivity and spatial resolution of hyperspectral imaging instruments. The first experiment tested the hypothesis that a near-infrared (IR) tunable diode-based remote sensing system is capable of monitoring degradation of hard gelatin capsules at a relatively long distance (0.5 km). Spectra from the capsules were used to differentiate among capsules exposed to an atmosphere containing 150 ppb formaldehyde for 0, 2, 4, and 8 h. Robust median-based principal component regression with Bayesian inference was employed for outlier detection. The second experiment tested the hypothesis that near-IR imaging spectrometry of tablets permits the identification and composition of multiple individual tablets to be determined simultaneously. A near-IR camera was used to collect thousands of spectra simultaneously from a field of blister-packaged tablets. The number of tablets that a typical near-IR camera can currently analyze simultaneously was estimated to be approximately 1300. The bootstrap error-adjusted single-sample technique chemometric-imaging algorithm was used to draw probability-density contour plots that revealed tablet composition. The single-capsule analysis provides an indication of how far apart the sample and instrumentation can be and still maintain adequate signal-to-noise ratio (S/N), while the multiple-tablet imaging experiment gives an indication of how many samples can be analyzed simultaneously while maintaining an adequate S/N and pixel coverage on each sample.

Administration, Oral↗

Uncertainty analysis in matched-field geoacoustic inversions.

Quantifying uncertainty for parameter estimates obtained from matched-field geoacoustic inversions using a Bayesian approach requires estimation of the uncertainties in the data due to ambient noise as well as modeling errors. In this study, the variance parameter of the Gaussian error model, hereafter called error variance, is assumed to describe the data uncertainty. In practice, this parameter is not known a priori, and choosing a particular value is often difficult. Hence, to account for the uncertainty in error variance, several methods are introduced for implementing both the full and empirical Bayesian approaches. A full Bayesian approach that permits uncertainty of the error variance to propagate through the parameter estimation processes is a natural way of incorporating the uncertainty of error variance. Due to the large number of unknown parameters in the full Bayesian uncertainty analysis, an alternative, the empirical Bayesian approach, is developed, in which the posterior distributions of model parameters are conditioned on a point estimate of the error variance. Comparisons between the full and empirical Bayesian inferences of model parameters are presented using both synthetic and experimental data.

Journal Article↗