Advice on comparing two independent samples of circular data in biology.
Ontology highlight
ABSTRACT:
SUBMITTER: Landler L
PROVIDER: S-EPMC8514454 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Circular data are common in biological studies. The most fundamental question that can be asked of a sample of circular data is whether it suggests that the underlying population is uniformly distributed around the circle, or whether it is concentrated around at least one preferred direction (e.g. a migratory goal or activity phase). We compared the statistical power of five commonly used tests (the Rayleigh test, the V-test, Watson's test, Kuiper's test and Rao's spacing test) across a range of different unimodal scenarios. The V-test showed higher power for symmetrical distributions, Rao's spacing performed worst for all explored unimodal distributions tested and the remaining three tests showed very similar performance. However, the V-test only applies if the hypothesis is restricted to one (pre-specified) direction of interest. In all other unimodal cases, we recommend using the Rayleigh test. Much less explored is the multimodal case with data concentrated around several directions. We performed power simulations for a variety of multimodal situations, testing the performance of the widely used Rayleigh, Rao's, Watson, and Kuiper's tests as well as the more recent Bogdan and Hermans-Rasson tests. Our analyses of alternative statistical methods show that the commonly used tests lack statistical power in many of multimodal cases. Transformation of the raw data (e.g. doubling the angles) can overcome some of the issues, but only in the case of perfect f-fold symmetry. However, the Hermans-Rasson method, which is not yet implemented in any software package, outcompetes the alternative tests (often by substantial margins) in most of the multimodal situations explored. We recommend the wider uptake of the powerful but hitherto neglected Hermans-Rasson method. In summary, we provide guidance for biologists helping them to make decisions when testing circular data for single or multiple departures from uniformity.
| S-EPMC6060829 | biostudies-other
Project description:AbstractThe most common statistical procedure with a sample of circular data is to test the null hypothesis that points are spread uniformly around the circle without a preferred direction. An array of tests for this has been developed. However, these tests were designed for continuously distributed data, whereas often (e.g. due to limited precision of measurement techniques) collected data is aggregated into a set of discrete values (e.g. rounded to the nearest degree). This disparity can cause an uncontrolled increase in type I error rate, an effect that is particularly problematic for tests that are based on the distribution of arc lengths between adjacent points (such as the Rao spacing test). Here, we demonstrate that an easy-to-apply modification can correct this problem, and we recommend this modification when using any test, other than the Rayleigh test, of circular uniformity on aggregated data. We provide R functions for this modification for several commonly used tests. In addition, we tested the power of a recently proposed test, the Gini test. However, we concluded that it lacks sufficient increase in power to replace any of the tests already in common use. In conclusion, using any of the standard circular tests (except the Rayleigh test) without modifications on rounded/aggregated data, especially with larger sample sizes, will increase the proportion of false-positive results-but we demonstrate that a simple and general modification avoids this problem.Significance statementCircular data are widespread across biological disciplines, e.g. in orientation studies or circadian rhythms. Often these data are rounded to the nearest 1-10 degrees. We have shown previously that this leads to an inflation of false-positive results when testing whether the data is significantly different from a random distribution using the Rao test. Here we present a modification that avoids this increase in false-positives for rounded data while retaining statistical power for a variety of tests. In sum, we provide comprehensive guidance on how best to test for departure from uniformity in non-continuous data.
| S-EPMC7373216 | biostudies-literature
Project description:BackgroundStudent's two-sample t test is generally used for comparing the means of two independent samples, for example, two treatment arms. Under the null hypothesis, the t test assumes that the two samples arise from the same normally distributed population with unknown variance. Adequate control of the Type I error requires that the normality assumption holds, which is often examined by means of a preliminary Shapiro-Wilk test. The following two-stage procedure is widely accepted: If the preliminary test for normality is not significant, the t test is used; if the preliminary test rejects the null hypothesis of normality, a nonparametric test is applied in the main analysis.MethodsEqually sized samples were drawn from exponential, uniform, and normal distributions. The two-sample t test was conducted if either both samples (Strategy I) or the collapsed set of residuals from both samples (Strategy II) had passed the preliminary Shapiro-Wilk test for normality; otherwise, Mann-Whitney's U test was conducted. By simulation, we separately estimated the conditional Type I error probabilities for the parametric and nonparametric part of the two-stage procedure. Finally, we assessed the overall Type I error rate and the power of the two-stage procedure as a whole.ResultsPreliminary testing for normality seriously altered the conditional Type I error rates of the subsequent main analysis for both parametric and nonparametric tests. We discuss possible explanations for the observed results, the most important one being the selection mechanism due to the preliminary test. Interestingly, the overall Type I error rate and power of the entire two-stage procedure remained within acceptable limits.ConclusionThe two-stage procedure might be considered incorrect from a formal perspective; nevertheless, in the investigated examples, this procedure seemed to satisfactorily maintain the nominal significance level and had acceptable power properties.
| S-EPMC3444333 | biostudies-literature
Project description:College students' performance in introductory-level biology course work is an important predictor of ongoing persistence in the major. This study reports on a researcher-educator partnership that designed and compared two cocurricular workshops. Seventeen laboratory sections of an undergraduate biology course were randomly assigned to one of two educational interventions during the regularly scheduled lab class section after students had completed and received the results for the first exam. The baseline Metacognition intervention was an hourlong workshop focused on effective learning strategies and self-awareness in the learning process; the extended Metacognition plus Time Management (Metacognition+TM) intervention included the aforementioned workshop plus a second hourlong workshop on time management and procrastination. Based on three exams and self-report surveys administered before the intervention and at the end of the semester, students who participated in the Metacognition+TM intervention experienced greater increases in their exam scores and degree commitment than those in the baseline intervention. Additionally, group status moderated the effect of the intervention, as the Metacognition+TM intervention was especially effective in increasing use of time management tools by students from minoritized groups.
| S-EPMC8715789 | biostudies-literature
Project description:Epidemiological data indicate a potential relationship between milk and dairy product consumption and the incidence of breast cancer, as well as neurodegenerative diseases. We report the isolation of two novel circular DNA molecules isolated from commercially available milk.
| S-EPMC5408121 | biostudies-literature
Project description:BACKGROUND: A number of statistical models has been proposed for studying the association between gene expression and copy number data in integrated analysis. The next step is to compare association patterns between different groups of samples. RESULTS: We propose a method, named dSIM, to find differences in association between copy number and gene expression, when comparing two groups of samples. Firstly, we use ridge regression to correct for the baseline associations between copy number and gene expression. Secondly, the global test is applied to the corrected data in order to find differences in association patterns between two groups of samples. We show that dSIM detects differences even in small genomic regions in a simulation study. We also apply dSIM to two publicly available breast cancer datasets and identify chromosome arms where copy number led gene expression regulation differs between positive and negative estrogen receptor samples. In spite of differing genomic coverage, some selected arms are identified in both datasets. CONCLUSION: We developed a flexible and robust method for studying association differences between two groups of samples while integrating genomic data from different platforms. dSIM can be used with most types of microarray/sequencing data, including methylation and microRNA expression. The method is implemented in R and will be made part of the BioConductor package SIM.
| S-EPMC4227098 | biostudies-literature
Project description:Using standard systems biology methodologies a 14-compartment dynamic model was developed for the Corona virus epidemic. The model predicts that: (i) it will be impossible to limit lockdown intensity such that sufficient herd immunity develops for this epidemic to die down, (ii) the death toll from the SARS-CoV-2 virus decreases very strongly with increasing intensity of the lockdown, but (iii) the duration of the epidemic increases at first with that intensity and then decreases again, such that (iv) it may be best to begin with selecting a lockdown intensity beyond the intensity that leads to the maximum duration, (v) an intermittent lockdown strategy should also work and might be more acceptable socially and economically, (vi) an initially intensive but adaptive lockdown strategy should be most efficient, both in terms of its low number of casualties and shorter duration, (vii) such an adaptive lockdown strategy offers the advantage of being robust to unexpected imports of the virus, e.g. due to international travel, (viii) the eradication strategy may still be superior as it leads to even fewer deaths and a shorter period of economic downturn, but should have the adaptive strategy as backup in case of unexpected infection imports, (ix) earlier detection of infections is the most effective way in which the epidemic can be controlled, whilst waiting for vaccines.
| S-EPMC7293297 | biostudies-literature
Project description:Single-cell transcriptomics has revolutionized our understanding of basic biology and disease. Since transcript levels often do not correlate with protein expression, it is crucial to complement transcriptomics approaches with proteome analyses at single-cell resolution. Despite continuous technological improvements in sensitivity, mass-spectrometry-based single-cell proteomics ultimately faces the challenge of reproducibly comparing the protein expression profiles of thousands of individual cells. Here, we combine two hitherto opposing analytical strategies, DIA and Tandem-Mass-Tag (TMT)-multiplexing, to generate highly reproducible, quantitative proteome signatures from ultralow input samples. We developed a novel, identification-independent proteomics data-analysis pipeline that allows to quantitatively compare DIA-TMT proteome signatures across hundreds of samples independent of their biological origin to identify cell types and single protein knockouts. These proteome signatures overcome the need to impute quantitative data due to accumulating detrimental amounts of missing data in standard multibatch TMT experiments. We validate our approach using integrative data analysis of different human cell lines and standard database searches for knockouts of defined proteins. Our data establish a novel and reproducible approach to markedly expand the numbers of proteins one detects from ultralow input samples.
| S-EPMC8717550 | biostudies-literature
Project description:Mass spectrometry-based metaproteomics is a relatively new field of research that enables the characterization of the functionality of microbiota. Recently, we demonstrated the applicability of data-independent acquisition (DIA) mass spectrometry to the analysis of complex metaproteomic samples. This allowed us to circumvent many of the drawbacks of the previously used data-dependent acquisition (DDA) mass spectrometry, mainly the limited reproducibility when analyzing samples with complex microbial composition. However, the DDA-assisted DIA approach still required additional DDA data on the samples to assist the analysis. Here, we introduce, for the first time, an untargeted DIA metaproteomics tool that does not require any DDA data, but instead generates a pseudospectral library directly from the DIA data. This reduces the amount of required mass spectrometry data to a single DIA run per sample. The new DIA-only metaproteomics approach is implemented as a new open-source software package named glaDIAtor, including a modern web-based graphical user interface to facilitate wide use of the tool by the community.
| S-EPMC9723653 | biostudies-literature
Project description:Protein ubiquitination is involved in virtually all cellular processes. Enrichment strategies employing antibodies targeting ubiquitin-derived diGly remnants combined with mass spectrometry (MS) have enabled investigations of ubiquitin signaling at a large scale. However, so far the power of data independent acquisition (DIA) with regards to sensitivity in single run analysis and data completeness have not yet been explored. Here, we develop a sensitive workflow combining diGly antibody-based enrichment and optimized Orbitrap-based DIA with comprehensive spectral libraries together containing more than 90,000 diGly peptides. This approach identifies 35,000 diGly peptides in single measurements of proteasome inhibitor-treated cells - double the number and quantitative accuracy of data dependent acquisition. Applied to TNF signaling, the workflow comprehensively captures known sites while adding many novel ones. An in-depth, systems-wide investigation of ubiquitination across the circadian cycle uncovers hundreds of cycling ubiquitination sites and dozens of cycling ubiquitin clusters within individual membrane protein receptors and transporters, highlighting new connections between metabolism and circadian regulation.
| S-EPMC7801436 | biostudies-literature