A copula-based method to build diffusion models with prescribed marginal and serial dependence

This paper investigates the probabilistic properties that determine the existence of space-time transformations between diffusion processes. We prove that two diffusions are related by a monotone space-time transformation if and only if they share the same serial dependence. The serial dependence of a diffusion process is studied by means of its copula density and the effect of monotone and non-monotone space-time transformations on the copula density is discussed. This provides us a methodology to build diffusion models by freely combining prescribed marginal behaviors and temporal dependence structures. Explicit expressions of copula densities are provided for tractable models. A possible application in neuroscience is sketched as a proof of concept

Estimation in discretely observed diffusions killed at a threshold

Parameter estimation in diffusion processes from discrete observations up to a first-hitting time is clearly of practical relevance, but does not seem to have been studied so far. In neuroscience, many models for the membrane potential evolution involve the presence of an upper threshold. Data are modeled as discretely observed diffusions which are killed when the threshold is reached. Statistical inference is often based on the misspecified likelihood ignoring the presence of the threshold causing severe bias, e.g. the bias incurred in the drift parameters of the Ornstein-Uhlenbeck model for biological relevant parameters can be up to 25-100%. We calculate or approximate the likelihood function of the killed process. When estimating from a single trajectory, considerable bias may still be present, and the distribution of the estimates can be heavily skewed and with a huge variance. Parametric bootstrap is effective in correcting the bias. Standard asymptotic results do not apply, but consistency and asymptotic normality may be recovered when multiple trajectories are observed, if the mean first-passage time through the threshold is finite. Numerical examples illustrate the results and an experimental data set of intracellular recordings of the membrane potential of a motoneuron is analyzed.Comment: 29 pages, 5 figure

Estimating input parameters from intracellular recordings in the Feller neuronal model

We study the estimation of the input parameters in a Feller neuronal model from a trajectory of the membrane potential sampled at discrete times. These input parameters are identified with the drift and the infinitesimal variance of the underlying stochastic diffusion process with multiplicative noise. The state space of the process is restricted from below by an inaccessible boundary. Further, the model is characterized by the presence of an absorbing threshold, the first hitting of which determines the length of each trajectory and which constrains the state space from above. We compare, both in the presence and in the absence of the absorbing threshold, the efficiency of different known estimators. In addition, we propose an estimator for the drift term, which is proved to be more efficient than the others, at least in the explored range of the parameters. The presence of the threshold makes the estimates of the drift term biased, and two methods to correct it are proposed

Applications of The Reflected Ornstein-Uhlenbeck Process

An Ornstein-Uhlenbeck process is the most basic mean-reversion model and has been used in various fields such as finance and biology. In some instances, reflecting boundary conditions are needed to restrict the state space of this process. We study an Ornstein-Uhlenbeck diffusion process with a reflecting boundary and its application to finance and neuroscience. In the financial application, the Vasicek model which is an Ornstein-Uhlenbeck process has been used to capture the stochastic movement of the short term interest rate in the market. The shortcoming of applying this model is that it allows a negative interest rate theoretically. Thus we use a reflected Ornstein-Uhlenbeck process as an interest rate model to get around this problem. Then we price zero-coupon bond and European options with respect to our model. In the application to neuroscience, we study integrate-and-fire (I-F) neuron models. We assume that the membrane voltage follows a reflected Ornstein-Uhlenbeck process and fires when it reaches a threshold. In this case, the interspike intervals (ISIs) are the same as the first hitting times of the process to a certain barrier. We find the first passage time density given ISIs using numerical inversion integration of the Laplace transform of the first passage time pdf. Then we estimate the unknown identifiable parameters in our model

Motoneuron membrane potentials follow a time inhomogeneous jump diffusion process

Stochastic leaky integrate-and-fire models are popular due to their simplicity and statistical tractability. They have been widely applied to gain understanding of the underlying mechanisms for spike timing in neurons, and have served as building blocks for more elaborate models. Especially the Ornstein–Uhlenbeck process is popular to describe the stochastic fluctuations in the membrane potential of a neuron, but also other models like the square-root model or models with a non-linear drift are sometimes applied. Data that can be described by such models have to be stationary and thus, the simple models can only be applied over short time windows. However, experimental data show varying time constants, state dependent noise, a graded firing threshold and time-inhomogeneous input. In the present study we build a jump diffusion model that incorporates these features, and introduce a firing mechanism with a state dependent intensity. In addition, we suggest statistical methods to estimate all unknown quantities and apply these to analyze turtle motoneuron membrane potentials. Finally, simulated and real data are compared and discussed. We find that a square-root diffusion describes the data much better than an Ornstein–Uhlenbeck process with constant diffusion coefficient. Further, the membrane time constant decreases with increasing depolarization, as expected from the increase in synaptic conductance. The network activity, which the neuron is exposed to, can be reasonably estimated to be a threshold version of the nerve output from the network. Moreover, the spiking characteristics are well described by a Poisson spike train with an intensity depending exponentially on the membrane potential