Bias correction in ARMA models

We give a general matrix formula for computing the bias of the exact unconditional maximum likelihood estimate in ARMA models, with known and unknown mean, up to order 1/n, where n is the length of the series. Some illustrative examples are presented.