We consider the problem of reconstructing the cross-power spectrum of an unobservable multivariate stochastic process from indirect measurements of a second multivariate stochastic process, related to the first one through a linear operator. In the two-step approach, one would first compute a regularized reconstruction of the unobservable signal, and then compute an estimate of its cross-power spectrum from the regularized solution. We investigate whether the optimal regularization parameter for reconstruction of the signal also gives the best estimate of the cross-power spectrum. We show that the answer depends on the regularization method, and specifically we prove that, under a white Gaussian assumption: (i) when regularizing with truncated SVD the optimal parameter is the same; (ii) when regularizing with the Tikhonov method, the optimal parameter for the cross-power spectrum is lower than half the optimal parameter for the signal. We also provide evidence that a one-step approach would likely have better mathematical properties than the two-step approach. Our results apply particularly to the brain connectivity estimation from magneto/electro-encephalographic recordings and provide a formal interpretation of recent empirical results.

On the two-step estimation of the cross-power spectrum for dynamical linear inverse problems

Vallarino E.;Sommariva S.;Piana M.;Sorrentino A.
2020-01-01

Abstract

We consider the problem of reconstructing the cross-power spectrum of an unobservable multivariate stochastic process from indirect measurements of a second multivariate stochastic process, related to the first one through a linear operator. In the two-step approach, one would first compute a regularized reconstruction of the unobservable signal, and then compute an estimate of its cross-power spectrum from the regularized solution. We investigate whether the optimal regularization parameter for reconstruction of the signal also gives the best estimate of the cross-power spectrum. We show that the answer depends on the regularization method, and specifically we prove that, under a white Gaussian assumption: (i) when regularizing with truncated SVD the optimal parameter is the same; (ii) when regularizing with the Tikhonov method, the optimal parameter for the cross-power spectrum is lower than half the optimal parameter for the signal. We also provide evidence that a one-step approach would likely have better mathematical properties than the two-step approach. Our results apply particularly to the brain connectivity estimation from magneto/electro-encephalographic recordings and provide a formal interpretation of recent empirical results.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11567/1005060
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 5
  • ???jsp.display-item.citation.isi??? 4
social impact