Fitting heteroscedastic regression models to individual pharmacokinetic data using standard statistical software

In the analysis of individual pharmacokinetic data by nonlinear regression it is important to allow for possible heterogeneity of variance in the response. Two common methods of doing this are weighted least squares with appropriate weights or data transformation using a suitable transform. With either approach it is appealing to let the data determine the appropriate choice of weighting scheme or transformation. This article describes two methods of doing this which are easy to compute using standard statistical software. The first method is a generalized least squares scheme for the case where the variance is assumed proportional to an unknown power of the mean. The second involves applying a power transformation to both sides of the regression equation. It is shown that both techniques may be implemented using only nonlinear regression routines. Sample code is provided for their implementation using the SAS software package. However, the proposed methods are feasible using any software package that incorporates a nonlinear least squares routine, and are thus well suited to routine use.

[1]  D. Ruppert,et al.  A comparison between maximum likelihood and generalized least squares in a heteroscedastic linear model , 1982 .

[2]  D. Cox,et al.  An Analysis of Transformations , 1964 .

[3]  D. J. Finney,et al.  The Form and Estimation of a Variance Function, with Particular Reference to Radioimmunoassay , 1977 .

[4]  Noel A Cressie,et al.  A Transformation/Weighting Model for Estimating Michaelis-Menten Parameters, , 1989 .

[5]  D. Ruppert,et al.  Diagnostics and robust estimation when tranforming the regression model and the response , 1987 .

[6]  M. Lewis,et al.  Pharmacokinetic profile of sulfisoxazole following intravenous, intramuscular, and oral administration to man. , 1972, Journal of pharmaceutical sciences.

[7]  J C van Houwelingen.  Use and abuse of variance models in regression. , 1988, Biometrics.

[8]  J. C. van Houwelingen,et al.  Use and abuse of variance models in regression. , 1988 .

[9]  D. W. Bacon,et al.  Further Consideration of Heteroscedasticity in Fitting Kinetic Models , 1977 .

[10]  Lewis B. Sheiner,et al.  Pharmacokinetic parameter estimates from several least squares procedures: Superiority of extended least squares , 1985, Journal of Pharmacokinetics and Biopharmaceutics.

[11]  D. Ruppert,et al.  Transformation and Weighting in Regression , 1988 .

[12]  Lewis B. Sheiner,et al.  Heteroscedastic nonlinear regression , 1988 .

[13]  Difficulties in obtaining consistent estimators of variance parameters , 1981 .

[14]  David Ruppert,et al.  Power Transformations When Fitting Theoretical Models to Data , 1984 .

[15]  J. S. Hunter,et al.  Statistics for Experimenters: An Introduction to Design, Data Analysis, and Model Building. , 1979 .

[16]  George E. P. Box,et al.  Correcting Inhomogeneity of Variance with Power Transformation Weighting , 1974 .

[17]  Lewis B. Sheiner,et al.  Extended least squares nonlinear regression: A possible solution to the “choice of weights” problem in analysis of individual pharmacokinetic data , 1984, Journal of Pharmacokinetics and Biopharmaceutics.