我们可以使用normCounts(RUVr)输出来计算TPM吗?

dsekswqp  于 2023-02-01  发布在  其他
关注(0)|答案(1)|浏览(193)

我已经使用了RUVSeq R软件包中的RUVr()来校正数据中的批次,现在我想使用校正后的表格来计算TPM。
是否可以使用normCounts()提取TPM的批校正数据?

irlmq6kh

irlmq6kh1#

RUVSeq包中的normCounts函数不是提取批量校正的数据,而是设计用于归一化读取计数。
每个样本的归一化因子可以通过calcNormFactors函数获得,然后您可以使用该函数计算数据的TPM(每百万个转录)值。

library(RUVSeq)

# Your batch-corrected data
bc_data <- # ...

# Obtain normalization factors
factors <- calcNormFactors(bc_data)

# Calculate TPM values
tpm_values <- tpm(bc_data, factors)

相关问题