predict-methods: Predict values from PCA.

Description Usage Arguments Details Value Author(s) Examples

Description

Predict data using PCA model

Usage

1
2
3
4
5
6
7
## S3 method for class 'pcaRes'
predict(object, newdata, pcs = nP(object), pre = TRUE,
  post = TRUE, ...)

## S4 method for signature 'pcaRes'
predict(object, newdata, pcs = nP(object),
  pre = TRUE, post = TRUE, ...)

Arguments

object

pcaRes the pcaRes object of interest.

newdata

matrix new data with same number of columns as the used to compute object.

pcs

numeric The number of PC's to consider

pre

pre-process newdata based on the pre-processing chosen for the PCA model

post

unpre-process the final data (add the center back etc)

...

Not passed on anywhere, included for S3 consistency.

Details

This function extracts the predict values from a pcaRes object for the PCA methods SVD, Nipals, PPCA and BPCA. Newdata is first centered if the PCA model was and then scores (T) and data (X) is 'predicted' according to : That=XnewP Xhat=ThatP'. Missing values are set to zero before matrix multiplication to achieve NIPALS like treatment of missing values.

Value

A list with the following components:

scores

The predicted scores

x

The predicted data

Author(s)

Henning Redestig

Examples

1
2
3
4
5
6
data(iris)
hidden <- sample(nrow(iris), 50)
pcIr <- pca(iris[-hidden,1:4])
pcFull <- pca(iris[,1:4])
irisHat <- predict(pcIr, iris[hidden,1:4])
cor(irisHat$scores[,1], scores(pcFull)[hidden,1])

Example output

Loading required package: Biobase
Loading required package: BiocGenerics
Loading required package: parallel

Attaching package: 'BiocGenerics'

The following objects are masked from 'package:parallel':

    clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
    clusterExport, clusterMap, parApply, parCapply, parLapply,
    parLapplyLB, parRapply, parSapply, parSapplyLB

The following objects are masked from 'package:stats':

    IQR, mad, sd, var, xtabs

The following objects are masked from 'package:base':

    Filter, Find, Map, Position, Reduce, anyDuplicated, append,
    as.data.frame, basename, cbind, colMeans, colSums, colnames,
    dirname, do.call, duplicated, eval, evalq, get, grep, grepl,
    intersect, is.unsorted, lapply, lengths, mapply, match, mget,
    order, paste, pmax, pmax.int, pmin, pmin.int, rank, rbind,
    rowMeans, rowSums, rownames, sapply, setdiff, sort, table, tapply,
    union, unique, unsplit, which, which.max, which.min

Welcome to Bioconductor

    Vignettes contain introductory material; view with
    'browseVignettes()'. To cite Bioconductor, see
    'citation("Biobase")', and for packages 'citation("pkgname")'.


Attaching package: 'pcaMethods'

The following object is masked from 'package:stats':

    loadings

[1] 0.9999995

pcaMethods documentation built on Nov. 8, 2020, 6:19 p.m.