powered by
Get eigen molecule
get.eigen.molecule(data, groups, whichgroups = NULL, methods = "svd", n = 10)
the resulting vector.
a data matrix ([data.frame object] row: molecules, col: samples or replicates)
a vector of group memberships as returned by cutree
a vector of group numbers to examine
c("svd", "nipals", "rnipals", "bpca", "ppca"). See also pca() function in pcaMethods package
top n principal components
Atsushi Fukushima
library(pcaMethods) data(golub, package = "multtest") hc.mol1 <- cluster.molecule(golub[1:100, 1:27], "pearson", "average") g1 <- cutree(hc.mol1, h=0.6) res1 <- get.eigen.molecule(golub[1:100,], g1)
Run the code above in your browser using DataLab