In statistics, Whittle likelihood is an approximation to the likelihood function of a stationary Gaussian time series. It is named after the mathematician and statistician Peter Whittle, who introduced it in his PhD thesis in 1951. It is commonly used in time series analysis and signal processing for parameter estimation and signal detection.
Context In a stationary Gaussian time series model, the likelihood function is (as usual in Gaussian models) a function of the associated mean and covariance parameters. With a large number ( N {\displaystyle N} ) of observations, the ( N × N {\displaystyle N\times N} ) covariance matrix may become very large, making computations very costly in practice. However, due to stationarity, the covariance matrix has a rather simple structure, and by using an approximation, computations may be simplified considerably (from O ( N 2 ) {\displaystyle O(N^{2})} to O ( N log ( N ) ) {\displaystyle O(N\log(N))} ). The idea effectively boils down to assuming a heteroscedastic zero-mean Gaussian model in Fourier domain; the model formulation is based on the time series' discrete Fourier transform and its power spectral density.
Definition Let X 1 , … , X N {\displaystyle X_{1},\ldots ,X_{N}} be a stationary Gaussian time series with (one-sided) power spectral density S 1 ( f ) {\displaystyle S_{1}(f)} , where N {\displaystyle N} is even and samples are taken at constant sampling intervals Δ t {\displaystyle \Delta _{t}} . Let X ~ 1 , … , X ~ N / 2 + 1 {\displaystyle {\tilde {X}}_{1},\ldots ,{\tilde {X}}_{N/2+1}} be the (complex-valued) discrete Fourier transform (DFT) of the time series. Then for the Whittle likelihood one effectively assumes independent zero-mean Gaussian distributions for all X ~ j {\displaystyle {\tilde {X}}_{j}} with variances for the real and imaginary parts given by
Var ( Re ( X ~ j ) ) = Var ( Im ( X ~ j ) ) = S 1 ( f j ) {\displaystyle \operatorname {Var} \left(\operatorname {Re} ({\tilde {X}}_{j})\right)=\operatorname {Var} \left(\operatorname {Im} ({\tilde {X}}_{j})\right)=S_{1}(f_{j})}
where f j = j N Δ t {\displaystyle f_{j}={\frac {j}{N\,\Delta _{t}}}} is the j {\displaystyle j} th Fourier frequency. This approximate model immediately leads to the (logarithmic) likelihood function
… excerpt ends here. Continue reading the full article.
