mixOmicsTeam / mixOmics

Development repository for the Bioconductor package 'mixOmics '
http://mixomics.org/
150 stars 51 forks source link

How many datasets should be listed in the new version? #295

Open chunkaili217 opened 1 year ago

chunkaili217 commented 1 year ago

test.keepX = list ( nitrogen = c(1:5), enzyme = c(1:4), soil = c(1:3), otu = c(seq(100, 700, 10)) ) t1 = proc.time() BPPARAM <- BiocParallel::SnowParam(workers = parallel::detectCores()-1) tune.TCGA = tune.block.splsda(X = data, Y = Y, ncomp = ncomp, test.keepX = test.keepX, design = design, validation = 'Mfold', folds = 5, nrepeat = 1, BPPARAM = BPPARAM, dist = "centroids.dist")

list.keepX = tune.TCGA$choice.keepX

final.diablo.model = block.splsda(X = data, Y = Y, ncomp = ncomp, keepX = list.keepX, design = design)

network(final.diablo.model, blocks = c(1,2,3,4), color.node = c('orchid','orange','tomato','sky blue'), layout.auto, cutoff = 0.6)

Error: 'comp' must be either NULL or a list of length 4.