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.