powered by
O2P(x)
B2P
X <- tcrossprod(matrix(rnorm(100),ncol=10)) # Orthogonal projector based on the first three eigenvectors of X P <- O2P(eigen(X)$vectors[,1:3])
Run the code above in your browser using DataLab