project_pca: Project new data in existing PCA space

Description Usage Arguments Value Examples

View source: R/project_pca.R

Description

inputs:

Usage

1
project_pca(Xnew = NULL, pc = NULL)

Arguments

Xnew:

new data matrix to project.

pc:

prcomp object used for projection generated by prcomp. Must have the same number of columns as X.

Value

Xpc: a matrix that corresponds to the projection of X in the PCA space. Xpc[,1] is the first PC, Xpc[,2] the second PC, etc.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
set.seed(1)
data(iris)

#randomly select 100 data points from the iris data set
#to perform PCA
ix=rep(FALSE,nrow(iris))
ix[sample(nrow(iris),100)]=TRUE

# perform PCA
pc=stats::prcomp(iris[ix,1:4])

# plot the PCA
plot_pca(pc,col=iris$Species[ix])

# project the remaining 50 data points (given by !ix) in the existing PCA space
Xpc=project_pca(iris[!ix,1:4],pc)

#plot the new data points in the existing plot
points(Xpc[,1],Xpc[,2])

nchlis/pca.utils documentation built on May 23, 2017, 11:33 a.m.