Abstract. Discrete Phase-Type (DPH) distributions have one property that is not shared by Continuous Phase-Type (CPH) distributions, i.e., representing a deterministic value as a DPH random variable. This property distinguishes the application of DPH in stochastic modeling of real-life problems, such as stochastic scheduling, in which service time random variables should be compared with a deadline that is usually a constant value. In this paper, we consider a restricted class of DPH distributions, called Mixed Shifted Negative Binomial (MSNB), and show its exibility in producing a wide range of variances as well as its adequacy in tting fat-tailed distributions. These properties render MSNB applicable to represent data on certain types of service time. Therefore, we adapt an ExpectationMaximization (EM) algorithm to estimate the parameters of MSNB distributions that accurately t trace data. To present the applicability of the proposed algorithm, we use it to t real operating room times and a set of benchmark traces generated from continuous distributions as case studies. Finally, we illustrate the e ciency of the proposed algorithm by comparing its results with those of two existing algorithms in the literature. We conclude that our proposed algorithm outperforms other DPH algorithms in tting trace data and distributions.