if (FALSE) {
CICA_data <- Sim_CICA(Nr = 15, Q = 5, R = 4, voxels = 100, timepoints = 10,
E = 0.4, overlap = .25, externalscore = TRUE)
#Compute single subject ICAs (nClus equals length(ExampleData))
output <- CICA(DataList = CICA_data$X, nStarts = 1,
nComp = 5, nClus = 9, verbose = FALSE)
RV <- computeRVmat(DataList = output$Sr, dist = TRUE,
verbose = FALSE)
# apply hierarchical clustering on RV output
hcl <- hclust(RV)
plot(hcl)
# low dimensional visualisation using Classical Multidimensional Scaling
mds <- cmdscale(RV)
plot(mds)
}
Run the code above in your browser using DataLab