ACCURACY OF NUMERICAL INVERSION OF LAPLACE TRANSFORMS FOR PHARMACOKINETIC PARAMETER-ESTIMATION

被引:23
作者
PURVES, RD
机构
[1] Department of Pharmacology, University of Otago Medical School, Dunedin
关键词
D O I
10.1002/jps.2600840117
中图分类号
R914 [药物化学];
学科分类号
100701 ;
摘要
Numerical inversion of the Laplace transform is a useful technique for pharmacokinetic modeling and parameter estimation when the model equations can be solved in the Laplace domain but the solutions cannot be inverted back to the time domain. The accuracy of numerical inversion of the Laplace transform using an infinite series approximation due to Hosono was systematically studied by reference to 17 widely differing functions having known inverse transforms. The error of inversion was found to be very sensitive to the details of the computer implementation of the method; for example, double-precision arithmetic is essential. The method used to sum the series in the least-squares program MULTI(FILT) was often unable to achieve a relative error of less than 10(-4), and a Monte Carlo simulation showed that this method is insufficiently accurate for reliable least-squares parameter estimation. Improvements to the algorithm are described whereby a better method of applying Euler's transformation is used and the number of terms summed is determined automatically by the rate of convergence of the series. The improved algorithm is more efficient in inverting easy functions and more reliable in inverting difficult functions, especially those involving a time lag. With its use, pharmacokinetic parameter estimation can be performed with essentially the same accuracy as when the function is defined in the time domain.
引用
收藏
页码:71 / 74
页数:4
相关论文
共 20 条
[1]   NUMERICAL INVERSION OF THE LAPLACE TRANSFORM - SURVEY AND COMPARISON OF METHODS [J].
DAVIES, B ;
MARTIN, B .
JOURNAL OF COMPUTATIONAL PHYSICS, 1979, 33 (01) :1-32
[2]   APPLICATION OF THE AXIAL-DISPERSION MODEL OF HEPATIC DRUG ELIMINATION TO THE KINETICS OF DIAZEPAM IN THE ISOLATED PERFUSED-RAT-LIVER [J].
DIAZGARCIA, JM ;
EVANS, AM ;
ROWLAND, M .
JOURNAL OF PHARMACOKINETICS AND BIOPHARMACEUTICS, 1992, 20 (02) :171-193
[3]   INFLUENCE OF ALBUMIN ON THE DISTRIBUTION AND ELIMINATION KINETICS OF DICLOFENAC IN THE ISOLATED PERFUSED-RAT-LIVER - ANALYSIS BY THE IMPULSE-RESPONSE TECHNIQUE AND THE DISPERSION MODEL [J].
EVANS, AM ;
HUSSEIN, Z ;
ROWLAND, M .
JOURNAL OF PHARMACEUTICAL SCIENCES, 1993, 82 (04) :421-428
[4]   A 2-COMPARTMENT DISPERSION MODEL DESCRIBES THE HEPATIC OUTFLOW PROFILE OF DICLOFENAC IN THE PRESENCE OF ITS BINDING-PROTEIN [J].
EVANS, AM ;
HUSSEIN, Z ;
ROWLAND, M .
JOURNAL OF PHARMACY AND PHARMACOLOGY, 1991, 43 (10) :709-714
[5]   NUMERICAL INVERSION OF LAPLACE TRANSFORM AND SOME APPLICATIONS TO WAVE OPTICS [J].
HOSONO, T .
RADIO SCIENCE, 1981, 16 (06) :1015-1019
[6]  
KOMURA H, 1992, DRUG METAB DISPOS, V20, P585
[7]   ANALYSIS OF DRUG PENETRATION THROUGH THE SKIN BY THE 2-LAYER SKIN MODEL [J].
OKAMOTO, H ;
YAMASHITA, F ;
SAITO, K ;
HASHIDA, M .
PHARMACEUTICAL RESEARCH, 1989, 6 (11) :931-937
[8]  
Press W., 1992, NUMERICAL RECIPES FO
[9]   ANOMALOUS PARAMETER ESTIMATES IN THE ONE-COMPARTMENT MODEL WITH 1ST-ORDER ABSORPTION [J].
PURVES, RD .
JOURNAL OF PHARMACY AND PHARMACOLOGY, 1993, 45 (10) :934-936