reconstruct {EESPCA}

Calculates the reduced rank reconstruction

Description

Utility function for computing the reduced rank reconstruction of X using the PC loadings in V.

Usage

    reconstruct(X,V,center=TRUE)

Arguments

X

An n-by-p data matrix whose top k principal components are contained the p-by-k matrix V.

V

A p-by-k matrix containing the loadings for the top k principal components of X.

center

If true (the default), X will be mean-centered before the reconstruction is computed. If the PCs in V were computed via SVD on a mean-centered matrix or via eigen-decomposition of the sample covariance matrix, this should be set to true.

Value

Reduced rank reconstruction of X.

Examples

    set.seed(1)
    # Simulate 10x5 MVN data matrix
    X=matrix(rnorm(50), nrow=10)
    # Perform PCA
    prcomp.out = prcomp(X)
    # Get rank 2 reconstruction
    reconstruct(X, prcomp.out$rotation[,1:2])