Abstract
The main purpose in the paper is to discuss algorithm of computing the best prediction in MA(q) time series model. Except in Gaussian process, the best prediction in MA(q) model is nonlinear, generally there is no exact solution, it could compute numeric solution by numeric method. But the exact solution of best linear prediction could be evaluate via characteristic of Gaussian process. Breidt and Hsu(2005) evaluate numeric solution effectively by Monte Carlo method. This paper will use control variate to improve the Breidt and Hsu's method, we hope that we could get better numeric solution in same computation quantity.