svdvis (version 0.1)

svd.scatter: Visualizing Singular Vectors or Principal Components by Scatterplot Matrices

Description

Creates a set of multiple scatter plots from all pairs of selected singular vectors or principal components. Principal components can be plotted by setting weights = "sv". Since it largely uses ggpairs from the GGally package, optional arguments for ggpairs can be specified.

Usage

svd.scatter(svd.obj, r = NULL, group = NULL, weights = NULL, alpha = 0.7, axisLabels = "none", ...)

Arguments

svd.obj

A list, resulted from applying svd to a dataset, with u, d, and v corresponding to left singular vector, singular values, and right singular vectors, respectively. Alternatively, supply singular vectors, v.

r

A positive integer to use only the first r vectors in visualization. If not specified, all vectors available in svd.obj$v are visualized.

group

A vector of length n, specifying groups (e.g., phenotypes or conditions for n samples).

weights

A vector of length r. If "sv", singular values contained in svd.obj$d[1:r] are used.

alpha

A numeric value for transparency.

axisLabels

Set to either "none" (default), "show", or "internal".

...

Additional arguments to pass onto ggpair.

Value

svd.scatter creates and draws a figure, which is a ggpair object.

See Also

ggpairs

Examples

set.seed(1234)
dat = matrix(rnorm(1000), 100, 10)
svd.obj = svd(dat)
colnames(svd.obj$v) = paste0("V",1:10)
svd.scatter(svd.obj, r=3, group=c(rep("Group1",5), rep("Group2",5)))