The pharmacokinetic behavior of furosemide concentrations in performance horses is of great interest to the equine industry and equine researchers. Specifically, such information is useful for the evaluation of the existing guidelines in several racing jurisdictions and for the possible development of new guidelines for varying time periods after administration. We studied several approaches within the framework of nonlinear mixed effects models to increase the accuracy of evaluating these guidelines. Theoretical properties of the proposed methods were examined and the variances of the resulting estimators compared. Their numerical performances were further evaluated through simulations. Finally, we applied these methods to a furosemide concentration profile data set and used our findings to address certain important practical issues.