locmulti {coloredICA}R Documentation

Local polynomial spectral density estimation


This function implements a local polinomial estimation for the log spectral density at a point x0 \in R^2.


locmulti(x0, l_period, n, freq, h)



Point \in R^2 at which the spectral density estimate is evaluated.


Vector of length n with the log-periodogram evaluations at the n Fourier frequencies.


Number of points in the analyzed lattice.


n \times 2 matrix with the n Fourier frequencies.


Kernel bandwidth.


locmulti function is auxiliary for the nonparametric estimation of the sources spectral density step of the scICA function. locmulti function implements the initial estimates for the local maximum likelihood estimator of the log spectral density m(\code{x0}) at a point x0 \in R^2. To obtain an estimate of m(\code{x0}) the local likelihood function

\sum_{k}\left(Y_k - a - \textbf{b}'(\boldsymbol{\omega}_k - x0) - e^{Y_k - a - \textbf{b}'(\boldsymbol{\omega}_k - x0)} \right)K_H(\boldsymbol{\omega}_k - x0)

is constructed, where Y_k denotes the log-periodogram value at the Fourier frequency \boldsymbol{\omega}_k, K_H a surface kernel and H=(h,h). The local maximum estimator \widehat{m}_{LK}(\code{x0}) is \widehat{a} in the maximizer (\widehat{a},\widehat{\textbf{b}}). The estimate is implemented directly in the scICA function through a Newton-Rapshon algorithm. The initialization for the Newton-Rapshon algorithm is derived through a local polynomial approximation implemented in this locmulti function. In particular the following function is minimized to find a local polynomial approximation for m(\code{x0})

\sum_{k}\left(Y_k - a - \textbf{b}'(\boldsymbol{\omega}_k - x0) \right)^{2}K_H(\boldsymbol{\omega}_k - x0)

and the minimizer \widehat{a} is used as an initial value in order to obtain the local maximum likelihood estimator \widehat{m}_{LK}(\code{x0}).


It returns a list containing the following component:


local polynomial estimate of the log spectral density at x0.


It is auxiliary for scICA function.


Lee, S., Shen, H., Truong, Y. and Zanini, P.


Shen, H., Truong, Y., Zanini, P. (2014). Independent Component Analysis for Spatial Processes on a Lattice. MOX report 38/2014, Department of Mathematics, Politecnico di Milano.

Fan, J., Kreutzberger, E. (1998). Automatic Local Smoothing for Spectral Density Estimation. Scandinavian Journal of Statistics, 25, 359–369.

See Also

scICA, kern

[Package coloredICA version 1.0.0 Index]