Given an OTS of length \(T\) with range \(\mathcal{S}=\{s_0, s_1, s_2, \ldots, s_n\}\) (\(s_0 < s_1 < s_2 < \ldots < s_n\)),
\(\overline{X}_t=\{\overline{X}_1,\ldots, \overline{X}_T\}\), the function computes the standard
estimated dispersion given by \(\widehat{disp}_{loc, d}=\frac{1}{T}\sum_{t=1}^Td\big(\overline{X}_t, \widehat{x}_{loc, d}\big)\),
where \(\widehat{x}_{loc, d}\) is the standard estimate of the location and \(d(\cdot, \cdot)\) is a distance between ordinal states.
If normalize = TRUE, then the normalized dispersion is computed, namely
\(\widehat{disp}_{loc, d}/\)max\(_{s_i, s_j \in \mathcal{S}}d(s_i, s_j)\).