Description Usage Arguments Details Value References See Also Examples

`pdSpecEst2D`

calculates a *(d,d)*-dimensional HPD wavelet-denoised time-varying spectral matrix estimator
by applying the following steps to an initial noisy HPD time-varying spectral estimate (obtained with e.g., `pdPgram2D`

):

a forward intrinsic AI wavelet transform, with

`WavTransf2D`

,(tree-structured) thresholding of the wavelet coefficients, with

`pdCART`

,an inverse intrinsic AI wavelet transform, with

`InvWavTransf2D`

.

The complete estimation procedure is described in more detail in Chapter 5 of \insertCiteC18pdSpecEst.

1 2 | ```
pdSpecEst2D(P, order = c(3, 3), metric = "Riemannian", alpha = 1,
return_val = "f", ...)
``` |

`P` |
a ( |

`order` |
a 2-dimensional numeric vector |

`metric` |
the metric that the space of HPD matrices is equipped with. The default choice is |

`alpha` |
an optional tuning parameter in the wavelet thresholding procedure. The penalty (or sparsity)
parameter in the tree-structured wavelet thresholding procedure in |

`return_val` |
an optional argument that specifies whether the denoised spectral estimator is returned or not. See the Details section below. |

`...` |
additional arguments for internal use. |

The input array `P`

corresponds to an initial noisy HPD time-varying spectral estimate of the (*d, d*)-dimensional
spectral matrix at a time-frequency grid of size *m_1 \times m_2*, with *m_1, m_2* dyadic numbers. This can be e.g.,
a multitaper HPD time-varying periodogram given as output by the function `pdPgram2D`

.

`P`

is transformed to the wavelet domain by the function `WavTransf2D`

, which applies an intrinsic
2D AI wavelet transform based on a metric specified by the user. The noise is removed by tree-structured
thresholding of the wavelet coefficients based on the trace of the whitened coefficients with `pdCART`

by
minimization of a *complexity penalized residual sum of squares* (CPRESS) criterion via the fast tree-pruning algorithm
in \insertCiteD97pdSpecEst. The penalty (i.e., sparsity) parameter in the optimization procedure is set equal to `alpha`

times the universal threshold, where the noise variance of the traces of the whitened wavelet
coefficients are determined from the finest wavelet scale. See Chapter 5 of \insertCiteC18pdSpecEst
for further details.

The function computes the forward and inverse intrinsic 2D AI wavelet transform in the space of HPD matrices equipped with
one of the following metrics: (i) the affine-invariant Riemannian metric (default) as detailed in e.g., \insertCiteB09pdSpecEst[Chapter 6]
or \insertCitePFA05pdSpecEst; (ii) the log-Euclidean metric, the Euclidean inner product between matrix logarithms;
(iii) the Cholesky metric, the Euclidean inner product between Cholesky decompositions; (iv) the Euclidean metric; or
(v) the root-Euclidean metric. The default choice of metric (affine-invariant Riemannian) satisfies several useful properties
not shared by the other metrics, see \insertCiteCvS17pdSpecEst or \insertCiteC18pdSpecEst for more details. Note that this comes
at the cost of increased computation time in comparison to one of the other metrics.

If `return_val = 'f'`

the thresholded wavelet coefficients are transformed back to the time-frequency domain by
the inverse intrinsic 2D AI wavelet transform via `InvWavTransf2D`

, returning the wavelet-denoised
HPD time-varying spectral estimate.

The function returns a list with the following five components:

`f ` |
a ( |

`D ` |
the 2D pyramid of threshold wavelet coefficients. This is a list of arrays, where each array contains the rectangular grid
( |

`M0 ` |
a numeric array containing the midpoint(s) at the coarsest scale |

`tree.weights ` |
a list of logical values specifying which coefficients to keep, with each list component
corresponding to an individual wavelet scale starting from the coarsest wavelet scale |

`D.raw ` |
the 2D pyramid of non-thresholded wavelet coefficients in the same format as the component |

`pdPgram2D`

, `WavTransf2D`

, `InvWavTransf2D`

, `pdCART`

1 2 3 4 5 | ```
## Not run:
P <- rExamples2D(c(2^6, 2^6), 2, example = "tvar")$P
f <- pdSpecEst2D(P)
## End(Not run)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.