Unbiased estimates for moments and cumulants in linear regression

We consider the linear regression model Y = X β + e with independent residuals having zero mean and unknown distribution F. A method is given for obtaining unbiased estimates of the moments and cumulants of F and β ^ , the least squares estimate. Explicit formulas are given for those of order less than six.