Description Usage Arguments Details Value Author(s) References See Also Examples
Function to combine multiple independent studies measured on the same variables or predictors (P-integration) using variants of multi-group PLS-DA for supervised classification.
1 2 3 4 5 6 7 8 9 10 11 | mint.plsda(
X,
Y,
ncomp = 2,
study,
scale = TRUE,
tol = 1e-06,
max.iter = 100,
near.zero.var = FALSE,
all.outputs = TRUE
)
|
X |
numeric matrix of predictors combining multiple independent studies
on the same set of predictors. |
Y |
A factor or a class vector indicating the discrete outcome of each sample. |
ncomp |
Positive Integer. The number of components to include in the model. Default to 2. |
study |
Factor, indicating the membership of each sample to each of the studies being combined |
scale |
Logical. If scale = TRUE, each block is standardized to zero means and unit variances (default: TRUE) |
tol |
Positive numeric used as convergence criteria/tolerance during the
iterative process. Default to |
max.iter |
Integer, the maximum number of iterations. Default to 100. |
near.zero.var |
Logical, see the internal |
all.outputs |
Logical. Computation can be faster when some specific
(and non-essential) outputs are not calculated. Default = |
mint.plsda
function fits a vertical PLS-DA models with ncomp
components in which several independent studies measured on the same
variables are integrated. The aim is to classify the discrete outcome
Y
. The study
factor indicates the membership of each sample in
each study. We advise to only combine studies with more than 3 samples as
the function performs internal scaling per study, and where all outcome
categories are represented.
X
can contain missing values. Missing values are handled by being
disregarded during the cross product computations in the algorithm
mint.plsda
without having to delete rows with missing data.
Alternatively, missing data can be imputed prior using the
impute.nipals
function.
The type of deflation used is 'regression'
for discriminant algorithms.
i.e. no deflation is performed on Y.
Useful graphical outputs are available, e.g. plotIndiv
,
plotLoadings
, plotVar
.
mint.plsda
returns an object of class "mint.plsda",
"plsda"
, a list that contains the following components:
X |
the centered and standardized original predictor matrix. |
Y |
original factor |
ind.mat |
the centered and standardized original response vector or matrix. |
ncomp |
the number of components included in the model. |
study |
The study grouping factor |
mode |
the algorithm used to fit the model. |
variates |
list containing the variates of X - global variates. |
loadings |
list containing the estimated loadings for the variates - global loadings. |
variates.partial |
list containing the variates of X relative to each study - partial variates. |
loadings.partial |
list containing the estimated loadings for the partial variates - partial loadings. |
names |
list containing the names to be used for individuals and variables. |
nzv |
list containing the zero- or near-zero predictors information. |
iter |
Number of iterations of the algorithm for each component |
prop_expl_var |
Percentage of explained variance for each component and each study after setting possible missing values to zero (note that contrary to PCA, this amount may not decrease as the aim of the method is not to maximise the variance, but the covariance between X and the dummy matrix Y). |
Florian Rohart, Kim-Anh Lê Cao, Al J Abadi
Rohart F, Eslami A, Matigian, N, Bougeard S, Lê Cao K-A (2017). MINT: A multivariate integrative approach to identify a reproducible biomarker signature across multiple experiments and platforms. BMC Bioinformatics 18:128.
Eslami, A., Qannari, E. M., Kohler, A., and Bougeard, S. (2014). Algorithms for multi-group PLS. J. Chemometrics, 28(3), 192-201.
mixOmics article:
Rohart F, Gautier B, Singh A, Lê Cao K-A. mixOmics: an R package for 'omics feature selection and multiple data integration. PLoS Comput Biol 13(11): e1005752
spls
, summary
, plotIndiv
,
plotVar
, predict
, perf
,
mint.pls
, mint.spls
, mint.splsda
and http://www.mixOmics.org/mixMINT for more details.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | data(stemcells)
res = mint.plsda(X = stemcells$gene, Y = stemcells$celltype, ncomp = 3,
study = stemcells$study)
plotIndiv(res)
#plot study-specific outputs for all studies
plotIndiv(res, study = "all.partial")
## Not run:
#plot study-specific outputs for study "2"
plotIndiv(res, study = "2", col = 1:3, legend = TRUE)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.