细胞基因完整矩阵转10xGenomics稀疏矩阵文件

经常分析单细胞数据的同学都知道,10x Genomics 的细胞基因矩阵通常存储在三个关键文件中,这些文件共同描述了单细胞测序数据。这三个文件分别是:

Barcodes.tsv,通常是一列多行,每一行代表一个细胞。

Features.tsv,通常是两列多行,每一行代表一个基因,每行第一个是基因ID,第二个是对应的基因名称。

Features.tsv。通常是三列多行,前两行可以理解为标题。从第三行开始,第一个数n1代表基因(即Features.tsv中第n1行对应的基因),第二个数代表细胞(即Barcodes.tsv中第n2行对应的细胞),第三列是表达量。

通常来说,将这三个文件放在一个文件夹,并且命名正确,就可以直接使用Read10X函数来创建细胞矩阵了。

有的组由于原始数据处理的代码写的不够优雅,搞出的细胞-基因矩阵是tsv格式的也就算了,矩阵的行列都搞反了,搞成了每一列代表一个基因,每一行代表一个细胞,与seurat的RNA矩阵刚好相反。

如果你遇到的矩阵也是以列代表基因,行代表细胞,那么为了将完整矩阵改成10xGenomics的稀疏矩阵格式,你可以采用使用以下R代码将它改成符合10x genomics标准的稀疏矩阵三文件:

library(Matrix)
setwd("修改为输出文件夹路径")
mat1 <- read.table("修改为细胞-基因名矩阵文件,包含文件路径")
mat2 <- read.table("修改为细胞-基因ID矩阵文件,包含文件路径")
mat1 <- as.matrix(mat1)
mat2 <- as.matrix(mat2)

genename <- data.frame(Column1 = colnames(mat2),Column2 = colnames(mat1))
write.table(rownames(mat2), file = 'barcodes.tsv', sep = '\t', row.names = FALSE, col.names = FALSE, quote=F)
write.table(genename, file = "genes.tsv", sep = "\t", row.names = FALSE, col.names = FALSE, quote=F)

mtx1<-Matrix(t(mat1))
writeMM(obj = mtx1, file='matrix.mtx')

相关推荐

  1. 细胞基因完整矩阵10xGenomics稀疏矩阵文件

    2024-01-19 21:48:03       51 阅读
  2. BC107 矩阵

    2024-01-19 21:48:03       46 阅读
  3. Scipy 高级教程——稀疏矩阵

    2024-01-19 21:48:03       49 阅读

最近更新

  1. docker php8.1+nginx base 镜像 dockerfile 配置

    2024-01-19 21:48:03       94 阅读
  2. Could not load dynamic library ‘cudart64_100.dll‘

    2024-01-19 21:48:03       100 阅读
  3. 在Django里面运行非项目文件

    2024-01-19 21:48:03       82 阅读
  4. Python语言-面向对象

    2024-01-19 21:48:03       91 阅读

热门阅读

  1. C语言栈实现就近匹配原则

    2024-01-19 21:48:03       56 阅读
  2. C语言之三子棋游戏(棋盘)

    2024-01-19 21:48:03       53 阅读
  3. SQL注入常用命令详解

    2024-01-19 21:48:03       42 阅读
  4. GoLang刷题之leetcode

    2024-01-19 21:48:03       61 阅读
  5. 【go】依赖倒置demo

    2024-01-19 21:48:03       65 阅读
  6. 网络安全产品之认识入侵检测系统

    2024-01-19 21:48:03       60 阅读
  7. 对docker的一些粗略理解

    2024-01-19 21:48:03       56 阅读
  8. 免费chartGPT网站汇总--

    2024-01-19 21:48:03       47 阅读
  9. Kafka系列(三)

    2024-01-19 21:48:03       52 阅读
  10. Python之归并排序

    2024-01-19 21:48:03       55 阅读
  11. Mysql的in与exits

    2024-01-19 21:48:03       52 阅读