In forensic science, DNA mixture interpretation is traditionally based on a binary model, which does not account for peak-height information in DNA profiles. In recent years, some countries have adopted a continuous model in which peak heights are used and stochastic effects are considered to enable rigorous calculation of likelihood ratios. However, this model requires certain biological parameters which affect the expected allelic and stutter peak heights. In this paper, we focused on estimating the distribution of the stutter ratio (SR) in 15 short tandem repeat loci in relation to the allele repeat number. We estimated the SR values of 234 single-source DNA samples by using a commercially available kit. In all loci except for D8S1179, D21S11, and D2S1338, a simple log-normal distribution model was fitted to the variability of SR. For D21S11, we developed a new distribution model in which distinct log-normal distributions between complete and incomplete repeat units are used (a separate log-normal distribution model). For D8S1179 and D2S1338, we developed another new distribution model that mixes two log-normal distributions to explain two types of repeat structures appearing within the same number of allele repeats. These two models were fitted to the observed SR values more accurately than the simple log-normal distribution model. We expected these new SR models to be applied to DNA mixture interpretation based on a continuous model.