Dir Sirs,
Today, I tried this.
http://mixomics.org/mixdiablo/case-study-tcga/
It did not work now any more (see below).
Last time I downloaded from CRAN.
It did work.
Now from bioconductor. Then id does not work.
When I did
plot(perf.diablo)
It say that no finite numbers of xlim.
I am glad if you can tell me why.
I used R-3.6.1 x64 on windows 10.
Yours, tag
===============
if (!requireNamespace(“BiocManager”, quietly = TRUE))
install.packages(“BiocManager”)
BiocManager::install(“mixOmics”)
library(mixOmics)
data(‘breast.TCGA’)
extract training data
data = list(mRNA = breast.TCGA$data.train$mrna,
miRNA = breast.TCGA$data.train$mirna,
proteomics = breast.TCGA$data.train$protein)
check dimension
lapply(data, dim)
Y = breast.TCGA$data.train$subtype
summary(Y)
design = matrix(0.1, ncol = length(data), nrow = length(data),
dimnames = list(names(data), names(data)))
diag(design) = 0
design
sgccda.res = block.splsda(X = data, Y = Y, ncomp = 5,
design = design)
set.seed(123) # for reproducibility, only when the `cpus’ argument is not used
this code takes a couple of min to run
perf.diablo = perf(sgccda.res, validation = ‘Mfold’, folds = 10, nrepeat = 10)
#perf.diablo # lists the different outputs
plot(perf.diablo)