To correctly evaluate the glucose control system, it is crucial to account for both insulin sensitivity and secretion. The disposition index (DI) is the most widely accepted method to do so. The original paradigm (hyperbolic law) consists of the multiplicative product of indices related to insulin sensitivity and secretion, but more recently, an alternative formula has been proposed with the exponent α (power function law). Traditionally, curve-fitting approaches have been used to evaluate the DI in a population: the algorithmic implementations often introduce some critical issues, such as the assumption that one of the two indices is error free or the effects of the log transformation on the measurement errors. In this work, we review the commonly used approaches and show that they provide biased estimates. Then we propose a novel nonlinear total least square (NLTLS) approach, which does not need to use the approximations built in the previously proposed alternatives, and show its superiority. All of the traditional fit procedures, including NLTLS, account only for uncertainty affecting insulin sensitivity and secretion indices when they are estimated from noisy data. Thus, they fail when part of the observed variability is due to inherent differences in DI values between individuals. To handle this inevitable source of variability, we propose a nonlinear mixed-effects approach that describes the DI using population hyperparameters such as the population typical values and covariance matrix. On simulated data, this novel technique is much more reliable than the curve-fitting approaches, and it proves robust even when no or small population variability is present in the DI values. Applying this new approach to the analysis of real IVGTT data suggests a value of α significantly smaller than 1, supporting the importance of testing the power function law as an alternative to the simpler hyperbolic law.