Journal List > Korean J Physiol Pharmacol > v.13(5) > 1025628

Ahn and Yim: Comparison of Parametric and Bootstrap Method in Bioequivalence Test

Abstract

The estimation of 90% parametric confidence intervals (CIs) of mean AUC and Cmax ratios in bioequivalence (BE) tests are based upon the assumption that formulation effects in log-transformed data are normally distributed. To compare the parametric CIs with those obtained from nonparametric methods we performed repeated estimation of bootstrap-resampled datasets. The AUC and Cmax values from 3 archived datasets were used. BE tests on 1,000 resampled datasets from each archived dataset were performed using SAS (Enterprise Guide Ver.3). Bootstrap nonparametric 90% CIs of formulation effects were then compared with the parametric 90% CIs of the original datasets. The 90% CIs of formulation effects estimated from the 3 archived datasets were slightly different from nonparametric 90% CIs obtained from BE tests on resampled datasets. Histograms and density curves of formulation effects obtained from resampled datasets were similar to those of normal distribution. However, in 2 of 3 resampled log (AUC) datasets, the estimates of formulation effects did not follow the Gaussian distribution. Bias-corrected and accelerated (BCa) CIs, one of the nonparametric CIs of formulation effects, shifted outside the parametric 90% CIs of the archived datasets in these 2 non-normally distributed resampled log (AUC) datasets. Currently, the 80∼125% rule based upon the parametric 90% CIs is widely accepted under the assumption of normally distributed formulation effects in log-transformed data. However, nonparametric CIs may be a better choice when data do not follow this assumption.

References

Bonate PL. Pharmacokinetic-pharmacodynamic modeling and simulation. 1st ed.Springer;New York: p. p. 355–363. 2005.
Chow SC. Encyclopedia of biopharmaceutical statistics. 2nd ed.Marcel Dekker;New York: p. p. 83–88. 2003.
Efron B., Tibshirani R. An introduction to the bootstrap. Monographs on statistics and applied probability. 1st ed.Chapman & Hall;New York: p. p. 179–201. 372–391,. 1993.
EMEA. Questions & Answers on the Bioavailability and Bioequivalence Guideline (EMEA/CHMP/EWP/40326/2006). [EMEA Web site],. http://www.emea.europa.eu/pdfs/human/ewp/4032606en.pdf. European Medicines Agency, Evaluation of Medicines for Human Use. Accessed August 10,. 2009.
FDA. Guidance for industry: Statistical Approaches to Establishing Bioequivalence January 2001. [CDER Web site],. http://www.fda.gov/downloads/Drugs/GuidanceComplianceRegulatoryInformation/Guidances/ucm070244.pdf. US Food and Drug Administration, Center for Drug Evaluation and Research. Accessed August 10,. 2009.
Henderson AR. The bootstrap: a technique for data-driven statistics. Using computer-intensive analyses to explore experimental data. Clin Chim Acta. 359:1–26. 2005.
crossref
Lacey LF., Keene ON., Pritchard JF., Bye A. Common noncompartmental pharmacokinetic variables: are they normally or log-normally distributed? J Biopharm Stat. 7:171–178. 1997.
crossref
Midha KK., Ormsby ED., Hubbard JW., McKay G., Hawes EM., Gavalas L., McGilveray IJ. Logarithmic transformation in bioequivalence: application with two formulations of perphenazine. J Pharm Sci. 82:138–144. 1993.
crossref
Pabst G., Jaeger H. Review of methods and criteria for the evaluation of bioequivalence studies. Eur J Clin Pharmacol. 38:5–10. 1990.
crossref
Patterson SD., Jones B. Bioequivalence and the pharmaceutical industry. Pharmaceut Statist. 1:83–95. 2002.
crossref
Sprent P., Smeeton NC. Applied nonparametric statistical methods. 3rd ed.Chapman & Hall/CRC;Boca Raton: p. p. 404–437. 2001.
Steyn HS., Koeleman HA., Gouws E., Ritschel WA. An approach to select the appropriate statistical method for testing bioequivalence. Int J Clin Pharmacol Ther Toxicol. 29:156–160. 1991.

Fig. 1.
Histograms and normal Q-Q plots with skewness, kurtosis and Shapiro-Wilk test results showing the distribution of formulation effects (mean differences in log (AUC)s between test and reference formulation) estimates in bootstrap-resampled datasets.
kjpp-13-367f1.tif
Fig. 2.
Histograms and normal Q-Q plots with skewness, kurtosis and Shapiro-Wilk test results showing the distribution of formulation effects (mean differences in log (Cmax) between test and reference formulation) estimates in bootstrap-resampled datasets.
kjpp-13-367f2.tif
Fig. 3.
Comparison of BCa 90% CIs with parametric CIs (solid lines, parametric 5% and 95% points; dotted lines, BCa 5th and 95th percentile points; density curve, anticipated parametric distribution of formulation effects from the archived datasets). BE1_AUC and BE2_AUC were found to have non-normally distributed formulation effects by Shapro-Wilk tests.
kjpp-13-367f3.tif
Table 1.
The 90% confidence intervals (CI) of formulation effects in the archived and bootstrap-resampled datasets and percent coverages of nonparametric CIs in contrast with parametric CIs
    Log (AUC) Log (Cmax)
Parametric Nonparametric Parametric Nonparametric
Dataset Parameter Archived-t Percentile Bootstrap-t BC BCa Archived-t Percentile Bootstrap-t BC BCa
BE1 5% -0.059 -0.056 -0.064 -0.059 -0.06 -0.083 -0.077 -0.097 -0.07 -0.071
  95% 0.033 0.028 0.031 0.026 0.027 0.162 0.157 0.171 0.161 0.162
  90% Interval 0.092 0.084 0.095 0.085 0.086 0.245 0.234 0.268 0.23 0.233
  % coverage 90 82.297 92.828 82.687 84.15 90 86.031 98.379 84.598 85.59
BE2 5% 0.061 0.07 0.071 0.075 0.077 0.259 0.267 0.256 0.272 0.272
  95% 0.168 0.168 0.205 0.18 0.173 0.43 0.422 0.434 0.426 0.424
  90% Interval 0.107 0.098 0.134 0.105 0.096 0.171 0.155 0.178 0.154 0.152
  % coverage 90 82.444 112.332 88.489 80.177 90 81.484 93.522 80.748 79.854
BE3 5% -0.202 -0.191 -0.203 -0.197 -0.197 -0.203 -0.19 -0.197 -0.194 -0.19
  95% -0.079 -0.079 -0.083 -0.086 -0.086 -0.005 -0.013 -0.005 -0.014 -0.016
  90% Interval 0.123 0.112 0.12 0.111 0.111 0.198 0.177 0.192 0.18 0.175
  % coverage 90 81.686 87.374 81.248 81.175 90 80.292 87.142 81.653 79.158

Archived-t, parameter estimation from t distribution using the REML method in the archived datasets; BC, bias-corrected CI; BCa, bias correctedand accelerated CI; 90% Interval, the interval between 5% (or percentile) and 95% (or percentile); % (percent) coverage, the length of nonparametric 90% CIs measured in relation to the parametric 90% CI. (i.e., when the 2 lengths are the same, the % coverage is 90) percent of the CI which bootstrap 90% CI covers in contrast to the parametric 90% CI; %, percent for the archived datasets and percentile for the bootstrap-resampled datasets.

TOOLS
Similar articles