There seems to be some missing units in the example calculation of the kurtosis coefficient, since the expected value of v can contain units, but the standard deviation function in the denominator strips the units out, so presumably, we need to apply UnitsOf(v) in the denominator:

Note also, that the calculation of the skewness coefficient has units, but shouldn't.
TTFN,
Eden