An inequality-constrained least-squares deconvolution method

The output-function (Y) of a linear system is the convolution of the input function (I) to the system with the disposition-function (H) of the system. Given Y and H deconvolution yields I. A non-parametric method for numerical deconvolution is described. The method is based on an inequality-constrained least-squares criterion and approximates I by a discontinuous function. No assumptions are made about the form of H or Y. Numerical stability and physical realism are obtained by constraining the estimated I to be nonnegative and piecewise-monotonic (nonincreasing, nondecreasing, or alternating segments of both). When I is constrained to be monotonic, the deconvolution yields a staircase function. The method can be used to calculate drug input rates. It is compared to previously published deconvolution methods for this purpose, using simulated data and real theophylline and pentobarbital data.

[1]  B. Hunt.  The inverse problem of radiography , 1970 .

[2]  Peter Veng-Pedersen,et al.  An algorithm and computer program for deconvolution in linear pharmacokinetics , 1980, Journal of Pharmacokinetics and Biopharmaceutics.

[3]  L. Dittert,et al.  Pharmacokinetics of pentobarbital after intravenous and oral administration , 1973, Journal of Pharmacokinetics and Biopharmaceutics.

[4]  M. E. Muller,et al.  A Note on the Generation of Random Normal Deviates , 1958 .

[5]  D P Vaughan,et al.  Mathematical basis of point-area deconvolution method for determining in vivo input functions. , 1978, Journal of pharmaceutical sciences.

[6]  S. Riegelman,et al.  New method for calculating the intrinsic absorption rate of drugs. , 1968, Journal of pharmaceutical sciences.

[7]  S. Van Huffel,et al.  Reliable and efficient deconvolution technique based on total linear least squares for calculating the renal retention function , 2006, Medical and Biological Engineering and Computing.

[8]  M. J. Berg,et al.  Linear systems approach to the analysis of an induced drug removal process. Phenobarbital removal by oral activated charcoal , 1986, Journal of Pharmacokinetics and Biopharmaceutics.

[9]  Yasuaki Ogawa,et al.  Estimation of drug absorption rates using a deconvolution method with nonequal sampling times , 1986, Journal of Pharmacokinetics and Biopharmaceutics.

[10]  A. Staib,et al.  A crossover study after oral and intravenous administration of theophylline in male volunteers (absolute bioavailability of Afonilum tablets). , 1983, International journal of clinical pharmacology, therapy and toxicology.

[11]  Grace Wahba,et al.  Constrained Regularization for Ill Posed Linear Operator Equations, with Applications in Meteorology and Medicine. , 1982 .

[12]  D. Cutler,et al.  Numerical deconvolution by least squares: Use of prescribed input functions , 1978, Journal of Pharmacokinetics and Biopharmaceutics.

[13]  A note on appropriate constraints on the initial input response when applying deconvolution , 1986, Journal of Pharmacokinetics and Biopharmaceutics.

[14]  B. R. Frieden.  Image enhancement and restoration , 1979 .

[15]  Numerical deconvolution by least squares: Use of polynomials to represent the input function , 1978, Journal of Pharmacokinetics and Biopharmaceutics.

[16]  L B Sheiner,et al.  Semiparametric approach to pharmacokinetic-pharmacodynamic data. , 1989, The American journal of physiology.

[17]  P V Pedersen,et al.  Novel deconvolution method for linear pharmacokinetic systems with polyexponential impulse response. , 1980, Journal of pharmaceutical sciences.

[18]  P. Pedersen,et al.  Model-independent method of analyzing input in linear pharmacokinetic systems having polyexponential impulse response I: Theoretical analysis. , 1980, Journal of pharmaceutical sciences.

[19]  B. R. Hunt,et al.  Biased estimation for nonparametric identification of linear systems , 1971 .

[20]  Peter Veng-Pedersen,et al.  A polyexponential deconvolution method. Evaluation of the “gastrointestinal bioavailability” and meanin vivo dissolution time of some ibuprofen dosage forms , 1985, Journal of Pharmacokinetics and Biopharmaceutics.

[21]  P. Jansson,et al.  Resolution Enhancement of Spectra , 1970 .

[22]  J. M. Bremner,et al.  Statistical Inference under Restrictions , 1973 .

[23]  Donald Goldfarb,et al.  A numerically stable dual method for solving strictly convex quadratic programs , 1983, Math. Program..

[24]  H. D. Brunk,et al.  Statistical inference under order restrictions : the theory and application of isotonic regression , 1973 .

[25]  S. Gull,et al.  Maximum entropy and its application to the calculation of drug absorption rates , 1987, Journal of Pharmacokinetics and Biopharmaceutics.