PCA of variable features vs all genes in integrated analysis
1
0
Entering edit mode
10 months ago
aa123 • 0

In this vignette (ctrl vs stim), does PCA use the 2000 variable features or use all the genes from the integrated analysis?

https://satijalab.org/seurat/v3.2/immune_alignment.html

Thank you.

single cell RNA seq integrated analysis • 358 views
ADD COMMENT
0
Entering edit mode
10 months ago
ATpoint 55k

It is the variable features identified in this command:

ifnb.list <- lapply(X = ifnb.list, FUN = function(x) {
    x <- NormalizeData(x)
    x <- FindVariableFeatures(x, selection.method = "vst", nfeatures = 2000)
})

The idea is to first identify genes with large variance assuming that these will drive separation between cells, then reduce them into principle components and then perform integration in PCA space.

ADD COMMENT
0
Entering edit mode

Thank you.

Is this because the default assay is “integrated”?

The RunPCA command does not specify using variable features through features=VariableFeatures(object=)).

ADD REPLY

Login before adding your answer.

Traffic: 1469 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6