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.066298 0.0000000 0.000000
#> [2,] 0.000000 0.9823207 0.000000
#> [3,] 0.000000 0.0000000 1.075064