The CUSUM statistics of change-point models based on dependent sequences.
Ontology highlight
ABSTRACT:
SUBMITTER: Ding S
PROVIDER: S-EPMC9225615 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:In a natural environment, sensory systems are faced with ever-changing stimuli that can occur, disappear or change their properties at any time. For the animal to react adequately the sensory systems must be able to detect changes in external stimuli based on its neuronal responses. Since the nervous system has no prior knowledge of the stimulus timing, changes in stimulus need to be inferred from the changes in neuronal activity, in particular increase or decrease of the spike rate, its variability, and shifted response latencies. From a mathematical point of view, this problem can be rephrased as detecting changes of statistical properties in a time series. In neuroscience, the CUSUM (cumulative sum) method has been applied to recorded neuronal responses for detecting a single stimulus change. Here, we investigate the applicability of the CUSUM approach for detecting single as well as multiple stimulus changes that induce increases or decreases in neuronal activity. Like the nervous system, our algorithm relies exclusively on previous neuronal population activities, without using knowledge about the timing or number of external stimulus changes. We apply our change point detection methods to experimental data obtained by multi-electrode recordings from turtle retinal ganglion cells, which react to changes in light stimulation with a range of typical neuronal activity patterns. We systematically examine how variations of mathematical assumptions (Poisson, Gaussian, and Gamma distributions) used for the algorithms may affect the detection of an unknown number of stimulus changes in our data and compare these CUSUM methods with the standard Rate Change method. Our results suggest which versions of the CUSUM algorithm could be useful for different types of specific data sets.
| S-EPMC4916211 | biostudies-literature
Project description:This paper investigates a change-point estimation problem in the context of high-dimensional Markov random field models. Change-points represent a key feature in many dynamically evolving network structures. The change-point estimate is obtained by maximizing a profile penalized pseudo-likelihood function under a sparsity assumption. We also derive a tight bound for the estimate, up to a logarithmic factor, even in settings where the number of possible edges in the network far exceeds the sample size. The performance of the proposed estimator is evaluated on synthetic data sets and is also used to explore voting patterns in the US Senate in the 1979-2012 period.
| S-EPMC5571889 | biostudies-literature
Project description:A threshold effect takes place in situations where the relationship between an outcome variable and a predictor variable changes as the predictor value crosses a certain threshold/change point. Threshold effects are often plausible in a complex biological system, especially in defining immune responses that are protective against infections such as HIV-1, which motivates the current work. We study two hypothesis testing problems in change point models. We first compare three different approaches to obtaining a p-value for the maximum of scores test in a logistic regression model with change point variable as a main effect. Next, we study the testing problem in a logistic regression model with the change point variable both as a main effect and as part of an interaction term. We propose a test based on the maximum of likelihood ratios test statistic and obtain its reference distribution through a Monte Carlo method. We also propose a maximum of weighted scores test that can be more powerful than the maximum of likelihood ratios test when we know the direction of the interaction effect. In simulation studies, we show that the proposed tests have a correct type I error and higher power than several existing methods. We illustrate the application of change point model-based testing methods in a recent study of immune responses that are associated with the risk of mother to child transmission of HIV-1.
| S-EPMC4390452 | biostudies-literature
Project description:Plant development studies often generate data in the form of multivariate time series, each variable corresponding to a count of newly emerged organs for a given development process. These phenological data often exhibit highly structured patterns, and the aim of this study was to identify such patterns in cultivated strawberry. Six strawberry genotypes were observed weekly for their course of emergence of flowers, leaves, and stolons during 7 months. We assumed that these phenological series take the form of successive phases, synchronous between individuals. We applied univariate multiple change-point models for the identification of flowering, vegetative development, and runnering phases, and multivariate multiple change-point models for the identification of consensus phases for these three development processes. We showed that the flowering and the runnering processes are the main determinants of the phenological pattern. On this basis, we propose a typology of the six genotypes in the form of a hierarchical classification. This study introduces a new longitudinal data modeling approach for the identification of phenological phases in plant development. The focus was on development variables but the approach can be directly extended to growth variables and to multivariate series combining growth and development variables.
| S-EPMC6812722 | biostudies-literature
Project description:Change-point models are widely used by statisticians to model drastic changes in the pattern of observed data. Least squares/maximum likelihood based estimation of change-points leads to curious asymptotic phenomena. When the change-point model is correctly specified, such estimates generally converge at a fast rate (n) and are asymptotically described by minimizers of a jump process. Under complete mis-specification by a smooth curve, i.e. when a change-point model is fitted to data described by a smooth curve, the rate of convergence slows down to n1/3 and the limit distribution changes to that of the minimizer of a continuous Gaussian process. In this paper we provide a bridge between these two extreme scenarios by studying the limit behavior of change-point estimates under varying degrees of model mis-specification by smooth curves, which can be viewed as local alternatives. We find that the limiting regime depends on how quickly the alternatives approach a change-point model. We unravel a family of 'intermediate' limits that can transition, at least qualitatively, to the limits in the two extreme scenarios. The theoretical results are illustrated via a set of carefully designed simulations. We also demonstrate how inference for the change-point parameter can be performed in absence of knowledge of the underlying scenario by resorting to subsampling techniques that involve estimation of the convergence rate.
| S-EPMC4678008 | biostudies-literature
Project description:There is a growing research interest in the detection of changes in hydrologic and climatic time series. Stationarity can be assessed using the autocorrelation function, but this is not yet common practice in hydrology and climate. Here, we use a global land-based gridded annual precipitation (hereafter P) database (1940-2009) and find that the lag 1 autocorrelation coefficient is statistically significant at around 14% of the global land surface, implying nonstationary behavior (90% confidence). In contrast, around 76% of the global land surface shows little or no change, implying stationary behavior. We use these results to assess change in the observed P over the most recent decade of the database. We find that the changes for most (84%) grid boxes are within the plausible bounds of no significant change at the 90% CI. The results emphasize the importance of adequately accounting for natural variability when assessing change.
| S-EPMC5878000 | biostudies-literature
Project description:In standard treatments of atomic multipole models, interaction energies, total molecular forces, and total molecular torques are given for multipolar interactions between rigid molecules. However, if the molecules are assumed to be flexible, two additional multipolar atomic forces arise because of (1) the transfer of torque between neighboring atoms and (2) the dependence of multipole moment on internal geometry (bond lengths, bond angles, etc.) for geometry-dependent multipole models. In this study, atomic force expressions for geometry-dependent multipoles are presented for use in simulations of flexible molecules. The atomic forces are derived by first proposing a new general expression for Wigner function derivatives partial derivative D(m'm)(l)/partial derivative Omega. The force equations can be applied to electrostatic models based on atomic point multipoles or gaussian multipole charge density. Hydrogen-bonded dimers are used to test the intermolecular electrostatic energies and atomic forces calculated by geometry-dependent multipoles fit to the ab initio electrostatic potential. The electrostatic energies and forces are compared with their reference ab initio values. It is shown that both static and geometry-dependent multipole models are able to reproduce total molecular forces and torques with respect to ab initio, whereas geometry-dependent multipoles are needed to reproduce ab initio atomic forces. The expressions for atomic force can be used in simulations of flexible molecules with atomic multipoles. In addition, the results presented in this work should lead to further development of next generation force fields composed of geometry-dependent multipole models.
| S-EPMC2941241 | biostudies-literature
Project description:We examine the ability of observers to extract summary statistics (such as the mean and the relative-variance) from rapid numerical sequences of two digit numbers presented at a rate of 4/s. In four experiments (total N = 100), we find that the participants show a remarkable ability to extract such summary statistics and that their precision in the estimation of the sequence-mean improves with the sequence-length (subject to individual differences). Using model selection for individual participants we find that, when only the sequence-average is estimated, most participants rely on a holistic process of frequency based estimation with a minority who rely on a (rule-based and capacity limited) mid-range strategy. When both the sequence-average and the relative variance are estimated, about half of the participants rely on these two strategies. Importantly, the holistic strategy appears more efficient in terms of its precision. We discuss implications for the domains of two pathways numerical processing and decision-making.
| S-EPMC8517333 | biostudies-literature
Project description:Compressed nano-pillars crackle from moving dislocations, which reduces plastic stability. Crackling noise is characterized by stress drops or strain bursts, which scale over a large region of sizes leading to power law statistics. Here we report that this "classic" behaviour is not valid in Ti-based nanopillars for a counterintuitive reason: we tailor precipitates inside the nano-pillar, which "regulate" the flux of dislocations. It is not because the nano-pillars become too small to sustain large dislocation movements, the effect is hence independent of size. Our precipitates act as "rotors": local stress initiates the rotation of inclusions, which reduces the stress amplitudes dramatically. The size distribution of stress drops simultaneously changes from power law to exponential. Rotors act like revolving doors limiting the number of passing dislocations. Hence each collapse becomes weak. We present experimental evidence for Ti-based nano-pillars (diameters between 300 nm and 2 μm) with power law distributions of crackling noise P(s) ∼ s-τ with τ ∼ 2 in the defect free or non-rotatable precipitate states. Rotors change the size distribution to P(s) ∼ exp(-s/s0). Rotors are inclusions of ω-phase that aligns under stress along slip planes and limit dislocation glide to small distances with high nucleation rates. This opens new ways to make nano-pillars more stable.
| S-EPMC6405840 | biostudies-other
Project description:The identification of the limits between the cell division, elongation and mature zones in the root apex is still a matter of controversy when methods based on cellular features, molecular markers or kinematics are compared while methods based on cell length profiles have been comparatively underexplored. Segmentation models were developed to identify developmental zones within a root apex on the basis of epidermal cell length profiles. Heteroscedastic piecewise linear models were estimated for maize lateral roots of various lengths of both wild type and two mutants affected in auxin signaling (rtcs and rum-1). The outputs of these individual root analyses combined with morphological features (first root hair position and root diameter) were then globally analyzed using principal component analysis. Three zones corresponding to the division zone, the elongation zone and the mature zone were identified in most lateral roots while division zone and sometimes elongation zone were missing in arrested roots. Our results are consistent with an auxin-dependent coordination between cell flux, cell elongation and cell differentiation. The proposed segmentation models could extend our knowledge of developmental regulations in longitudinally organized plant organs such as roots, monocot leaves or internodes.
| S-EPMC5662930 | biostudies-literature