lines 1616 to 1621 look redundant to me
why do we need mu3?
code below goes through all the mvsk cases and fills in the missing ones.
(code for moments from stats is somewhere else)