Hard-Thresholding Opreator on Covariance Matrix

Share:

Description

Apply hard-thresholding operator on a covariance matrix with a hard-thresholding parameter.

Usage

1
hard.thresholding(sigma, threshold = 0.5)

Arguments

sigma

a p*p covariance matrix

threshold

hard-thresholding parameter

Value

a regularized covariance matrix after hard-thresholding operation

References

"High-Dimensional Covariance Estimation" by Mohsen Pourahmadi

Examples

1
2
3
data(m.excess.c10sp9003)
cov.SAM <- cov(m.excess.c10sp9003)
hard.thresholding(cov.SAM, threshold = 0.001)