R语言 | R语言分析琐碎

---

首先，文件夹下应该存在三个文件，分别命名为：

• barcodes.tsv.gz
• matrix.mtx.gz
• features.tsv.gz
  然后使用Seurat包中的，Read10X方法读取。

library(Seurat)
data<-CreateSeuratObject(Read10X("./"),"ACC")

如上，即读取成功！

---

如何跳过注释行（如#），读取.txt文件。使用参数：comment.char

data<-read.table("GSE25219-GPL5175_series_matrix.txt",comment.char = "!",header = T)

---

当需要==多个值,应该使用%in%代替，多个值用子集c()表示。

meta.data[,which(meta.data[2,]%in%c("TC","VF","URL","PC","OC"))]

---

• 为什么要进行数据归一化？

为了数据组间可以相互比较，消除量纲上的差异。