spikedOperatorShrinkEst: Operator Norm Shrinkage Estimator, Spiked Covariance Model

In cvCovEst: Cross-Validated Covariance Matrix Estimation

Operator Norm Shrinkage Estimator, Spiked Covariance Model

Description

spikedOperatorShrinkEst() implements the asymptotically optimal shrinkage estimator with respect to the operator loss in a spiked covariance matrix model. Informally, this model admits Gaussian data-generating processes whose covariance matrix is a scalar multiple of the identity, save for a few number of large "spikes". A thorough review of this estimator, or more generally spiked covariance matrix estimation, is provided in \insertCitedonoho2018;textualcvCovEst.

Usage

spikedOperatorShrinkEst(dat, p_n_ratio, num_spikes = NULL, noise = NULL)

Arguments

dat	A numeric data.frame, matrix, or similar object.
p_n_ratio	A numeric between 0 and 1 representing the asymptotic ratio of the number of features, p, and the number of observations, n.
num_spikes	A numeric integer equal to or larger than one which providing the known number of spikes in the population covariance matrix. Defaults to NULL, indicating that this value is not known and must be estimated.
noise	A numeric representing the known scalar multiple of the identity matrix giving the approximate population covariance matrix. Defaults to NULL, indicating that this values is not known and must be estimated.

Value

A matrix corresponding to the covariance matrix estimate.

References

\insertAllCited

Examples

spikedOperatorShrinkEst(dat = mtcars, p_n_ratio = 0.1, num_spikes = 2L)

cvCovEst documentation built on May 29, 2024, 5:51 a.m.