Parameter estimation in linear models based on outage probability minimization

Author(s):  
Sergiy A. Vorobyov ◽  
Yonina C. Eldar ◽  
Alex B. Gershman
Author(s):  
Tim Loossens ◽  
Kristof Meers ◽  
Niels Vanhasbroeck ◽  
Nil Anarat ◽  
Stijn Verdonck ◽  
...  

AbstractComputational modeling plays an important role in a gamut of research fields. In affect research, continuous-time stochastic models are becoming increasingly popular. Recently, a non-linear, continuous-time, stochastic model has been introduced for affect dynamics, called the Affective Ising Model (AIM). The drawback of non-linear models like the AIM is that they generally come with serious computational challenges for parameter estimation and related statistical analyses. The likelihood function of the AIM does not have a closed form expression. Consequently, simulation based or numerical methods have to be considered in order to evaluate the likelihood function. Additionally, the likelihood function can have multiple local minima. Consequently, a global optimization heuristic is required and such heuristics generally require a large number of likelihood function evaluations. In this paper, a Julia software package is introduced that is dedicated to fitting the AIM. The package includes an implementation of a numeric algorithm for fast computations of the likelihood function, which can be run both on graphics processing units (GPU) and central processing units (CPU). The numerical method introduced in this paper is compared to the more traditional Euler-Maruyama method for solving stochastic differential equations. Furthermore, the estimation software is tested by means of a recovery study and estimation times are reported for benchmarks that were run on several computing devices (two different GPUs and three different CPUs). According to these results, a single parameter estimation can be obtained in less than thirty seconds using a mainstream NVIDIA GPU.


2019 ◽  
Vol 7 (1) ◽  
pp. 78-91
Author(s):  
Stephen Haslett

Abstract When sample survey data with complex design (stratification, clustering, unequal selection or inclusion probabilities, and weighting) are used for linear models, estimation of model parameters and their covariance matrices becomes complicated. Standard fitting techniques for sample surveys either model conditional on survey design variables, or use only design weights based on inclusion probabilities essentially assuming zero error covariance between all pairs of population elements. Design properties that link two units are not used. However, if population error structure is correlated, an unbiased estimate of the linear model error covariance matrix for the sample is needed for efficient parameter estimation. By making simultaneous use of sampling structure and design-unbiased estimates of the population error covariance matrix, the paper develops best linear unbiased estimation (BLUE) type extensions to standard design-based and joint design and model based estimation methods for linear models. The analysis covers both with and without replacement sample designs. It recognises that estimation for with replacement designs requires generalized inverses when any unit is selected more than once. This and the use of Hadamard products to link sampling and population error covariance matrix properties are central topics of the paper. Model-based linear model parameter estimation is also discussed.


2015 ◽  
Vol 3 (1-2) ◽  
pp. 52-87 ◽  
Author(s):  
Nori Jacoby ◽  
Naftali Tishby ◽  
Bruno H. Repp ◽  
Merav Ahissar ◽  
Peter E. Keller

Linear models have been used in several contexts to study the mechanisms that underpin sensorimotor synchronization. Given that their parameters are often linked to psychological processes such as phase correction and period correction, the fit of the parameters to experimental data is an important practical question. We present a unified method for parameter estimation of linear sensorimotor synchronization models that extends available techniques and enhances their usability. This method enables reliable and efficient analysis of experimental data for single subject and multi-person synchronization. In a previous paper (Jacoby et al., 2015), we showed how to significantly reduce the estimation error and eliminate the bias of parameter estimation methods by adding a simple and empirically justified constraint on the parameter space. By applying this constraint in conjunction with the tools of matrix algebra, we here develop a novel method for estimating the parameters of most linear models described in the literature. Through extensive simulations, we demonstrate that our method reliably and efficiently recovers the parameters of two influential linear models: Vorberg and Wing (1996), and Schulze et al. (2005), together with their multi-person generalization to ensemble synchronization. We discuss how our method can be applied to include the study of individual differences in sensorimotor synchronization ability, for example, in clinical populations and ensemble musicians.


2015 ◽  
Vol 3 (1-2) ◽  
pp. 32-51 ◽  
Author(s):  
Nori Jacoby ◽  
Peter E. Keller ◽  
Bruno H. Repp ◽  
Merav Ahissar ◽  
Naftali Tishby

The mechanisms that support sensorimotor synchronization — that is, the temporal coordination of movement with an external rhythm — are often investigated using linear computational models. The main method used for estimating the parameters of this type of model was established in the seminal work of Vorberg and Schulze (2002), and is based on fitting the model to the observed auto-covariance function of asynchronies between movements and pacing events. Vorberg and Schulze also identified the problem of parameter interdependence, namely, that different sets of parameters might yield almost identical fits, and therefore the estimation method cannot determine the parameters uniquely. This problem results in a large estimation error and bias, thereby limiting the explanatory power of existing linear models of sensorimotor synchronization. We present a mathematical analysis of the parameter interdependence problem. By applying the Cramér–Rao lower bound, a general lower bound limiting the accuracy of any parameter estimation procedure, we prove that the mathematical structure of the linear models used in the literature determines that this problem cannot be resolved by any unbiased estimation method without adopting further assumptions. We then show that adding a simple and empirically justified constraint on the parameter space — assuming a relationship between the variances of the noise terms in the model — resolves the problem. In a follow-up paper in this volume, we present a novel estimation technique that uses this constraint in conjunction with matrix algebra to reliably estimate the parameters of almost all linear models used in the literature.


2014 ◽  
Vol 43 (3) ◽  
pp. 181-193 ◽  
Author(s):  
Roland Fried ◽  
Tobias Liboschik ◽  
Hanan Elsaied ◽  
Stella Kitromilidou ◽  
Konstantinos Fokianos

We discuss the analysis of count time series following generalised linear models in the presence of outliers and intervention effects. Different modifications of such models are formulated which allow to incorporate, detect and to a certain degree distinguish extraordinary events (interventions) of different types in count time series retrospectively. An outlook on extensions to the problem of robust parameter estimation, identification of the model orders by robust estimation of autocorrelations and partial autocorrelations, and online surveillance by sequential testing for outlyingness is provided. 


Sign in / Sign up

Export Citation Format

Share Document