Compute diagonal covariance matrix

Returns a diagonal covariance matrix from the column-wise variances of Y.

Usage

computeCovDiag(Y)

Arguments

Y

Numeric matrix (samples x conditions).

Value

A diagonal covariance matrix of dimension ncol(Y) x ncol(Y).

Examples

data(eqtlRegionExample)
X <- eqtlRegionExample$X[, 1:30]
Y <- matrix(rnorm(nrow(X) * 3), nrow(X), 3)
computeCovDiag(Y = Y)
#>         [,1]     [,2]     [,3]
#> [1,] 1.02831 0.000000 0.000000
#> [2,] 0.00000 1.063156 0.000000
#> [3,] 0.00000 0.000000 0.863185