Biostatistics Research Today is a free monthly online journal that collates and summarizes the latest research about Biostatistics, including details on statistics, uncertainty, probability, modeling.

Gamma generalized linear models for pharmacokinetic data.

Salway R, Wakefield J

Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, UK.

This article considers the modeling of single-dose pharmacokinetic data. Traditionally, so-called compartmental models have been used to analyze such data. Unfortunately, the mean function of such models are sums of exponentials for which inference and computation may not be straightforward. We present an alternative to these models based on generalized linear models, for which desirable statistical properties exist, with a logarithmic link and gamma distribution. The latter has a constant coefficient of variation, which is often appropriate for pharmacokinetic data. Inference is convenient from either a likelihood or a Bayesian perspective. We consider models for both single and multiple individuals, the latter via generalized linear mixed models. For single individuals, Bayesian computation may be carried out with recourse to simulation. We describe a rejection algorithm that, unlike Markov chain Monte Carlo, produces independent samples from the posterior and allows straightforward calculation of Bayes factors for model comparison. We also illustrate how prior distributions may be specified in terms of model-free pharmacokinetic parameters of interest. The methods are applied to data from 12 individuals following administration of the antiasthmatic agent theophylline.

Published 16 May 2008 in Biometrics, 64(2): 620-6.
Full-text of this article is available online (may require subscription).

© 2005-2008 Biostatistics Research Today. All Rights Reserved.