# band.chol: Computes the banded covariance estimator by banding the... In PDSCE: Positive Definite Sparse Covariance Estimators

 band.chol R Documentation

## Computes the banded covariance estimator by banding the covariance Cholesky factor

### Description

Computes the k-banded covariance estimator by k-banding the covariance Cholesky factor as described by Rothman, Levina, and Zhu (2010).

### Usage

band.chol(x, k, centered = FALSE, method = c("fast", "safe"))


### Arguments

 x A data matrix with n rows and p columns. The rows are assumed to be a realization of n independent copies of a p-variate random vector. k The banding parameter (the number of sub-diagonals to keep as non-zero). Should be a non-negative integer. centered Logical: centered=TRUE should be used if the columns of x have already been centered. method The method to use. The default is method="fast", which uses the Grahm-Schmidt style algorithm and must have k ≤q \min(n-2, p-1). Alternatively, method="safe" uses an inverse or generalized inverse to compute estimates of the regression coefficients and is more numerically stable (and capable of handling k \in \{0,…,p-1\} regardless of n).

### Details

method="fast" is much faster than method="safe". See Rothman, Levina, and Zhu (2010).

### Value

The banded covariance estimate (a p by p matrix).

### References

Rothman, A. J., Levina, E., and Zhu, J. (2010). A new approach to Cholesky-based covariance regularization in high dimensions. Biometrika 97(3): 539-550.

band.chol.cv

### Examples

set.seed(1)
n=50
p=20
true.cov=diag(p)
true.cov[cbind(1:(p-1), 2:p)]=0.4
true.cov[cbind(2:p, 1:(p-1))]=0.4
eo=eigen(true.cov, symmetric=TRUE)
z=matrix(rnorm(n*p), nrow=n, ncol=p)
x=z%*% tcrossprod(eo$vec*rep(eo$val^(0.5), each=p),eo\$vec)
sigma=band.chol(x=x, k=1)
sigma


PDSCE documentation built on June 23, 2022, 9:12 a.m.