Logo: to the web site of Uppsala University

The smorgasbord of methods that we use within the field of pharmacometrics has developed steadily over several decades and is now a well-laid-out buffet. This thesis adds some garnish to the table in the form of small improvements to the handling of certain problems.

The first problem tackled by the thesis was the challenge of saddle points and local non-identifiability when estimating pharmacometric model parameters. Substituting the common method of randomly perturbing the initial parameter estimates with one saddle-reset step enhances the accuracy of maximum likelihood estimates by overcoming saddle points parameter values, a common issue in nonlinear mixed-effects models. This algorithm, as implemented in the NONMEM software, was applied to various identifiable and nonidentifiable pharmacometric models, showing improved performance over traditional methods.

Part of the thesis was dedicated to the development of a paediatric pharmacokinetic model for ethionamide, a drug used in treating multidrug-resistant tuberculosis. The resulting model was then used to simulate drug exposure under different dosing regimens, a new dosing regimen for children was proposed. The developed model, and therefore the proposed paediatric dosing regimen, considers factors like maturation of pharmacokinetic pathways and, administration by nasogastric tube, and concurrent rifampicin treatment. The regimen, with some modifications, was adopted in the 2022 update to the World Health Organization operational handbook on tuberculosis.

Finally, the thesis explored novel model-integrated evidence (MIE) approaches for bioequivalence (BE) determination. Such methods could offer more robust alternatives to standard BE approached using non-compartmental analysis (NCA). Model-based methods have been shown to be advantageous in sparse data situations, such as is found in studies of ophthalmic formulations, but have suffered from inflated type I error rates. MIE BE approaches using a single model or using model averaging were presented and shown to control type I error at the nominal level while demonstrating increased power in bioequivalence determination.