Discrete approximation of multivariate densities with application to Bayesian estimation

doi:10.1016/0167-9473(84)90030-6

Computational Statistics & Data Analysis

Volume 2, Issue 1, June 1984, Pages 27-36

https://doi.org/10.1016/0167-9473(84)90030-6 Get rights and content

Abstract

A procedure is presented for finding a discrete approximation to a continuous multivariate density function. It is based on a previously developed algorithm [2] for determining the L₁ optimal discrete approximation to a univariate density. Results of approximating continuous bivariate density functions, which represent distributions of the parameters of a pharmacokinetic model, show good agreement between the mean and covariance matrix of the approximated and approximating densities. The distribution of a predicted drug conceptration was also calculated using a continuous density and discrete approximations with both 25 and 81 points. The expected values of the predicted concentration, as well as selected percentile points, obtained using each density are in close agreement.

Reference (6)

KatzD.
Discrete approximations to continuous density functions that are L₁ optimal
Computat. Statist. Data Anal.
(1983)
KatzD. et al.
Reduction of dimensionality in Bayesian nonlinear regression with a pharmacokinetic application
Math. Biosci.
(1982)
KatzD. et al.
Bayesian approach to the analysis of nonlinear models: Implementation and evaluation
Biometrics
(1981)

There are more references available in the full text version of this article.

Cited by (14)

Sampling schedule design towards optimal drug monitoring for individualizing therapy
2005, Computer Methods and Programs in Biomedicine
We study the individualization of therapy by simultaneously taking into account the design of sampling schedule and optimal therapeutic drug monitoring. The sampling schedule design in this work is to determine the number of samples, the sampling times, the switching time from the loading to the maintenance period, and the drug dosages. A closed-loop control policy is employed to determine the sampling schedule, and an advanced stochastic global optimization algorithm, which integrates the stochastic approximation and simulated annealing techniques, is implemented to search the optimal sampling schedule. A simulated one-compartment model of intravenous theophylline therapy is used to illustrate our method. This method can be readily extended to multiple compartment systems and allow incorporating other criteria of drug control. While currently the method is mainly of theoretical interest, it offers a starting point for practical applications and thus is hopefully of great value for the clinically individualizing therapy in the future.
Uncertain pharmacokinetic/pharmacodynamic systems: Design, estimation and control
1997, Control Engineering Practice
In the study of many biological systems, measurement of some process variables occurs only infrequently and at irregular intervals relative to system time constants, while others are completely unobservable. The quantitative study of such sparse data systems (common in the fields of biochemistry, endocrinology, immunology, metabolism, pharmacology, pharmaceutical sciences, toxicology, and other areas), requires modeling methodologies developed expressly to handle the challenges of modeling and data analysis under the constraints of limited data. This presentation reviews current methods for design, estimation and control of sparse data systems, focusing on methods that formally incorporate important sources of uncertainty (both biological and experimental) into the modeling and analysis processes. The methods are illustrated using examples from pharmacokinetics and pharmacodynamics.
Stochastic optimization algorithms of a Bayesian design criterion for Bayesian parameter estimation of nonlinear regression models: Application in pharmacokinetics
1997, Mathematical Biosciences
This article proposes three stochastic algorithms to optimize a Bayesian design criterion for Bayesian estimation of the parameters of nonlinear regression models; this criterion is the information expected from an experiment. The first algorithm is based on a stochastic version of the simplex with an adaptive sampling procedure. The others are stochastic approximation algorithms: the Kiefer-Wolfowitz and the pseudogradient algorithms. We first present the information criterion and the optimization algorithms. The efficiency of each algorithm for optimizing this Bayesian design criterion is then assessed by a simulation study for a nonlinear model assuming a discrete prior distribution. An application for designing an experiment to estimate the kinetics of radioiodine thyroid uptake is then proposed.
Nonparametric EM Algorithms for estimating prior distributions
1991, Applied Mathematics and Computation
Consider a sequence of repetitive experiments E₁, E₂,..., E_N. Each E_i is represented by a pair of random vectors (Y_i, θ_i) such that Y_i is observed, θ_i is not observed, and the conditional density of Y_i given θ_i is known. The θ_i's are independent and identically distributed with common but unknown distribution function F^*(θ). The problem is to estimate F^* given the data Y₁, Y₂,..., Y_N. Two nonparametric maximum likelihood methods based on the EM algorithm are presented. Motivation for this problem comes from the analysis of the distribution of parameters in pharmacokinetic models based on population data.
Incorporating prior parameter uncertainty in the design of sampling schedules for pharmacokinetic parameter estimation experiments
1990, Mathematical Biosciences
An experiment design procedure is proposed for nonlinear parameter estimation studies that formally incorporates prior parameter uncertainty. The design criterion derives from information theory considerations and involves an asymptotic interpretation of the expected posterior information provided by an experiment. A pharmacokinetic sample schedule design problem is used to illustrate and evaluate this information theoretic design strategy. The model considered is commonly used to describe the plasma concentration of a drug following its oral administration. The limitations and advantages of the proposed design procedure are discussed in relation to other previously reported design techniques for incorporating parameter uncertainty.
From Personalized to Precision Medicine in Oncology: A Model-Based Dosing Approach to Optimize Achievement of Imatinib Target Exposure
2023, Pharmaceutics

View all citing articles on Scopus

: This work was supported in part by NIH grant RR01629 made to The Laboratory of Applied Pharmacokinetics, University of Southern California School of Medicine.

View full text

Discrete approximation of multivariate densities with application to Bayesian estimation

Abstract

Computat. Statist. Data Anal.

Math. Biosci.

Bayesian approach to the analysis of nonlinear models: Implementation and evaluation

Biometrics