| big_pls_cox | R Documentation |
Compute Partial Least Squares (PLS) components tailored for
Cox proportional hazards models when predictors are stored as a
big.matrix from the bigmemory package.
big_pls_cox(
X,
time,
status,
ncomp = 2L,
control = survival::coxph.control(),
keepX = NULL
)
X |
A numeric matrix or a |
time |
Numeric vector of survival times. |
status |
Integer (0/1) vector of event indicators. |
ncomp |
Number of latent components to compute. |
control |
Optional list passed to |
keepX |
Optional integer vector specifying the number of variables to retain (naive sparsity) in each component. A value of zero keeps all predictors. If a single integer is supplied it is recycled across components. |
The function standardises each predictor column, iteratively builds latent scores using martingale residuals from Cox fits, and deflates the predictors without materialising the full design matrix in memory. Both in-memory and file-backed bigmemory matrices are supported.
A list with the computed scores, loadings, weights, scaling information and the
fitted Cox model returned by survival::coxph.fit.
Maumy, M., Bertrand, F. (2023). PLS models and their extension for big data. Joint Statistical Meetings (JSM 2023), Toronto, ON, Canada.
Maumy, M., Bertrand, F. (2023). bigPLS: Fitting and cross-validating PLS-based Cox models to censored big data. BioC2023 — The Bioconductor Annual Conference, Dana-Farber Cancer Institute, Boston, MA, USA. Poster. https://doi.org/10.7490/f1000research.1119546.1
Bastien, P., Bertrand, F., Meyer, N., & Maumy-Bertrand, M. (2015). Deviance residuals-based sparse PLS and sparse kernel PLS for censored data. Bioinformatics, 31(3), 397–404. doi:10.1093/bioinformatics/btu660
Bertrand, F., Bastien, P., Meyer, N., & Maumy-Bertrand, M. (2014). PLS models for censored data. In Proceedings of UseR! 2014 (p. 152).
big_pls_cox_gd(), predict.big_pls_cox(), select_ncomp(),
computeDR().
if (requireNamespace("survival", quietly = TRUE)) {
set.seed(1)
X <- matrix(rnorm(100), nrow = 20)
time <- rexp(20)
status <- rbinom(20, 1, 0.5)
fit <- big_pls_cox(X, time, status, ncomp = 2)
str(fit)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.