Unbiased Estimation of Population Moments;
Moments of Moments

mathStatica can find unbiased estimators of population moments. For instance, it offers h-statistics (unbiased estimators of population central moments), k-statistics (unbiased estimators of population cumulants), multivariate varieties of the same, polykays (unbiased estimators of products of cumulants) and more. Consider the k-statistic [Graphics:Images/index_gr_1.gif] which is an unbiased estimator of the [Graphics:Images/index_gr_2.gif] cumulant [Graphics:Images/index_gr_3.gif]; that is, [Graphics:Images/index_gr_4.gif], for [Graphics:Images/index_gr_5.gif] . Here are the [Graphics:Images/index_gr_6.gif] are [Graphics:Images/index_gr_7.gif] k-statistics:

[Graphics:Images/index_gr_8.gif]
[Graphics:Images/index_gr_9.gif]
[Graphics:Images/index_gr_10.gif]

As per convention, the solution is expressed in terms of power sums [Graphics:Images/index_gr_11.gif].

Moments of moments: Because the above expressions (sample moments) are functions of random variables [Graphics:Images/index_gr_12.gif], we might want to calculate population moments of them. With mathStatica, we can find any moment (raw, central, or cumulant) of the above expressions. For instance, [Graphics:Images/index_gr_13.gif] is meant to have the property that [Graphics:Images/index_gr_14.gif]. We test this by calculating the first raw moment of [Graphics:Images/index_gr_15.gif], and express the answer in terms of cumulants:

[Graphics:Images/index_gr_16.gif]
[Graphics:Images/index_gr_17.gif]

In 1928, Fisher published the product cumulants of the k-statistics, which are now listed in reference bibles such as Stuart and Ord (1994). Here is the solution to [Graphics:Images/index_gr_18.gif]:

[Graphics:Images/index_gr_19.gif]
[Graphics:Images/index_gr_20.gif]

This is the correct solution. Unfortunately, the solutions given in Stuart and Ord (1994, equation (12.70)) and Fisher (1928) are actually incorrect.