Bayesian structural equation modeling in multiple omics data with application to circadian genes.
Ontology highlight
ABSTRACT: Motivation
It is well known that the integration among different data-sources is reliable because of its potential of unveiling new functionalities of the genomic expressions, which might be dormant in a single-source analysis. Moreover, different studies have justified the more powerful analyses of multi-platform data. Toward this, in this study, we consider the circadian genes' omics profile, such as copy number changes and RNA-sequence data along with their survival response. We develop a Bayesian structural equation modeling coupled with linear regressions and log normal accelerated failure-time regression to integrate the information between these two platforms to predict the survival of the subjects. We place conjugate priors on the regression parameters and derive the Gibbs sampler using the conditional distributions of them.Results
Our extensive simulation study shows that the integrative model provides a better fit to the data than its closest competitor. The analyses of glioblastoma cancer data and the breast cancer data from TCGA, the largest genomics and transcriptomics database, support our findings.Availability and implementation
The developed method is wrapped in R package available at https://github.com/MAITYA02/semmcmc.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Maity AK
PROVIDER: S-EPMC7332567 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:ObjectivesWhen examining the relationship between smoking intensity and toxicant exposure biomarkers in an effort to understand the potential risk for smoking-related disease, individual biomarkers may not be strongly associated with smoking intensity because of the inherent variability in biomarkers. Structural equation modeling (SEM) offers a powerful solution by modeling the relationship between smoking intensity and multiple biomarkers through a latent variable.MethodsBaseline data from a randomized trial (N = 1250) were used to estimate the relationship between smoking intensity and a latent toxicant exposure variable summarizing five volatile organic compound biomarkers. Two variables of smoking intensity were analyzed: the self-report cigarettes smoked per day and total nicotine equivalents in urine. SEM was compared with linear regression with each biomarker analyzed individually or with the sum score of the five biomarkers.ResultsSEM models showed strong relationships between smoking intensity and the latent toxicant exposure variable, and the relationship was stronger than its counterparts in linear regression with each biomarker analyzed separately or with the sum score.ConclusionsSEM is a powerful multivariate statistical method for studying multiple biomarkers assessing the same class of harmful constituents. This method could be used to evaluate exposure from different combusted tobacco products.
| S-EPMC9075702 | biostudies-literature
Project description:Learning is an intentional activity, with several factors affecting students' intention to use new learning technology. Researchers have investigated technology acceptance in different contexts by developing various theories/models and testing them by a number of means. Although most theories/models developed have been examined through regression or structural equation modeling, Bayesian analysis offers more accurate data analysis results. To address this gap, the unified theory of acceptance and technology use in the context of e-learning via Facebook are re-examined in this study using Bayesian analysis. The data (S1 Data) were collected from 170 students enrolled in a business statistics course at University of Malaya, Malaysia, and tested with the maximum likelihood and Bayesian approaches. The difference between the two methods' results indicates that performance expectancy and hedonic motivation are the strongest factors influencing the intention to use e-learning via Facebook. The Bayesian estimation model exhibited better data fit than the maximum likelihood estimator model. The results of the Bayesian and maximum likelihood estimator approaches are compared and the reasons for the result discrepancy are deliberated.
| S-EPMC5590745 | biostudies-literature
Project description:Bayesian structural equation modeling (BSEM) is a flexible tool for the exploration and estimation of sparse factor loading structures; that is, most cross-loading entries are zero and only a few important cross-loadings are nonzero. The current investigation was focused on the BSEM with small-variance normal distribution priors (BSEM-N) for both variable selection and model estimation. The prior sensitivity in BSEM-N was explored in factor analysis models with sparse loading structures through a simulation study (Study 1) and an empirical example (Study 2). Study 1 examined the prior sensitivity in BSEM-N based on the model fit, population model recovery, true and false positive rates, and parameter estimation. Seven shrinkage priors on cross-loadings and five noninformative/vague priors on other model parameters were examined. Study 2 provided a real data example to illustrate the impact of various priors on model fit and parameter selection and estimation. Results indicated that when the 95% credible intervals of shrinkage priors barely covered the population cross-loading values, it resulted in the best balance between true and false positives. If the goal is to perform variable selection, a sparse cross-loading structure is required, preferably with a minimal number of nontrivial cross-loadings and relatively high primary loading values. To improve parameter estimates, a relatively large prior variance is preferred. When cross-loadings are relatively large, BSEM-N with zero-mean priors is not recommended for the estimation of cross-loadings and factor correlations.
| S-EPMC7565120 | biostudies-literature
Project description:Optimized breeding goals and management practices for the improvement of target traits requires knowledge regarding any potential functional relationships between them. Fitting a structural equation model (SEM) allows for inferences about the magnitude of causal effects between traits to be made. In recent years, an adaptation of SEM was proposed in the context of quantitative genetics and mixed models. Several studies have since applied the SEM in the context of animal breeding. However, fitting the SEM requires choosing a causal structure with prior biological or temporal knowledge. The inductive causation (IC) algorithm can be used to recover an underlying causal structure from observed associations between traits. The results of the papers, which are introduced in this review, showed that using the IC algorithm to infer a causal structure is a helpful tool for detecting a causal structure without proper prior knowledge or with uncertain relationships between traits. The reports also presented that fitting the SEM could infer the effects of interventions, which are not given by correlations. Hence, information from the SEM provides more insights into and suggestions on breeding strategy than that from a multiple-trait model, which is the conventional model used for multitrait analysis.
| S-EPMC7187322 | biostudies-literature
Project description:In classifications of vocabulary knowledge, vocabulary size and depth have often been separately conceptualized (Schmitt, 2014). Although size and depth are known to be substantially correlated, it is not clear whether they are a single construct or two separate components of vocabulary knowledge (Yanagisawa and Webb, 2020). This issue has not been addressed extensively in the literature and can be better examined using structural equation modeling (SEM), with measurement error modeled separately from the construct of interest. The current study reports on conventional and Bayesian SEM approaches (e.g., Muthén and Asparouhov, 2012) to examine the factor structure of the size and depth of second language vocabulary knowledge of Japanese adult learners of English. A total of 255 participants took five vocabulary tests. One test was designed to measure vocabulary size in terms of the number of words known, while the remaining four were designed to measure vocabulary depth in terms of word association, polysemy, and collocation. All tests used a multiple-choice format. The size test was divided into three subtests according to word frequency. Results from conventional and Bayesian SEM show that a correlated two-factor model of size and depth with three and four indicators, respectively, fit better than a single-factor model of size and depth. In the two-factor model, vocabulary size and depth were strongly correlated (r = 0.945 for conventional SEM and 0.943 for Bayesian SEM with cross-loadings), but they were distinct. The implications of these findings are discussed.
| S-EPMC7187790 | biostudies-literature
Project description:We meta-analyzed the relationship between team planning and performance moderated by task, team, context, and methodological factors. For testing our hypothesized model, we used a meta-analytic structural equation modeling approach. Based on K = 33 independent samples (N = 1,885 teams), a mixed-effects model indicated a non-zero moderate positive effect size (ρ = .31, 95% CI [.20, .42]). Methodological quality, generally rated as adequate, was unrelated to effect size. Sensitivity analyses suggest that effect sizes were robust to exclusion of any individual study and publication bias. The statistical power of the studies was generally low and significantly moderated the relationship, with a large positive relationship for studies with high-powered (k = 42, ρ = .40, 95% CI [.27, .54]) and a smaller, significant relationship for low-powered studies (k = 54, ρ = .16, 95% CI [.01, .30]). The effect size was robust and generally not qualified by a large number of moderators, but was more pronounced for less interdependent tasks, less specialized team members, and assessment of quality rather than quantity of planning. Latent class analysis revealed no qualitatively different subgroups within populations. We recommend large-scale collaboration to overcome several methodological weaknesses of the current literature, which is severely underpowered, potentially biased by self-reporting data, and lacks long-term follow-ups.
| S-EPMC9838875 | biostudies-literature
Project description:MotivationCircadian rhythms are approximately 24-h endogenous cycles that control many biological functions. To identify these rhythms, biological samples are taken over circadian time and analyzed using a single omics type, such as transcriptomics or proteomics. By comparing data from these single omics approaches, it has been shown that transcriptional rhythms are not necessarily conserved at the protein level, implying extensive circadian post-transcriptional regulation. However, as proteomics methods are known to be noisier than transcriptomic methods, this suggests that previously identified arrhythmic proteins with rhythmic transcripts could have been missed due to noise and may not be due to post-transcriptional regulation.ResultsTo determine if one can use information from less-noisy transcriptomic data to inform rhythms in more-noisy proteomic data, and thus more accurately identify rhythms in the proteome, we have created the Multi-Omics Selection with Amplitude Independent Criteria (MOSAIC) application. MOSAIC combines model selection and joint modeling of multiple omics types to recover significant circadian and non-circadian trends. Using both synthetic data and proteomic data from Neurospora crassa, we showed that MOSAIC accurately recovers circadian rhythms at higher rates in not only the proteome but the transcriptome as well, outperforming existing methods for rhythm identification. In addition, by quantifying non-circadian trends in addition to circadian trends in data, our methodology allowed for the recognition of the diversity of circadian regulation as compared to non-circadian regulation.Availability and implementationMOSAIC's full interface is available at https://github.com/delosh653/MOSAIC. An R package for this functionality, mosaic.find, can be downloaded at https://CRAN.R-project.org/package=mosaic.find.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC8098022 | biostudies-literature
Project description:BackgroundThe use of Bayesian Structural Equation Model (BSEM) to evaluate the impact of TB on self-reported health related quality of life (HRQoL) of TB patients has been not studied.ObjectiveTo identify the factors that contribute to the HRQoL of TB patients using BSEM.MethodsThis is a latent variable modeling with Bayesian approach using secondary data. HRQoL data collected after one year from newly diagnosed 436 TB patients who were registered and successfully completed treatment at Government health facilities in Tiruvallur district, south India under the National TB Elimination Programme (NTEP) were used for this analysis. In this study, the four independent latent variables such as physical well-being (PW = PW1-7), mental well-being (MW = MW1-7), social well-being (SW = SW1-4) and habits were considered. The BSEM was constructed using Markov Chain Monte Carlo algorithm for identifying the factors that contribute to the HRQoL of TB patients who completed treatment.ResultsBayesian estimates were obtained using 46,300 observations after convergence and the standardized structural regression estimate of PW, MW, SW on HRQoL were 0.377 (p<0.001), 0.543 (p<0.001) and 0.208 (p<0.001) respectively. The latent variables PW, MW and SW were significantly associated with HRQoL of TB patients. The age was found to be significantly negatively associated with HRQoL of TB patients.ConclusionsThe current study demonstrated the application of BSEM in evaluating HRQoL. This methodology may be used to study precise estimates of HRQoL of TB patients in different time points.
| S-EPMC8162623 | biostudies-literature
Project description:There has been great interest in developing nonlinear structural equation models and associated statistical inference procedures, including estimation and model selection methods. In this paper a general semiparametric structural equation model (SSEM) is developed in which the structural equation is composed of nonparametric functions of exogenous latent variables and fixed covariates on a set of latent endogenous variables. A basis representation is used to approximate these nonparametric functions in the structural equation and the Bayesian Lasso method coupled with a Markov Chain Monte Carlo (MCMC) algorithm is used for simultaneous estimation and model selection. The proposed method is illustrated using a simulation study and data from the Affective Dynamics and Individual Differences (ADID) study. Results demonstrate that our method can accurately estimate the unknown parameters and correctly identify the true underlying model.
| S-EPMC3371320 | biostudies-literature
Project description:IntroductionStructural equation modeling (SEM) is a multivariate data analytic technique used in many domains of addictive behaviors research. SEM results are usually summarized and communicated through statistical tables and path diagrams, which emphasize path coefficients and global fit without showing specific quantitative values of data points that underlie the model results. Data visualization methods are often absent in SEM research, which may limit the quality and impact of SEM research by reducing data transparency, obscuring unexpected data anomalies and unmodeled heterogeneity, and inhibiting the communication of SEM research findings to research stakeholders who do not have advanced statistical training in SEM.Methods and resultsIn this report, we show how data visualization methods can address these limitations and improve the quality of SEM-based addictive behaviors research. We first introduce SEM and data visualization methodologies and differentiate data visualizations from model visualizations that are commonly used in SEM, such as path diagrams. We then discuss ways researchers may utilize data visualization in SEM research, including by obtaining estimates of latent variables and by visualizing multivariate relations in two-dimensional figures. R syntax is provided to help others generate data visualizations for several types of effects commonly modeled in SEM, including correlation, regression, moderation, and simple mediation.DiscussionThe techniques outlined here may help spur the use of data visualization in SEM-based addictive behaviors research. Using data visualization in SEM may enhance methodological transparency and improve communication of research findings.
| S-EPMC6760248 | biostudies-literature