prcomp.gds.class {bigmelon} | R Documentation

## Principal Component Analysis for high-dimensional data

### Description

Performs principal components analysis on the given gds object and returns the results as an object of class "prcomp".

### Usage

    ## S3 method for class 'gds.class'
    prcomp(x, node.name, center = FALSE, scale. = FALSE,
    rank. = NULL, retx = FALSE, tol = NULL, perc = 0.01,
    npcs = NULL, parallel = NULL, method = c('quick', 'sorted'), verbose = FALSE, ...)

### Arguments

`x` | A gds.class object.
`node.name` | Name of the gdsn.class node to learn principal components from
`center` | Logical value indicating whether variables should be shifted to be zero centered.
`scale.` | Logical value indicating whether the variables should be scaled to have unit variance
`tol` | a value indicating the magnitude below which components should be omitted.
`rank.` | Old. (Still functional) Number of principal components to be returned
`retx` | a logical value indicating whether the rotated variables should be returned.
`perc` | Percentage of data to be used.
`npcs` | Number of principal components to be returned
`parallel` | Can supply either a cluster object (made from makeCluster) or a integer depicting number of cores to be used. This is only used if method="sorted".
`method` | Indicates what method to use out of "quick" and "sorted". "quick" stochastically selects number of rows according to perc. And the supplies them to svd. "sorted" determines the interquartile range for each row then selects the top percentage (according to perc) of probes with the largest interquartile range and supplies selected rows to svd.
`verbose` | A logical value indicating whether message outputs are displayed.
`...` | arguments passed to or from other methods. If "x" is a formula one might specify "scale." or "tol".

### Details

The calculation is done by a singular value decomposition of the (centered and possibly scaled) data matrix, not by using "eigen" on the covariance matrix. This is generally the preferred method for numerical accuracy. The "print" method for these objects prints the results in a nice format and the "plot" method produces a scree plot.

### Value

An object of prcomp class

### Examples

    data(melon)
    e <- es2gds(melon, "melon.gds")
    prcomp(e, node.name="betas", perc=0.01, method='quick')
    closefn.gds(e)
    unlink("melon.gds")

[Package bigmelon version 1.8.0 Index]