Seurat: NormalizeData()

2023-08-30  本文已影响0人  LET149

1. NormalizeData()

method = LogNormalize :

    1. Feature counts for each cell are divided by the total counts for that cell
    1. Multiplied by the scale.factor (10000)
    1. Log1p()

结果存放在 @assays$RNA@data

2. NormalizeData(Seurat_object, assay = "HTO", normalization.method = "CLR")

https://zhuanlan.zhihu.com/p/477173992

HTO 数据进行 CLR 归一化处理

Seurat_object@assays$HTO@counts中的每一行进行CLR处理,即对同一hashtagUMI值在细胞之间进行归一化

归一化后的值放在 Seurat_object@assays$HTO@data

处理步骤:

kk <- as.numeric(kk@assays$HTO@counts[1,])
data_0 <- log1p(kk[kk > 0])
data_1 <- sum(data_0, na.rm=TRUE)
data_2 <- length(kk)
data_3 <- exp(data_1/data_2)
data_4 <- kk/data_3
log1p(data_4) : CLR 归一化后的值

上一篇下一篇

猜你喜欢

热点阅读