Volume limitations and low yield thresholds of biological fluids have led to widespread use of passive microparticle rheology. The mean-squared-displacement (MSD) statistics of bead position time series (bead paths) are either applied directly to determine the creep compliance [Xu et al., Rheol. Acta 37, 387-398 (1998)] or transformed to determine dynamic storage and loss moduli [Mason and Weitz, Phys. Rev. Lett. 74, 1250-1253 (1995)]. A prevalent hurdle arises when there is a nondiffusive experimental drift in the data. Commensurate with the magnitude of drift relative to diffusive mobility, quantified by a Peclet number, the MSD statistics are distorted, and thus the path data must be "corrected" for drift. The standard approach is to estimate and subtract the drift from particle paths, and then calculate MSD statistics. We present an alternative, parametric approach using maximum likelihood estimation that simultaneously fits drift and diffusive model parameters from the path data; the MSD statistics (and consequently the compliance and dynamic moduli) then follow directly from the best-fit model. We illustrate and compare both methods on simulated path data over a range of Peclet numbers, where exact answers are known. We choose fractional Brownian motion as the numerical model, because it affords tunable, subdiffusive MSD statistics consistent with typical 30 s long, experimental observations of microbeads in several biological fluids. Finally, we apply and compare both methods on data from human bronchial epithelial cell culture mucus. (C) 2016 The Society of Rheology.