The objective of this project is to compare the results of the same study carried out on NONMEM and nlmixr2. This analysis consists of evaluating previously published population pharmacokinetic models of gentamicin and tobramycin in our population of interest with sparse concentrations. A literature review was performed to determine the gentamicin and tobramycin models in critically ill adult patients. In parallel, gentamicin and tobramycin dosing data, information on the treatment, the patient, and the bacteria were collected retrospectively in 2 Quebec establishments. The external evaluations were previously performed using NONMEM Version 7.5. Model equations were rewritten with R, and external evaluations were performed using nlmixr2. Predictive performance was assessed based on the estimation of bias and imprecision of the prediction error for maximum a posteriori (MAP) Bayesian PK parameter and observed concentrations. Comparison between nlmixr2 and NONMEM was performed on 4 gentamicin and 3 tobramycin population pharmacokinetic models. Compared to NONMEM, for gentamicin and tobramycin clearance and central volume of distribution, nlmixr2 produced individual pharmacokinetic parameters with bias values ranging from -32.5% to 5.67% and imprecision values ranging from 6.33% to 32.5%. Despite these differences, population bias and imprecision for sparse concentrations were low and ranged from 0% to 5.3% and 0.2% to 6.5%, respectively. The external evaluations performed with both software packages resulted in the same interpretation in terms of population predictive performance for all 7 models. Nlmxir2 showed comparable predictive performance with NONMEM with sparse concentrations that are, at most, sampled twice within a single dose administration (peak and trough).