In population pharmacokinetic studies, one of the main objectives is to estimate population pharmacokinetic parameters specifying the population distributions of pharmacokinetic parameters. Confidence intervals for population pharmacokinetic parameters are generally estimated by assuming the asymptotic normality, which is a large-sample property, that is, a property which holds for the cases where sample sizes are large enough. In actual clinical trials, however, sample sizes are limited and not so large in general. Likelihood functions in population pharmacokinetic modelling include a multiple integral and are quite complicated. We hence suspect that the sample sizes of actual trials are often not large enough for assuming the asymptotic normality and that the asymptotic confidence intervals underestimate the uncertainties of the estimates of population pharmacokinetic parameters. As an alternative to the asymptotic normality approach, we can employ a bootstrap approach. This paper proposes a bootstrap standard error approach for constructing confidence intervals for population pharmacokinetic parameters. Comparisons between the asymptotic and bootstrap confidence intervals are made through applications to a simulated data set and an actual phase I trial.