Learn R Programming

SplitKnockoff (version 2.1)

sk.decompose: make SVD as well as orthogonal complements

Description

make SVD as well as orthogonal complements

Usage

sk.decompose(A, D)

Value

U

S

V

U_perp: orthogonal complement for U

Arguments

A

the input matrix

D

the linear transform

Examples

Run this code
library(mvtnorm)
n = 350
p = 100
D <- diag(p)
Sigma = matrix(0, p, p)
X <- rmvnorm(n,matrix(0, p, 1), Sigma)
decompose.result <- sk.decompose(X, D)
U_perp <- decompose.result$U_perp

Run the code above in your browser using DataLab