A comparison of gene region simulation methods.
Ontology highlight
ABSTRACT: Background
Accurately modeling LD in simulations is essential to correctly evaluate new and existing association methods. At present, there has been minimal research comparing the quality of existing gene region simulation methods to produce LD structures similar to an existing gene region. Here we compare the ability of three approaches to accurately simulate the LD within a gene region: HapSim (2005), Hapgen (2009), and a minor extension to simple haplotype resampling.Methodology/principal findings
In order to observe the variation and bias for each method, we compare the simulated pairwise LD measures and minor allele frequencies to the original HapMap data in an extensive simulation study. When possible, we also evaluate the effects of changing parameters. HapSim produces samples of haplotypes with lower LD, on average, compared to the original haplotype set while both our resampling method and Hapgen do not introduce this bias. The variation introduced across the replicates by our resampling method is quite small and may not provide enough sampling variability to make a generalizable simulation study.Conclusion
We recommend using Hapgen to simulate replicate haplotypes from a gene region. Hapgen produces moderate sampling variation between the replicates while retaining the overall unique LD structure of the gene region.
SUBMITTER: Hendricks AE
PROVIDER: S-EPMC3399793 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundWe present a series of simulation studies that explore the relative performance of several phylogenetic network approaches (statistical parsimony, split decomposition, union of maximum parsimony trees, neighbor-net, simulated history recombination upper bound, median-joining, reduced median joining and minimum spanning network) compared to standard tree approaches, (neighbor-joining and maximum parsimony) in the presence and absence of recombination.Principal findingsIn the absence of recombination, all methods recovered the correct topology and branch lengths nearly all of the time when the substitution rate was low, except for minimum spanning networks, which did considerably worse. At a higher substitution rate, maximum parsimony and union of maximum parsimony trees were the most accurate. With recombination, the ability to infer the correct topology was halved for all methods and no method could accurately estimate branch lengths.ConclusionsOur results highlight the need for more accurate phylogenetic network methods and the importance of detecting and accounting for recombination in phylogenetic studies. Furthermore, we provide useful information for choosing a network algorithm and a framework in which to evaluate improvements to existing methods and novel algorithms developed in the future.
| S-EPMC2275308 | biostudies-literature
Project description:Six empirical force fields were tested for applicability to calculations for automated carbohydrate database filling. They were probed on eleven disaccharide molecules containing representative structural features from widespread classes of carbohydrates. The accuracy of each method was queried by predictions of nuclear Overhauser effects (NOEs) from conformational ensembles obtained from 50 to 100 ns molecular dynamics (MD) trajectories and their comparison to the published experimental data. Using various ranking schemes, it was concluded that explicit solvent MM3 MD yielded non-inferior NOE accuracy with newer GLYCAM-06, and ultimately PBE0-D3/def2-TZVP (Triple-Zeta Valence Polarized) Density Functional Theory (DFT) simulations. For seven of eleven molecules, at least one empirical force field with explicit solvent outperformed DFT in NOE prediction. The aggregate of characteristics (accuracy, speed, and compatibility) made MM3 dynamics with explicit solvent at 300 K the most favorable method for bulk generation of disaccharide conformation maps for massive database filling.
| S-EPMC7589101 | biostudies-literature
Project description:BackgroundMultiple imputation (MI) was developed as a method to enable valid inferences to be obtained in the presence of missing data rather than to re-create the missing values. Within the applied setting, it remains unclear how important it is that imputed values should be plausible for individual observations. One variable type for which MI may lead to implausible values is a limited-range variable, where imputed values may fall outside the observable range. The aim of this work was to compare methods for imputing limited-range variables, with a focus on those that restrict the range of the imputed values.MethodsUsing data from a study of adolescent health, we consider three variables based on responses to the General Health Questionnaire (GHQ), a tool for detecting minor psychiatric illness. These variables, based on different scoring methods for the GHQ, resulted in three continuous distributions with mild, moderate and severe positive skewness. In an otherwise complete dataset, we set 33% of the GHQ observations to missing completely at random or missing at random; repeating this process to create 1000 datasets with incomplete data for each scenario.For each dataset, we imputed values on the raw scale and following a zero-skewness log transformation using: univariate regression with no rounding; post-imputation rounding; truncated normal regression; and predictive mean matching. We estimated the marginal mean of the GHQ and the association between the GHQ and a fully observed binary outcome, comparing the results with complete data statistics.ResultsImputation with no rounding performed well when applied to data on the raw scale. Post-imputation rounding and imputation using truncated normal regression produced higher marginal means than the complete data estimate when data had a moderate or severe skew, and this was associated with under-coverage of the complete data estimate. Predictive mean matching also produced under-coverage of the complete data estimate. For the estimate of association, all methods produced similar estimates to the complete data.ConclusionsFor data with a limited range, multiple imputation using techniques that restrict the range of imputed values can result in biased estimates for the marginal mean when data are highly skewed.
| S-EPMC4021274 | biostudies-literature
Project description:In response time (RT) research, RT outliers are typically excluded from statistical analysis to improve the signal-to-noise ratio. Nevertheless, there exist several methods for outlier exclusion. This poses the question, how these methods differ with respect to recovering the uncontaminated RT distribution. In the present simulation study, two RT distributions with a given population difference were simulated in each iteration. RTs were replaced by outliers following two different approaches. The first approach generated outliers at the tails of the distribution, the second one inserted outliers overlapping with the genuine RT distribution. We applied ten different outlier exclusion methods and tested, how many pairs of distributions significantly differed. Outlier exclusion methods were compared in terms of bias. Bias was defined as the deviation of the proportion of significant differences after outlier exclusion from the proportion of significant differences in the uncontaminated samples (before introducing outliers). Our results showed large differences in bias between the exclusion methods. Some methods showed a high rate of Type-I errors and should therefore clearly not be used. Overall, our results showed that applying an exclusion method based on z-scores / standard deviations introduced only small biases, while the absence of outlier exclusion showed the largest absolute bias.
| S-EPMC8238084 | biostudies-literature
Project description:We have evaluated various density functional theory (DFT) methods to simulate geometric, energetic, electronic, and hydrogen adsorption properties of metal-nanoparticles supported on metal surfaces. We used Pt and Pd nanoislands on Au(111) as model systems. The evaluated DFT methods include GGA (PW91, PBE, RPBE, revPBE, and PBESol), GGA with van der Waals (vdW) corrected (PBE-D3), GGA with optimized vdW functionals (revPBE-vdW), meta-GGA (SCAN and MS2), and the machine learning-based method BEEF-vdW. The results show that the various DFT methods yield similar geometric and electronic properties for Pt (or Pd) nanoislands on Au(111). The DFT methods also produce similar relative energetics for small Pt (or Pd) clusters with different conformations on Au(111). The results show that a triatomic cluster of Pt on Au(111) is more stable with a linear conformation. In contrast, a triatomic cluster of Pd is more stable with a triangular conformation. For clusters with four or more atoms, Pt and Pd clusters on Au(111) prefer non-linear conformation. We found that the various DFT methods yield different results only for the adsorption energy of hydrogen.
| S-EPMC8224827 | biostudies-literature
Project description:PurposePrivacy-protecting analytic and data-sharing methods that minimize the disclosure risk of sensitive information are increasingly important due to the growing interest in utilizing data across multiple sources. We conducted a simulation study to examine how avoiding sharing individual-level data in a distributed data network can affect analytic results.MethodsThe base scenario had four sites of varying sizes with 5% outcome incidence, 50% treatment prevalence, and seven confounders. We varied treatment prevalence, outcome incidence, treatment effect, site size, number of sites, and covariate distribution. Confounding adjustment was conducted using propensity score or disease risk score. We compared analyses of three types of aggregate-level data requested from sites: risk-set, summary-table, or effect-estimate data (meta-analysis) with benchmark results of analysis of pooled individual-level data. We assessed bias and precision of hazard ratio estimates as well as the accuracy of standard error estimates.ResultsAll the aggregate-level data-sharing approaches, regardless of confounding adjustment methods, successfully approximated pooled individual-level data analysis in most simulation scenarios. Meta-analysis showed minor bias when using inverse probability of treatment weights (IPTW) in infrequent exposure (5%), rare outcome (0.01%), and small site (5,000 patients) settings. SE estimates became less accurate for IPTW risk-set approach with less frequent exposure and for propensity score-matching meta-analysis approach with rare outcomes.ConclusionsOverall, we found that we can avoid sharing individual-level data and obtain valid results in many settings, although care must be taken with meta-analysis approach in infrequent exposure and rare outcome scenarios, particularly when confounding adjustment is performed with IPTW.
| S-EPMC6135666 | biostudies-literature
Project description:BackgroundInterprofessional collaborative practice is essential for meeting patients' needs and improving their health outcomes; thus, the effectiveness of interprofessional education (IPE) should be clearly identified. There is insufficient evidence in the literature to determine the outcomes of IPE compared to traditional single-profession education (SPE). This study aimed to compare the outcomes of IPE and SPE during a simulation training course.MethodsThe study design was a mixed-methods, incorporated cross-over design and a qualitative survey. A total of 54 students including 18 medical students and 36 nursing students were recruited from March to April 2019. The 4-week simulation course was designed based on Kolb's experimental learning theory and Bandura's social learning theory. Participants were evenly divided into group 1 (received IPE-learning followed by SPE-learning), and group 2 (received SPE-learning followed by IPE-learning). Students' medical task performance, team behavior performance, teamwork attitude, and patient safety attitude were collected at pretest, mid-test, and posttest. Descriptive statistics and repeated measures analysis of variance were used. End-of-study qualitative feedback was collected, and content analysis was performed.ResultsBoth groups demonstrated moderate-to-large within-group improvements for multiple learning outcomes at mid-test. Group 1 students' medical task performance (F = 97.25; P < 0.001) and team behavior performance (F = 31.17; P < 0.001) improved significantly. Group 2 students' medical task performance (F = 77.77; P < 0.001), team behavior performance (F = 40.14; P < 0.001), and patient safety attitude (F = 6.82; P < 0.01) improved significantly. Outcome differences between groups were nonsignificant. Qualitative themes identified included: personal factor, professional factor, interprofessional relationship, and learning. The IPE program provided students with exposure to other professions and revealed differences in expertise and responsibilities.ConclusionIPE-simulation and SPE-simulation were effective interventions that enabled medical and nursing students to develop critical medical management and team behavior performance. IPE-simulation provided more opportunities for improving competencies in interprofessional collaborative practice. In circumstances with limited teaching resources, SPE-simulation can be an acceptable alternative to IPE-simulation.
| S-EPMC9429663 | biostudies-literature
Project description:Comparison of 16S rRNA regions for ecological studies
| PRJEB7317 | ENA
Project description:Generalized structured component analysis (GSCA) is a theoretically well-founded approach to component-based structural equation modeling (SEM). This approach utilizes the bootstrap method to estimate the confidence intervals of its parameter estimates without recourse to distributional assumptions, such as multivariate normality. It currently provides the bootstrap percentile confidence intervals only. Recently, the potential usefulness of the bias-corrected and accelerated bootstrap (BCa) confidence intervals (CIs) over the percentile method has attracted attention for another component-based SEM approach-partial least squares path modeling. Thus, in this study, we implemented the BCa CI method into GSCA and conducted a rigorous simulation to evaluate the performance of three bootstrap CI methods, including percentile, BCa, and Student's t methods, in terms of coverage and balance. We found that the percentile method produced CIs closer to the desired level of coverage than the other methods, while the BCa method was less prone to imbalance than the other two methods. Study findings and implications are discussed, as well as limitations and directions for future research.
| S-EPMC6797821 | biostudies-literature
Project description:Background:Traditional simulation-based education prioritizes participation in simulated scenarios. The educational impact of observation in simulation-based education compared with participation remains uncertain. Our objective was to compare the performances of observers and participants in a standardized simulation scenario. Methods:We assessed learning differences between simulation-based scenario participation and observation using a convergent, parallel, quasi-experimental, mixed-methods study of 15 participants and 15 observers (N = 30). Fifteen first-year residents from six medical specialties were evaluated during a simulated scenario (cardiac arrest due to critical hyperkalemia). Evaluation included predefined critical actions and performance assessments. In the first exposure to the simulation scenario, participants and observers underwent a shared postevent debriefing with predetermined learning objectives. Three months later, a follow-up assessment using the same case scenario evaluated all 30 learners as participants. Wilcoxon signed rank and Wilcoxon rank sum tests were used to compare participants and observers at 3-month follow-up. In addition, we used case study methodology to explore the nature of learning for participants and observers. Data were triangulated using direct observations, reflective field notes, and a focus group. Results:Quantitative data analysis comparing the learners' first and second exposure to the investigation scenario demonstrated participants' time to calcium administration as the only statistically significant difference between participant and observer roles (316 seconds vs. 200 seconds, p = 0.0004). Qualitative analysis revealed that both participation and observation improved learning, debriefing was an important component to learning, and debriefing closed the learning gap between observers and participants. Conclusions:Participants and observers had similar performances in simulation-based learning in an isolated scenario of cardiac arrest due to hyperkalemia. Findings support current limited literature that observation should not be underestimated as an important opportunity to enhance simulation-based education. When paired with postevent debriefing, scenario observers and participants may reap similar educational benefits.
| S-EPMC6339532 | biostudies-literature