## The Annals of Statistics

- Ann. Statist.
- Volume 28, Number 5 (2000), 1279-1301.

### Parametric deconvolution of positive spike trains

#### Abstract

This paper describes a parametric deconvolution method (PDPS) appropriate for a particular class of signals which we call spike-convolution models. These models arise when a sparse spike train —Dirac deltas according to our mathematical treatment —is convolved with a fixed point-spread function, and additive noise or measurement error is superimposed.We view deconvolution as an estimation problem, regarding the locations and heights of the underlying spikes, as well as the baseline and the measurement error variance as unknown parameters.Our estimation scheme consists of two parts: model fitting and model selection.To fit a spike-convolution model of a specific order, we estimate peak locations by trigonometric moments, and heights and the baseline by least squares. The model selection procedure has two stages. Its first stage is so designed that we expect a model of a somewhat larger order than the truth to be selected. In the second stage, the final model is obtained using backwards deletion. This results in not only an estimate of the model order, but also an estimate of peak locations and heights with much smaller bias and variation than that found in a direct trigonometric moment estimate. A more efficient maximum likelihood estimate can be calculated from these estimates using a Gauss–Newton algorithm. We also present some relevant results concerning the spectral structure of Toeplitz matrices which play a key role in the estimation. Finally, we illustrate the behavior of these estimates using simulated and real DNA sequencing data.

#### Article information

**Source**

Ann. Statist., Volume 28, Number 5 (2000), 1279-1301.

**Dates**

First available in Project Euclid: 12 March 2002

**Permanent link to this document**

https://projecteuclid.org/euclid.aos/1015957394

**Digital Object Identifier**

doi:10.1214/aos/1015957394

**Mathematical Reviews number (MathSciNet)**

MR1805784

**Zentralblatt MATH identifier**

1105.62382

**Subjects**

Primary: 62F10: Point estimation

Secondary: 62F12: Asymptotic properties of estimators 86A22: Inverse problems [See also 35R30]

**Keywords**

Deconvolution spike train model selection DNA sequencing Toeplitz matrix

#### Citation

Li, Lei; Speed, Terence P. Parametric deconvolution of positive spike trains. Ann. Statist. 28 (2000), no. 5, 1279--1301. doi:10.1214/aos/1015957394. https://projecteuclid.org/euclid.aos/1015957394