#load iris data
data(iris)
#Define grouping variable, grouping the observations by whether their species is
#Setosa or not
iris$group <- "setosa"
iris$group[iris$Species != "setosa"] <- "non-setosa"
iris$Species <- NULL
## Not run: ------------------------------------
# #make a partial PCADSC object, splitting the data by "group"
# irisPCADSC <- PCADSC(iris, "group", doChroma = FALSE)
#
# #No chromaInfo available
# irisPCADSC$chromaInfo
#
# #Add and show chromaInfo
# irisPCADSC <- doChroma(irisPCADSC)
# irisPCADSC$chromaInfo
## ---------------------------------------------
#Make a partial PCADSC object and only add chroma information for a
#faster runtime
irisPCADSC_fast <- PCADSC(iris, "group", doAngle = FALSE,
doChroma = FALSE, doCE = FALSE)
irisPCADSC_fast <- doChroma(irisPCADSC_fast)
irisPCADSC_fast$chromaInfo
Run the code above in your browser using DataLab