Identification of pharmacokinetic models via population smoothing splines