#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", doCE = FALSE)
#
# #No CEInfo available
# irisPCADSC$CEInfo
#
# #Add and show CEInfo
# irisPCADSC <- doCE(irisPCADSC)
# irisPCADSC$CEInfo
## ---------------------------------------------
#Make a partial PCADSC object and only add CE information with no
#bootstrapping (and thus no test)
irisPCADSC_fast <- PCADSC(iris, "group", doAngle = FALSE,
doChroma = FALSE, doCE = FALSE)
irisPCADSC_fast <- doCE(irisPCADSC_fast, B = 100)
irisPCADSC_fast$CEInfo
Run the code above in your browser using DataLab