SDA model trained by tumor gene expression

This function trains a sparse discriminant analysis model using tumor gene expression data. The parameters can be self-determined or determined by the cross-validation. If cross-validation is prefered, users can refer to the function sda_model_cv.

Usage

sda_model(object, stop = NULL, lambda = NULL)

Arguments

object

A CASCAM object.

stop

A negative value with its absolute value representing the number of genes selected for model training.

lambda

A parameter for the L2-norm for elastic net regression.

Value

A CASCAM object with sda_model, camod_norm_data, tumor_norm_data, tumor_sda_project, and camod_sda_project slots.

References

Clemmensen, Line, Trevor Hastie, Daniela Witten, and Bjarne Ersbøll. 2011. “Sparse Discriminant Analysis.” Technometrics: A Journal of Statistics for the Physical, Chemical, and Engineering Sciences 53 (4): 406–13. https://doi.org/10.1198/tech.2011.08118.

Examples

if (FALSE) {
gene_info <- create_InformativeGenes(tumor_ct, tumor_label2, "ILC")
CASCAM_eg <- create_CASCAM(tumor_aligned, tumor_label, camod_aligned, gene_info)
CASCAM_eg <- sda_model(CASCAM_eg)
}