Properties of two of the stochastic circulatory models theoretically introduced by Smith et al., 1997, Bull. Math. Biol. 59, 1-22 were investigated. The models assumed the gamma distribution of the cycle time under either the geometric or Poisson elimination scheme. The reason for selecting these models was the fact that the probability density functions of the residence time of these models are formally similar to those of the Bateman and gamma-like function models, i.e., the two common deterministic models. Using published data, the analytical forms of the probability density functions of the residence time and the distributions of the simulated values of the residence time were determined on the basis of the deterministic models and the stochastic circulatory models, respectively. The Kolmogorov-Smirnov test revealed that even for 1000 xenobiotic particles, i.e., a relatively small number if the particles imply drug molecules, the probability density functions of the residence time based on the deterministic models closely matched the distributions of the simulated values of the residence time obtained on the basis of the stochastic circulatory models, provided that parameters of the latter models fulfilled selected conditions.