The Cluster Newton method (CN method) has proved to be very efficient at finding multiple solutions to underdetermined inverse problems. In the case of pharmacokinetics, underdetermined inverse problems are often given extra constraints to restrain the variety of solutions. In this paper, we propose an algorithm based on the two parameters of the Beta distribution to find families of solution near a solution of interest. This allows for a much greater control of the variety of solutions that can be obtained with the CN method. In addition, this algorithm facilitates the task of obtaining pharmacologically feasible parameters. Moreover, we also make some improvements to the original CN method including an adaptive margin of error for the perturbation of the target values and the use of an analytical Jacobian in the resolution of the forward problem.
Gaudreau et al. (Wed,) studied this question.