## Metabolomics starvation data obtained from http://genomics-pubs.princeton.edu/StarvationMetabolomics/Download.shtml
data(StarvationData)
## Computes the eigensystem for the actual data
eigensystem <- compute(StarvationData)
## Exclude the eigenfeatures representing steady-state intensity
eigensystem <- exclude(eigensystem)
## Computes the eigensystem on the variance in the data after filtering out stead-state intensity
eigensystem <- compute(eigensystem, apply="variance")
## No exclusion of eigenfeatures representing steady-scale variance
eigensystem <- exclude(eigensystem, excludeEigenfeatures=0)
## Generate report for eigenfeatures 1 and 2
params <- new("EigensystemPlotParam")
if (.Platform$OS.type %in% "windows") path(params) <- getwd()
report(eigensystem, params)
## Generate report for eigenfeatures 2 and 3
whichPolarAxes(params) <- c(3,2)
report(eigensystem, params)
Run the code above in your browser using DataLab