Semiparametric modelling of two-component mixtures with stochastic dominance.
Ontology highlight
ABSTRACT: Supplementary information
The online version contains supplementary material available at 10.1007/s10463-022-00835-5.
SUBMITTER: Wu J
PROVIDER: S-EPMC9127045 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:In longitudinal biomedical studies, there is often interest in the rate functions, which describe the functional rates of change of biomarker profiles. This paper proposes a semiparametric approach to model these functions as the realizations of stochastic processes defined by stochastic differential equations. These processes are dependent on the covariates of interest and vary around a specified parametric function. An efficient Markov chain Monte Carlo algorithm is developed for inference. The proposed method is compared with several existing methods in terms of goodness-of-fit and more importantly the ability to forecast future functional data in a simulation study. The proposed methodology is applied to prostate-specific antigen profiles for illustration. Supplementary materials for this paper are available online.
| S-EPMC3298426 | biostudies-literature
Project description:Studying time-dependent exposure mixtures has gained increasing attentions in environmental health research. When a scalar outcome is of interest, distributed lag (DL) models have been employed to characterize the exposures effects distributed over time on the mean of final outcome. However, there is a methodological gap on investigating time-dependent exposure mixtures with different quantiles of outcome. In this paper, we introduce semiparametric partial-linear single-index (PLSI) DL quantile regression, which can describe the DL effects of time-dependent exposure mixtures on different quantiles of outcome and identify susceptible periods of exposures. We consider two time-dependent exposure settings: discrete and functional, when exposures are measured in a small number of time points and at dense time grids, respectively. Spline techniques are used to approximate the nonparametric DL function and single-index link function, and a profile estimation algorithm is proposed. Through extensive simulations, we demonstrate the performance and value of our proposed models and inference procedures. We further apply the proposed methods to study the effects of maternal exposures to ambient air pollutants of fine particulate and nitrogen dioxide on birth weight in New York University Children's Health and Environment Study (NYU CHES).
| S-EPMC10718172 | biostudies-literature
Project description:Membrane fouling during particle filtration occurs through a variety of mechanisms, including internal pore clogging by contaminants, coverage of pore entrances and deposition on the membrane surface. In this paper, we present an efficient method for modelling the behaviour of a filter, which accounts for different retention mechanisms, particle sizes and membrane geometries. The membrane is assumed to be composed of a series of, possibly interconnected, pores. The central feature is a conductivity function, which describes the blockage of each individual pore as particles arrive, which is coupled with a mechanism to account for the stochastic nature of the arrival times of particles at the pore. The result is a system of ordinary differential equations based on the pore-level interactions. We demonstrate how our model can accurately describe a wide range of filtration scenarios. Specifically, we consider a case where blocking via multiple mechanisms can occur simultaneously, which have previously required the study through individual models; the filtration of a combination of small and large particles by a track-etched membrane and particle separation using interconnected pore networks. The model is significantly faster than comparable stochastic simulations for small networks, enabling its use as a tool for efficient future simulations.
| S-EPMC5415697 | biostudies-literature
Project description:Several researchers have described two-part models with patient-specific stochastic processes for analysing longitudinal semicontinuous data. In theory, such models can offer greater flexibility than the standard two-part model with patient-specific random effects. However, in practice, the high dimensional integrations involved in the marginal likelihood (i.e. integrated over the stochastic processes) significantly complicates model fitting. Thus, non-standard computationally intensive procedures based on simulating the marginal likelihood have so far only been proposed. In this paper, we describe an efficient method of implementation by demonstrating how the high dimensional integrations involved in the marginal likelihood can be computed efficiently. Specifically, by using a property of the multivariate normal distribution and the standard marginal cumulative distribution function identity, we transform the marginal likelihood so that the high dimensional integrations are contained in the cumulative distribution function of a multivariate normal distribution, which can then be efficiently evaluated. Hence, maximum likelihood estimation can be used to obtain parameter estimates and asymptotic standard errors (from the observed information matrix) of model parameters. We describe our proposed efficient implementation procedure for the standard two-part model parameterisation and when it is of interest to directly model the overall marginal mean. The methodology is applied on a psoriatic arthritis data set concerning functional disability.
| S-EPMC5723155 | biostudies-literature
Project description:BackgroundThe aim of this paper is to develop a flexible model for analysis of quantitative trait loci (QTL) in outbred line crosses, which includes both additive and dominance effects. Our flexible intercross analysis (FIA) model accounts for QTL that are not fixed within founder lines and is based on the variance component framework. Genome scans with FIA are performed using a score statistic, which does not require variance component estimation.ResultsSimulations of a pedigree with 800 F2 individuals showed that the power of FIA including both additive and dominance effects was almost 50% for a QTL with equal allele frequencies in both lines with complete dominance and a moderate effect, whereas the power of a traditional regression model was equal to the chosen significance value of 5%. The power of FIA without dominance effects included in the model was close to those obtained for FIA with dominance for all simulated cases except for QTL with overdominant effects. A genome-wide linkage analysis of experimental data from an F2 intercross between Red Jungle Fowl and White Leghorn was performed with both additive and dominance effects included in FIA. The score values for chicken body weight at 200 days of age were similar to those obtained in FIA analysis without dominance.ConclusionWe have extended FIA to include QTL dominance effects. The power of FIA was superior, or similar, to standard regression methods for QTL effects with dominance. The difference in power for FIA with or without dominance is expected to be small as long as the QTL effects are not overdominant. We suggest that FIA with only additive effects should be the standard model to be used, especially since it is more computationally efficient.
| S-EPMC2716366 | biostudies-literature
Project description:Natural olfactory stimuli are volatile-chemical mixtures in which relative perceptual saliencies determine which odor-components are identified. Odor identification also depends on rapid selective adaptation, as shown for 4 odor stimuli in an earlier experimental simulation of natural conditions. Adapt-test pairs of mixtures of water-soluble, distinct odor stimuli with chemical features in common were studied. Identification decreased for adapted components but increased for unadapted mixture-suppressed components, showing compound identities were retained, not degraded to individual molecular features. Four additional odor stimuli, 1 with 2 perceptible odor notes, and an added "water-adapted" control tested whether this finding would generalize to other 4-compound sets. Selective adaptation of mixtures of the compounds (odors): 3 mM benzaldehyde (cherry), 5 mM maltol (caramel), 1 mM guaiacol (smoke), and 4 mM methyl anthranilate (grape-smoke) again reciprocally unmasked odors of mixture-suppressed components in 2-, 3-, and 4-component mixtures with 2 exceptions. The cherry note of "benzaldehyde" (itself) and the shared note of "methyl anthranilate and guaiacol" (together) were more readily identified. The pervasive mixture-component dominance and dynamic perceptual salience may be mediated through peripheral adaptation and central mutual inhibition of neural responses. Originating in individual olfactory receptor variants, it limits odor identification and provides analytic properties for momentary recognition of a few remaining mixture-components.
| S-EPMC5863551 | biostudies-literature
Project description:In the field of asset allocation, how to balance the returns of an investment portfolio and its fluctuations is the core issue. Capital asset pricing model, arbitrage pricing theory, and Fama-French three-factor model were used to quantify the price of individual stocks and portfolios. Based on the second-order stochastic dominance rule, the higher moments of return series, the Shannon entropy, and some other actual investment constraints, we construct a multiconstraint portfolio optimization model, aiming at comprehensively weighting the returns and risk of portfolios rather than blindly maximizing its returns. Furthermore, the whale optimization algorithm based on FTSE100 index data is used to optimize the above multiconstraint portfolio optimization model, which significantly improves the rate of return of the simple diversified buy-and-hold strategy or the FTSE100 index. Furthermore, extensive experiments validate the superiority of the whale optimization algorithm over the other four swarm intelligence optimization algorithms (gray wolf optimizer, fruit fly optimization algorithm, particle swarm optimization, and firefly algorithm) through various indicators of the results, especially under harsh constraints.
| S-EPMC7474746 | biostudies-literature
Project description:In this work, we develop a data-driven modelling framework to reproduce the locomotion of fish in a confined environment. Specifically, we highlight the primary characteristics of the motion of individual zebrafish (Danio rerio), and study how these can be suitably encapsulated within a mathematical framework utilising a limited number of calibrated model parameters. Using data captured from individual zebrafish via automated visual tracking, we develop a model using stochastic differential equations and describe fish as a self propelled particle moving in a plane. Based on recent experimental evidence of the importance of speed regulation in social behaviour, we extend stochastic models of fish locomotion by introducing experimentally-derived processes describing dynamic speed regulation. Salient metrics are defined which are then used to calibrate key parameters of coupled stochastic differential equations, describing both speed and angular speed of swimming fish. The effects of external constraints are also included, based on experimentally observed responses. Understanding the spontaneous dynamics of zebrafish using a bottom-up, purely data-driven approach is expected to yield a modelling framework for quantitative investigation of individual behaviour in the presence of various external constraints or biological assays.
| S-EPMC4598355 | biostudies-literature
Project description:Phenotypic variation is a hallmark of cellular physiology. Metabolic heterogeneity, in particular, underpins single-cell phenomena such as microbial drug tolerance and growth variability. Much research has focussed on transcriptomic and proteomic heterogeneity, yet it remains unclear if such variation permeates to the metabolic state of a cell. Here we propose a stochastic model to show that complex forms of metabolic heterogeneity emerge from fluctuations in enzyme expression and catalysis. The analysis predicts clonal populations to split into two or more metabolically distinct subpopulations. We reveal mechanisms not seen in deterministic models, in which enzymes with unimodal expression distributions lead to metabolites with a bimodal or multimodal distribution across the population. Based on published data, the results suggest that metabolite heterogeneity may be more pervasive than previously thought. Our work casts light on links between gene expression and metabolism, and provides a theory to probe the sources of metabolite heterogeneity.
| S-EPMC6428880 | biostudies-literature
Project description:MALDI fingerprinting was first described two decades ago as a technique to identify microbial cell lines. Microbial fingerprinting has since evolved into an automated platform for microorganism identification and classification, which is now routinely used in clinical and environmental sectors. The extension of fingerprinting to mammalian cells has yet to progress partly due to compartmentalization of eukaryotic cells and overall higher cellular complexity. A number of publications on mammalian whole cell fingerprinting suggest that the method could be useful for classification of different cell types, cell states, and monitoring cell differentiation. We report the optimization of MALDI fingerprinting workflow parameters for mammalian cells and its application for differential profiling of mammalian cell lines and two-component cell line mixtures. Murine fallopian tube cells and high-grade ovarian carcinoma cell lines and their mixtures are used as model mammalian cell lines. Two-component cell mixtures serve to determine the method's feasibility for complex biological samples as the ability to detect cancer cells in a mixed cell population. The level of detection of cancer cells in the two-component mixture by principle component analysis (PCA) starts to deteriorate at 5% but with application of a different statistical approach, Wilcoxon rank sum test, the level of detection was determined to be 1%. The ability to differentiate heterogeneous cell mixtures will help further extend whole cell MALDI fingerprinting to complex biological systems. Graphical Abstract.
| S-EPMC6347503 | biostudies-literature