PubMed HealthSearch

PubMed · 8733949

A nonparametric subject-specific population method for deconvolution: II. External validation.

Abstract

A lot of attention has been given in the past to deconvolution and in particular to its nonparametric variants. In a companion paper (1), we present a fully nonparametric deconvolution method in which subject specificity is explicitly taken into account. To do so we use so-called "longitudinal splines." A longitudinal spline is a nonparametric function composed of a template spline, in common to all subjects, and of a distortion spline representing the difference of the subject's function from the template. In this paper we concentrate on testing and documenting the performance of this nonparametric methodology in terms of the approximation of unknown functions. We simulate population data using parametric functions, and use longitudinal splines to recover the unknown functions. We consider different estimation methods including (1) parametric nonlinear mixed effect, (2) least squares, and (3) two-stage. Methods 2-3 are more robust than Method 1, and obtain reliable estimates of the unknown functions. The lack of robustness of Method 1 appears to be due to the misspecifications of the distribution of the subjects' parameters. Results also suggest that in a data-rich situation nonparametric nonlinear mixed-effect models should be preferred.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

K E Fattinger, D Verotta. 1995. A nonparametric subject-specific population method for deconvolution: II. External validation.. https://doi.org/10.1007/bf02353464

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Stochastic resonance in non-dynamical systems without response thresholds.

The addition of noise to a system can sometimes improve its ability to transfer information reliably. This phenomenon--known as stochastic resonance--was originally proposed to account for periodicity in the Earth's ice ages, but has now been shown to occur in many systems in physics and biology. Recent experimental and theoretical work has shown that the simplest system exhibiting 'stochastic resonance' consists of nothing more than signal and noise with a threshold-triggered device (when the signal plus noise exceeds the threshold, the system responds momentarily, then relaxes to equilibrium to await the next triggering event). Here we introduce a class of non-dynamical and threshold-free systems that also exhibit stochastic resonance. We present and analyse a general mathematical model for such systems, in which a sequence of pulses is generated randomly with a probability (per unit time) that depends exponentially on an input. When this input is a sine-wave masked by additive noise, we observe an increase in the output signal-to-noise ratio as the level of noise increases. This result shows that stochastic resonance can occur in a broad class of thermally driven physico-chemical systems, such as semiconductor p-n junctions, mesoscopic electronic devices and voltage-dependent ion channels, in which reaction rates are controlled by activation barriers.

Models, Theoretical

Equivalent methods to analyse dynamic experiments in which the input function is noisy.

A comparison is made between two methods of parameter estimation for analysis of dynamic experiments in which the input function is noisy. Noise in the input function leads to uncertainties in the calculated model-predicted values, and therefore the covariance matrix of the residuals is a function of the model parameters. Statistical uncertainties in the model-predicted values significantly change the nature of the fitting process and the quality of the results. The initial method uses a weighted least-squares criterion where the weighting matrix is the inverse of the full covariance matrix of the residuals, incorporating both the noise in the output data and the noise in the input function. The methodology was applied to dynamic emission tomography studies of the heart, where the blood (input) and tissue (output) tracer concentrations at each time are derived from two regions of interest in the same tomographic section. The second method introduces additional parameters to describe the input function, and adds terms to the weighted sum of squares which comprise the criterion. Instead of only summing the weighted terms to account for differences between the model and the output function, there is a second set of terms to account for the differences between the model and the input function. The two methods have different theoretical bases and appear to optimize different criteria, but it is shown here that they are equivalent to one another. The criterion which they minimize is the same under certain matrix invertibility constraints, which must be satisfied to ensure the stability of either method.

Models, Theoretical