SUMMARY For a signal plus noise model, where both signal and noise are generated by nonstation- ary ARIMA, autoregressive integrated moving average, models, we use the transformation approach of Ansley & Kohn (1985) to construct an estimate of the signal and obtain its mean squared error given a finite number of observations. A method for efficient computa- tion of the signal estimate and its mean squared error is also presented. Consider observations z(t) generated by the signal plus noise nmodel z(t) = s(t)+ n(t), (1s1) where both signal s(t) and noise n(t) follow nonstationary ARIMA processes. We give a solution to the finite sample signal extraction problem by defining an estimate of the signal s(t) given observations z(1),... , z(N) and finding its mean squared error. The major difficulty in handling the model with nonstationary components is that the distribution of the observations is not defined unless assumptions are made about the initial values of the series. Moreover, it is usually very difficult to formulate reasonable assumptions about the initial values because they affect the correlation structure of the whole series. See, for example, the discussion of Bell (1984, p. 652). Kohn & Ansley (1986) use a transformation approach to overcome the initial value problem in defining the likelihood and obtaining predittors and interpolators for an ARIMA model which may have missing values. This approach does not require any assumptions about the initial values and generalizes the usual differencing approach of Box & Jenkins (1976) for the case where there are no missing observations. In the present paper we apply the transformation approach of Kohn & Ansley (1986) and Ansley & Kohn (1985) to solve the signal extraction problem. We define a linear estimator A( t) of s( t) so that s( t) - SA( t) is invariant to the initial values and, furthermore, is unbiased in the sense that E{s(t) - s(t)} = 0. Moreover, it has minimum mean squared error amongst all such estimators. Our approach allows missing observations and extends to signal extraction for several additive signals when each is an ARIMA process. We also show how to compute the signal estimate and its mean squared error efficiently by writing the model in state space form, placing a diffuse prior on the initial values of the process and then using the modified Kalman filter and the modified smoothing algorithm of Ansley & Kohn (1985) and Kohn & Ansley (1986). Gersch & Kitagawa (1983), Kitagawa & Gersch (1984) and Harvey & Todd (1983) use
Read full abstract