This function performs unsupervised PCA on each mixscape class separately and projects each subspace onto all cells in the data. Finally, it uses the first 10 principle components from each projection as input to lda in MASS package together with mixscape class labels.
MixscapeLDA( object, assay = NULL, ndims.print = 1:5, nfeatures.print = 30, reduction.key = "LDA_", seed = 42, pc.assay = "PRTB", labels = "gene", nt.label = "NT", npcs = 10, verbose = TRUE, logfc.threshold = 0.25 )
An object of class Seurat.
Assay to use for performing Linear Discriminant Analysis (LDA).
Number of LDA dimensions to print.
Number of features to print for each LDA component.
Reduction key name.
Value for random seed
Assay to use for running Principle components analysis.
Meta data column with target gene class labels.
Name of non-targeting cell class.
Number of principle components to use.
Print progress bar.
Limit testing to genes which show, on average, at least X-fold difference (log-scale) between the two groups of cells. Default is 0.25 Increasing logfc.threshold speeds up the function, but can miss weaker signals.
Returns a Seurat object with LDA added in the reduction slot.