Mlegp: statistical analysis for computer models of biological systems using R.
Ontology highlight
ABSTRACT:
SUBMITTER: Dancik GM
PROVIDER: S-EPMC2732217 | biostudies-literature | 2008 Sep
REPOSITORIES: biostudies-literature
Publications
mlegp: statistical analysis for computer models of biological systems using R.
Dancik Garrett M GM Dorman Karin S KS
Bioinformatics (Oxford, England) 20080717 17
<h4>Unlabelled</h4>Gaussian processes (GPs) are flexible statistical models commonly used for predicting output from complex computer codes. As such, GPs are well suited for the analysis of computer models of biological systems, which have been traditionally difficult to analyze due to their high-dimensional, non-linear and resource-intensive nature. We describe an R package, mlegp, that fits GPs to computer model outputs and performs sensitivity analysis to identify and characterize the effects ...[more]
Similar Datasets
Project description:Standard techniques for studying biological systems largely focus on their dynamical or, more recently, their informational properties, usually taking either a reductionist or holistic perspective. Yet, studying only individual system elements or the dynamics of the system as a whole disregards the organizational structure of the system-whether there are subsets of elements with joint causes or effects, and whether the system is strongly integrated or composed of several loosely interacting components. Integrated information theory offers a theoretical framework to (1) investigate the compositional cause-effect structure of a system and to (2) identify causal borders of highly integrated elements comprising local maxima of intrinsic cause-effect power. Here we apply this comprehensive causal analysis to a Boolean network model of the fission yeast (Schizosaccharomyces pombe) cell cycle. We demonstrate that this biological model features a non-trivial causal architecture, whose discovery may provide insights about the real cell cycle that could not be gained from holistic or reductionist approaches. We also show how some specific properties of this underlying causal architecture relate to the biological notion of autonomy. Ultimately, we suggest that analysing the causal organization of a system, including key features like intrinsic control and stable causal borders, should prove relevant for distinguishing life from non-life, and thus could also illuminate the origin of life problem.This article is part of the themed issue 'Reconceptualizing the origins of life'.
| S-EPMC5686412 | biostudies-literature
Project description:In this paper, we propose to use a linear system with switching methodology for description and analysis of complex biological systems. We show advantages of the proposed approach over the one usually used, which is based on ODE. We propose the detailed methodology of a full analysis of developed models, including analytical determination of the location and type of equilibrium points, finding an analytical solution, stability and bifurcation analysis. We illustrate the above with the example of the well-known p53 signalling pathway comparing the results with the results of a nonlinear, ODE-based version of the proposed model. The complex methodology proposed by us, especially due to the definition of model structure, which is easy to understand for biologists and medics, may be a bridge for closer cooperation between them and engineers in the future.
| S-EPMC7743982 | biostudies-literature
Project description:The linear "Ising" model, which has been around for nearly a century, treats the behavior of linear arrays of repetitive, interacting subunits. Linear "repeat-proteins" have only been described in the last decade or so, and their folding energies have only been characterized very recently. Owing to their repetitive structures, linear repeat-proteins are particularly well suited for analysis by the nearest-neighbor Ising formalism. After briefly describing the historical origins and applications of the Ising model to biopolymers, and introducing repeat protein structure, this chapter will focus on the application of the linear Ising model to repeat proteins. When applied to homopolymers, the model can be represented and applied in a fairly simplified form. When applied to heteropolymers, where differences in energies among individual subunits (i.e. repeats) must be included, some (but not all) of this simplicity is lost. Derivations of the linear Ising model for both homopolymer and heteropolymer repeat-proteins will be presented. With the increased complexity required for analysis of heteropolymeric repeat proteins, the ability to resolve different energy terms from experimental data can be compromised. Thus, a simple matrix approach will be developed to help inform on the degree to which different thermodynamic parameters can be extracted from a particular set of unfolding curves. Finally, we will describe the application of these models to analyze repeat-protein folding equilibria, focusing on simplified repeat proteins based on "consensus" sequence information.
| S-EPMC2967780 | biostudies-literature
Project description:Translational Statistics proposes to promote the use of Statistics within research and improve the communication of statistical findings in an accurate and accessible manner to diverse audiences. When statistical models become more complex, it becomes harder to evaluate the role of explanatory variables on the response. For example, the interpretation and communication of the effect of predictors in regression models where interactions or smoothing splines are included can be challenging. Informative graphical representations of statistical models play a critical translational role; static nomograms are one such useful tool to visualise statistical models. In this paper, we propose the use of dynamic nomogram as a translational tool which can accommodate models of increased complexity. In theory, all models appearing in the literature could be accompanied by the corresponding dynamic nomogram to translate models in an informative manner. The R package presented will facilitate this communication for a variety of linear and non-linear models.
| S-EPMC6857916 | biostudies-literature
Project description:Distribution networks-from vasculature to urban transportation pathways-are spatially embedded networks that must route resources efficiently in the face of pressures induced by the costs of building and maintaining network infrastructure. Such requirements are thought to constrain the topological and spatial organization of these systems, but at the same time, different kinds of distribution networks may exhibit variable architectural features within those general constraints. In this study, we use methods from network science to compare and contrast two classes of biological transport networks: mycelial fungi and vasculature from the surface of rodent brains. These systems differ in terms of their growth and transport mechanisms, as well as the environments in which they typically exist. Though both types of networks have been studied independently, the goal of this study is to quantify similarities and differences in their network designs. We begin by characterizing the structural backbone of these systems with a collection of measures that assess various kinds of network organization across topological and spatial scales, ranging from measures of loop density, to those that quantify connected pathways between different network regions, and hierarchical organization. Most importantly, we next carry out a network analysis that directly considers the spatial embedding and properties especially relevant to the function of distribution systems. We find that although both the vasculature and mycelia are highly constrained planar networks, there are clear distinctions in how they balance tradeoffs in network measures of wiring length, efficiency, and robustness. While the vasculature appears well organized for low cost, but relatively high efficiency, the mycelia tend to form more expensive but in turn more robust networks. As a whole, this work demonstrates the utility of network-based methods to identify both common features and variations in the network structure of different classes of biological transport systems.
| S-EPMC6145589 | biostudies-literature
Project description:The purpose of this work is to present a computer assisted diagnostic tool for radiologists in their diagnosis of Alzheimer's disease. A statistical likelihood-ratio procedure from signal detection theory was implemented in the detection of Alzheimer's disease. The probability density functions of the likelihood ratio were constructed by using medial temporal lobe (MTL) volumes of patients with Alzheimer's disease (AD) and normal controls (NC). The volumes of MTL as well as other anatomical regions of the brains were calculated by the FreeSurfer software using T1 weighted MRI images. The MRI images of AD and NC were downloaded from the database of Alzheimer's disease neuroimaging initiative (ADNI). A separate dataset of minimal interval resonance imaging in Alzheimer's disease (MIRIAD) was used for diagnostic testing. A sensitivity of 89.1% and specificity of 87.0% were achieved for the MIRIAD dataset which are better than the 85% sensitivity and specificity achieved by the best radiologists without input of other patient information.
| S-EPMC9937475 | biostudies-literature
Project description:Rare events such as genetic mutations or cell-cell interactions are important contributors to dynamics in complex biological systems, eg, in drug-resistant infections. Computational approaches can help analyze rare events that are difficult to study experimentally. However, analyzing the frequency and dynamics of rare events in computational models can also be challenging due to high computational resource demands, especially for high-fidelity stochastic computational models. To facilitate analysis of rare events in complex biological systems, we present a multifidelity analysis approach that uses medium-fidelity analysis (Monte Carlo simulations) and/or low-fidelity analysis (Markov chain models) to analyze high-fidelity stochastic model results. Medium-fidelity analysis can produce large numbers of possible rare event trajectories for a single high-fidelity model simulation. This allows prediction of both rare event dynamics and probability distributions at much lower frequencies than high-fidelity models. Low-fidelity analysis can calculate probability distributions for rare events over time for any frequency by updating the probabilities of the rare event state space after each discrete event of the high-fidelity model. To validate the approach, we apply multifidelity analysis to a high-fidelity model of tuberculosis disease. We validate the method against high-fidelity model results and illustrate the application of multifidelity analysis in predicting rare event trajectories, performing sensitivity analyses and extrapolating predictions to very low frequencies in complex systems. We believe that our approach will complement ongoing efforts to enable accurate prediction of rare event dynamics in high-fidelity computational models.
| S-EPMC6077899 | biostudies-literature
Project description:four biological sludge systems Raw sequence reads
| PRJNA532830 | ENA
Project description:BackgroundGrading of prostatic adenocarcinoma is based on the Gleason scoring system and the more recently established prognostic grade groups. Typically, prostate cancer grading is performed by pathologists based on the morphology of the tumor on hematoxylin and eosin (H and E) slides. In this study, we investigated the histopathological image features with various survival models and attempted to study their correlations.MethodsThree texture methods (speeded-up robust features, histogram of oriented gradient, and local binary pattern) and two convolutional neural network (CNN)-based methods were applied to quantify histopathological image features. Five survival models were assessed on those image features in the context with other prostate clinical prognostic factors, including primary and secondary Gleason patterns, prostate-specific antigen levels, age, and clinical tumor stages.ResultsBased on statistical comparisons among different image features with survival models, image features from CNN-based method with a recurrent neural network called CNN-long-short-term memory provided the highest hazard ratio of prostate cancer recurrence under Cox regression with an elastic net penalty.ConclusionsThis approach outperformed the other image quantification methods listed above. Using this approach, patient outcomes were highly correlated with the histopathological image features of the tissue samples. In future studies, we plan to investigate the potential use of this approach for predicting recurrence in a wider range of cancer types.
| S-EPMC6788183 | biostudies-literature
Project description:Dengue virus (DENV) infects hundreds of millions of people annually, yet there is only a limited knowledge of the host immune response to dengue. Here, we used a systems biological approach to perform a detailed analysis of the innate immune response to DENV infection in the whole blood samples of acutely infected humans in Bangkok, Thailand. Transcriptomic analysis revealed that genes encoding pro-inflammatory mediators and type I IFN related proteins, were associated with high levels of virus during the first few days of infection. Individuals with low or negative viremia at the late stage of fever were enriched with genes associated with pathways involved in cell cycle, proliferation, cell metabolism and translational control. Meta-analysis showed significant enrichment in genes specific for innate cells (monocytes, macrophages and DCs) in the specimens with high VL and enrichment in genes specific for NK cells, CD4+ and CD8+ T cells as well as B cells in specimens with low VL. Furthermore, flow cytometric analysis revealed an expansion in the numbers of CD14+CD16+ monocytes and depletion of CD14dimCD16++ cells and BDCA-1+ myeloid DC in blood. Consistent with this, in a non-human primate model, infection with DENV boosted the numbers of CD14+CD16+ monocytes in the blood and in secondary lymphoid organs. In vitro, freshly isolated blood monocytes infected with DENV up regulated CD16 and mediated robust differentiation of resting B cells to CD27++CD38++ plasmablasts and IgG and IgM secretion. Taken together, these data provide a detailed picture of the innate response to dengue infection in humans, and highlight an unappreciated role for CD14+CD16+ monocytes in promoting the differentiation of plasmablasts and mediating antibody response to DENV.
2014-06-27 | GSE51808 | GEO