cum3(a, b=a, c=a, unbiased=TRUE)a.
If supplied then it must be the same length as a.a.
If supplied then it must be the same length as a.n/((n-1)*(n-2)) where n is
the sample size, if unbiased is FALSE then a multiplier of 1/n is used.
This is multiplied by sum((a-mean(a))*(b-mean(b))*(c-mean(c))) to give the
required estimate.