Wiener Kolmogorov Estimators
Arguments
- ucm
An UCARIMA model returned by
ucarima_model().- cmp
Index of the component for which we want to compute the filter
- signal
TRUE for the signal (component), FALSE for the noise (complement)
- nspectrum
Number of points used to compute the (pseudo-) spectrum of the estimator
- nwk
Number of weights of the wiener-kolmogorov filter returned in the result
Value
A list with the (pseudo-)spectrum, the weights of the filter and the squared-gain function (with the same number of points as the spectrum)
Examples
mod1 <- arima_model("trend", delta = c(1,-2,1))
mod2 <- arima_model("noise", var = 1600)
hp<-ucarima_model(components=list(mod1, mod2))
wk1<-ucarima_wk(hp, 1, nwk=50)
wk2<-ucarima_wk(hp, 2)
plot(wk1$filter, type='h')
