scholarly journals Discrete Time Rescaling Theorem: Determining Goodness of Fit for Discrete Time Statistical Models of Neural Spiking

2010 ◽  
Vol 22 (10) ◽  
pp. 2477-2506 ◽  
Author(s):  
Robert Haslinger ◽  
Gordon Pipa ◽  
Emery Brown

One approach for understanding the encoding of information by spike trains is to fit statistical models and then test their goodness of fit. The time-rescaling theorem provides a goodness-of-fit test consistent with the point process nature of spike trains. The interspike intervals (ISIs) are rescaled (as a function of the model's spike probability) to be independent and exponentially distributed if the model is accurate. A Kolmogorov-Smirnov (KS) test between the rescaled ISIs and the exponential distribution is then used to check goodness of fit. This rescaling relies on assumptions of continuously defined time and instantaneous events. However, spikes have finite width, and statistical models of spike trains almost always discretize time into bins. Here we demonstrate that finite temporal resolution of discrete time models prevents their rescaled ISIs from being exponentially distributed. Poor goodness of fit may be erroneously indicated even if the model is exactly correct. We present two adaptations of the time-rescaling theorem to discrete time models. In the first we propose that instead of assuming the rescaled times to be exponential, the reference distribution be estimated through direct simulation by the fitted model. In the second, we prove a discrete time version of the time-rescaling theorem that analytically corrects for the effects of finite resolution. This allows us to define a rescaled time that is exponentially distributed, even at arbitrary temporal discretizations. We demonstrate the efficacy of both techniques by fitting generalized linear models to both simulated spike trains and spike trains recorded experimentally in monkey V1 cortex. Both techniques give nearly identical results, reducing the false-positive rate of the KS test and greatly increasing the reliability of model evaluation based on the time-rescaling theorem.

2011 ◽  
Vol 23 (6) ◽  
pp. 1452-1483 ◽  
Author(s):  
Felipe Gerhard ◽  
Robert Haslinger ◽  
Gordon Pipa

Statistical models of neural activity are integral to modern neuroscience. Recently interest has grown in modeling the spiking activity of populations of simultaneously recorded neurons to study the effects of correlations and functional connectivity on neural information processing. However, any statistical model must be validated by an appropriate goodness-of-fit test. Kolmogorov-Smirnov tests based on the time-rescaling theorem have proven to be useful for evaluating point-process-based statistical models of single-neuron spike trains. Here we discuss the extension of the time-rescaling theorem to the multivariate (neural population) case. We show that even in the presence of strong correlations between spike trains, models that neglect couplings between neurons can be erroneously passed by the univariate time-rescaling test. We present the multivariate version of the time-rescaling theorem and provide a practical step-by-step procedure for applying it to testing the sufficiency of neural population models. Using several simple analytically tractable models and more complex simulated and real data sets, we demonstrate that important features of the population activity can be detected only using the multivariate extension of the test.


Author(s):  
Shannon Wagner ◽  
John B. Ferris

Terrain topology is the principal source of vertical excitation into the vehicle system and must be accurately represented in order to correctly predict the vehicle response. It is desirable to evaluate vehicle models over a wide range of terrain, but it is computationally impractical to simulate long distances of every terrain type. A method to parsimoniously characterize terrain topology is developed in this work so that terrain can be grouped into meaningful sets with similar topological characteristics. Specifically, measured terrain profiles are considered realizations of an underlying stochastic process; an autoregressive model and a residual process provide the mathematical framework to describe this process. A statistical test is developed to determine if the residual process is independent and identically distributed (IID) and, therefore, stationary. A reference joint probability distribution of the residuals is constructed based on the assumption that the data are realizations of an IID stochastic process. The distribution of the residuals is then compared to this reference distribution via the Kolmogorov–Smirnov “goodness of fit” test to determine whether the IID assumption is valid. If the residual process is IID, a single probability distribution can be used to generate residuals and synthetic terrain of any desired length. This modeling method and statistical test are applied to a set of U.S. highway profile data and show that the residual process can be assumed to be IID in virtually all of these cases of nondeformable terrain surfaces.


2009 ◽  
Vol 24 (S1) ◽  
pp. 1-1 ◽  
Author(s):  
N. Wu ◽  
X. Zhang ◽  
L. Ye ◽  
Q. Xu ◽  
S. Jin ◽  
...  

Background:An increasing number of studies have described the relationship between velo-cardio-facial syndrome (VCFS) and schizophrenia. in a family-based study, we found that rs10314, a single nucleotide polymorphism (SNP) present in the 3’-flanking region of the CLDN5 gene, was associated with schizophrenia among a Chinese population. High false positive rate is a common problem with the association study of human diseases. It is very important to replicate an initial finding with different samples and experimental designs.Methods:A total of 749 patients with schizophrenia and 383 age and sex matched healthy control subjects in Chinese population were recruited. PCR-based RFLP protocol was applied to genotype rs10314 to see its disease association.Results:The χ2 goodness-of-fit test showed that the genotypic distributions of rs10314 were in Hardy-Weinberg equilibrium in both the patient group (χ2=1.12, P=0.289) and the control group (χ2=0.22, P=0.639). rs10314 was associated with schizophrenia with an odds ratio (OR) of 1.32 in the male subjects (χ2=5.45, P=0.02, 95% CI 1.05-1.67) but not in the female subjects (χ2=0.64, P=0.425, OR=1.14, 95% CI 0.83-1.57). the χ2 test showed a genotypic association only for combined samples (χ2=7.80, df=2, P=0.02). SNP rs10314 is a G to C base change. Frequency of the genotypes containing the C allele was significantly higher in the patient group than in the control group.Conclusions:The present work shows that the CLDN5 gene polymorphism is more likely to be involved in schizophrenic men than women, suggesting that this gene may contribute to the gender differences in schizophrenia.


Author(s):  
Xie Y ◽  
◽  
Dong H ◽  
Liao Y ◽  
Zhang J ◽  
...  

Background: COVID-19 nucleic acid swab tests have a high false positive rate; therefore, diagnosing COVID-19 pneumonia and predicting prognosis by CT scan are very important. Methods: In this retrospective single-centre study, we included consecutive suspected critical COVID-19 pneumonia cases in the intensive care unit of Wuhan Third Hospital from January 31, 2020, to March 16, 2020. 204 cases were confirmed by real-time RT-PCR, and all patients were evaluated with CT, cut-off values were obtained according to the Youden index and were divided into a high CT score group and a low CT score group. Epidemiological, demographic, clinical, and laboratory data were collected. Finally, Through multi-factor logistic regression model, a prediction model based on multiple prediction indicators was formed, and new joint predictive factors were calculated. The prediction model of mortality in COVID-19 pneumonia based on CT score and lymphocyte count was constructed through data processing analysis. Results: The major imaging feature of COVID-19 pneumonia is Ground Glass Opacities (GGOs). Multivariate regression analysis found that the CT score and absolute lymphocyte count were independent risk factors for death and that the CT score predicted mortality (AUC-ROC =0.7, cut-off=1.45). When the absolute lymphocyte count was lower, the patient’s CT score was also lower. Based on this, a prediction model was established. The prediction model was: In [P/(1-P)]=0.667*gender+0.057*age-0.086CT score-0.831 lymphocyte count-3.91, the goodness of fit test of the model was P=0.041, and the area under the curve of the ROC curve of the model was 0.779. Conclusion: CT score and absolute lymphocyte count are independent risk factors for mortality, and patients with a high CT score may have a worse prognosis. A lower absolute lymphocyte count may indicate that the patient’s CT score is also reduced. The model established by combining CT scores and lymphocyte count showed a good degree of calibration and differentiation.


Automatica ◽  
2021 ◽  
Vol 131 ◽  
pp. 109742
Author(s):  
Alexis J. Vallarella ◽  
Paula Cardone ◽  
Hernan Haimovich

Modelling ◽  
2021 ◽  
Vol 2 (1) ◽  
pp. 78-104
Author(s):  
Vasili B. V. Nagarjuna ◽  
R. Vishnu Vardhan ◽  
Christophe Chesneau

Every day, new data must be analysed as well as possible in all areas of applied science, which requires the development of attractive statistical models, that is to say adapted to the context, easy to use and efficient. In this article, we innovate in this direction by proposing a new statistical model based on the functionalities of the sinusoidal transformation and power Lomax distribution. We thus introduce a new three-parameter survival distribution called sine power Lomax distribution. In a first approach, we present it theoretically and provide some of its significant properties. Then the practicality, utility and flexibility of the sine power Lomax model are demonstrated through a comprehensive simulation study, and the analysis of nine real datasets mainly from medicine and engineering. Based on relevant goodness of fit criteria, it is shown that the sine power Lomax model has a better fit to some of the existing Lomax-like distributions.


Sign in / Sign up

Export Citation Format

Share Document