Comparison of Bayesian and numerical optimization-based diet estimation on herbivorous zooplankton.
Ontology highlight
ABSTRACT:
SUBMITTER: Litmanen JJ
PROVIDER: S-EPMC7333965 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Terrestrial organic matter inputs have long been thought to play an important role in aquatic food web dynamics. Results from recent whole lake (13)C addition experiments suggest terrestrial particulate organic carbon (t-POC) inputs account for a disproportionate portion of zooplankton production. For example, several studies concluded that although t-POC only represented approximately 20% of the flux of particulate carbon available to herbivorous zooplankton, this food source accounted for approximately 50% of the C incorporated by zooplankton. We tested the direct dietary impact of t-POC (from the leaves of riparian vegetation) and various phytoplankton on Daphnia magna somatic growth, reproduction, growth efficiency, and lipid composition. By itself, t-POC was a very poor quality resource compared to cryptophytes, diatoms, and chlorophytes, but t-POC had similar food quality compared to cyanobacteria. Small additions of high quality Cryptomonas ozolinii to t-POC-dominated diets greatly increased Daphnia growth and reproduction. When offered alone, t-POC resulted in a Daphnia growth efficiency of 5 +/- 1%, whereas 100% Cryptomonas and Scenedesmus obliquus diets resulted in growth efficiencies of 46 +/- 8% (+/- SD) and 36 +/- 3%, respectively. When offered in a 50:50 mixed diet with Cryptomonas or Scenedesmus, the t-POC fraction resulted in a partial growth efficiency of 22 +/- 9% and 15 +/- 6%, respectively. Daphnia that obtained 80% of their available food from t-POC assimilated 84% of their fatty acids from the phytoplankton component of their diet. Overall, our results suggest Daphnia selectively allocate phytoplankton-derived POC and lipids to enhance somatic growth and reproduction, while t-POC makes a minor contribution to zooplankton production.
| S-EPMC2795543 | biostudies-literature
Project description:Due to the high occupational pressure suffered by intensive care units (ICUs), a correct estimation of the patients' length of stay (LoS) in the ICU is of great interest to predict possible situations of collapse, to help healthcare personnel to select appropriate treatment options and to predict patients' conditions. There has been a high amount of data collected by biomedical sensors during the continuous monitoring process of patients in the ICU, so the use of artificial intelligence techniques in automatic LoS estimation would improve patients' care and facilitate the work of healthcare personnel. In this work, a novel methodology to estimate the LoS using data of the first 24 h in the ICU is presented. To achieve this, XGBoost, one of the most popular and efficient state-of-the-art algorithms, is used as an estimator model, and its performance is optimized both from computational and precision viewpoints using Bayesian techniques. For this optimization, a novel two-step approach is presented. The methodology was carefully designed to execute codes on a high-performance computing system based on graphics processing units, which considerably reduces the execution time. The algorithm scalability is analyzed. With the proposed methodology, the best set of XGBoost hyperparameters are identified, estimating LoS with a MAE of 2.529 days, improving the results reported in the current state of the art and probing the validity and utility of the proposed approach.
| S-EPMC9919941 | biostudies-literature
Project description:A two-layer model based on the integrated form of Richards' equation (RE) was recently developed to simulate the soil water movement in the roots layer and the vadose zone with a relatively shallow and dynamic water table. The model simulates thickness-averaged volumetric water content and matric suction as opposed to point values and was numerically verified for three soil textures using HYDRUS as a benchmark. However, the strengths and limitations of the two-layer model and its performance in stratified soils and under actual field conditions have not been tested. This study further examined the two-layer model using two numerical verification experiments and, most importantly, tested its performance at site-level under actual, highly variable hydroclimate conditions. Moreover, model parameters were estimated and uncertainty and sources of errors were quantified using a Bayesian framework. First, the two-layer model was evaluated for 231 soil textures under varying soil layer thicknesses with a uniform soil profile. Second, the two-layer model was assessed for stratified conditions where the top and bottom soil layers have contrasting hydraulic conductivities. The model was evaluated by comparing soil moisture and flux estimates to those from the HYDRUS model. Last, a case study of model application using data from a Soil Climate Analysis Network (SCAN) site was presented. Bayesian Monte Carlo (BMC) method was implemented for model calibration and quantifying sources of uncertainty under real hydroclimate and soil conditions. For a homogeneous soil profile, the two-layer model generally had excellent performance in estimating volumetric water content and fluxes, while the model performance slightly declined with increasing layer thickness and coarser textured soils. The model configurations regarding layer thicknesses and soil textures that generate accurate soil moisture and flux estimations were further suggested. With the two layers of contrasting permeability, model-simulated soil moisture contents and fluxes agreed well with those computed by HYDRUS, indicating that the two-layer model accurately handles the water flow dynamics around the layer interface. In the field application, given the highly variable hydroclimate conditions, the two-layer model combined with the BMC method showed good agreement with the observed average soil moisture of the root zone and the vadose zone below (RMSE <0.021 during calibration and <0.023 during validation periods). The contribution of parametric uncertainty to the total model uncertainty was too small compared to other sources. The numerical tests and the site level application showed that the two-layer model can reliably simulate thickness-averaged soil moisture and estimate fluxes in the vadose zone under various soil and hydroclimate conditions. Results also indicated that the BMC method could be a robust framework for vadose zone hydraulic parameters identification and model uncertainty estimation.
| S-EPMC10266505 | biostudies-literature
Project description:Quantitative kinetic parameters derived from dynamic contrast-enhanced (DCE) data are dependent on signal measurement quality and choice of pharmacokinetic model. However, the fundamental optimization analysis method is equally important and its impact on pharmacokinetic parameters has been mostly overlooked. We examine the effects of those choices on accuracy and performance of parameter estimation using both computer processing unit and graphical processing unit (GPU) numerical optimization implementations and evaluate the improvements offered by a novel optimization approach. A test framework was developed where experimentally derived population-average arterial input function and randomly sampled parameter sets {Ktrans , Kep , Vb , ?} were used to generate known tissue curves. Five numerical optimization algorithms were evaluated: sequential quadratic programming, downhill simplex (Nelder-Mead), pattern search, simulated annealing, and differential evolution. This was combined with various objective function implementation details: delay approximation, discretization and varying sampling rates. Then, impact of noise and CPU/GPU implementation was tested for speed and accuracy. Finally, the optimal method was compared to conventional implementation as applied to clinical DCE computed tomography. Nelder-Mead, differential evolution and sequential quadratic programming produced good results on clean and noisy input data outperforming simulated annealing and pattern search in terms of speed and accuracy in the respective order of 10-8%, 10-7%, and ×10-6%). A novel approach for DCE numerical optimization (infinite impulse response with fractional delay approximation) was implemented on GPU for speed increase of at least 2 orders of magnitude. Applied to clinical data, the magnitude of overall parameter error was <10%.
| S-EPMC6403032 | biostudies-literature
Project description:We compare three state-of-the-art Bayesian inference methods for the estimation of the unknown parameters in a stochastic model of a genetic network. In particular, we introduce a stochastic version of the paradigmatic synthetic multicellular clock model proposed by Ullner et al., 2007. By introducing dynamical noise in the model and assuming that the partial observations of the system are contaminated by additive noise, we enable a principled mechanism to represent experimental uncertainties in the synthesis of the multicellular system and pave the way for the design of probabilistic methods for the estimation of any unknowns in the model. Within this setup, we tackle the Bayesian estimation of a subset of the model parameters. Specifically, we compare three Monte Carlo based numerical methods for the approximation of the posterior probability density function of the unknown parameters given a set of partial and noisy observations of the system. The schemes we assess are the particle Metropolis-Hastings (PMH) algorithm, the nonlinear population Monte Carlo (NPMC) method and the approximate Bayesian computation sequential Monte Carlo (ABC-SMC) scheme. We present an extensive numerical simulation study, which shows that while the three techniques can effectively solve the problem there are significant differences both in estimation accuracy and computational efficiency.
| S-EPMC5552360 | biostudies-literature
Project description:Quantitative analysis of stable isotopes (SI) and, more recently, fatty acid profiles (FAP) are useful and complementary tools for estimating the relative contribution of different prey items in the diet of a predator. The combination of these two approaches, however, has thus far been limited and qualitative. We propose a mixing model for FAP that follows the Bayesian machinery employed in state-of-the-art mixing models for SI. This framework provides both point estimates and probability distributions for individual and population level diet proportions. Where fat content and conversion coefficients are available, they can be used to improve diet estimates. This model can be explicitly integrated with analogous models for SI to increase resolution and clarify predator-prey relationships. We apply our model to simulated data and an experimental dataset that allows us to illustrate modeling strategies and demonstrate model performance. Our methods are provided as an open source software package for the statistical computing environment R.
| S-EPMC4411522 | biostudies-literature
Project description:Benchmark dose (BMD) modeling is an important step in human health risk assessment and is used as the default approach to identify the point of departure for risk assessment. A probabilistic framework for dose-response assessment has been proposed and advocated by various institutions and organizations; therefore, a reliable tool is needed to provide distributional estimates for BMD and other important quantities in dose-response assessment.We developed an online system for Bayesian BMD (BBMD) estimation and compared results from this software with U.S. Environmental Protection Agency's (EPA's) Benchmark Dose Software (BMDS).The system is built on a Bayesian framework featuring the application of Markov chain Monte Carlo (MCMC) sampling for model parameter estimation and BMD calculation, which makes the BBMD system fundamentally different from the currently prevailing BMD software packages. In addition to estimating the traditional BMDs for dichotomous and continuous data, the developed system is also capable of computing model-averaged BMD estimates.A total of 518 dichotomous and 108 continuous data sets extracted from the U.S. EPA's Integrated Risk Information System (IRIS) database (and similar databases) were used as testing data to compare the estimates from the BBMD and BMDS programs. The results suggest that the BBMD system may outperform the BMDS program in a number of aspects, including fewer failed BMD and BMDL calculations and estimates.The BBMD system is a useful alternative tool for estimating BMD with additional functionalities for BMD analysis based on most recent research. Most importantly, the BBMD has the potential to incorporate prior information to make dose-response modeling more reliable and can provide distributional estimates for important quantities in dose-response assessment, which greatly facilitates the current trend for probabilistic risk assessment. https://doi.org/10.1289/EHP1289.
| S-EPMC6014690 | biostudies-other
Project description:BackgroundVarious methods exist for statistical inference about a prevalence that consider misclassifications due to an imperfect diagnostic test. However, traditional methods are known to suffer from truncation of the prevalence estimate and the confidence intervals constructed around the point estimate, as well as from under-performance of the confidence intervals' coverage.MethodsIn this study, we used simulated data sets to validate a Bayesian prevalence estimation method and compare its performance to frequentist methods, i.e. the Rogan-Gladen estimate for prevalence, RGE, in combination with several methods of confidence interval construction. Our performance measures are (i) error distribution of the point estimate against the simulated true prevalence and (ii) coverage and length of the confidence interval, or credible interval in the case of the Bayesian method.ResultsAcross all data sets, the Bayesian point estimate and the RGE produced similar error distributions with slight advantages of the former over the latter. In addition, the Bayesian estimate did not suffer from the RGE's truncation problem at zero or unity. With respect to coverage performance of the confidence and credible intervals, all of the traditional frequentist methods exhibited strong under-coverage, whereas the Bayesian credible interval as well as a newly developed frequentist method by Lang and Reiczigel performed as desired, with the Bayesian method having a very slight advantage in terms of interval length.ConclusionThe Bayesian prevalence estimation method should be prefered over traditional frequentist methods. An acceptable alternative is to combine the Rogan-Gladen point estimate with the Lang-Reiczigel confidence interval.
| S-EPMC7370479 | biostudies-literature
Project description:Number magnitude estimation has been investigated over the last decades using different tasks including non-symbolic numerosity but also number line estimation tasks. Recently, a bi-directional mapping process was suggested for numerosity estimation accounting for underestimation in a perception version of the task (i.e., indicating the number of non-symbolic dots in a set) and overestimation in the corresponding production task (i.e., produce the number of dots indicated by a symbolic number). In the present study, we evaluated the generalizability of these estimation biases in perception and production tasks to bounded and unbounded number line estimation. Importantly, target numbers were underestimated/overestimated by participants in the perception/production version of numerosity estimation as well as unbounded number line estimation. However, this pattern was reversed for bounded number line estimation. Thereby, the present data indicate a conceptual similarity of unbounded number line estimation and the established non-symbolic numerosity estimation task as a measure of numerical estimation. Accordingly, this corroborates the notion that unbounded number line estimation may reflect a purer measure of number magnitude representation than the bounded task version. Furthermore, our findings strengthen the bi-directional mapping hypothesis for numerical estimation by providing evidence for its generalizability to unbounded number line estimation for the first time.
| S-EPMC6417782 | biostudies-literature
Project description:BackgroundLaurasiatheria contains taxa with diverse diets, while the molecular basis and evolutionary history underlying their dietary diversification are less clear.ResultsIn this study, we used the recently developed molecular phyloecological approach to examine the adaptive evolution of digestive system-related genes across both carnivorous and herbivorous mammals within Laurasiatheria. Our results show an intensified selection of fat and/or protein utilization across all examined carnivorous lineages, which is consistent with their high-protein and high-fat diets. Intriguingly, for herbivorous lineages (ungulates), which have a high-carbohydrate diet, they show a similar selection pattern as that of carnivorous lineages. Our results suggest that for the ungulates, which have a specialized digestive system, the selection intensity of their digestive system-related genes does not necessarily reflect loads of the nutrient components in their diets but appears to be positively related to the loads of the nutrient components that are capable of being directly utilized by the herbivores themselves. Based on these findings, we reconstructed the dietary evolution within Laurasiatheria, and our results reveal the dominant carnivory during the early diversification of Laurasiatheria. In particular, our results suggest that the ancestral bats and the common ancestor of ruminants and cetaceans may be carnivorous as well. We also found evidence of the convergent evolution of one fat utilization-related gene, APOB, across carnivorous taxa.ConclusionsOur molecular phyloecological results suggest that digestive system-related genes can be used to determine the molecular basis of diet differentiations and to reconstruct ancestral diets.
| S-EPMC9210794 | biostudies-literature